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: when holes are introduced into the copper oxide planes of high t@xmath3 cuprates , spin and charge correlations change dramatically . the _ local _ magnetization @xmath0 , measured by @xmath4sr @xcite on e.g. @xmath5 , reveals a qualitative difference between the insulating and superconducting phases : @xmath0 is rather insensitive to doping in the poorly conducting regime @xmath6 , but drops precipitously above the onset of superconductivity at @xmath7 , becoming undetectable at optimal doping @xmath8 . theoretically , holes can cause _ dilution _ and _ frustration _ in the heisenberg antiferromagnet , which create spin textures : either random ( `` spin glass '' ) or with ordering wavevector away from @xmath9 ( sometimes called `` stripes'')@xcite . however , the apparent reduction of local magnetization by the onset of superonductivity , is a novel and poorly understood effect . theory must go beyond purely magnetic models , and involve the superconducting degrees of freedom . we find that this problem is amenable to a variational approach , using hole - doped resonating valence bonds ( rvb ) states . the rvb states were originally suggested by anderson to describe the spin and charge correlations in the high @xmath10 cuprates@xcite . they are excellent trial wave functions for the doped mott insulators , with large hubbard repulsion @xmath11 since:(i ) configurations with doubly occupied sites are excluded.(ii ) marshall s sign criterion for the magnetic energy@xcite is satisfied , and heisenberg ground state energy and antiferromagnetism at zero doping is accurately recovered@xcite . the hole - doped rvb state is a new class of variational states , in which spin and charge correlations are parameterized independently , without explicit spin nor gauge symmetry breaking . such parameterization allows states with magnetic and independently d or s - wave superconducting ( off - diagonal ) order or disorder , thus permit an _ unbiased _ determination of ground state spin and charge correlations appropriate for the cuprates . these are important advantages over commonly used spin density wave , hartree - fock and bcs wavefunctions . a comprehensive study of the state is performed using monte carlo and mean field calculations . phenomenological low energy effective hamiltonian is proposed , with two major components : heisenberg interaction for spins , and single or cooper pair hopping kinetic energy for fermion holes .    regarding this model our key results are:(i ) for the magnetic energy alone , the local magnetization @xmath0 is _ weakly dependent _ on doping concentration . this holds independently of inter - hole correlations for either randomly localized or extended states.(ii ) in contrast to ( i ) , @xmath0 is strongly reduced by the kinetic energy of _ cooper pair hopping _ , which correlates the reduction of @xmath0 with the rise of superconducting stiffness , and hence@xcite the transition temperature @xmath10 . these results agree with the experimentally reported correlation between @xmath0 and @xmath10@xcite . this relation appears to be _ weakly _ dependent on the precise hole density . we also find that rvb states have the following properties:(i ) the magnetic energy is correlated with the average loop density:@xmath12 , where @xmath13 is the linear size off the lattice.(ii ) the gutzwiller mean field approximation ( ga ) for magnetic correlations is in good agreement with the monte carlo results.(iii ) long range magnetic correlations in rvb states are extremely sensitive to changes in the singlet bond amplitude @xmath14 . for example with @xmath15 the spin - spin correlation function decays exponentially with correlation length @xmath16 , where @xmath17 is the hole concentration . the paper is organized as followed : sec.[srv ] introduces the hole - doped rvb state , and discusses the numerical procedure . sec.[uuuu ] defines our variational parameters . sec.[orpar ] deals with the antiferromagnetic and superconducting order parameters . sec.[efh ] deals with the components of the effective hamiltonian . sec.[tcm0 ] correlates between superconducting t@xmath3 and local magnetization . sec.[diss ] is a summary and discussion . the paper has 3 appendices . [ ferpart ] reduces the hole part of the doped rvb to a numerically convenient format . [ esr ] derives expressions for expectation values . particularly , an alternative procedure to calculate magnetic correlation is derived and used to check the computer program . in app.[largen ] , the ga is performed analytically . a valence bond ( vb ) state is @xmath18 where @xmath19 is a pair covering of the lattice , @xmath20 are schwinger bosons , and @xmath21 is a site index on a square lattice . rvb states are superposition of vb states . we restrict the disscussion to @xmath22\rangle=\sum_{\alpha}\prod_{(i , j)\in \alpha}u_{ij }   ( a^{\dagger}_i b^{\dagger}_j - b^{\dagger}_ia^{\dagger}_j)|0\rangle . \label{rvb1}\ ] ] where @xmath23 is a variational singlet bond amplitude , which connects sites of different sublattices @xmath24 and @xmath25 only . this ensures marshall s sign @xcite . the hole doped rvb state is defined by : @xmath26\rangle   & = &     { \cal p}_g(x ) |\bar{\psi}[u , v]\rangle \nonumber\\    |\bar{\psi } [ u , v]\rangle & \equiv & \exp\left [ \sum_{i\in a , j\in b}\left ( v_{ij } f^\dagger_i f^\dagger_j + u_{ij } ( a^{\dagger}_i b^{\dagger}_j - b^{\dagger}_ia^{\dagger}_j)\right)\right]|0\rangle \label{states2}\end{aligned}\ ] ] where @xmath27 are spinless hole fermions , @xmath28 , and @xmath29 is an _ independent _ hole bond parameter . the gutzwiller projector @xmath30 imposes two constraints . a constraint of no double occupancy : @xmath31 and a global constraint on the total number of holes : @xmath32 due to @xmath30 , @xmath33 can be written as a sum over bond configurations of singlets and hole pairs which cover the lattice as depicted in fig . [ rvbfig ] . an overlap of two vb states , @xmath34 , is expressed in terms of a directed loop covering of the lattice ( dlc ) @xcite , hence : @xmath35|\psi[u]\rangle = \sum_{\lambda}\omega_{\lambda }   \label{n0}\ ] ] where @xmath36 is a dlc , @xmath37 and @xmath38 is a directed loop .    with the results of app.[ferpart ] , the norm of the doped rvb state is @xmath39 | \psi [ u , v;x]\rangle= \sum_{\gamma , \lambda(\gamma ) }   w\left ( \gamma , \lambda(\gamma)\right ) , \label{exnorm}\ ] ] where @xmath40 is a _ distinct _ configuration of @xmath41 holes sites : @xmath42 @xmath43 is a dlc which coveres the lattice but the hole sites , @xmath44 and @xmath45 is an @xmath46 matrix with @xmath47    expectation value of an operator @xmath48 is expressed as a weighted sum @xmath49 where @xmath50 is defined by eqs . ( [ wwf ] ) and ( [ ooooo ] ) . we use standard metropolis algorithm@xcite for the evaluation of sum ( [ ooooo ] ) . the basic monte carlo step for updating the dlcs is the one used by ref . @xcite : choose at random a site and one of its next - nearest neighbors and exchange , with transition probability that satisfy detailed balance , the bonds connecting each of them , either to the next site ( forward - bond ) , or the previous site in their loops . in ref . @xcite we show , that for @xmath51 , these steps are _ ergodic _ , that is , any dlc can be reached from any other by a sequence of monte carlo steps . for the fermion holes our update scheme is a simple generalization of the `` inverse - update '' algorithm of ceperley , chester and kalos @xcite . according to eq . ( [ calv ] ) , changing a position of an @xmath24 ( @xmath25 ) sublattice hole amounts to changing one row ( column ) in the matrix @xmath45 . in the our calculation boundary conditions are periodic .    for dimer doped rvb state , where @xmath52 we obtained exact results using transfer matrix technique . for a @xmath53 undoped lattice @xcite the magnetic energy is @xmath54 . the monte carlo result is @xmath55 . @xmath56 is exponentially decaying with correlation length of @xmath57 , the monte carlo result is @xmath58 . exact and monte carlo results for the doped @xmath59 ladder@xcite appear in fig.([evmc ] ) . our program successfully reproduced existing data for rvb states @xcite . other tests of the program appear below . we also use the gutzwiller approximation ( ga ) to evaluate expectation values in the doped rvb state . the ga is discussed in app.[largen ] . in the undoped , @xmath60 , case we treat three classes of singlet bond amplitude @xmath14 . @xmath61 with @xmath62 and @xmath63 where for @xmath64 ( @xmath65 ) @xmath66 and @xmath67 . @xmath68 determines the short range decay of @xmath64 and @xmath65@xcite . we also use @xmath69 . @xmath70 is derived from the schwinger - boson mean field theory of the heisenberg model @xcite . for @xmath71 we use @xmath72 , eq . ( [ uppl ] ) , and @xmath64 , eq . ( [ uex ] ) .    for the function @xmath73 the following cases of inter - hole correlations are treated : @xmath74 where @xmath75 , @xmath76 are nearest neighbor vectors on the square lattice , @xmath77 and @xmath78 . @xmath79 puts the @xmath41 holes on random sites . this state describes an insulator with disordered localized charges .    @xmath80 describes weakly interacting holes in a `` metallic '' state : @xmath81 where the product is over @xmath41 states , @xmath82 and @xmath83 . here we check @xmath84 which is centered at @xmath85 . results for @xmath84 centered at @xmath86 are _ not _ qualitatively different @xcite . @xmath80 obey @xmath87 hence @xmath88 only connects @xmath89 to @xmath90 . correlations in a state with @xmath91 were previously computed by bonesteel and wilkins@xcite . @xmath92 and @xmath93 describe tightly bound hole pairs in relative @xmath94 and @xmath95-wave symmetry respectively . the local magnetization is @xmath96 where e.g. @xmath97 . with respect to eq . ( [ ooooo ] ) , @xmath98 is calculated using@xcite @xmath99 where the sign is + if @xmath100 and @xmath101 are on the same sublattice . to check our program we also used an alternative procedure to calculate magnetic cerrelations . see app.[c222 ] .    in fig . [ sppvsp](a ) @xmath102 is plotted for @xmath103 $ ] and @xmath104 $ ] for various choices of @xmath73 . finite size scaling in fig.[sppvsp](b ) for @xmath105 indicates vanishing long range order , @xmath106 , at @xmath107 . it lowers the bound given previously by ref . @xcite : at @xmath108 . in fig . [ sppvsp](c ) @xmath109 is plotted for @xmath110 $ ] and @xmath111 $ ] . finite size scaling in fig.[sppvsp](d ) indicates @xmath106 , at @xmath112 . in all the cases the ga ( lines ) works well . good agreement between ga and monte carlo is also seen in fig.[sofr](a ) and ( b ) , where @xmath113 is plotted for @xmath114 , and @xmath64 respectively . note how slow @xmath113 decays for @xmath115 . by fig.[sppvsp](b ) @xmath116 in this state . exponentially decaying spin correlations are seen , both by monte carlo and ga , for @xmath114 with @xmath117 and @xmath64 with @xmath118 @xcite . details of @xmath68 , eq.([usr ] ) , have strong effect on long range spin correlations . we use ga to extrapolate monte carlo calculations for @xmath113 . in app.[exuyy ] we find for exponential bond amplitude , @xmath15 and @xmath119 , that @xmath113 decays exponentially with correlation length @xmath120 for gaussian bond amplitude , @xmath121 with @xmath122 , we find in app.[andel ] that @xmath113 decays exponentially with correlation length @xmath123 for @xmath114 , app.[pluan ] suggests vanishing long range order , @xmath106 , at @xmath124 . correlation lengths ( [ xiex ] ) and ( [ xigus ] ) explain the slow decay of @xmath113 in fig.[sofr](b ) . it also indicate , that in the @xmath125 system , a small change in the ground state parameters brings an extremely sharp change in long range magnetic correlations . the superconducting singlet order parameters are @xmath126 where @xmath127 the expressions of @xmath128 s matrix elements are discussed in app.[resddd ] . by gauge invariance imposed by the gutzwiller projector , @xmath129 . however , @xmath130 $ ] describes true @xmath94 ( @xmath95)-wave superconductors as seen by the singlet pair correlation function @xmath131 , @xmath132 , in fig.[odlro ] . for @xmath133 , @xmath134 and @xmath33 has ( off - diagonal ) long range order in @xmath135 . in contrast , the insulator states @xmath136 $ ] and the `` metallic '' states @xmath137 $ ] have no long range superconducting order of either symmetry . magnetic order is driven by the diluted antiferromagnetic quantum heisenberg model @xmath138    _ magnetic energy for @xmath60 _ : in fig . [ gammas](a ) @xmath139 , @xmath140 and @xmath141 are plotted as a function of @xmath142 , @xmath143 and @xmath144 for @xmath64 , @xmath114 and @xmath65 , eqs . ( [ uppl ] ) , ( [ uex ] ) and ( [ ug ] ) respectively . in @xmath60 , _ all _ the three bond amplitude yield lowest magnetic energy of @xmath145 for @xmath114 , the optimal value of @xmath146 is @xmath147 , and @xmath148 . the ground state parameters of the heisenberg model on an @xmath149 lattice are : e(ground state)= 0.3347 j / bond and @xmath150(ground state)=0.109@xcite . table [ paraatmin ] contains a summary of results for the optimal choice of parameters in all the classes .    _ magnetic energy for @xmath105 _ : in fig.[hemag ] @xmath139 and @xmath140 are plotted as a function of @xmath142 and @xmath143 , for @xmath105 and various choices of @xmath73 from eq . ( [ uv ] ) . within numerical errors , all states minimize @xmath151 at the same optimal parameters as for @xmath60 ( table [ paraatmin ] ) . for @xmath114 , by fig . [ sppvsp](a ) it yields local magnetization of @xmath152 . for @xmath64 , fig . [ sppvsp](c ) shows @xmath153 . thus we conclude that aside from the trivial kinematical constraints , _ the hole density and correlations have little effect on the magnetic energy at low doping . _ a better understanding of the properties of the optimal bond amplitude for @xmath154 is gained by the _ average loop density _ defined below . from eq . ( [ csisj ] ) , a dlc contributes to @xmath150 , eq . ( [ sp ] ) , it s number of pairs of sites , which share the same loop hence @xmath155 where @xmath156 is the loop length , and @xmath157 . thus @xmath158 is proportional to the average loop length per site . the average radius of gyration of a loop is : @xmath159 where @xmath160 with @xmath161 , and @xmath162 is the number of loops in the dlc @xmath36 . with eqs . ( [ allps ] ) and ( [ rgg ] ) we define the average density of a loop per site @xmath163    the average loop density , @xmath164 , is plotted in fig.[gammas](b ) , in the undoped case for all the bond amplitudes ( [ uppl ] ) , ( [ uex ] ) , and ( [ ug ] ) ; and in the doped case for @xmath165 $ ] . comparison with fig.[gammas](a ) shows that @xmath164 is correlated with the magnetic energy . for vanishing @xmath0 , @xmath164 converge to its value in the dimer rvb state , eq . ( [ di ] ) , where @xmath166 . this value of @xmath164 is only slightly larger than @xmath164 s value for an ensemble of dlcs , which include only configurations with two ( or four ) sites loops with dimer bonds . for such loops @xmath167 ( or @xmath168 ) and @xmath169 . the occurrence of loop lengths ( @xmath170 ) is interesting . in fig.[histogram ] we plot an histogram of the number of loops ( @xmath171 ) , versus the number of sites on a loop ( @xmath170 ) . the size of the lattice is @xmath172 , and @xmath69 , which is derive from the schwinger bosons mean field theory of @xmath154 @xcite . for all the bond amplitudes and lattice sizes we have checked @xmath173 decays either algebraically or exponentially . a single hole hopping in the antiferromagnetic background has been shown by semiclassical arguments@xcite , to be effectively restricted at low energies to hopping between sites on the same sublattice : @xmath174 where @xmath175 are removed by two adjacent lattice steps , and @xmath176 . unconstrained , the single hole ground state of @xmath177 has momentum on the edge of the magnetic brillouin zone , in agreement with exact diagonalization of @xmath178 clusters@xcite . previous investigations have found that _ inter_-sublattice hopping ( the @xmath179-term in the @xmath178 model ) , is a high energy processes in the afm correlated state@xcite . we thus expect the same to hold even in rvb spin liquids with strong short range afm correlations but no long range order . the primary effects at low doping may be to shift the ordering wavevector . we denote by @xmath180 ( @xmath181 ) the coefficients of second ( third ) nearest neighbor hopping terms . for @xmath182 the single hole bend minimum is at @xmath183 , otherwise @xmath184 . here we put @xmath185 , @xmath186 . results for the expectation value of @xmath177 are plotted in fig.[ehop ] . the single holes hopping , eq . ( [ tprime ] ) , prefers the metallic states @xmath91 over states with @xmath187@xcite . it also prefers longer range @xmath188 and thus actually _ enhances _ magnetic order at low doping . this is a type of a nagaoka effect , where mobile holes separately polarize each of the sublattices ferromagnetically . we consider cooper pairs hopping terms @xmath189 calculation of @xmath190 matrix elements is discussed in app.[resddd ] . the first term in @xmath190 is derived from the large @xmath11 hubbard model to order @xmath191@xcite . it includes terms ( a ) and ( b ) in fig.[duforms ] . term ( a ) is a _ rotation _ of the singlet pair . it is positive for @xmath92 @xcite and hence prefers @xmath93 over @xmath92 . term ( c ) in fig.[duforms ] is a parallel _ translation _ of singlets . it prefers superconductivity with @xmath192 or @xmath193 over metallic states with @xmath91 @xcite . for @xmath105 , @xmath190 is minimize by @xmath93 .    in fig . [ eph ] the ground state energy @xmath194 of ( [ jprime ] ) is plotted for @xmath195 , @xmath196 and @xmath197 . @xmath17=0.1 , and the size of the lattice is @xmath13=40 . the variational energy is minimize at @xmath198 and @xmath199 , for @xmath64 and @xmath114 , respectively . in both cases , by the finite size scaling of fig.[sppvsp](b ) and ( d ) , it indicates vanishing @xmath0 at @xmath200 . _ thus , cooper pair hopping drives the groundstate toward a spin liquid phase ! _    the gutzwiller approximation fails to predict this effect . according to the ga , the minimum of the double hopping energy roughly coincides with the minimum of the magnetic energy ( @xmath201 ) . this is understood by ( see app . [ largen ] ) : @xmath202 where @xmath89 , @xmath90 . the ga agrees with monte carlo results for matrix elements of _ long range _ pair hopping . the matrix element of ( d ) in fig.[duforms ] , and @xmath203 also drives the groundstate toward a spin liquid , and prefer superconducting over metallic states@xcite . these terms are excluded due to relatively large thermal noise . since @xmath190 is the effective model which drives superconductivity it produces phase stiffness , which in the continuum approximation is given by @xmath204 the stiffness constant @xmath205 can be determined variationally from the doped rvb states . imposing a uniform gauge field twist on @xmath128 , @xmath206 , @xmath207 becomes , to second order in @xmath208 , @xmath209 following ref . @xcite , at low doping for the square lattice @xmath205 is roughly equal to @xmath10 .    in fig.([jproverj ] ) we show our main result : the staggered magnetization @xmath210 for @xmath211 is plotted against the superconducting to magnetic stiffness ratio @xmath212 for different doping concentrations @xmath213 , @xmath192 , and @xmath196 . the actual free parameter in the graph is @xmath214 , from which @xmath0 and @xmath205 are determined variationally . two primary observations are made : ( i ) the local magnetization is sharply reduced at relatively low superconducting stiffness ( and @xmath215 ) . ( ii ) the relation between @xmath0 and @xmath216 appears to be independent of @xmath17 .    for @xmath64 , eq . ( [ uex ] ) , it requires @xmath217 for @xmath218 to be minimized at @xmath219 . by fig.[sppvsp](b ) this leads to @xmath220 . in this paper we used extensive monte carlo calculations to study properties of hole doped rvb states . we found that an effective model which include heisenberg and pair hopping terms is consistent with the experimental connection between superconductivity and reduction of local magnetic moment . within checked variational options we showed that the properties of the model are independent of particular choice of parameters for the state . gutzwiller mean field approximation for magnetic correlations was found to agree with monte carlo calculation , and used for analytical extrapolation of numerical results . we showed that long range magnetic correlations in rvb states are extremely sensitive to variational parameters . we found that the average loop density is well correlated with the magnetic energy . we conclude this paper in several arguments and insights regarding our results .    _ magnetic energy and long range magnetic correlation : _ note the contrast between correlation lengths ( [ xiex ] ) and ( [ xigus ] ) , and the `` shallowness '' of the minima of the magnetic energy in fig.[hemag ] . it imply that a very weak pair hopping term in the hamiltonian causes a dramatic change in long range magnetic correlations .    _ magnetic energy and loop density : _ a comparison between loops ( a ) and ( b ) in fig.[tkol ] shows that large amplitude ( @xmath221 ) of dlcs with `` denser '' loops enhance the probability to find nearest neighbor sites on the same loop and reduce the magnetic energy . the loop density shows that the optimal bond amplitude is determined by an intricate balance between @xmath0 and @xmath222 . this relates _ quantum spin fluctuations _ to the average loop density of the ensemble .    _ effective model for doped system : _ @xmath223 describes the low energy physics of the lightly doped cuprates . as the lattice is doped , its variational ground state is a d - wave superconductor , with a sharply reduced local magnetic moment . the model includes built - in pairing . such a model is supported by the existence of a pseudogap in the normal state of the high t@xmath224 materials . _ relation between phase stiffness and local magnetization : _ because of finite size uncertainty , @xmath225 in fig.([jproverj ] ) is an _ upper _ bound on the thermodynamic local magnetization . a sharper reduction of the local magnetization occurs if:(a ) the ga result of app.[pluan ] , @xmath226 , is correct to the discrete lattice . in that case @xmath0 vanishs already at @xmath227.(b ) in finite doping the optimal bond amplitude for @xmath223 decays exponentially . in that case @xmath0 vanishes for @xmath228 . variationally , we can not rule out this possibility.in both of these cases there is a qualitative agreement with the doping dependence of the local magnetization and t@xmath3 , as measured by refs . @xcite .    useful conversations with c. henley , s. kivelson and s - c . zhang , are gratefully acknowledged . mh thanks taub computing center for support . aa is supported by the israel science foundation and the fund for promotion of research at technion . the fermion part of @xmath229\rangle$ ] is @xmath230 } |0\rangle = \frac{1}{\frac{{n_h}}{2 } ! } \left[\sum_{i\in a , j\in b } v_{ij } f^{\dagger}_i f^{\dagger}_j\right]^{\frac{{n_h}}{2 } } |0\rangle   \label{fs1}\end{aligned}\ ] ] where @xmath231 . we write this state as @xmath232 where @xmath40 is a _ distinct _ configurations of @xmath41 holes sites : @xmath233 from eq.([fs1 ] ) @xmath234 where @xmath235 . for each @xmath40 we commute _ pairs _ of operators , without any affect of sign , and order @xmath24 holes operators in an increasing order of their site index @xmath236 where @xmath237 . commuting the @xmath25 operators we get @xmath238 where @xmath45 is an @xmath46 matrix with @xmath239 when @xmath73 connects same - sublattice sites the determinant is replaced by a pffian@xcite . we use the operator identity @xcite @xmath240(1-f^{\dagger}_i f_i)(1-f^{\dagger}_j f_j ) , \label{ida}\end{aligned}\ ] ] where @xmath241 , to express @xmath242 .    show that for @xmath243 and @xmath244 @xmath245 @xmath246 is the set of forward bonds in @xmath43 , of the sites on sub - lattice @xmath24 . we demonstrate eq . ( [ csisj222 ] ) for half filled lattice ( @xmath247 ) . with @xmath248 , @xmath249 @xmath250 and hence @xmath251\rangle =   \sum_{\alpha_{ij } }   \left[\left(\sum\limits _ { \scriptstyle ( l , n)\in \alpha_{ij } \atop   \scriptstyle l\neq i } \frac{u_{lj}u_{in}}{u_{ln}u_{ij}}\right ) + 2\right ] \left(\prod_{(l , n)\in \alpha_{ij}}u_{ln}\right )    @xmath252 is a valence bond state , with @xmath253 . the term in the square brackets requires further explanation . from eq . ( [ z2 ] ) , for any pair @xmath254 , @xmath255 , where @xmath256 , @xmath257 , and otherwise @xmath258 , see fig.([tconf ] ) . in @xmath259\rangle$ ] , each @xmath260 carries a factor @xmath261 . ( [ z1 ] ) indicates an additional option to get @xmath262 , from @xmath262 .    taking the overlap with @xmath263|$ ] , we get the matrix element which is expressed in terms of eq.([csisj222 ] ) . @xmath262 represents the ket . a possible definition of the bonds of the ket is the forward bonds of the sites on sublattice @xmath24 . fig.[phopp ] describes the effect of a single hole hoping term on a hole - pair configuration . using definition ( [ ooooo ] ) , for @xmath264 @xmath265 where @xmath266 , @xmath267 , and otherwise @xmath268 ; @xmath269 is the forward bond of @xmath270 ( i.e. originated in the ket ) ; and @xmath271 comes from reordering the fermion operators . relation ( [ hme ] ) is simplified using@xcite @xmath272 with @xmath273      for @xmath274 we express @xmath275 , where @xmath128 is given in eq . ( [ dalte ] ) . @xmath276 creates a singlet bond . @xmath277 creates a pair of holes . with the results of app.[c222 ] , for @xmath278 @xmath279 where @xmath280 if @xmath281 exclusive - or @xmath282 and 1 otherwise , @xmath283 is defined like eq . ( [ ahtv ] ) , with a possible replacement of a row _ and _ a column , and @xmath284 is the set of forward bonds of a sub - lattice sites in @xmath285 . the gutzwiller approximation amounts to dropping the projector @xmath286 in definition ( [ states2 ] ) and setting @xmath287\rangle \rightarrow |\bar{\psi}[yu , zv]\rangle=|yu\rangle \otimes |zv\rangle$ ] . the constants @xmath288 and @xmath289 are determined by global constraint equations @xmath290 for @xmath291 , @xmath292 respectively . in this section @xmath293 . @xmath294 is a schwinger bosons mean field wave function@xcite , on which we preform the marshall transformation@xcite : @xmath295 . hence @xmath296 operators are transformed accordingly , for example @xmath297 for @xmath298 .    from eqs . ( [ mfeq ] ) and ( [ ffmfeq ] ) @xmath299 whereas @xmath300 . thus we use @xmath301 as the ga for the long range magnetic correlations and @xmath0 in the doped rvb state . _ empirically _ we omit the @xmath302 factor in the estimates of magnetic energy . we calculate spin - spin correlation function for [ exuyy ] @xmath312 with @xmath313 . @xmath314 where mbz= magnetic brillouin zone . eq.([kcoeq ] ) becomes @xmath315 where we multiplied the left side in @xmath316 to account for the integration over the complete brillouin zone . in all our calculations we took the continuum limit ( lattice constant@xmath3170 ) , where the upper bound of the integration @xmath318 this approximation works very well for slow decaying bond amplitude @xcite . eq.([cex ] ) gives @xcite @xmath319 }    \right . \nonumber \\    & + & \left . \ln{\left [ 4 \left(\frac{\pi y } { \xi } \right)^{\frac{4}{3 } } + 2 \left(\frac{2\pi^2 y^2}{\xi^8}\right)^{\frac{1}{3 } } +   \frac{2^{\frac{2}{3 } } } { \xi^4 } \right ] } - 2 \ln{\left [ \frac{2^{\frac{1}{3 } } } { \xi^2 } - 2 \left(\frac{\pi y } { \xi } \right)^{\frac{2}{3}}\right ] }   \right\ } = \frac{1 -   x}{2 } \label{ix}\end{aligned}\ ] ] the argument of the last logarithm has to be sufficiently close to zero for eq.([ix ] ) to be satisfied . therefore @xmath320 hence we can neglect in the left side of eq.([ix ] ) all terms but the last log . consequently @xmath321}\right\ } \equiv \frac{1}{(2\pi)^2\xi^{4}}\cal{d } \label{x2}\ ] ]    eq.([deqd ] ) becomes @xmath322 where @xmath323 is the bessel function@xcite . since the integrand in eq.([sq2 ] ) vanishes as @xmath324 , we can replace @xmath323 with its approximation for @xmath325 . expanding the denominator to first order in @xmath326 @xmath327 \equiv \frac { { \cal d}}{6 \pi \sqrt{r\pi } \xi^{2}}y_0 \label{sq4}\end{aligned}\ ] ] where @xmath328 . in the definition @xmath329 @xmath330 . let us consider the integral @xmath331 where the close contour encircles the upper half of the complex plane . the part of the contour along the negative real axis is @xmath332 where we substituted @xmath333 . hence @xmath334 . using the residue method for @xmath335 : @xmath336}{\sqrt{r } } , \label{sir}\ ] ] and with eq . ( [ x2 ] ) we find for the correlation length of the spin - spin correlation function , @xmath337 : @xmath338}\ ] ]      we calculate spin - spin correlation function for @xmath339 for @xmath340 mbz @xmath341 eq.([kcoeq ] ) is @xmath342 with the solution @xmath343}\right\ } \equiv \frac{\mu^2}{\pi^2}{\cal g } \label{xgxx}\ ] ] calculation of @xmath344 is identical to the exponential case . substituting in eq . ( [ sq4 ] ) @xmath345 @xmath346 and hence the spin - spin correlation function decays exponentially with correlation length @xmath347}\ ] ]        calculations of the ga on lattices of size @xmath352 , show that for any @xmath353 , the spin - spin correlation function calculated with function ( [ plu ] ) , decayes slower than with @xmath354 . this suggests that @xmath351 for @xmath354 .      where we used eq.([deqd ] ) . for @xmath340 mbz @xcite @xmath356 from eqs.([deq ] ) and ( [ sigq ] ) @xmath357 , and @xmath358 . hence @xmath350 might diverges only if there is @xmath359 , such that @xmath360 , where @xmath361 . therefore if @xmath362 , @xmath350 is finite . eq.([kcoeq ] ) for @xmath291 is @xmath363 which becomes @xcite @xmath364 the right side of eq.([sumpl ] ) is @xmath291 independent , and increases with @xmath353 , hence @xmath365 increases with @xmath353 . therefore @xmath366 . for @xmath367 , the left side of eq.([sumpl ] ) is @xmath368 and @xmath369 and for @xmath370 , @xmath350 is finite and hence @xmath351 . .minimal magnetic energy ( @xmath371 ) , optimal choice of parameters , and square of the ground state magnetization ( @xmath372 ) , for various bond amplitudes in the undoped , @xmath60 , case . the size of the lattice is @xmath149 . @xmath70 is derived from the schwinger bosons mean field theory of @xmath154 . the ground state parameters of the heisenberg model on an @xmath149 lattice are : e(ground state)= 0.3347 j / bond and @xmath150(ground state)=0.109@xcite . [ cols="^,^,^,^",options="header " , ]                                          m. inui . s. doniach and m. gabay phys . rev . * b 38 * , 6631 ( 1988 ) ; b. i. shraiman , and e. d. siggia , phys . * 62 * , 1564 ( 1989 ) ; v. cherepanov , i. ya . korenblit , amnon aharony , o. entin - wohlman , cond - mat/9808235 . pair hopping terms with _ weak _ hubbard interactions can drive superconductivity even at half filling , as shown by f. f. assaad , m. imada , and d. j. scalapino , phys . lett . * 77 * , 4592 ( 1996 ) , and s - c . zhang , science * 275 * , 1089 ( 1997 ) . here we treat the large - u system .        for simplicity , we ignore frustration and limit @xmath373 $ ] to have purely antiferromagnetic correlations and uniform charge distributions . although experiments observe small doping dependent shifts in ordering wavevectors ( stripes ) , these effects should be small on the quantum fluctuations of the _ local _ magnetization .
Abstract:

