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: as has been shown recently in refs . @xcite , an antiproton implanted in a heavy nucleus serves as an attractor for surrounding nucleons that can lead to a sizable increase of the central nucleon density . this effect is caused by the strong attractive scalar and vector potentials acting on the antiproton , as follows from the @xmath1-parity transformation of nuclear potentials @xcite . correspondingly , the antiproton also creates an attractive potential acting on nucleons . this leads to the concentration of nucleons around the antiproton and , as result , to a considerable increase of the nucleon density .    within the relativistic mean field ( rmf ) model , the @xmath1-parity transformed nuclear optical potential is about @xmath2 mev at the normal nuclear matter density @xmath3 @xmath4 , while a phenomenological value of an antiproton optical potential is limited within the range of @xmath5 mev @xcite . therefore , in order to fit the empirical optical potential , the antiproton coupling constants with @xmath6- , @xmath7- and @xmath8-meson fields should be reduced with respect to the values given by the @xmath1-parity transformation . the rmf calculations with reduced coupling constants @xcite still show quite strong compressional effects for light and medium nuclei . an important question , which arises here is whether the compression process is fast enough to develop before the @xmath9-annihilation . the total @xmath10-annihilation cross section in vacuum can be parameterized at low relative velocities @xmath11 as @xmath12 where @xmath13 mb and @xmath14 mb@xmath15c @xcite . using these numbers we can estimate the life time of an antiproton inside the nuclear matter at normal density : @xmath16 this is , of course , a very short time in nuclear scale . however , as argued in ref . @xcite , this time can become much longer , up to 20 fm / c , for deeply bound antiprotons due to the phase space suppression factors . therefore , the compression effects can , in - principle , show up in @xmath9-nuclear interactions . in the present work , we apply a dynamical transport model in order to study the formation and decay of the compressed @xmath9-nuclear system . our calculations are based on the giessen boltzmann - uehling - uhlenbeck ( gibuu ) model @xcite , which has been recently supplemented by the relativistic mean fields @xcite . apart from collision terms , the gibuu model solves the coupled ( through the mean fields ) vlasov equations for nucleon and antiproton phase space distribution functions . as well known @xcite , the vlasov equation provides a semiclassical limit of the time - dependent hartree - fock calculations . thus , the compressional effects found in refs . @xcite should also be reproduced as a static solution of the coupled vlasov equations . it will be demonstrated that the compression process is characterized by the time scale which is comparable with the @xmath9 life time in nuclear environment . thus , @xmath9 has , indeed , a chance to annihilate inside the compressed nucleus . we will show , that the @xmath9-annihilation in a compressed nucleus should lead to the collective expansion of the residual nuclear system . the appearance of the high - energy tails in the kinetic energy spectra of the emitted nucleons is predicted . the distributions in the total invariant mass of produced mesons reveal a noticeable shift toward lower invariant masses , when the annihilation takes place inside the compressed nucleus . the annihilation of slow antiprotons inside heavy nuclei was , first , proposed by rafelski @xcite as a unique opportunity to study nuclear matter in unusual conditions . later , cahay et al . @xcite studied the @xmath9 annihilation inside nuclei within an intranuclear cascade model . in ref . @xcite , antiproton annihilation events into pions at the center of @xmath17ca and @xmath18ag nuclei were simulated . the mean field effects were , however , completely neglected in @xcite .    in sect . ii , we describe the theoretical model applied in calculations . sect . iii contains the results of the time evolution study for the compression and explosion dynamics . in sect . iv , we propose several observable signals sensitive to the @xmath9-annihilation in the compressed nucleus . the summary and outlook are given in sect . in calculations , we apply the gibuu model developed in giessen university . for the detailed description and related references , we refer the reader to the web page @xcite , where the new version of the model is presented . below , we mostly describe the new features implemented in the present work .      below we consider a system composed of an antinucleon interacting with baryons . this system is described by the rmf lagrangian of the following form @xcite : @xmath19 \psi_j                                            \nonumber \\                 & + & { 1 \over 2}\partial_\mu\sigma\partial^\mu\sigma - u(\sigma )                   -   { 1 \over 4 } f_{\mu\nu } f^{\mu\nu } + { 1 \over 2 } m_\omega^2 \omega_\mu\omega^\mu~ ,                                                                                \label{lagr}\end{aligned}\ ] ] where @xmath20 are the baryon ( @xmath21 ) and antinucleon ( @xmath22 ) fields , respectively ; @xmath6 is the isoscalar - scalar meson field ( @xmath23 , @xmath24 ) ; @xmath25 is the isoscalar - vector meson field ( @xmath26 , @xmath27 ) ; and @xmath28 . here @xmath29 and @xmath30 denotes , respectively , the isospin 1/2 and 3/2 nonstrange baryonic resonances , and @xmath31 stands for the @xmath32 baryons explicitly propagated in the gibuu model @xcite . in the case of the spin 3/2 , 5/2 and 7/2 baryonic resonances , their fields @xmath20 carry also one or more vector indices , which are dropped in eq.([lagr ] ) and below for brevity . when appropriate , the covariant summation is assumed over these indices . for simplicity , the isovector and electromagnetic terms are disregarded in ( [ lagr ] ) . the selfinteractions of the @xmath6-field are included in ( [ lagr ] ) via the term @xmath33 in order to avoid an unrealistically high compressibility coefficient of the nuclear matter @xcite : @xmath34    some comments are in order to gain more insight into eq . ( [ lagr ] ) . following ref . @xcite , the antinucleon field @xmath35 in the lagrangian density ( [ lagr ] ) is represented in terms of wave functions of _ physical antinucleons_. these wave functions can be obtained by the @xmath1-parity transformation acting on the wave functions of the dirac sea nucleons ( see ref . @xcite for details ) , which appear in the relativistic description of the nucleon @xcite . by applying the same transformation , the nonlinear rmf lagrangian of refs . @xcite ( neglecting terms responsible for the baryon - antibaryon annihilation ) can be expressed as ( [ lagr ] ) with the following relations between coupling constants : @xmath36 the relations ( [ gpar_relation ] ) are satisfied if the physical system would be exactly symmetric with respect to the @xmath1-parity transformation . however , this is not necessary to be true in a many - body system @xcite . the reason is that the concept of the @xmath1-parity symmetry is strictly applicable on the level of the elementary processes only . however , the rmf lagrangian ( [ lagr ] ) is dealing with the effective interactions , which are usually tuned to describe the bulk properties of the nuclear medium and/or the properties of some selected nuclei . due to the many - body effects , such as the pauli blocking or mixed scalar - vector terms in the scattering amplitudes , these effective interactions may not obey the exact @xmath1-parity symmetry anymore . to take into account possible deviations from the @xmath1-parity symmetry , we introduce an overall scaling of the antinucleon - meson coupling constants with respect to the values given by ( [ gpar_relation ] ) ( see ref . @xcite ) : @xmath37 where @xmath38 is a scaling factor .    throughout the paper , we consider two options for the scaling factor of the antinucleon - meson coupling constants : @xmath39 , motivated by the @xmath1-parity , and @xmath40 , which is in a better agreement with the empirical @xmath9a optical potential . for other baryonic fields we put in the present work , for simplicity , the same coupling constants as for the nucleon  : @xmath41 , @xmath42 . all calculations have been performed emloying the nl3 parameterization @xcite of the rmf model . this parameterization provides quite reasonable nuclear matter properties : the binding energy 16.299 mev / nucleon , the compressibility coefficient @xmath43 mev and the nucleon effective mass @xmath44 at @xmath45 . moreover , the nl3 parameterization reproduces the ground state properties of spherical and deformed nuclei very well @xcite . the dirac equations of motion for baryons have the following form : @xmath46 where @xmath47 is the effective ( dirac ) mass .    within the mean field approximation the @xmath6- and @xmath7-fields are treated classically . they satisfy the ( nonlinear ) klein - gordon - like equations with the source terms due to coupling to baryons and an antinucleon : @xmath48 where @xmath49 is the partial scalar density and @xmath50 is the partial baryon current . equation ( [ kgeqomega ] ) has to be supplemented by the four - transversality condition @xmath51      instead of solving the dirac equations ( [ diraceqb ] ) , we will describe the baryons and antinucleon dynamics by the coupled set of the semiclassical kinetic equations @xcite : @xmath52 f_j(x,{\bf p^\star } )     = i_j[f_b , f_m]~ ,                                                                                \label{kineq}\ ] ] where @xmath53 ; @xmath54 ; @xmath55 ; and @xmath56 is the distribution function ( df ) in a six - dimensional phase space @xmath57 with @xmath58 being the spatial components of the kinetic four - momentum @xmath59 the baryons and antinucleon are assumed to be on the respective effective mass shells : @xmath60 the l.h.s . of eq . ( [ kineq ] ) describes the propagation of the @xmath61-th type particles in the classical @xmath6- and @xmath7-fields . the r.h.s . ( [ kineq ] ) is a collision integral , which represents the ( in)elastic two - body collisions with corresponding vacuum cross sections as well as the resonance decays . the complete description of the collision integral structure , in - particular , the differential elementary cross sections included into the gibuu model can be found in @xcite and in refs . therein . the in - medium modification of the baryon - baryon and baryon - meson cross sections is neglected in the present work . we will apply the full kinetic equations , including collision terms , only to describe the _ post - annihilation _ evolution of a system . by this reason , the antiproton df is excluded from the collision integral . instead , we enforce @xmath9 to annihilate into mesons at some preselected time ( see sect . iie ) . thus , in the present work the collision integral includes the nucleon , @xmath62 and higher baryon resonances up to the mass of 2 gev , which can be excited in the meson - baryon and baryon - baryon collisions . a possible hyperon formation in the processes @xmath63 and @xmath64 is included too . the `` valence mesons '' @xmath65 are explicitly taken into account . they are assumed to propagate freely between collisions , i.e. we neglect the mean field potentials acting on these mesons . the scalar density and the baryon current of the @xmath61-th type baryons are expressed in terms of df as follows : @xmath66 where @xmath67 is the spin - isospin degeneracy factor ( @xmath68 etc . ) .    one can show @xcite that the kinetic equations ( [ kineq ] ) with the @xmath6- and @xmath7-fields evolving according to eqs . ( [ kgeqsigma]),([kgeqomega ] ) lead to the continuity equations @xmath69 and the energy - momentum conservation @xmath70 where the energy - momentum tensor is written as @xmath71 here we have also included possible contributions of the `` valence '' mesons @xmath72 , which can be produced at the annihilation . it is assumed that @xmath73 for the valence mesons .      in order to solve eq . ( [ kineq ] ) numerically , df is represented by the set of point - like test particles : @xmath74 where @xmath75 is the number of physical particles of the type @xmath61 and @xmath76 is the number of test particles per physical particle ( the same for all types @xmath61 ) . the test particle positions @xmath77 and kinetic momenta @xmath78 are evolving in time according to the following equations : @xmath79 with @xmath53 and @xmath54 . it is easy to check that df ( [ tp ] ) with @xmath77 and @xmath80 satisfying eqs . ( [ rdot]),([pstardot ] ) gives a formal solution of the vlasov equation in the case when the collision integral in ( [ kineq ] ) is equal to zero . equations ( [ rdot]),([pstardot ] ) are equivalent to the hamiltonian equations of motion for the test particle positions @xmath77 and canonical momenta @xmath81 : @xmath82 where @xmath83 is the single - particle energy ( see @xcite ) . however , it is more convenient to propagate in time the test particle kinetic momenta rather than the canonical ones , since then eq .  ( [ kgeqsigma ] ) for the @xmath6-field decouples from eq . ( [ kgeqomega ] ) for the @xmath7-field .    when the collision integral in eq . ( [ kineq ] ) is taken into account , the test particles are propagated between the two - body collisions using eqs . ( [ rdot]),([pstardot ] ) . all calculations have been performed in the parallel ensemble mode . in this mode , the two - body collisions are permitted between the test particles belonging to the same parallel ensemble only , while the mean field is averaged over @xmath76 parallel ensembles of the test particles propagated simultaneously ( see eq . ( [ tp ] ) ) . therefore , a single parallel ensemble can be considered as a physical event .    in actual calculations , we have neglected the time derivatives of the meson fields in eqs . ( [ kgeqsigma]),([kgeqomega ] ) . however , the spatial derivatives were treated without any simplifying assumptions . the reason for such a strategy is that we are dealing with nuclear systems which have large density gradients , but evolving slowly , as compared with the spatial and temporal scales involved in the mesonic equations of motion . indeed , including the temporal gradients would lead to the frequent oscillations of the mesonic fields with a period of less than @xmath84 , where @xmath85 is the meson mass this gives the period of 2.5 fm / c ( 1.5 fm / c ) for the @xmath6- ( @xmath7- ) field . by taking into account the finite wave lenghts of these oscillations would further reduce the periods . the treatment of such oscillations would strongly complicate the numerical calculations , in - particular , due to the classical meson field radiation . on the other hand , the characteristic periods of the oscillations are significantly smaller than the characteristic compression times ( @xmath0 fm / c , see sect . iii a below ) . therefore , one can approximately average - out the mesonic fields with respect to these oscillations , that is actually assumed in our model . the @xmath6- and @xmath7-fields are , therefore , calculated from the equations @xmath86 within the same approximation , the energy - momentum tensor has the following form : @xmath87 the factor @xmath88 in eq . ( [ twotd ] ) reflects the fact , that due to the omission of the time derivatives of the meson fields in the lagrangian density , only the first term in the r.h.s . contributes to the three - momentum density @xmath89 ( @xmath90 ) . although eqs . ( [ eqsigma]),([eqomega ] ) are not covariant , they provide a better description of the nuclear surface than pure local fields @xcite . this improves the stability of a nuclear ground state and is more appropriate for studying nuclear response to external hadronic and electromagnetic probes . equations ( [ eqsigma]),([eqomega ] ) have been solved numerically by applying the alternating direction implicit iterative method of douglas described in ref . @xcite . due to the scalar density dependence on the effective mass ( see eqs.([mstar]),([rhosj ] ) ) , additional iterations are needed to solve eq.([eqsigma ] ) . in other words , the scalar density has to be computed selfconsistently . to evaluate the meson fields , we used a uniform grid in coordinate space with steps @xmath91 . for the systems @xmath9@xmath92o , @xmath9@xmath17ca and @xmath9@xmath93pb considered below , the grid covered a cubic volume with the side of 10 , 20 and 30 fm , respectively , centered at the center - of - mass ( c.m . ) of a @xmath9a system . by numerical reasons , the @xmath94-functions in coordinate space , introduced in eq . ( [ tp ] ) , have been replaced by the gaussians of the width @xmath95 : @xmath96 the width of the gaussian and the grid step sizes are pure numerical parameters which should resolve the coordinate space nonuniformities of the system . in our case the characteristic space scale is given by the radius of the smallest considered nucleus @xmath92o , i.e. @xmath97 fm . on the other hand , in order to have smooth density distributions , the number @xmath76 of test particles per physical particle ( see eq . ( [ tp ] ) ) should be correlated to the width of the gaussian as @xmath98 . this puts a restriction on too small width due to cpu time increase . as an optimum choice , we fixed in the present work @xmath99 fm . we have realised , however , that there is a rather moderate tendency of increasing maximum compression ( see discussion in sect . iii a below ) with decreasing gaussian width . the number of test particles per nucleon was set to @xmath100 in the most of calculations . the equations of motion ( [ rdot]),([pstardot ] ) have been solved by applying the second - order in time predictor - corrector method @xcite with the time step of 0.1 fm / c . this value is small enough to resolve the time scale of a few fm / c for the compression processes ( see figs . [ fig : rho_u_vs_z ] and [ fig : rho_vs_t ] below ) . we have checked , that taking smaller time step does not influence the results . the full numerical scheme conserves the total energy with the accuracy of about 5% of the initial total binding energy of the @xmath9a system . the nucleons were distributed in coordinate space according to the woods - saxon density profile . the momenta of nucleons were sampled according to the local fermi distribution . the initial antiproton df was chosen as a gaussian wave packet in coordinate and momentum space @xcite located at the center of a nucleus ( @xmath101 ) : @xmath102 where @xmath103 is the width in coordinate space . equation ( [ cohstate ] ) implies , that the antiproton is at rest . the width of the initial antiproton distribution in momentum space is @xmath104 , which follows from the uncertainty relation . if not mentioned explicitly , the calculation is done with the choice @xmath105 fm . this value agrees with results of the static rmf calculations of ref . @xcite . as for nucleons , the antiproton df ( [ cohstate ] ) is projected onto test particles according to eq . ( [ tp ] ) with the @xmath94-functions in coordinate space replaced by gaussians . to avoid misunderstanding , we note that one should distinguish the width @xmath103 of the physical antiproton spatial distribution in eq . ( [ cohstate ] ) and the width of the test particle gaussian .      after the initialization , the system of nucleons and antiproton was propagated in time according to eqs . ( [ rdot]),([pstardot ] ) . the meson fields have been calculated by eqs . ( [ eqsigma]),([eqomega ] ) with the source terms given by the scalar densities ( [ rhosj ] ) and the baryon currents ( [ jbj ] ) . in such a way , the evolution of the system toward compressed state has been followed . in this calculation , the collision term in the r.h.s . of the kinetic equation ( [ kineq ] ) has been set to zero , i.e. we considered a pure mean - field vlasov dynamics . this was done to see most clearly the role of the mean fields . an introduction of the @xmath106 and @xmath107 elastic collisions would mainly lead to a dissipation of the collective energy into heat . as pointed out in ref . @xcite , this effect is rather small and , therefore , can not change significantly the compression dynamics . the reason is , that the elastic collisions are not frequent on the time scale of compression ( see figs . [ fig : rho_u_vs_z ] , [ fig : rho_vs_z ] and [ fig : rho_vs_t ] ) . indeed , the mean time @xmath108 between nucleon - nucleon collisions can be estimated as @xmath109 , where @xmath110 mb is the elastic nucleon - nucleon cross section ( c.f . @xcite ) and @xmath111 is the fermi velocity . this gives @xmath112 fm / c for the nucleon density @xmath113 . the pauli blocking effect will further increase @xmath108 . a similar estimate can also be done for @xmath106 elastic collisions .    at certain time moment @xmath114 , which is an external parameter to our model , we simulated the annihilation of an antiproton . this implies that annihilation occurs instantaneously , as a single quantum mechanical transition , in distinction to description of this process via the collision term in a kinetic equation . in the last case , the antiproton distribution function would gradually disappear on the way to the compressed state . the purpose of the present work is to look at the strongest possible effect of the nuclear compression on observables . therefore , we let the compressed system to be formed , and simulate the sudden annihilation afterwards . the ambiguity in the in - medium annihilation cross sections is taken into account by varying the parameter @xmath114 .    in the actual calculations , the annihilation was simulated as follows : for each antiproton test particle , the closest in coordinate space nucleon test particle was chosen to be the annihilation partner . at large enough values of the total in - medium c.m . energy @xmath115 of the annihilating @xmath116 pair ( see below ) , the annihilation event of the test particle pair into mesons was simulated using the quark model @xcite , which has been already implemented in the gibuu model @xcite earlier . a quark and an antiquark with the same flavour are assumed to annihilate and transfer their total four - momentum to the remaining ( anti)quarks . the remaining four ( anti)quarks form two orthogonal @xmath117 jets with equal energies in the c.m . frame . the jets were hadronized via the lund string fragmentation model @xcite in the jetset version included into the pythia 6.225 program package . the applied annihilation model corresponds to the r2 type diagram in classification of ref . @xcite , i.e. to the quark rearrangement with one @xmath118 annihilation vertex . in this sense , the model has some similarity with the two - meson doorway models of refs . @xcite . to illustrate how the model works we have performed simulations of the @xmath119 annihilation . [ fig : pimul_ppbar ] shows the pion multiplicity distribution for the @xmath119 annihilation at rest in vacuum compared to the data compilation from refs . the calculated distribution is somewhat shifted to smaller pion multiplicities with respect to the data : the calculated average pion multiplicity @xmath120 compared with the experimental value of @xmath121 . we would like to remark , that non - vanishing contribution of the @xmath122 channel in calculations is completely due to the final states with other particles : @xmath123 ( 78% ) , @xmath124 ( 14% ) , @xmath125 ( 6% ) and @xmath126 + photons ( 2% before @xmath127 decay ) . in calculations , we took into account @xmath127 decays into @xmath128 or into final states with pions and disregarded photons afterwards . however , it is not clear to us how photons were counted in the data ( see also ref . @xcite ) .    fig . [ fig : pimom_ppbar ] shows the calculated charged pion momentum distributions in the c.m . frame of the annihilating @xmath119 pair at rest in vacuum . from the partial contributions of the channels with various pion multiplicities we observe , as expected , that the hard(soft ) part of the total momentum distribution is populated mainly by the low(high ) pion multiplicity events . the experimental data are described reasonably well , except for the momenta @xmath129 gev / c , where the calculations significantly overestimate the data . we believe that the accuracy of the model in describing the data in figs . [ fig : pimul_ppbar ] and [ fig : pimom_ppbar ] is sufficient for the exploratory studies of global observables in the present work . certainly , the improvement of the annihilation model is needed to perform more detailed study of the mesonic final states in the annihilation . below we concentrate more on the in - medium effects on the annihilation .    due to the mean field , the invariant energy of the annihilating @xmath116 pair can be substantially below the vacuum threshold value of @xmath130 . this makes the direct application of the jetset model for the @xmath9-annihilation in nuclei physically and numerically problematic . to overcome this difficulty , we introduced the corrected invariant energy as follows @xcite : @xmath131 where @xmath132 . the quantity @xmath133 is a vacuum analog of the total in - medium invariant c.m . energy @xmath115 with @xmath134 . provided that @xmath135 , we have used @xmath133 in the jetset simulation in order to produce the mesonic final states . this lower limit of @xmath115 is due to the fact that the jetset model does not generate enough direct @xmath136 and @xmath137 annihilation final states .    in order to take into account the in - medium effects , in - particular , to ensure the correct in - medium threshold condition @xmath138 , where @xmath139 are the vacuum masses of the produced mesons , the annihilation event was accepted with the probability @xmath140 where @xmath141 is the invariant phase space volume , @xmath142 are the four - momenta of the produced mesons satisfying the vacuum mass shell conditions @xmath143 , @xmath144 . finally , the three - momenta of the produced mesons in the c.m . frame of the annihilating @xmath116 pair were multiplied by the common factor adjusted to get the correct in - medium total c.m . energy @xmath115 .    in - fact , the way we simulate the in - medium effects eq.([p ] ) implies using the vacuum matrix elements of the annihilation channels , which are given by the jetset model , while taking into account the in - medium effects in the phase space factors only . similar procedures have been applied earlier in refs . @xcite .    at @xmath145 , the final @xmath136 or @xmath137 channel was chosen by monte - carlo according to the probability ratio @xmath146 where @xmath147 is the ratio of the @xmath136 and @xmath137 final state probabilities for the @xmath119 annihilation at rest ( see table vi in ref . @xcite ) . for the zero total charge @xmath148 of the annihilating @xmath116 pair , the charge states of the outgoing pions were also determined from the data compilation of ref . since for @xmath149 the data are absent , the charges of the @xmath137 final states were determined by assuming that the @xmath150 and @xmath151 final channels have equal probabilities . the momenta of the outgoing pions were distributed microcanonically according to the available two- or three - body phase space . after the annihilation is simulated , the residual nucleons and produced mesons were propagated in time according to the full kinetic equations ( [ kineq ] ) , including both the baryonic mean fields and collision integrals . this takes into account the entropy production caused by the two - body collisions at the expansion stage . moreover , important processes of the meson rescattering and absorption , e.g. @xmath152 or @xmath153 , @xmath154 are included in the collision integral . these processes influence the observed particle spectra . as demonstrated in refs . @xcite by static rmf calculations , a deeply - bound antiproton - nucleus system can be significantly compressed as compared with a normal nucleus . now we want to study the real dynamics of such a system starting from the unperturbed nuclear ground state at @xmath155 .    fig . [ fig : rho_u_vs_z ] ( top panels ) shows the nucleon and antiproton density profiles calculated at different times along the axis @xmath156 drawn through the center of the @xmath9@xmath17ca system . [ fig : rho_u_vs_z ] ( bottom panels ) also shows the nucleon and antiproton potentials @xmath157 , @xmath158 , along the same axis . left and right panels present results for @xmath40 and @xmath159 , respectively . we see that the initial configuration is unstable and the system starts to shrink . both nucleon and antiproton central densities grow quite fast , reaching their maxima within several fm / c . in the course of the compression process , the nucleon potential becomes deeper in the case of the reduced antiproton coupling constants ( @xmath40 ) and does not , practically , change in the @xmath1-parity motivated case ( @xmath159 ) . the antiproton potential deepens quite strongly with time for the both sets of the antiproton coupling constants .    in the case of @xmath40 , the first maximum of the central nucleon density ( @xmath160 @xmath4 ) is reached at @xmath161 fm / c . at later time the system rebounds and oscillates approaching gradually a static configuration with the nucleon density @xmath162 @xmath4 at the center . since the annihilation is switched off in this calculation , the compressed configuration may exist , in principle , infinitely long time . however , due to numerical reasons , the stability is destroyed by a gradual test particle escape from a box in the coordinate space , where the mean field is computed . nevertheless , the numerical accuracy is good enough to trace the stable system up to at least @xmath163 fm / c .    in the case of @xmath159 , the compression process is much faster than in the case of @xmath40 . already at @xmath164 fm / c we observe the first maximum of the central nucleon density with @xmath165 @xmath4 . a smaller value @xmath166 @xmath4 is reached asymptotically after some oscillations .    in fig . [ fig : rho_vs_z ] we compare time evolution of the nucleon density distribution along the central axis @xmath156 for the light ( @xmath9@xmath92o ) and heavy ( @xmath9@xmath93pb ) systems . for @xmath9@xmath92o , the bell - like shape of the density distribution is reached quite fast . however , in the case of @xmath9@xmath93pb , we observe a quickly growing peak in the center , while peripheral nucleons still do not react on the compression . this leads to the delayed shape equilibration via a complicated compression - decompression cycle . [ fig : rho_vs_t ] shows the time evolution of the central nucleon density for the three systems : @xmath9@xmath92o , @xmath9@xmath17ca and @xmath9@xmath93pb . the case without annihilation is shown by the dotted lines . for a comparison , we also present the central density time evolution in the respective ground state nuclei without an antiproton inside ( dashed lines ) . we see , that at long enough times of the order of several tens fm / c , the static compressed configuration is indeed reached . the small oscillations of the central density in the compressed system visible at @xmath167 fm / c are approximately of the same amplitude as the oscillations of the respective ground state . thus , the reason for these small amplitude oscillations is the fermionic ground state instability due to the classical treatment of particles @xcite . in principle , this instability can be removed by either employing the pauli potential as in refs . @xcite or by adding a friction force to the hamiltonian equations of motions for the test particles @xcite . such modifications are beyond the scope of the present work . we expect that they would not essentially modify the density profiles of the compressed configurations , which are in an overall agreement with the previous static hartree calculations of refs . @xcite . one can also notice ( see lower panels of fig . [ fig : rho_vs_t ] ) a peculiar feature of the @xmath9@xmath93pb system : the dip in the central nucleon density at @xmath168 fm / c . this is mostly a consequence of the delayed shape equilibration mentioned above in discussing fig . [ fig : rho_vs_z ] . the especially strong density drop for @xmath40 is partly caused by the symmetry loss due to the finite number of test particles . the dip is absent in lighter systems , since the shape equilibration for them is much faster ( @xmath169 fm / c , see upper panels in figs . [ fig : rho_u_vs_z ] and [ fig : rho_vs_z ] ) . since the compression time is of primary importance , we have also studied the sensitivity of our results to the width @xmath103 of the initial antiproton df ( [ cohstate ] ) . we have found , that for a larger ( smaller ) width the compression time becomes somewhat longer ( shorter ) . in - particular , for the @xmath9@xmath92o system at @xmath159 the time needed to reach the first density maximum is 5 fm / c ( 2.5 fm / c ) for @xmath170 fm ( @xmath171 fm ) . for the same system at @xmath40 the first density maximum is reached at 10 fm / c ( 3 fm / c ) for @xmath170 fm ( @xmath171 fm ) . the value of the maximum density reached in the compression process is practically insensitive to @xmath103 . we have to admit also , that there is some numerical uncertainty in our calculations due to the choice of the width @xmath95 of the test particle gaussian and the grid step size . e.g. in the calculation with @xmath172 fm for the lightest system @xmath9@xmath92o the maximum and saturation densities are @xmath173 % higher than in calculation with @xmath99 fm . setting the smaller grid step is not feasible for technical reasons . overall , this numerical uncertainty is comparable to the one due to different choices of the scaling factor @xmath174 in our calculations . next , we study the dynamics of a residual nucleus after sudden annihilation of an antiproton . the annihilation was simulated as described in sect . iie . for each considered @xmath9a system and the scaling factor @xmath174 , three different annihilation times @xmath114 have been chosen : they correspond to ( i ) the early ( @xmath175 ) annihilation from a non - compressed ground - state nucleus , ( ii ) the annihilation at the time moment when the first maximum of the central density is reached , and ( iii ) the late annihilation from an asymptotic compressed configuration . one can see from fig . [ fig : rho_vs_t ] , that the annihilation from a non - compressed ground state nucleus ( thin solid lines ) does not lead to significant expansion of the residual nuclear system . the central nucleon density stays always below but close to @xmath45 in this case . in the case of annihilation from the compressed configurations ( thick solid and dash - dotted lines ) , we observe that , for the light systems @xmath9@xmath92o and @xmath9@xmath17ca , the central nucleon density decreases sharply after annihilation and reaches values well below @xmath45 . this is a clear indication of the collective expansion of a system from the initially compressed state . on the other hand , for the heavy @xmath9@xmath93pb system , the expansion is not very pronounced at any choice of the annihilation time . it is interesting , that if the annihilation is switched on at the first density maximum ( thick solid lines ) , then after an abrupt falling down the central density stays for some time @xmath176 fm / c close to @xmath45 before decreasing further . this is explained by an inertial compression : after annihilation , the periphery of a residual nucleus still continues to move to the center during some time until rebound . in the case of the @xmath9@xmath17ca system , this is demonstrated in figs . [ fig : vrad_ca40_xip03 ] and [ fig : vrad_ca40_xip1 ] , where we show the baryon density and the radial collective velocity , @xmath177 at several times as a function of the radius @xmath178 . the collective velocity has been determined as @xmath179 at @xmath180 fm , the evolution of the radial collective velocity field is noticeably influenced by the fermionic ground state instability discussed above . nevertheless , one can still observe the inertial compression . indeed , the radial collective velocity at @xmath181 fm / c for @xmath40 and at @xmath182 fm / c for @xmath159 reveals the fast outward motion at the center , while the peripheral nucleons still continue to move to the center . this explains plateaus in the central density evolution near @xmath183 fm / c in fig . [ fig : rho_vs_t ] . at @xmath184 fm / c the whole system starts to expand . this is reflected in the monotonically increasing radial collective velocity with radius . the especially strong rise of @xmath185 at large @xmath178 is due to emission of fast particles . at later times @xmath186 fm / c the expansion is replaced by the inward motion of the matter in the central zone . however , the fast particles are still continuing to escape from the dense region . we expect that in reality the system will break - up into fragments before the inward motion will start ( see discussion sect . it is presently well established ( see @xcite and refs . therein ) that nuclear matter at low densities ( @xmath187 ) becomes unstable with respect to small density perturbations , so called spinodal instability . however , in order these density perturbations to develop into nuclear fragments , the system must stay long enough time @xmath188 fm / c in the spinodal region . one can see from fig . [ fig : rho_vs_t ] , that the light systems @xmath9@xmath92o and @xmath9@xmath17ca spend a long time in this region . therefore , the residual nuclear systems can undergo a multifragment break - up , if the annihilation happens in the compressed configurations . in other words , the multifragment break - up of nuclei after the @xmath9-annihilation may serve as a signal of the compression prior the annihilation .    in table [ tab : fragmenting_source ] we collect the estimated parameters of fragmenting sources for the @xmath9@xmath92o and @xmath9@xmath17ca systems . the case of @xmath9-annihilation from the state of maximum central density is considered here ( see thick solid lines in fig . [ fig : rho_vs_t ] ) . the sources have been determined by selecting nucleons in the space region where the baryon density is larger than @xmath189 . they are characterized by the neutron ( @xmath190 ) and proton ( @xmath191 ) numbers , the collective kinetic energy per nucleon ( @xmath192 ) and the residual excitation energy per nucleon ( @xmath193 ) . these parameters are defined as @xmath194 where @xmath195 and @xmath196 are the neutron and proton densities , respectively ; @xmath197 ; and @xmath198 is the ground state energy per nucleon of a nucleus with neutron number @xmath190 and proton number @xmath191 computed within our model . the collective kinetic energy ( [ e_kin^coll ] ) is calculated neglecting pressure effects . due to the initially isospin - symmetric nuclei and the neglect of the isovector mesons in the lagrangian density ( [ lagr ] ) , we obtained in all cases @xmath199 in the source . it turned out also , that the coulomb interaction optionally included in some of the calculations does not change this result . thus , only the charge numbers @xmath191 are given in table [ tab : fragmenting_source ] . the time moment for determination of the source parameters has been chosen at @xmath200 fm / c , when the central nucleon density is about @xmath201 , i.e. inside the spinodal region . the calculated residual excitation energy is typically @xmath202 mev / nucleon that corresponds to the temperatures @xmath203 mev . according to the statistical multifragmentation model @xcite , in this energy domain the multifragment break - up is the dominating decay channel of residual nuclei . the collective kinetic energy is @xmath204 mev / nucleon in the case of @xmath9@xmath92o . this is well above the coulomb energy of the source , which is only about 0.4 mev / nucleon for @xmath205 . unfortunately , such a source is too small to experience the real multifragment break - up , rather a fermi break - up into small clusters @xcite . in the case of larger sources , produced in @xmath9@xmath17ca annihilation the collective kinetic energy is considerably smaller , @xmath206 mev / nucleon , but still significant with respect to the total coulomb energy @xmath207 mev / nucleon for @xmath208 . thus , we expect some signs of collective expansion to be visible in kinetic energy spectra of produced fragments . let us now consider other observable effects . [ fig : ekin_spectra ] shows the c.m . kinetic energy spectra of the nucleons emitted from the @xmath9@xmath92o , @xmath9@xmath17ca and @xmath9@xmath93pb systems after the @xmath9-annihilation . in order to separate emitted nucleons from the bound nucleons of a residual nucleus , we used a simple criterion : only those nucleons , both protons and neutrons , were included in the spectra which are separated by at least 3 fm from the other test particles of a given parallel ensemble at @xmath163 fm / c . one can see , that nucleons with the kinetic energy @xmath209 , where @xmath210 mev is the fermi energy of the nuclear matter at @xmath45 , are abundantly emitted . such nucleons are knocked - out from the nucleus by the mesons produced after the annihilation @xcite .    in table [ tab : slopes ] we list the slope parameters @xmath211 of the nucleon kinetic energy spectra obtained by the maxwell - boltzmann fit @xmath212 in the region of @xmath213 mev . we would like to mention , that the authors of ref . @xcite report the slope temperature of about 60 mev for the kinetic energy spectrum of the emitted protons in the case of @xmath9@xmath17ca system , which is not so far from our results for the annihilation in the ground - state nucleus at @xmath175 . we want to emphasize that the kinetic energy spectra of nucleons emitted after the @xmath9-annihilation from the compressed @xmath9@xmath92o , and @xmath9@xmath17ca systems are significantly harder than the spectra of nucleons from the annihilation at @xmath175 . this can be explained by two effects . first , the collective expansion of the outer shell will increase the slope temperature , typically , by several mev ( see table [ tab : fragmenting_source ] and figs .  [ fig : vrad_ca40_xip03 ] , [ fig : vrad_ca40_xip1 ] ) . second , just after the annihilation the nucleon potential at the center of a nucleus grows suddenly by @xmath214 mev . this creates an additional push for the fast nucleons emitted from the nucleus . although the hardening effect is most pronounced for the lightest system @xmath9@xmath92o , it is also quite visible for @xmath9@xmath17ca . for the heaviest system @xmath9@xmath93pb , we observe almost identical high energy tails of the nucleon spectra for the different annihilation times . the reason is that the collective expansion is practically absent in this system ( see fig . [ fig : rho_vs_t ] ) . also , the yield of fast nucleons is reduced by their subsequent rescatterings in the residual nucleus . the meson production from the @xmath9-annihilation inside nucleus is influenced both by the mean field via the potentials of the annihilating pair and by the final state interactions ( fsi ) , i.e. the two - body collisions and resonance decays . it is instructive to disentangle the contributions of the mean field effects from the rescattering and absorption effects . to this aim , we have performed additional calculations by subsequently switching off the fsi and the mean field . corresponding results are shown in figs . [ fig : pimul_pbaro16],[fig : pimom_pbaro16 ] and [ fig : minv_spectra_pbaro16 ] . figure  [ fig : pimul_pbaro16 ] shows pion multiplicity distributions for the @xmath9@xmath92o system . in the case of reduced antiproton coupling constants ( @xmath40 ) , the mean field alone does not produce any noticeable modification of the pion multiplicity distribution . on the other hand , fsi leads to a rather substantial shift toward smaller @xmath215 . for the case of @xmath159 , we observe a strong pion multiplicity reduction due to smaller @xmath115 for the @xmath116 annihilation , while the fsi effects are much weaker . one can also see a quite significant compressional effect for @xmath159 , which is only very weak for @xmath40 .    in fig . [ fig : pimom_pbaro16 ] , we present charged pion momentum spectra in the c.m . frame of the @xmath9@xmath92o system . fsi strongly modifies these spectra , mostly due to the @xmath152 processes , which effectively decelerate pions . as we have already observed earlier in fig . [ fig : ekin_spectra ] , emitted nucleons gain energy , correspondingly . the effect of the baryonic mean field on the pion momentum spectrum is relatively moderate : we observe some depletion of the high - momentum tail , which is more pronounced in the case of @xmath159 . the compressional effect is visible in the reduction of the pion yield in the momentum range @xmath216 gev / c . this can be understood from figs . [ fig : pimom_ppbar ] and [ fig : pimul_pbaro16 ] : in vacuum , the events with @xmath217 and @xmath218 contribute substantially to this momentum range , while the probability of such events is substantially suppressed in a compressed nucleus due to the reduced annihilation phase space @xcite . the pion momentum spectrum has clearly the two components : the slow pions , which have undergone rescatterings via the @xmath30-resonance excitation ( @xmath219 less than about 300 mev / c ) , and the high energy pions , which were emitted from the system almost without secondary interactions . the similar result has been obtained in earlier intranuclear cascade calculations @xcite . following ref . @xcite , we have also fitted the low energy part of the pion spectrum ( @xmath220 mev , @xmath221 ) by the maxwell - boltzmann distribution @xmath222 this fit has produced the following slope temperatures of the charged pion momentum spectrum for the @xmath9@xmath92o system : @xmath223 mev and 44 mev for @xmath40 ( @xmath224 mev and 36 mev for @xmath159 ) in the case of early annihilation ( @xmath175 ) and annihilation at the time of maximum compression , respectively . the extracted @xmath225 is smaller than the slope temperature of @xmath226 mev of the low energy pions for the @xmath9@xmath17ca system reported in ref . however , in the calculation without mean field , we obtain @xmath227 mev , which is in a good agreement with the result of ref . therefore , the difference between our results and those of ref . @xcite is caused by the mean field acting on the annihilating pair in medium . moreover , we see the softening of the pion spectrum in the case of larger antinucleon couplings due to smaller @xmath115 of the annihilating pair . figure  [ fig : minv_spectra_pbaro16 ] presents the distributions of the annihilation events in the total invariant mass of produced mesons from the @xmath9@xmath92o system . the invariant mass is defined as @xmath228 where @xmath229 is the sum of four - momenta of the mesons produced in a given annihilation event . the calculations were done for the case of @xmath175 .    in the absence of fsi , @xmath230 should be equal to the invariant energy @xmath115 of the annihilating @xmath9n pair . indeed , without fsi and without mean field , as expected , we get a quite sharp peak at @xmath130 only slightly smeared out due to the fermi motion of nucleons and momentum spread of the initial antiproton df ( [ cohstate ] ) . the baryonic mean field leads to the shift of a peak position toward smaller @xmath230 and to some broadening of the distribution . the broadening is due to the spatial spread of the initial antiproton df ( [ cohstate ] ) , which results in different mean - field potentials acting on different annihilating @xmath9n test particle pairs . additionally , the fsi leads to a very strong broadening of the invariant mass spectrum due to the deceleration and absorption of the annihilation mesons . this is clearly seen in fig . [ fig : minv_spectra_pbaro16 ] for the case when the rmf was switched off ( dotted line ) . nevertheless , the full calculation ( solid line ) shows quite strong softening of the @xmath230 distribution due to the mean - field effects . obviously , this effect is stronger for the case of @xmath159 as compared with @xmath40 . finally , in fig . [ fig : minv_spectra ] , we systematically study how the meson invariant mass spectra are affected by the nuclear compression effects . the results are shown for the different @xmath9-nucleus systems . due to the strong reduction of the nucleon effective mass with the scalar density , the meson invariant mass spectra become softer when the annihilation happens in the compressed configurations , as compared with the annihilation in the normal state at @xmath175 . the effect is , again , more pronounced for the light systems @xmath9@xmath92o and @xmath9@xmath17ca . in the @xmath9@xmath92o system the shift is almost 500 mev even in the case of @xmath40 . we believe that results presented in figs .  [ fig : ekin_spectra],[fig : pimul_pbaro16],[fig : pimom_pbaro16 ] and [ fig : minv_spectra ] constitute the set of observables sensitive to the compressional effects in nuclei induced by an antiproton . we have performed dynamical modeling of possible compression effects in nuclei due to the presence of an antiproton . the semiclassical transport gibuu model @xcite incorporating the relativistic mean fields for the baryons and antibaryons has been employed in calculations . the model reproduces reasonably well the earlier static calculations of bound @xmath9-nuclear systems @xcite .    in this work , we did not consider the stopping process of an incident antiproton in a target nucleus . this is a rather complicated problem due to the unknown in - medium cross sections of the @xmath9-scattering and annihilation . this problem will be addressed in a forthcoming paper . instead , we have assumed , that the antiproton has penetrated to the center of the nucleus , stopped there due to an inelastic collision , and then get captured to the lowest energy state . such events should be very rare , with a probability of the order of @xmath231 for the central collisions @xcite . as proposed in ref . @xcite , formation of bound antiproton - nucleus states can be triggered by the emission of fast nucleons , pions or kaons . we have shown , that during the time interval of @xmath0 fm / c after creation of the initial state the central density of the target nucleus grows up to the values of @xmath232 depending on somewhat uncertain values of the antiproton coupling constants . we expect that the life time of strongly bound antiprotons can be long enough to observe this cold compression effect . detailed kinetic simulations of the post - annihilation evolution of residual nuclei have been carried out at different assumptions on the annihilation time . it is shown , that the @xmath9-annihilation in compressed light systems , like @xmath9@xmath92o and @xmath9@xmath17ca , leads to the pronounced collective expansion of the residual nucleus , which may result in the multifragment break - up . another clear signature of the nuclear compression is the hardening of the kinetic energy spectra of emitted nucleons . on the other hand , the invariant mass distribution of produced mesons is shifted to smaller invariant masses due to the in - medium reduction of the nucleon effective mass at high scalar density . similar phenomena are expected also for the case of @xmath233-nucleus bound states which can be produced via the @xmath234 reaction on nuclei . another interesting possibility is that the compressed zone of the nucleus might undergo a deconfinement phase transition . then one can expect formation of a quark - antiquark droplet with a non - zero baryon number and relatively low temperature @xcite . our main assumption in the present study was that the annihilation takes place in the central region of a nucleus . the experimental selection of the central annihilation events is a difficult problem . no clear trigger condition for such events has been invented so far . one suggestion is that the central annihilation events , in average , will be characterised by isotropic emission of secondary particles and high fragment multiplicity @xcite . however , further theoretical and experimental efforts are needed to develop a good trigger condition for the central annihilation . despite of these difficulties , we propose to study the above predictions in antiproton - nucleus reactions at the future facility for antiproton and ion research ( fair ) at gsi ( darmstadt ) . the support by the frankfurt center for scientific computing is gratefully aknowledged . the authors are grateful to a.s . botvina , th.j . brvenich , i.a . pshenichnov , j. ritman and h. stcker for stimulating discussions . this work has been partially supported by the dfg grant 436 rus 113/711/0 - 2 ( germany ) and the grant ns-3004.2008.2 ( russia ) . .[tab : fragmenting_source ] fragmenting source parameters for the different annihilating systems and values of the scaling factor @xmath174 of the antiproton coupling constant . @xmath114 denotes the annihilation time moment . @xmath236 is the time moment when the source parameters have been determined . @xmath237 is the central nucleon density . @xmath191 , @xmath192 and @xmath193 are the charge number , collective kinetic energy per nucleon and residual excitation energy per nucleon , respectively ( see eqs . ( [ nz]),([e_kin^coll ] ) and ( [ e_res^coll ] ) ) . [ cols="^,^,^,^,^,^,^,^,^ " , ]
Abstract:

Answer: on the basis of the kinetic equation with selfconsistent relativistic mean fields acting on baryons and antibaryons , we study dynamical response of the nucleus to an antiproton implanted in its interior . by solving numerically the time - dependent vlasov equation , we show that the compressed state is formed on a rather short time scale of about @xmath0 fm / c . this justifies the assumption , that the antiproton annihilation may happen in the compressed nuclear environment . the evolution of the nucleus after antiproton annihilation is described by the same kinetic equation including collision terms . we show , that nucleon kinetic energy spectra and the total invariant mass distributions of produced mesons are quite sensitive observables to the antiproton annihilation in the compressed nucleus .