Answer: a comprehensive study of doped rvb states is performed . it reveals a fundamental connection between superconductivity and quantum spin fluctuations in underdoped cuprates : _ cooper pair hopping _ strongly reduces the local magnetization @xmath0 . this effect pertains to recent muon spin rotation measurements in which @xmath0 varies weakly with hole doping in the poorly conducting regime , but drops precipitously above the onset of superconductivity . gutzwiller mean field approximation ( ga ) is found to agree with numerical monte carlo calculation . ga shows for example that for a bond amplitude @xmath1 , spin spin correlations decay exponentially with a correlation length @xmath2 . expectation value of the heisenberg model is found to be correlated with average loop density .


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: vortex instabilities are observed throughout the universe . on the smallest scales in the laboratory , @xcite have experimentally investigated the 3d interaction of a strong shock with a spherical density inhomogeneity via laser - driven experiments . they found a double vortex ring structure and an azimuthal instability that ultimately results in the breakup of the ring . they also performed 3d computational simulations which were shown to be in remarkable agreement with their experimental investigations .    in planetary atmospheres , vortices can be regularly observed in von karman vortex streets often associated with atmospheric flow past stationary objects such as mountains . they have also been suggested as the explanation of the giant storm on jupiter known as the great red spot @xcite . atmospheric simulations have successfully reproduced vortices comparable to the great red spot @xcite and other vortex anomalies in the atmosphere of jupiter @xcite . @xcite have hydrodynamically simulated vortices at the stellar size scale in thin accretion disks . vortices at the interstellar size scale have been predicted since some of the first high resolution calculations by @xcite where kelvin - helmholtz instabilities begin the breakup of interstellar clouds . more recently , @xcite predicted the existence of such vortices forming from the interaction of strong shocks with the medium through which they propagate . a series of laboratory experiements as well as two- and three - dimensional simulations @xcite confirmed these earlier predictions . they modelled the interaction of a supernova shock and a cloud in the interstellar medium ( ism ) as the effect of a high density sphere embedded in a low density medium after the passage of a strong shock wave . they found vortex rings are shed from the sphere and , importantly , these vortex rings break up at high mach number due to three - dimensional bending - mode vortex - ring instabilities . since their simulations are isothermal , they can be scaled to any size , for the given mach number , and so are relevant to clouds in the ism ranging from the smallest resolved structures at @xmath2 pc to diffuse clouds up to 100 pc . vortices can act as a source of local angular momentum and turbulence , both of which are important in the evolution of molecular clouds . for a discussion of the importance of interstellar turbulence , we refer to the two - part review by elmegreen & scalo ( 2004a , b ) . such molecular clouds are themselves important for star formation @xcite . vortices in the ism can also affect anything expanding into this medium e.g. supernovae , novae and planetary nebulae .    in this paper , we present simulations that show under certain conditions von karman - like vortices are produced in the wakes of asymptotic giant branch ( agb ) stars with size scales from @xmath0 pc to @xmath1 pc . specifically , we model the interaction of a stellar wind with the ism where the star has a significant motion through the ism . the stellar wind in our simulation is produced by an agb star @xcite . agb stellar winds have typical speeds of @xmath3 @xmath4  and mass - loss rates @xmath5@xmath6 @xmath7 . this wind ejects up to a few solar masses of material into the surrounding medium over the course of @xmath8 years . in our two - wind model , we have found that the movement of the star through the ism ( modelling the ism as a second wind ) causes the ejected stellar material to form into a bow shock upstream of the star with a tail flowing downstream . it is in this tail that we find evidence of von karman - like vortices . observational evidence for interaction between the ism and agb winds has been found by @xcite and @xcite and a model of the interaction has been used to explain the shape of the structure around r hya as a bow shock @xcite .    in hindsight , previous 2d astrophysical simulations considering this phase of stellar evolution have produced such vortices ( see figure 18 of * ? ? ? * ) , but the authors did not discuss them . we consider the vortices in 3d allowing an investigation and description of their behaviour . a computational fluid dynamics ( cfd ) scheme , cubempi , has been used to hydrodynamically simulate the interaction of a star moving through the ism ejecting a stellar wind during the agb phase of evolution . the numerical scheme is based on a second - order godunov scheme due to @xcite . the riemann solver is due to @xcite . the scheme is eulerian , posed in three dimensions and cartesian coordinates and is fully parallel using the mpi library of subroutines . the scheme also includes the effect of radiative cooling via a parametric fit to the cooling curve of @xcite above @xmath9 k. the scheme has been fully tested using standard cfd tests and its performance in these tests compared to results obtained by @xcite for other cfd schemes . astrophysical tests have also been performed using the sedov and primakoff analytical supernova models . the full details of the scheme and its performance can be found in @xcite . we have used a two - wind model with a numerical domain of @xmath10 cells . this model is the same as that used to successfully model the structure around r hya @xcite . it is also the agb stage of our triple - wind model following the agb and post - agb phases of evolution used to successfully model the planetary nebula ( pn ) sh 2 - 188 @xcite . we have considered four values of the mass - loss rate on the agb : @xmath11 @xmath7 . we have considered the agb phase , and hence the simulation , to last 500000 years with an agb wind velocity of 15 @xmath4 . these wind parameters are typical of previous models in the literature @xcite . an unphysical temperature of @xmath9 k was used for the agb wind but this does not affect the overall result @xcite . in view of the still considerable uncertainties on the detailed properties and evolution of these winds , more detailed temporal variations were not modelled . relative velocities of the star with respect to the ism have been considered from 0 to 200 @xmath4  in steps of 25 @xmath4corresponding to the range of velocities found within the galaxy . with speeds of relative movement above 75 @xmath4 , we have not simulated mass - loss rates of @xmath12 @xmath7  due to cpu constraints . the ism has been modelled with densities n@xmath13 = @xmath14 and a temperature of @xmath15 k. these are the properties of the warm ism or wim , the larger constituent of the ism @xcite . the simulations begin at the start of the agb phase of evolution of the central star and are performed in the frame of reference of the star . the wind from the star drives a shock into the ism . in the case of a stationary star we have found the shock drives a spherical shell of agb wind material sweeping up ism material . inside the reverse shock is a bubble of undisturbed agb wind material . this bubble is of low density and temperature , increasing towards the star , and of high pressure . if the central star is moving through the ism , a bow shock forms upstream of the star . such bow shocks have been simulated in other cases of winds interacting with injected flows @xcite . eventually , the simulations show the bow shock reaches a maximum distance ahead of star which can be understood in terms of a ram pressure balance between the stellar wind and the ism . strong shock theory predicts the temperature of the shocked material at the head of the bow shock : t @xmath16 in general agreement with the simulations . ram - pressure - stripped material from the head of the bow shock forms a tail behind the nebula . as material moves down the tail , it cools and mixes with ism material . figure [ density ] shows two snapshots from the agb evolution of a central star moving at 75 @xmath4 . both panels show the gas state at the end of the agb phase , 500000 years into the simulation . in the left panel , the bow shock driven by the agb wind has reached a stationary position at the point of ram pressure balance , 0.02 pc ahead of the central star , with a smooth tail of shocked material extending downstream of the bow shock . in the right panel , where the mass - loss rate is 50 times higher , the tail contains vortices flowing downstream . figure [ evolution ] shows the evolution of the vortex in the right panel of figure [ density ] . in the top panel , the indentation near the head of bow shock indicates an instability , which in the second panel is shed downstream from the head of the bow shock . the remaining panels show the instability forming into a vortex , spiralling and flowing downstream . this episode of vortex shedding from the head of the bow shock is not the first such event . after 50000 years of evolution , the tail appears to become turbulent and the first episode of vortex shedding begins at 90000 years . the initial episode appears to be symmetric around the bow shock and the instability forms into a vortex ring flowing downstream . the next episode begins at 205000 years and continues with multiple episodes of vortex shedding with no clear periodicity or ring symmetry . the vortices flow out of the simulation domain in 100000 years having achieved a stable structure which suggests they will have an extended lifetime in the wake .    to show the breakdown of the later vortices azimuthally , figure [ 3d ] shows a visualisation of the density datacube 500000 years into the simulation . an azimuthally - unstable partial vortex ring can be seen mid - way down the tail and more partial vortex - rings are in the process of shedding off the bow shock .    in order to investigate the effect of grid resolution on the production of vortices , the same simulation which produces vortices with a movement through the ism of 75 @xmath4was rerun at a grid resolution of @xmath17 and @xmath18 . the low resolution simulation showed a smooth , steady state tail with no clear vortices produced from the head of the bow shock . undulations at the head of the bow shock do appear to move down the tails . we suggest the initial instabilities which lead to vortices downstream are strongly suppressed by this low resolution . at @xmath18 , many more vortices than at @xmath10 seem to be produced from similar instabilities at the head of the bow shock , but importantly the physical size of the vortices is similar and therefore not dependent on resolution . a 2d numerical simulation was also performed in cylindrical polar coordinates for further comparison . similar vortex - like structures form in the tail of agb material and flow downstream . the initial instability appears to be slower in developing and once developed , the instability causes a periodic movement of the tails more like von karman vortices . the vortices themselves appear as density enhancements moving down the tail structure . clearly in 3d simulations the vortices form fully and the periodicity is destroyed by the instability of the bow shock following the launching of the first vortex ring . we postulate that the vortices in the tail are the result of vortex shedding from the head of the bow shock . vortex shedding occurs when a fluid passes by an object and the shear layer near the object creates a velocity gradient . the kelvin - helmholtz instability often accompanies flows where shear is present , making the flow unstable to perturbations . the simulations show that the bow shock does not remain stationary in the frame of reference of the star after it has reached the ram pressure balance point . in fact , its distance ahead of the star is oscillating and it is the kelvin - helmholtz instability which is responsible for this oscillation . as the instabilities spread from the head of the bow shock , peel off and flow downstream to form vortices , the bow shock loses material in the vortices and is driven back towards the star by the pressure of the ism . as the stellar wind supplies more material to the bow shock , it recovers its ram - pressure - balance position and this instability - driven oscillation begins again . note that this shows that instabilities set in only after the bow shock has reached its balance position , and thus vortices only appear in the tail after this time . figure [ sh2 - 188 ] shows the pn sh 2 - 188 in the left panel and a detail of the nebula in the right panel . in a previous paper , @xcite modelled sh 2 - 188 as a nebula shaped by wind - ism interaction where the bright arc is a bow shock ahead of the central star . the structures highlighted in the right panel can be interpreted in our vortex - shedding scenario as instabilities being peeled off the bow - shock which will form vortices downstream . theoretically , vortices form over a certain range of reynolds number . the reynolds number of a flow is defined as the ratio of internal forces to viscous forces and can be formulated as @xmath19 where @xmath20 is the fluid density , @xmath21 is the free - stream fluid velocity , @xmath22 is the characteristic length and @xmath23 is the dynamic fluid viscosity . for a reynolds number of less than 50 , the wake behind an object is thin and the flow can be classed as laminar . with a reynolds number around 100 , features similar to the von karman vortex street develop as we see here . the wake also can become several times wider than the characteristic width of the object . at reynolds numbers of 1000 and greater , the wake behind the object becomes turbulent . it is difficult to estimate the reynolds numbers in our simulations as this involves the viscosity of the fluid and the simulations include an artificial viscosity to suppress the quirk instability @xcite . since this artificial viscosity is constant in all the simulations , it is possible to estimate _ relative _ reynolds numbers between the simulations . the characteristic length in a particular simulation is estimated as the diameter of the object in the flow , defined by twice the ram - pressure - balance distance . note that the point of ram pressure balance is inversely proportional to flow - velocity and thus so is the characteristic length . therefore , there is no dependency on flow - velocity in the estimation of relative reynolds numbers . the relative reynolds number increases with increasing mass - loss rate and increasing ism density . in the simulations , there is a general trend to see vortices when the mass - loss rates or ism densities are not extreme . these simulations correspond to the mid - range of relative reynolds numbers . however , in some cases of extreme mass - loss rate or ism density , the bow shock structure is either so large that it is possible we do not see vortices on the simulation timescale , or so small that it is possible the ism flow confines the bow shock enough to suppress the formation of vortices . thus , it is not clear whether the vortices are dependent on relative reynolds number although it seems there is a range where vortices are more apparent on the timescales of these agb simulations , which we might expect from theoretical indications . all the simulations are supersonic from a mach number of 1.8 at 25 @xmath4  to 14.4 at 200 @xmath4  suggesting no dependency on mach number . the vortex shedding lacks the regular period seen in the case of von karman vortices . we do note though that the first episode of vortex shedding is in the form of a symmetric ring shed downstream from the head of the bow shock . further episodes of vortex shedding do not shed such stable rings and vortices are shed from either side of the bow shock . this can be understood in terms of the regular shape of the bow shock at the time of the first vortex shedding instability ; following this , vortex shedding occurs from an irregular bow shock and results in irregular behaviour . agb winds are important sources of dust and light elements and thus for the evolution of the global ism , but their importance to the dynamics of the ism has not previously been considered . our results show that several vortices can be launched into the ism during the agb phase of evolution . these improve the mixing of material , and provide the local ism with a source of a angular momentum . our simulations suggest the vortices are long - lived since they appear to form stable structures which flow out of the simulation domain . the space density of mass - losing agb stars is low ( @xmath24 kpc@xmath25 ) , and assuming a typical life time of the phase of @xmath26 yrs , each region in the ism may be affected by the tail of an agb star every @xmath27 yrs . we estimate the power input into the ism from agb winds to be @xmath28 ergcm@xmath25s@xmath30 . this is several orders of magnitude less than the contribution of supernovae and winds from massive stars at around @xmath31 ergcm@xmath25s@xmath30 @xcite . we find the temperature of material in the vortices to be on the order of 5 - 10000 k , cooling as the vortex moves downstream . the global ( as opposed to the local ) importance of the vortices may therefore be limited and possibly the angular momentum returned to the local ism may be most important . it is also possible that the galactic distribution of agb stars provides a source of turbulence where supernovae and winds from massive stars do not . one may also expect similar vortices to form behind other types of mass - losing stars . we suggest that the swirls seen in the light echo of v838 mon @xcite could represent such vortices . in our simulations we have modelled the ism with a constant density . in reality , the ism is unlikely to be homogeneous and this will more than likely affect the structures forming around the star . it is possible that inhomogenities would fracture the bow shock , as seen in the case of the pn sh 2 - 188 @xcite , and seed more vortex shedding . further , a star travelling at 200 @xmath4  would travel approximately 100 pc during the agb phase and with a bow shock cross - section of up to @xmath32 pc@xmath33 , interact with a volume of ism of up to 200 pc@xmath34 , which could contain on the order of a few tens of stars which could interact with the vortices .    in our simulations , we have neglected viscosity and observed a variety of behaviour from low to high relative reynolds numbers . magnetic fields may suppress viscosity significantly implying a greater degree of turbulence . however , large amounts of viscosity in the ism would smooth the structures which develop in our simulations . it is for this reason that it would be premature to quantify when to expect vortices for particular stellar parameters . the authors acknowledge very useful comments from the anonymous referee and professor john scalo . cjw acknowledges the support of pparc . the numerical computations were carried out using the cobra supercomputer at jodrell bank observatory . binney , j. , & merrifield , m. 1998 , galactic astronomy ( princeton , nj : princeton university press ) , ch . 10 bond , h. 2006 , preprint(astro - ph/0608222 ) burton , w. b. 1988 , galactic and extragalactic radio astronomy ( new york : springer ) , 295 drew , j. e. , et al . , 2005 , mnras , 362 , 753 elmegreen , b. g. , & scalo , j. 2004 , , 42 , 211 elmegreen , b. g. , & scalo , j. 2004 , , 42 , 275 falle , s. a. e. g. 1991 , , 250 , 581 frank , a. , & mellema , g. 1994 , , 430 , 800 klein , r. i. , mckee , c. f. , colella , p. 1994 , , 420 , 213 klein , r. i. , budil , k. s. , perry , t. s. , bach , d. r. 2003 , , 583 , 245 li , h. , colgate , s. a. , wendroff , b. , liska , r. 2001 , apj , 551 , 874 liska , r. , & wendroff , b. 2003 , siam j. sci . comput . , 25 , 995 van leer , b. 1979 , , 32 , 101 mauron , n. , & huggins , p. j. 2006 , , 452 , 257 mellema , g. , eulderink , f. , & icke , v. 1991 , , 252 , 718 morales - juberias , r. , & dowling t. e. , 2005 , p&ss , 53 , 1221 - 1233 pittard , j. m. , et al . , 2005 , , 361 , 1077 quirk , j. j. 1994 , int . j. numer . meth . , 18 , 555 raymond , j. c. , cox , d. p. , & smith , b. w. 1976 , , 204 , 290 robey , h. f. , et al . , 2002 , phys . lett . , 89 , 085001 schier , f. l. , lindqvist , m. , & olofsson , h. 2005 , , 436 , 633 szentgyorgyi , a. , raymond , j. , franco , f. , villaver , e. & lopez - martin , l. 2003 , , 594 , 874 tilley , d. a. , & pudritz , r. e. 2004 , , 353 , 769 ueta , t. , et al . 2006 , apj , 648 , l39 villaver , e. , garcia - segura , g. , manchado , a. , 2003 , apjl , 585 , l49 wareing , c. j. 2005 , phd thesis , univ . of manchester wareing , c. j. , et al . , 2006a , , 366 , 387 wareing , c. j. , et al . , 2006b , , 372 , l63 williams , g. p. 1997 , jgr , 102 , 9303 woodward , p. r. 1976 , apj , 207 , 484 yamazaki , y. h. , skeet , d. r. , & read , p. l. 2004 , p&ss , 52 , 423 - 445 zijlstra , a. a. , & weinberger , r. 2002 , , 572 , 1006
Abstract:

Answer: vortices have been postulated at a range of size scales in the universe including at the stellar size - scale . whilst hydrodynamically simulating the wind from an asymptotic giant branch ( agb ) star moving through and sweeping up its surrounding interstellar medium ( ism ) , we have found vortices on the size scale of @xmath0 pc to @xmath1 pc in the wake of the star . these vortices appear to be the result of instabilities at the head of the bow shock formed upstream of the agb star . the instabilities peel off downstream and form vortices in the tail of agb material behind the bow shock , mixing with the surrounding ism . we suggest such structures are visible in the planetary nebula sh 2 - 188 .


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 early as 1934 , w.  elsasser noticed the existence of `` special numbers '' of neutrons and protons that confer to the corresponding nuclei a particularly stable configuration  @xcite . in analogy with atomic electrons , he correlated these numbers with closed shells in a model of non - interacting nucleons occupying energy levels generated by a potential well . it took 15 more years until mayer , haxel , suess and jensen  @xcite showed that this nuclear potential could be modeled by a one - body harmonic oscillator ( ho ) , an l@xmath2 term and a spin - orbit ( so ) potential to give rise to the ho shell gaps at 8 , 20 , 40 and so shell gaps at 28 , 50 , 82 and 126 . in 1954 , elliott wrote that this @xmath3 spin - orbit `` term may be a caricature of a more complicated force '' as it must be changed when particular shells are filled  @xcite . to give an example , using yet available experimental data , the n=14 , 28 and 50 so shell gaps increase by the _ same _ value of about 2.7 mev when the neutron orbits d@xmath4 , f@xmath5 and g@xmath6 are filled , respectively(see fig . 10 in  @xcite ) . this points to a universal mechanism that can not be accounted for by one - body potentials for which a single parametrization can not lead to such modifications in a given isotopic chain . indeed recent theoretical achievements point out the role of two and even three - body interactions to account for these behaviors  @xcite . in 1960 , talmi and unna  @xcite noticed that nuclear forces can modify the ordering of the shells , stressing that `` the order of the filling of the neutron shells may depend on the proton configuration '' . these combined works of elliott and talmi and unna introduced the possibility of shell evolutions due to specific two - body nuclear forces . it took however about 20 more years until the combined works of thibault  @xcite , hubert  @xcite , dtraz  @xcite and guillemaud - mueller  @xcite on atomic masses , nuclear radii and nuclear spectra demonstrated the first disappearance of a magic shell ( n=20 ) and the onset of collectivity far from stability . on the theoretical point of view the works of campi , poves and warburton  @xcite described this onset of collectivity as due to the combination of a shell gap reduction and the excitations of particles from the normally occupied orbital to the first orbital of the upper shell , often referred to as _ intruder orbit_. some nuclei , in which the ordering of the intruder and normal configurations is inverted belong to the so - called `` island of inversion ''  @xcite . these pioneering discoveries in the n=20 region opened a new branch of experimental and theoretical researches related to the structure of nuclei far from stability . this triggered the development of radioactive ion beam facilities and of highly sensitive detector arrays , that in turn bring new discoveries in other regions of the chart of nuclides . this territorial exploration of the chart of nuclides goes in concert with the inspection of various facets of the nuclear force as different combinations of protons and neutron orbits are filled , from the proton to the neutron drip lines .    combining all experimental discoveries , some remarkable and general shell evolutions can be observed in the chart of nuclides as the one evoked above about the n=14 , 28 and 50 shell gaps . in the present contribution , the striking analogy between the behavior of the ho gaps n=8 , 20 and 40 from the valley of stability to more neutron - rich regions will be presented in sect . the onset of deformation through intruder configurations at n=20 will be depicted in sect . [ island ] using two recent experimental studies on the @xmath7si and @xmath8 mg isotones . [ n17 ] proposes to study the evolution of the neutron single particle energies of the neutron d@xmath9 , f@xmath5 and p@xmath9 orbitals between which the n=20 and n=28 gaps are formed . the underlying nuclear forces leading to the disappearance of the n=20 shell gap and the swapping between the f@xmath5 and p@xmath9 orbit will be described . [ hierarchy ] shows that a hierarchy in the nuclear forces is responsible for these drastic shell evolutions . a generalization of the hierarchy of nuclear forces mechanism to other regions of the chart of nuclides ( around @xmath10ca , below @xmath11ni and below @xmath12sn where the r - process nucleosynthesis occurs ) is proposed . sect.[drip ] addresses the role of continuum in changing effective nuclear forces in close - to - drip - line nuclei . there are several means to test the existence of neutron shell gaps , based on ( i ) the search of sudden drops in neutron separation energy between adjacent nuclei , ( ii ) the study of evolution of neutron single particle energies , ( iii ) the discovery of rises in 2@xmath13 energy and ( iv ) reduction in the b(e2 ; 2@xmath140@xmath13 ) values in a given isotopic chain , to quote a few . as shown in fig . [ e2plus ] , the increase of the 2@xmath15 energies at the neutron numbers n=8,20 and 40 reveals the presence of large gaps between the last occupied neutron orbit and the valence one . by looking at the same systematics for neutron - rich nuclei in the @xmath16be , @xmath17ne and @xmath18 mg , @xmath19fe and @xmath20cr chains there is no more rise of the 2@xmath13 energies at any of the n=8 , 20 and 40 neutron numbers . this partly arises from the fact that the neutron shell gaps have been reduced there , but also from the fact the proton ( sub-)shells are no longer fully filled , bringing additional correlations provided by this open proton shell . from the similarity of these observations at n=8 , 20 and 40 , we can anticipate that the same process is occurring in all these regions . note for instance that , since the pioneering works of refs  @xcite , many recent experimental studies have been carried out to study the n=40 sub - shell closure worldwide this last decade  @xcite . we shall demonstrate in section . [ hierarchy ] that the quantum characteristics of the orbits in which protons are added do matter for the creation of shell gaps , by specific proton - neutron interaction , from the neutron drip line to the valley of stability . remarkably , far from stability , before these proton - neutron forces act , other local magic numbers at n=16 and n=32 exist , as proposed in refs . @xcite . in the next section we study how abrupt is the transition to deformation in the n=20 isotones by looking at the systematics of 0@xmath21 energies . energies in the n=8 ( be , c , o ) , 20 ( o , ne , mg , si , s and ca ) and 40 ( ni , fe , cr , ti , ca ) regions . the rise in 2@xmath13 energies at n=8 , 20 and 40 no longer exists for neutron - rich nuclei . rather rises in 2@xmath13 energies appear at n=16 and n=32.,width=313 ] the 0@xmath13 ground state of even - even magic nuclei can be described , in a first approximation , with all nucleons paired in the occupied orbitals , the valence orbit higher up in energy being unoccupied . it is known however that , for instance , correlations of pairing origin can slightly dilute the occupancy of pairs to the upper shell . the main configuration of the 0@xmath21 states corresponds to the promotion of a pair of nucleons to the upper shell ( labelled as a 2p2h configuration ) . the energy of this state is twice the energy to promote a nucleon across the shell gap minus the correlation energy gained by the particles and holes . the constancy of the energy of the 0@xmath21 at n=20 in the @xmath22ca , @xmath23ar and @xmath24s nuclei ( e(0@xmath253.5mev ) indicates that the size of the neutron shell gap ( formed between the d@xmath9 and f@xmath5 or p@xmath9 orbits ) and the amount of correlations are unchanged between z=20 and z=16 . with a reduction of the shell gap and the increase of correlations due to the open proton orbit d@xmath4 , an inversion of the 0@xmath13 and 0@xmath21 states should be occurring as predicted more than 30 years ago .     and 0@xmath21 energies in the n=20 isotones . an abrupt inversion between the 0p0h and 2p2h configurations occurs between the @xmath7si and @xmath8 mg nuclei.,width=313 ]    to ascertain this mechanism , the energies and configurations of the 0@xmath13 and 0@xmath21 in the @xmath7si@xmath22 and @xmath8mg@xmath22 , between which the transition happens , should be determined . the study of these missing links was achieved only these last three years . the @xmath26mg(t , p)@xmath8 mg reaction was used at the rex - isolde facility to identify the 0@xmath21 in @xmath8 mg at 1058 kev  @xcite . protons were detected in a highly segmented si detector array ( t - rex ) , while the @xmath27-rays following the 0@xmath282@xmath14 0@xmath13 transitions were observed in the miniball detector array . the 0@xmath21 being produced with a similar cross section than the ground state 0@xmath13 , one can think that being close in energy , the two 0@xmath15 interact and become largely mixed . the 0@xmath21 in @xmath7si was produced at the ganil facility through the @xmath29-decay of a newly discovered isomeric state j=1@xmath15 of 26(1 ) ms in the mother nucleus @xmath7al , whose neutron configuration 2p1h favors the feeding of a 2p2h configuration in @xmath7si through the conversion of a neutron d@xmath9 into a proton d@xmath4  @xcite . the selective decay of this 1@xmath15 isomer though a gamow teller @xmath30-decay ( @xmath31 j=0 @xmath32 1 ) was the key of success for populating the 0@xmath21 state at 2.72 mev that decays by pair creation with a half - life of 19.4(7 ) ns to the 0@xmath13 state . the 2@xmath13 state at 3327 kev was also fed in this @xmath30 decay . information on its competing decay branches to the 0@xmath21 and 0@xmath13 was used to infer that the mixing between the two 0@xmath15 states is of the order of 20%  @xcite . the 0@xmath13 ground state has 80% 0p0h and 20% 2p2h configuration while that of the 0@xmath21 at 2.72 mev is reversed . as shown in fig . [ 0plus ] an abrupt transition from spherical to deformed configurations is observed between @xmath7si and @xmath8 mg , both from the trend in 2@xmath13 and 0@xmath21 energies . concomitant information from the s@xmath33 , 2@xmath13 , 0@xmath21 energies and b(e2 ) systematics suggest an erosion of the n=20 shell gap . but these global properties do not prove that the neutron shell gap is reduced below z=14 . indeed as the z@xmath3412 nuclei that lie in the island of inversion are deformed , states are very mixed and their spherical single particle ( sp ) energies can not be determined there . we have to stand as much as possible @xmath35 the island of inversion to reveal the evolution of sp energies and probe the mechanism that is conducting the shell evolution in the n=20 region . very recent experimental studies using transfer reactions have brought significant insight upon this . they took advantage of the fact that the n=14 and n=16 shell gaps are relatively large in the ne isotopic chain to study the evolution of the neutron n=20 and n=28 gaps at z=10 .    , 7/2@xmath36 and 3/2@xmath36 states in the n=17 isotonic chain . the red and blue dashed lines corresponds to a simple linear extrapolations of the 7/2@xmath36 and 3/2@xmath36 energies towards the drip line , respectively . they show the swapping between the 7/2@xmath36 and 3/2@xmath36 states and the reduction of the n=20 gap but are not intended to have a significant predictive power.,width=309 ]    the @xmath37ne and @xmath38ne , produced by the spiral1 facility at an energy of about 10 @xmath39.mev , impinged into a cd@xmath40 target in which the ( d , p ) reaction took place to populate the 3/2@xmath15 , 7/2@xmath36 and 3/2@xmath36 states that arise mainly from neutron d@xmath9 , f@xmath5 and p@xmath9 configurations , respectively  @xcite . protons were detected by the tiara detector array , @xmath27-rays by the exogam detector , while the transfer - like products were identified at the focal plane of the vamos spectrometer . the energy , spin and spectroscopic factors ( sf ) of the 3/2@xmath15 , 7/2@xmath36 and 3/2@xmath36 states were obtained from the proton energy , angular distribution and cross section , respectively . sf values are in most of the cases larger than 0.6 , suggesting that these states carry a significant fraction of single particle strength . it is however important to bear in mind that the forthcoming discussion is semi - quantitative only and that detailed theoretical calculations should confirm the conclusions drawn here . [ n17fig ] shows the evolution of the 3/2@xmath15 , 7/2@xmath36 and 3/2@xmath36 states in the n=17 isotones based on the result of  @xcite for @xmath41ne . two different regimes can be distinguished in this figure , with a turning point at z=14 , as already found for the n=20 isotones . above z=14 , it is seen that the spacing between the 3/2@xmath15 and 7/2@xmath36 states forming the n=20 gap is almost constant after the addition of 6 protons in the s@xmath42 and d@xmath9 orbits . this means that nuclear forces between the protons in these orbits and neutrons in the two orbits forming the n=20 gap are almost the same . the spacing between the 7/2@xmath36 and 3/2@xmath36 states , which traditionally forms the n=28 gap ( very small here @xmath43 1 mev ) , decreases slightly between z=14 and z=16 and increases afterwards until z=20 . this trend is in agreement with what is predicted in the sdpf shell - model interactions  @xcite . below z=14 two major changes occur . the 7/2@xmath36 and 3/2@xmath36 states are swapped and the n=20 gap collapses . even if states are not of fully single particle origin , this drastic trend should carry valuable information upon the intensity of the proton - neutron forces involved here and their spin , angular momentum and radial overlap dependencies . to give an example on the spin dependance , while the 3/2@xmath15 and 7/2@xmath36 spacing is unchanged by adding four protons in the 1d@xmath9 orbit from z=16 to z=20 , it increases by about 1.2 mev by adding four protons in the 1d@xmath4 orbit from z=10 to z=14 . an even more dramatic effect occurs for the 7/2@xmath36 and 3/2@xmath36 spacing which increases by 1 mev ( 2 mev ) with the filling of the d@xmath9 ( d@xmath4 ) orbit . note that the results on the 3/2@xmath15 , 7/2@xmath36 and 3/2@xmath36 states in @xmath44ne obtained through ( d , p ) reaction also confirm that a reduction of the 3/2@xmath15  7/2@xmath36 spacing and a swapping between the 7/2@xmath36 and 3/2@xmath36 states take place between the @xmath45si and @xmath44ne isotones  @xcite , as between the @xmath46si and @xmath41ne isotones . we shall see in the next section that these modifications can be intuited using known properties of the nuclear force . we shall then propose that a universal mechanism , based on the hierarchy of few nuclear forces , is expected to provide these significant shell evolutions and re - ordering around n=20 , but also in other regions of the chart of nuclides . as shown above , a significant change in the energy of the neutron orbits 1d@xmath9 , 1f@xmath5 and 2p@xmath9 occurs when changing the occupancy of the proton 1d@xmath4 orbit . by using simple properties of the nuclear force , we can show that the corresponding two - body proton - neutron forces v@xmath47 at play here follow this hierarchy @xmath48v@xmath49v@xmath50v@xmath51 . indeed the nucleon wave functions are characterized by the number of nodes @xmath52 of the radial part , the angular momentum @xmath53 and spin orientation with respect to the angular momentum . when the proton and neutron wave functions have the same @xmath52 and @xmath53 values and opposite spin orientations with respect to @xmath53 , the amplitude of the corresponding interaction is maximized  @xcite . as @xmath48v@xmath54 @xmath48 contains all these properties , it has the largest intensity . on the other extreme , when @xmath52 differs , the protons and neutrons radial wave functions oscillate out of phase , hereby reducing their overlap considerably ( see the schematic view of fig .  [ hierarchyfig ] ) . added to the fact that the proton and neutron angular momenta differ by one unit and that the proton and neutron spin orientations are aligned , the @xmath48 v@xmath55@xmath48 interaction has all reasons to be very weak . the remaining interaction @xmath48v@xmath56@xmath48 falls in the middle , @xmath52 is similar for protons and neutrons , but their @xmath53 values differ by one unit , and their spin are aligned with respect to @xmath53 . we let the reader find that such a clean hierarchy can not be obtained when filling the proton 1d@xmath9 orbital during which no such drastic changes in the spacing between neutron orbits is observed . this hierarchy of proton - neutron forces derived at n=17 applies as well to the regions south to @xmath7si and @xmath57si , as shown in fig . [ hierarchyfig ] . below @xmath7si , this accounts for the reduction of the n=20 gap and a likely swapping between the neutron 7/2@xmath36 and 3/2@xmath36 states . below the deformed @xmath57si nucleus , the size of the n=28 gap between the f@xmath5 and p@xmath9 orbits is expected to be further reduced . it follows that n=28 isotones , as @xmath58 mg , will be also deformed as seems to be indicated by the recent work on the neutron - rich @xmath59 mg isotopes @xcite . then no sharp drop in s@xmath33 is expected at ( z@xmath6014,n=28 ) , but rather a smooth variation , bringing the drip - line further from stability as observed in ref .  @xcite . as there is no n=20 effect as well in the mg isotopic chain , it follows as suggested in ref.@xcite that a merging of the two islands of inversions n=20 and n=28 is occurring there . interesting extrapolations of this hierarchy of forces to other regions of the chart of nuclides can be seen in fig . [ hierarchyfig ] . by changing all involved proton and neutron angular momenta by one unit , other shell gaps are concerned . indeed , by applying @xmath61 , the proton orbit becomes 1f@xmath5 , while the neutron orbits are 1f@xmath4 , 1g@xmath6 and 2d@xmath4 . between these neutron orbits the n=40 and 50 are formed , the intensity of which depends on the action of proton - neutron forces involved . starting from the @xmath62ni nucleus , the size of the n=40 gap between 1f@xmath4 ( or 2p@xmath42 ) and 1g@xmath6 as well as of the n=50 gap between the 1g@xmath6 and 2d@xmath4 are expected to be reduced significantly when moving toward @xmath10ca ( see for instance the predictions of ref.@xcite ) . these latter orbits can eventually be inverted there , as proposed in the theoretical calculations of  @xcite that use realistic interactions . this hierarchy of forces will also impact the structure of nuclei south to @xmath11ni by reducing the n=50 gap there . the last extrapolation ( @xmath61 ) will bring us into the region below @xmath63sn . this is a region in which the r - process spend enough time to build the a@xmath43130 peak in the abundance curve of the elements . the location , height and shape of this peak could be traced back from the neutron separation energies ( s@xmath64 ) , the half - lives ( t@xmath42 ) , the neutron delayed emission probabilities ( @xmath65 ) , and the neutron - capture cross sections ( @xmath66 ) of the nuclei located around the @xmath67 shell closure . they vary with the structural evolution of the nuclei . using the hierarchy @xmath48v@xmath68v@xmath69 v@xmath70 both the size of the n=82 gap and the energy between the neutron 1h@xmath71 and 1g@xmath5 orbits are expected to be reduced far from stability ( see fig .  [ hierarchyfig ] ) . we can anticipate two important consequences for the r - process . first the reduction of the n=82 gap is expected to change the location of the r - process path , changing the shape of the r - process peak accordingly . instead of being all accumulated at n=82 , isotopes with lower neutron number will become r - process progenitors , filling the low - mass wing of the a=130 peak . second the @xmath29-decay half - lives t@xmath42 , that mainly depend on the gt transition @xmath721g@xmath731g@xmath6 , will be drastically shortened far from stability  @xcite , speeding up the r - process accordingly . indeed as t@xmath42 scales with the phase factor value 1/(q@xmath74-e*)@xmath75 , where e * is the energy to which the gt transition occurs , the reduction of the neutron 1g@xmath5  1h@xmath71 spacing is decreasing e * and shortening the half - lives t@xmath42 of the nuclei that are paving the r - process . therefore the proton - neutron interaction strongly influence the half - lives values in the n=82 isotones . to finish this section few remarks can be put forward . _ first _ , we proposed that a hierarchy of nuclear forces is applicable to several regions of the chart of nuclides , leading to significant shell evolutions . very recent calculations using realistic forces have been made for the o ( ca ) isotopic chain , showing that the inversion between the neutron fp ( gd ) shell has occurred when reaching the drip line . this simplified vision of the hierarchy of nuclear forces somehow finds echo with the properties of realistic forces . _ second _ , the same forces play a role while changing occupancies of the proton orbits 1d@xmath4 , 1f@xmath5 and 1g@xmath6 by up to 6 , 8 and 10 protons , respectively . a relatively good predictive power of the evolution of shells is therefore expected far from stability , provided that the aforementioned interactions are well constrained in a given region of mass . _ third _ , correlations have to be taken into account to fully model the nuclei along the regions shown in fig . [ hierarchyfig ] . _ fourth _ , in the n=82 region other shells are present between the major ones mentioned above . while they are not expected to change this global picture , they act to dilute the action of the aforementioned nuclear forces as other intermediate shells are filled . it is therefore expected that structural evolutions will be occurring more slowly in heavy nuclei . this could explain why recent isomer studies on the n=82 isotones of @xmath76cd ( z=48 ) and @xmath77pd ( z=46 ) @xcite do not indicate that a change in nuclear structure has yet occurred there . important is however to note that these two experiments are rather sensitive to proton configurations ( g@xmath78 and rather indirectly to neutron configurations . therefore , a more direct evolution of neutron single - particle energies should be probed there in the future in order to see if the n=82 gap is reduced below z=50 . _ fifth _ , the same hierarchy seems also applicable to lighter nuclei to account for the evolution of the neutron 2s@xmath42 , 1p@xmath42 and 1d@xmath4 orbits when the proton 1p@xmath9 is filled in the n=7 isotones between the recently studied @xmath79he @xcite and @xmath80c . in the previous sections , it was assumed that the proton - neutron nuclear forces do not change when approaching the neutron drip - line . there , a large asymmetry is present between the proton and neutron binding energies . while the proton is deeply bound and its wave function well confined inside the nucleus , the neutron is weakly bound and its wave function can extend out of the mean nuclear radius . in such a case , the overlap of their wave functions is reduced and the effective proton - neutron interaction is weakened . this change in effective interaction depends on the binding energy of the neutron and its angular momentum . indeed when the neutron centrifugal barrier is high as compared to the binding energy of the unbound neutron , it can be considered as quasi - bound and the effective interaction can be similar to that in the valley of stability . several experimental works are being carried out worldwide ( e.g.  @xcite ) to probe the structure of unbound states and theoretical works are treating the interaction with the continuum explicitly  @xcite .        as described in ref.@xcite the @xmath81f@xmath82 nucleus is a suitable choice for studying the proton - neutron d@xmath4-d@xmath9 interaction that has a decisive role for changing the n=20 gap when approaching the drip lines . the binding energies of all states of the multiplet @xmath83 arising from the @xmath84d@xmath85d@xmath9 coupling in @xmath81f@xmath82 were determined using four experimental techniques ( time of flight measurement for the j=1 ground state  @xcite ) , in - beam @xmath27-ray spectroscopy for the j=2 state  @xcite , m3 isomer decay study for the j=4 state  @xcite , and in - beam neutron spectroscopy for the unbound j=3 state  @xcite . the corresponding proton - neutron interaction energy , int(j ) , was determined in @xcite from the difference between the experimental binding energy of a state j in @xmath38f and that of the @xmath37o+1p+1n system , in which proton and neutron are free of interaction ( int=0 ) . [ int ] compares the experimental results to shell model ( sm ) calculations using the usda interaction fitted to nuclei close to the valley of stability and to coupled cluster ( cc ) calculations that use interactions from chiral effective field theory , the treatment of three - body forces and the interaction with the continuum . comparison with sm calculations shows that the experimental binding energy is globally weaker than expected by about 17% , and that the multiplet of j=1 - 4 states is compressed by about 25% . these discrepancies possibly point out the deficiencies of the sm that is treating nucleon interactions using a harmonic oscillator potential basis that is not suited for taking the proton - to - neutron binding energy asymmetry and continuum effects into account . on the other hand cc calculations , in which the weak binding energy of the neutron is considered explicitly , agree with experimental results . further studies on other nuclei should be carried out to confirm the present findings upon the proton - neutron interaction when approaching the drip - line . the field of nuclear structure far from stability is experiencing tremendous breakthroughs . benchmarking nuclei as @xmath37o , @xmath57si , @xmath86ca , @xmath87cr , @xmath88sn and @xmath12sn - to quote a few - have been or are on the way to be studied . extremes of isospin , ranging from the proton to neutron drip lines were probed , as for instance between @xmath89f@xmath90 and @xmath91f@xmath92 . when sorting out the data available to date , amazingly similar pictures arise throughout the chart of nuclides , which have been recently identified to take root from specific properties of the nuclear force and from nuclear correlations . the magic nuclei n= 8 , 20 and 40 ( leading to so - called ho shell gaps ) all disappear far from stability , leading to three islands of inversion through the same mechanism . it was proposed here that the reduction of the n=8 , 20 , and 40 shell gaps and the swapping between certain levels arise from a simple hierarchy between proton - neutron forces . this hierarchy was applied to yet unexplored regions of the chart of nuclides such as below @xmath57si , @xmath87cr , @xmath11ni or below @xmath12sn where the r - process nucleosynthesis is accumulating peaks in the r abundance curve of the elements . as for the the growing of the so - called so shell gaps at n=14 , 28 and 50 by about 2.7 mev , it is likely due to three - body neutron - neutron forces , which were only recently implemented in theoretical models . these remarkable successes and similarities in different regions of the chart of nuclides however reach limit when looking at the evolution of so shell gaps far from stability . for so far unknown reasons , the valence mirror nuclei built on proton and neutron so numbers such as @xmath93c@xmath94 , @xmath95si@xmath96 and @xmath63sn@xmath97 , behave totally differently . the first two are deformed , while the latter is a doubly magic nucleus . it will be therefore very important to study the @xmath98ni@xmath99 nucleus , which lies in between , to better understand why similar forces do not reach the same effect there . w. elsasser , j. de phys . * 5 * 625 ( 1934 ) m.  goeppert - mayer , phys . rev . * 75 * 1969 ( 1949 ) o. haxel , j.h.d . jensen and h.e . suess , phys . rev . * 75 * 1766 ( 1949 ) j. p. elliott and a. m. lane , phys . rev . * 96 * 1160 ( 1954 ) o. sorlin and m .- porquet , phys . t * 152 * 014003 ( 2013 ) j. d holt et al . , j. phys . g * 39 * , 085111 ( 2012 ) k. sieja and f. nowacki , phys . c * 85 * 051301 ( r ) ( 2012 ) i. talmi and i. unna , phys lett . * 4 * 469 ( 1960 ) c. thibault et al . c * 12 * 644 ( 1975 ) g. hubert et al . c * 18 * 2342 ( 1978 ) c. dtraz et al . c * 19 * 164 ( 1979 ) d. guillemaud - mueller et al . , nucl . phys . a * 426 * 37 ( 1984 ) x. campi et al . a * 251 * 193 ( 1975 ) a. poves and j. retamosa , phys . b * 184 * 311 ( 1987 ) w. warburton , j. a. becker and b. a. brown , phys . c * 41 * 1147 ( 1990 )    m. hannawald et al . lett . * * 82**1391 ( 1999 ) ; o. sorlin et al . j. a * 16 * 55 ( 2003 ) t. baugher et al . c * 86 * 011305(r ) ( 2012 ) ; a. gade et al . c * 81 * 051304(r ) ( 2010 ) ; p. adrich et al . c * 77 * 054306 ( 2008 ) ; j. ljungvall et al . , phys . c * 81 * 061301(r ) ( 2010 ) ; w. rother et al . lett . * 106 * 022502 ( 2011 ) ; a. brger et al . , phys . lett . b * 622 * , 29 ( 2005 ) ; n. aoi et al . * 102 * 012502 ( 2009 ) ; h. l. crawford et al . lett . * 110 * 242701 ( 2013 )    a. ozawa et al . , phys .  rev . * 84 * 5493 ( 2000 ) c. r. hoffmann et al . b * 672 * 17 ( 2009 ) k. tshoo et al . * 109 * 022501 ( 2012 ) a. huck et al . c  * 31 * 2226 ( 1985 ) f. wienholtz et al . , nature * 498 * 346 ( 2013 ) k. wimmer et al . , phys .  rev . lett .  * 105 * 252501 ( 2010 ) f. rotaru et al . phys .  rev . lett .  * 109 * 092503 ( 2012 ) w.n . catford et al . , phys . * 104 * 192501 ( 2010 ) s. m. brown et al . phys . rev . c  * 85 * 011302 ( r ) ( 2012 ) f. nowacki and a. poves , phys c * 79 * 014310 ( 2009 ) l. gaudefroy et al . , phys . * 97 * 092501 ( 2006 ) t. otsuka , phys . t * 152 * 014007 ( 2013 ) p. doornenbal et al . lett . * 111 * 212502 ( 2013 ) t. baumann et al . , nature * 449 * 1022 ( 2007 )      j.j . cuenca - garcia et al . , eur . j. a * 34 * 99 ( 2007 ) i. dillmann et al . * 91 * 162503 ( 2003 ) a. jungclaus et al . * 99 * 132501 ( 2007 ) h. watanabe et al . lett . * 111 * 152501 ( 2013 ) t. al kalanee et al . c * 88 * 034301 ( 2013 ) a. lepailleur et al . lett . * 110 * 082502 ( 2013 ) t. baumann , a. spyrou and m. thoennessen rep . . phys . * 75 * 036301 ( 2012 ) c.forss ' en et al . t * 152 * 014022 ( 2013 )
Abstract:

Answer: during the last 30 years , and more specifically during the last 10 years , many experiments have been carried out worldwide using different techniques to study the shell evolution of nuclei far from stability . what seemed not conceivable some decades ago became rather common : all known magic numbers that are present in the valley of stability disappear far from stability and are replaced by new ones at the drip line . by gathering selected experimental results , beautifully consistent pictures emerge , that very likely take root in the properties of the nuclear forces.the present manuscript describes some of these discoveries and proposes an intuitive understanding of these shell evolutions derived from observations . extrapolations to yet unstudied regions , as where the explosive r - process nucleosynthesis occurs , are proposed . some remaining challenges and puzzling questions are also addressed . @xmath0 : the present manuscript is a summary of a plenary talk given at the 25@xmath1 international nuclear physics conference ( inpc ) , firenze , italy , 2 - 7 june 2013 . some new remarkable experimental results released since then have been implemented here .