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: spatially inhomogeneous cosmologies constitute almost a `` terra incognita '' for general relativity . in particular , the nature of the singularity in generic cosmologies is unknown although belinskii , khalatnikov , and lifshitz ( bkl ) have conjectured that it is locally mixmaster - like @xcite . spatially homogeneous universes approach to the singularity can be either asymptotically velocity term dominated ( avtd ) @xcite or mixmaster - like @xcite . in general , avtd singularities occur when the influence of spatial derivatives can be neglected . this means that the solution to the full einstein equations eventually comes arbitrarily close to the solution to the truncated equations one obtains by eliminating the terms containing spatial derivatives @xcite . in spatially homogeneous models , spatial derivatives in the metric create the spatial scalar curvature which appears as a potential in minisuperspace @xcite . an avtd homogeneous cosmology eventually evolves toward the singularity as a kasner solution . ( the kasner universe is an exact solution of the vacuum einstein equations characterized by different expansion rates in different directions with flat spacelike hypersurfaces @xcite . ) mixmaster dynamics describes the evolution ( of non - flat space - like hypersurfaces ) toward the singularity as an infinite sequence of ( approximately ) kasner epochs . the kasner indices ( which determine the anisotropic expansion rates ) change whenever the influence of the spatial scalar curvature reappears . its recurring influence means that the mixmaster singularity is not avtd . bkl claimed to prove that the generic singularity in einstein s equations and , in particular , in spatially inhomogeneous cosmologies is mixmaster - like . while this result remains controversial @xcite , it provides a conjecture which can be tested by evolving spatially inhomogeneous collapsing universes numerically .    as part of such a study of this issue @xcite , we have considered the gowdy universes on @xmath0 as a test case . these models are one of the classes discovered by gowdy @xcite as a reinterpretation of einstein - rosen waves @xcite . while their plane symmetry precludes local mixmaster behavior , their nonlinear gravitational wave interactions lead to an interesting phenomenology which forms the subject of this paper . the gowdy singularity has been conjectured to be asymptotically velocity term dominated ( avtd ) @xcite and has been shown rigorously to be so for the polarized case @xcite . in addition , it has been shown that for the polarized case , cosmic censorship holds : there exists a natural foliation labeled by @xmath2 such that regular initial data evolve without singularity for all @xmath3 where ( in the collapse direction ) a curvature singularity occurs at @xmath4 @xcite . here we shall use numerical studies to provide support for the extension of the polarized gowdy results to generic gowdy models . gowdy cosmologies represent the simplest spatially inhomogeneous cosmological solutions to einstein s equations . the vacuum case may be interpreted as the two polarizations of gravitational waves propagating in an inhomogeneous background spacetime @xcite . the model is plane symmetric with the propagation direction of the waves orthogonal to the symmetry plane . near the singularity , the metric may be rewritten as locally kasner with spatially dependent kasner indices @xcite . ( this is a non - rigorous statement of avtd behavior . ) . matter fields compatible with the gowdy metric have been discussed by carmeli , charach , and malin @xcite , and ryan @xcite . here we shall consider only vacuum solutions . the main phenomena of interest in the gowdy universe are the growth of small scale spatial structure and the approach to the avtd regime . we shall demonstrate that both result from the nonlinear terms in the wave equation for @xmath5 , the @xmath6 polarization of the gravitational waves . these act as space and time dependent potentials whose effect is to drive the system to velocity dominance . as the avtd regime is reached , the potentials approach zero . their spatial dependence means that wave amplitudes grow at some spatial points and decrease at others , eventually producing small scale spatial dependence in the waveforms . non - generic behavior occurs at points where the coefficients of the potentials ( separately ) vanish identically . we note that this analysis in terms of potentials explains the observed episodic evolution of the system to the avtd regime and the fact that the conditions for avtd are reached at some spatial points before others . this should be considered together with the asymptotic expansion developed by grubisi and moncrief ( gm ) @xcite which displays an exponential decay to the avtd limit , applicable in the last episode . the primary advantage of the gowdy models as a test case for the study of inhomogeneous cosmologies is the existence of variables in which the dynamical equations for the wave amplitudes decouple from the constraints . initial data for the wave equations may be freely specified ( restricted only in that the total momentum along the direction of propagation must vanish ) while the constraints appear as prescriptions for the construction of the background metric from the known solution to the dynamical equations . thus the two principal difficulties of numerical relativity  solving the initial value problem and preserving the constraints during the evolution  become trivial . this decoupling between dynamical and constraint equations disappears in more complicated inhomogeneous models . however , studies of @xmath7 symmetric cosmologies ( with only one killing field ) indicate that many features of the gowdy phenomenology such as spiky small scale spatial structure persist in these more complicated models @xcite and , presumably , would be important in the dynamics of generic cosmologies .    in section 2 , the gowdy model on @xmath0 will be reviewed . section 3 will contain a brief description of the numerical method . in section 4 , the phenomenology will be described . section 5 will contain a unified explanation for the phenomenology . discussion will be given in section 6 . the gowdy model on @xmath0 is described by the metric @xcite @xmath8\end{aligned}\ ] ] where @xmath9 , @xmath5 , @xmath10 are functions of @xmath11 , @xmath12 . the sign of @xmath9 differs from that in @xcite as required for consistency with the form given there for the equations for @xmath9 . the arguments in @xcite did not depend on this error . we impose @xmath13 spatial topology by requiring @xmath14 and the metric functions to be periodic in @xmath11 . if we assume @xmath5 and @xmath10 to be small , we find them to be respectively the amplitudes of the @xmath6 and @xmath15 polarizations of the gravitational waves with @xmath9 describing the background in which they propagate . the time variable @xmath12 measures the area in the symmetry plane with @xmath4 a curvature singularity . einstein s equations split into two groups . the first is nonlinearly coupled wave equations for @xmath5 and @xmath10 ( where @xmath16 ) : @xmath17 @xmath18 the second contains the hamiltonian and @xmath11-momentum constraints which can be expressed as first order equations for @xmath9 in terms of @xmath5 and @xmath10 : @xmath19=0,\ ] ] @xmath20 this split into dynamical and constraint equations removes two of the most problematical areas of numerical relativity from this model : ( 1 ) the normally difficult initial value problem becomes trivial since @xmath5 , @xmath10 and their first time derivatives may be specified arbitrarily . the only restriction , that the total @xmath11 momentum in the waves vanishes , is a consequence of the requirement that metric variables be periodic in @xmath11 . ( 2 ) the constraints , while guaranteed to be preserved in an analytic evolution by the bianchi identities , are not automatically preserved in a numerical evolution with einstein s equations in differenced form . however , in the gowdy model , the constraints are trivial since @xmath9 may be constructed from the numerically determined @xmath5 and @xmath10 . for the special case of the polarized gowdy model ( @xmath21 ) , @xmath5 satisfies a linear wave equation whose exact solution is well - known @xcite . for this case , it has been proven that the approach to the singularity is avtd @xcite . this has also been conjectured to be true for generic gowdy models @xcite . we shall show in section 5 how a numerical study of this model provides strong support for this conjecture . the wave equations ( [ gowdywavep ] ) and ( [ gowdywaveq ] ) can be obtained by variation of the hamiltonian @xcite @xmath22}\nonumber \\    & + & { 1 \over 2}\int\limits_0^{2\pi } { d\theta \,\left [ { e^{- 2\tau } \left ( { p,_\theta ^2+\;e^{2p}q,_\theta ^2 } \right ) } \right]}=h_k+h_v .\end{aligned}\ ] ] the variation of the kinetic energy term , @xmath23 , yields equations of motion for the avtd solution which arises when spatial derivatives are neglected . these equations are exactly solvable with the solution @xmath24 \ ; \ ; \to v\tau \quad { \rm as}\;\tau \to \infty , \nonumber \\ q&=&q_0 \ ; + \ ; { { { e^{- { p_0 } } } \ , \sin \psi \ , \tanh v \tau } \over { 1 \ , + \ , \cos \psi \ , \tanh v \tau } } \quad \quad \quad \;\to q_\infty \quad { \rm as}\;\tau \to \infty , \nonumber \\ \pi _ p&=&v \ ; { { \tanh v \tau \ , + \ , \cos \psi } \over { 1 \ , + \ , \cos \psi \ , \tanh v \tau } } \quad \quad \quad \quad \quad \quad \to v\quad \;\;{\rm as}\;\tau \to \infty , \nonumber \\ \pi _ q&= & { e^{p_0 } } \ , v \ , \sin \psi \equiv \pi_q^0\end{aligned}\ ] ] given in terms of four functions of @xmath11 : @xmath25 , @xmath26 , @xmath27 , and @xmath28 . here the large @xmath30 limits given above are for @xmath31 . the special case @xmath32 will be discussed later . while ( [ avtdeq ] ) may be solved for all these functions , it is convenient to note that @xcite @xmath33 is the geodesic velocity in the target space ( with metric @xmath34 ) of the wavemap @xcite . if , as the singularity is approached , the behavior is avtd , the true solution will approach the avtd one . the exponential prefactor @xmath35 in @xmath36 of ( [ gowdywaveh ] ) makes plausible the conjectured avtd singularity . however , @xmath37 ( for @xmath38 ) as @xmath39 . if @xmath40 , the term @xmath41 in @xmath36 can grow rather than decay as @xmath39 . in fact , if one assumes the avtd behavior as @xmath42 in ( [ avtdeq ] ) , the wave equations become @xcite @xmath43\tau}[q'_\infty ] ^2+e^{-2\tau}v '' \tau + \pi_q^2 e^{-2 v \tau}\nonumber\\ \pi_q,_\tau & = & e^{-2[1-v]\tau}[q_\infty '' + 2 v ' q_\infty ' \tau]\end{aligned}\ ] ] where @xmath44 . these equations are inconsistent if @xmath40 unless @xmath45 , which does not occur in a generic ( unpolarized ) gowdy model . however , in generic gowdy models , one also expects isolated points ( a set of measure zero ) with @xmath46 but @xmath47 . in section 5 , we shall see that @xmath40 can persist at such isolated values of @xmath11 where the avtd conditions may not be satisfied . gm thus conjecture that the avtd limit as @xmath42 requires @xmath48 everywhere except , perhaps , at a set of measure zero ( isolated values of @xmath11 ) @xcite . we note that the polarized gowdy model ( @xmath49 ) has the avtd solution @xmath50 , @xmath51 with @xmath52 unrestricted in sign or magnitude @xcite . the fact that this behavior is not a limit of the generic avtd solution ( [ avtdeq ] ) illustrates further the qualitative distinction between the polarized case and @xmath53 but arbitrarily small . the numerical method is discussed in detail elsewhere @xcite . however , we shall summarize it here . the work reported here was performed by using a symplectic pde solver@xcite . consider a system with one degree of freedom described by @xmath54 and its canonically conjugate momentum @xmath55 with a hamiltonian @xmath56 note that the subhamiltonians @xmath23 and @xmath36 separately yield equations of motion which are exactly solvable no matter the form of @xmath57 . variation of @xmath23 yields @xmath58 , @xmath59 with solution @xmath60 variation of @xmath36 yields @xmath61 , @xmath62 with solution @xmath63 note that the absence of momenta in @xmath36 makes ( [ hvsoln ] ) exact for any @xmath64 . one can then demonstrate that to evolve from @xmath65 to @xmath66 an evolution operator @xmath67 can be constructed from the evolution sub - operators @xmath68 and @xmath69 obtained from ( [ hksoln ] ) and ( [ hvsoln ] ) . one can show that @xcite @xmath70 reproduces the true evolution operator through order @xmath71 . suzuki has developed a prescription to represent the full evolution operator to arbitrary order @xcite . for example @xmath72\,{\cal u}_{(2)}(s\delta t)\ ] ] where @xmath73 . the advantage of suzuki s approach is that one only needs to construct @xmath74 explicitly . @xmath75 is then constructed from appropriate combinations of @xmath76 . the generalization of this method to @xmath77 degrees of freedom and to fields is straightforward . in the latter case , @xmath78 \to v[\vec q(\vec x , t)]$ ] so that @xmath79 becomes the functional derivative @xmath80 . on the computational spatial lattice , the derivatives that are obtained in the expression for the functional derivative must be represented in differenced form . we note that , to preserve @xmath81th order accuracy in time , @xmath81th order accurate spatial differencing is required . some discussion of this has been given elsewhere @xcite . convergence studies have been performed to test the algorithm . we have already shown @xcite that at any fixed value of @xmath12 , there exists a spatial resolution at which every spiky feature is resolved . any finer resolution will yield identical results . however , as we shall see , generic spiky features grow narrower as the singularity is approached . this means ( see @xcite ) that a spatial resolution which was adequate at some @xmath82 will fail to resolve all features at some @xmath83 . this type of `` resolution dependence '' is well understood to be a consequence of the physical behavior of the model . if the physical system has features which become arbitrarily narrow as @xmath39 , no resolution will be adequate everywhere for all time . here we report the results of a numerical study of the evolution of generic gowdy models toward the singularity ( @xmath39 ) . rather than perform simulations for hundreds of gowdy models , we consider a single class of models and argue that the observed behavior is generic . for the most part , we shall use the initial data @xmath84 , @xmath85 , @xmath86 , and @xmath87 . this model is actually generic for the following reasons : the @xmath88 dependence is the smoothest nontrivial possibility . since we shall discuss the growth of small scale spatial structure , we expect the cleanest effects to arise from the smoothest initial data . a more complicated initial state will not yield any qualitatively new phenomena . for example , choosing @xmath89 instead yields the same solution repeated @xmath81 times on the grid thus yielding the same result with poorer resolution . in addition , the amplitude of @xmath10 is irrelevant since the hamiltonian ( [ gowdywaveh ] ) is invariant under @xmath90 , @xmath91 for any constant @xmath92 . this also means that any unpolarized model is qualitatively different from a polarized ( @xmath93 ) one no matter how small @xmath10 . traveling waves ( subject to @xmath94 where @xmath95 is the total @xmath11-momentum ) yield qualitatively similar behavior @xcite . we are interested in the spatial structure that can develop from smooth initial data . for this reason , we have used @xmath88 spatial dependence . starting with a more complicated wave form will yield an eventually more complicated spatial structure but nothing qualitatively different will appear . similarly , since rescaling @xmath10 does not change the solution qualitatively , the amplitude of @xmath10 is kept fixed . in this discussion , we shall consider only @xmath5 and @xmath10 and construct @xmath9 later if desired .    from the chosen class of initial data , with high spatial resolution , one obtains after some time ( @xmath96 ) the profiles for @xmath5 and @xmath10 shown in fig . [ gowdypq ] . the numbered features are typical for any generic solution . the peaks labeled # 1 , # 2 , and # 3 arise from the same features of the nonlinear interactions . on a finer spatial scale in fig . [ pqfine ] , we notice that ( a ) the peak in @xmath5 is associated with an extremum of @xmath10 ( which may be a maximum or minimum ) and ( b ) @xmath97 becomes large with @xmath98 for a maximum ( minimum ) in @xmath10 . the apparent discontinuity in @xmath10 ( labeled # 4 in fig .  [ gowdypq ] ) is shown on a finer scale in fig . [ qdisc ] . we see that this feature occurs as @xmath99 goes through zero where @xmath100 . in fig . [ pqfine18 ] , the same features # 1 and # 4 are shown at @xmath101 . we see that the `` spikes '' have narrowed in @xmath11 .    the evolution of @xmath5 and @xmath10 is shown in fig . [ pqevol ] ( see also @xcite ) where it is clearly seen that ever shorter wavelength modes develop until the point at which the avtd behavior sets in . [ scaling ] shows an interesting self - similarity in the early evolution . the parameter @xmath102 in the initial data is varied in a sequence of simulations . in each case , the number of peaks in @xmath5 is counted after each time step . we find @xmath103 where @xmath104 is the time at which @xmath77 peaks are first present in the waveform and @xmath105 is a constant . from ( [ selfsim ] ) , @xmath106 denotes the value of @xmath102 such that the @xmath77th peak appears only at @xmath4 . the scaling shown in fig . [ scaling ] is for @xmath107 . similar ( but less accurate ) straight line fits are found for @xmath77 equal 3 and 7 . ( the wave equations are even in @xmath11 as are the initial data so that @xmath5 and @xmath10 remain even functions of @xmath11 . thus , except for a central peak , other peaks will form two at a time . ) thus , one may summarize that spatial structure forms more rapidly and evolves a more complicated form , the greater the initial `` energy '' in the waves as denoted by , e.g. , @xmath108 from ( [ vdef ] ) . eventually , the waveforms of @xmath5 and @xmath10 `` freeze '' to the avtd behavior of @xmath5 increasing linearly and @xmath10 constant in @xmath12 . as has been emphasized by hern and stuart ( hs ) @xcite , the duration in @xmath12 of the observed non - avtd regime depends on the spatial resolution of the simulation .    in fig . [ vmax ] , the maximum value over the spatial grid , @xmath109 , of the avtd parameter @xmath110 is shown as a function of time @xmath12 for two simulations with different spatial resolutions . the steps in @xmath109 occur when the location of the spatial point with the maximum value of @xmath110 changes . note that , for the lower resolution simulation , @xmath109 eventually falls below unity in support of the gm conjecture @xcite . we also see that the higher resolution simulation starts to diverge from the lower resolution one , signaling the presence of narrow features that are unresolved by the coarser resolution . ( the former was not run longer primarily because the courant condition requires small time steps for high resolution studies . ) this suggests ( as confirmed by hs ) that the higher spatial resolution simulation will take longer before it becomes avtd everywhere . this is shown explicitly in fig . [ pipres ] where @xmath111 vs @xmath12 is plotted for three different spatial resolutions . note that for two of the spikes , @xmath112 is larger for the finest resolution . in fig . [ pbyt ] , the difference between @xmath113 and @xmath110 is plotted vs.  @xmath12 . since @xmath114 as @xmath39 , we see the influence of the next ( constant in @xmath12 ) term @xmath115 $ ] , in the avtd solution for @xmath5 ( [ avtdeq ] ) . to understand the gowdy phenomenology of section iii , we must explore the structure of the wave equations for @xmath5 and @xmath10 . let us consider ( [ gowdywaveq ] ) first . the equation for @xmath10 may be rewritten in first order form as @xmath116 this immediately provides the explantion for both the large central peak which quickly forms in @xmath10 and for the apparent discontinuity in @xmath10 seen in figs .  [ gowdypq ] and [ qdisc ] . from ( [ qequations ] ) , @xmath10 will grow rapidly where @xmath5 is negative . since , early in the evolution , @xmath117 , one expects a large peak ( as is observed in figs . [ gowdypq ] and [ pqevol ] ) where @xmath118 . we shall see later that @xmath5 is driven to positive values eventually stopping the growth of @xmath10 . however , we shall also see that if @xmath119 , @xmath5 can remain negative longer . say @xmath120 for @xmath121 at a value of @xmath122 where @xmath123 . thus we find that @xmath124 so that @xmath125 is large and @xmath10 increases ( decreases ) exponentially depending ( say for @xmath126 ) on whether @xmath127 . note that we assume @xmath128 . however , ( [ qequations]b ) implies this for @xmath5 large and negative . more interesting is the equation for @xmath5 . consider two approximate forms of ( [ gowdywavep ] ) depending on which nonlinear term dominates . since the nonlinear terms depend exponentially on @xmath5 , they will dominate the linear spatial derivative term when the exponential is large . for example , if @xmath5 is large and negative , @xmath129 which has the first integral @xmath130 this allows one to define an effective potential @xmath131 which is important when @xmath132 and ( especially if ) @xmath133 . this allows immediate explanation of the following phenomena : ( a ) in the avtd regime , @xmath134 and its asymptotic time derivative @xmath38 since if @xmath132 , a bounce off @xmath135 will occur at which @xmath136 so that @xmath5 will eventually reach positive values . ( b ) however , if @xmath119 , @xmath137 so that both @xmath5 and @xmath112 can remain negative for a long time . non - generic behavior occurs at the isolated point @xmath122 where @xmath135 vanishes since @xmath5 and @xmath112 will never be driven to positive values there .    on the other hand , if ( a ) @xmath5 is large and positive and ( b ) @xmath138 , the wave equation for @xmath5 is given by @xmath139 which is easy to integrate ( for @xmath140 ) as @xmath141 thus if @xmath138 , the dynamics of @xmath5 is dominated by an effective potential @xmath142 again , from ( [ qequations]a ) , our implicit assumption that @xmath143 is consistent . one sees immediately , that if @xmath138 a bounce off @xmath144 will occur with the result that @xmath145 or @xmath146 thus the existence of @xmath144 provides a mechanism to drive @xmath112 and thus asymptotically @xmath110 to values below unity . in the actual simulation , the asymptotic regions for the scattering are not reached so that ( [ zbounce ] ) can not be used to compute the numerical value of the change in @xmath112 with @xmath12 . again one sees that non - generic behavior can arise , in this case , at @xmath147 where @xmath148 . at @xmath149 the mechanism to drive @xmath110 below unity is absent so that larger values of @xmath110 are allowed . if @xmath144 is non - vanishing but small , it will take a long time for the bounce to occur so that the avtd limit will take a long time to arise . this also explains the peaks in @xmath5 seen in fig . [ gowdypq ] . where @xmath150 , @xmath144 is small so that @xmath5 can increase without hindrance with a large positive value of @xmath112 .    for a more quantitative understanding of the small scale spatial structure , note that both equations ( 18 ) and ( 21 ) have closed form solutions . equation ( 18 ) holds in the avtd regime , where the solution is given in equation ( 7 ) . the exceptional point @xmath151 is the point where @xmath152 . near this point we have @xmath153 where @xmath154 . it then follows that for large @xmath30 @xmath155 \nonumber \\ q \approx { q_0 } \ ; - \ ; { { { e^ { - { p_0 } } } \ , { \psi ' } \ , ( \theta - { \theta _ 1 } ) } \over { 2 \ , { e^ { - 2 v \tau } } \ , + \ , { { ( \psi ' ) } ^2 } { { ( \theta - { \theta _ 1})}^2}/2}}\end{aligned}\ ] ] thus there is a spiky feature in both @xmath5 and @xmath10 at @xmath156 as is shown in fig .  [ analytic ] . furthermore , these features steepen exponentially with time . because of this steepening , the spatial derivatives of @xmath5 and @xmath10 become large at late times . however , the avtd equations are formed by dropping spatial derivatives from the exact equations ( 2 ) and ( 3 ) . therefore , one might worry that the terms that have been neglected are not actually negligible near the peaks . however , using equation ( [ peqgm ] ) one can show that near @xmath157 the neglected terms go like @xmath158 and are therefore exponentially damped for @xmath159 . the solution of equation ( 21 ) is @xmath160 \ ; \ ; \ ; .\ ] ] here , @xmath161 , and @xmath162 are functions of @xmath163 , and @xmath164 . note that @xmath165 for large @xmath30 where @xmath166 . the exceptional point @xmath167 is the point where @xmath168 . near this point we have @xmath169 where @xmath170 . it then follows that for large @xmath30 @xmath171.\ ] ] thus the feature in @xmath5 narrows exponentially with time as can be seen in fig . [ analytic ] . the spatial derivatives of @xmath5 become very large . however , the neglected terms in the equations go like @xmath172 . since @xmath173 for @xmath174 asymptotically , it follows that these terms become negligible after the last bounce .    the differences in figs . [ vmax ] and [ pipres ] seen by increasing the spatial resolution are also understood in terms of our analysis of structure growth . at higher spatial resolution it is more likely that a spatial grid point of the numerical simulation will be very close to a non - generic point where @xmath148 . thus the non - generic behavior will become more apparent and the time to reach the avtd limit characterized by @xmath159 everywhere will be longer . clearly , the observed resolution dependence may be interpreted to be evidence for the existence of non - generic points . the apparent discontinuity in @xmath10 seen in figs . [ gowdypq ] and [ qdisc ] is evidence for the non - generic point in @xmath11 where @xmath87 . its effect may be seen as a resolution dependence of the shape of this feature . it becomes steeper and narrower with increasing spatial resolution for sufficiently large @xmath12 ( see fig .  [ qspikeres ] ) or for increasing @xmath12 with fixed spatial resolution ( see fig .  [ pqfine18]b ) . as in fig . [ pqfine ] , we note that the spikes in fig . [ pipres ] where @xmath40 are associated with extrema in @xmath10 . this is shown in fig . [ pipvsq ] . note that the alignment is not perfect between the peak in @xmath112 and the extremum in @xmath10 where the spike is well resolved ( see also fig . [ qdisc]a ) . this is understood as evidence that the avtd regime has not yet been reached . in particular , the site of the extremum may be distorted by a nearby larger feature in @xmath10 . running the simulation twice as long ( fig . [ pipreslater ] ) shows much closer alignment and is evidence that the predicted non - generic behavior will eventually dominate the spikes .    the effective potentials @xmath135 and @xmath144 are displayed in fig .  [ potentials ] . we see that the approach to the avtd regime is characterized by repeated alternate interactions of @xmath5 with the two potentials . this is displayed explicitly for a typical value of @xmath11 in fig . [ pv1v2 ] where @xmath5 , @xmath135 , and @xmath144 are shown as functions of @xmath12 . after each interaction with @xmath144 , the quantity @xmath175 changes sign . if this yields @xmath176 , no further bounces occur . if @xmath132 , there will be an interaction with @xmath135 . eventually , @xmath177 falls below unity so that @xmath144 permanently disappears . @xmath5 then continues to increase for all subsequent @xmath12 with constant slope @xmath178 as is characteristic of the avtd limit . ( note that @xmath135 will also decrease exponentially as @xmath5 increases even though it is part of the avtd solution . ) examination of ( [ qequations ] ) shows that , when @xmath177 falls below unity and @xmath39 , @xmath179 as required by the avtd limit . once @xmath99 becomes constant , as @xmath5 increases , @xmath180 , again as required . the development of @xmath10 with time is shown in fig . [ qv1 ] for the same value of @xmath11 as in the previous figure . here we see clearly how @xmath10 grows fastest when @xmath135 dominates . this becomes even more dramatic in fig . [ qminusp ] which is the same graph for @xmath181 , the site of the apparent discontinuity in @xmath10 . it is important to emphasize that the magnitude and time evolution of @xmath135 and @xmath144 as well as the initial values of the wave amplitude time derivatives are different at different values of @xmath11 . this explains the development of complicated small scale spatial structure . from ( [ k1eq ] ) or ( [ k2eq ] ) , one computes the time to the first bounce off @xmath135 or @xmath144 as @xmath182 where @xmath183 or @xmath184 and @xmath185 at @xmath186 . here @xmath187 and @xmath188 . evaluation of ( [ k1eq ] ) and ( [ k2eq ] ) at @xmath189 yields an approximate inverse scaling with @xmath102 seen in fig . [ scaling ] . the central peak appears first and is easy to understand . since , initially , @xmath190 , the highest negative velocity occurs at @xmath191 causing the bounce off @xmath135 to occur there first . @xmath5 then begins to increase at @xmath191 both earlier and faster than at neighboring values of @xmath11 . more complicated is the next pair of peaks which occur near @xmath192 . initially , the profile of @xmath5 vs @xmath11 steepens . however , bounces off @xmath144 will begin to occur . this causes the tangent to @xmath5 vs @xmath11 to evolve from negative ( positive ) to positive ( negative ) slope for @xmath193 ( @xmath194 ) eventually creating a pair of peaks . the higher the original value of @xmath102 , the greater will be the number of bounces before @xmath144 disappears for good . ever finer small scale structures are created at @xmath11 values where many bounces occur since @xmath5 will increase for some spatial regions and decrease for others in a manner that changes with time . as @xmath110 falls below unity , at a given @xmath11 , the bounces , and thus structure formation , will cease there but structure formation will continue in the ever smaller spatial regions where @xmath110 continues to exceed unity . we see that both potentials @xmath135 and @xmath144 drive the cosmology toward the avtd regime characterized by ( [ avtdeq ] ) as @xmath39 with @xmath195 . this is because @xmath135 appears and influences the dynamics only when and where @xmath132 while @xmath144 does so only for @xmath138 . in the avtd regime , with @xmath196 , both corresponding terms in the wave equation for @xmath5 are exponentially small . we note that the wave equations ( [ gowdywavep ] ) and ( [ gowdywaveq ] ) provide no other mechanism to drive @xmath112 into the required range if it is not there initially . recall that the polarized gowdy model can have any value for @xmath112 . thus we argue that at isolated points ( a set of measure zero in @xmath11 ) where either @xmath99 or @xmath197 is precisely zero , @xmath112 could remain outside the allowed range . in the former case ( say at @xmath198 ) , the consistency equations ( [ peqgm ] ) are satisfied so that one expects the solution at @xmath199 to be like an avtd polarized solution . however , in the latter case ( say at @xmath200 ) , the consistency conditions require @xmath201 as well as @xmath148something one does not expect generically . if @xmath201 at @xmath149 , then the solution is of the avtd large @xmath110 type there . if @xmath202 , the behavior can not be analyzed either numerically or by these analytic methods . we may therefore conclude that the simulations provide strong support for the conjecture that the gowdy models have a singularity which is avtd except , perhaps , at a set of measure zero . the existence of this set of measure zero of `` non - generic '' spatial points can be inferred from the simulations since near them the potentials are very flat . thus , the closer one is to e.g. @xmath199 or @xmath149 defined as above , the longer it will take to drive @xmath112 to the range @xmath203 there . this manifests itself as a resolution dependence of the simulations of a very specific type . at any value of @xmath12 ( say @xmath204 ) , there is a spatial resolution sufficient to resolve all spiky features . ( this was demonstrated in @xcite . ) but the `` center '' of each spiky feature is the site of a non - generic point . as one moves away from this point in either direction in @xmath11 , the coefficient of the potential will be larger . thus the interaction with the potential will occur first at regions away from the center of the feature and move closer to the center as @xmath12 increases . since the effect of the interaction is to change the sign of either @xmath112 or @xmath205 , the interactions will tend to eliminate the outer part of the spiky feature . thus the spiky feature gets narrower in time until for some @xmath206 , the original resolution is no longer adequate . although , in a peak near @xmath148 , @xmath112 should be @xmath207 , inadequate resolution may act as an effective averaging to give a measured smaller value of @xmath112 . thus the peak in @xmath110 appears to decrease in amplitude only because , at most @xmath11 values , the interaction with @xmath144 has already occurred . we note that this effective averaging will give @xmath48 everywhere for some @xmath208 for any fixed spatial resolution . as the spatial resolution is increased , @xmath208 will also increase since accurate resolution of spiky features can be maintained farther into the simulation and the effective averaging will take longer to come into play . thus the apparent resolution dependence is consistent with the behavior of the potentials @xmath135 and @xmath144 expected near the non - generic points . several examples are shown in figs . [ ptau2res][pneg ] . in fig . [ ptau2res ] , a spiky feature in @xmath112 is shown for three different spatial resolutions at @xmath209 . while the coarsest resolution shows @xmath159 , it also clearly fails to resolve the spiky feature . the highest resolution both resolves the spiky feature and shows that it has a core where @xmath40 as expected . a similar spiky feature in @xmath112 is shown in fig . [ threetaus ] for three different times for the same spatial resolution . one expects that the clearly resolved spike observed early in the simulation will narrow as @xmath12 increases remaining always @xmath207 at the center of the peak . this is not observed ( rather one sees the constant @xmath112 behavior expected in the avtd limit ) . however , one again sees clearly that the spatial resolution which was adequate at @xmath210 has become inadequate at later times and that the properties of the spike are not correctly represented . [ qnarrow ] shows that the apparent discontinuity in @xmath10 becomes narrower and steeper with time and eventually is not adequately resolved at a given spatial resolution . one expects @xmath5 to be increasingly negative at the center of this feature since @xmath132 is maintained at that non - generic point . at fixed spatial resolution , we see in fig . [ pneg ] that the negative region of @xmath5 both narrows and deepens as @xmath12 increases from @xmath211 to @xmath212 as expected . however , at @xmath213 , the expected behavior is not seen since the depth of the negative region fails to increase . however , we also see that the spatial resolution has become inadequate . presumably , adequate spatial resolution would show the narrowness and depth of the negative region suggested by the curvature of @xmath5 vs @xmath11 in the region just outside the core .    in conclusion , the behavior of generic gowdy cosmologies on @xmath0 can be completely understood in terms of the nonlinear interactions between the two polarizations of the gravitational waves . these act as effective potentials which drive the system to the prediced avtd behavior and then cease to play a role . at the set of measure zero where the potentials are absent , @xmath112 may lie outside its allowed avtd range of @xmath203 for all @xmath12 if it does so at any time since the mechanism to correct its value is absent . the existence of this set of measure zero is inferred from the details of the shape and time dependence of the spiky features and from the inability to resolve them with a given spatial resolution beyond a certain time . the authors would like to thank vincent moncrief and boro grubisi for useful discussions . bkb would like to thank the institute for geophysics and planetary physics at lawrence livermore national laboratory and the albert einstein institute at potsdam for hospitality . this work was supported in part by national science foundation grants phy9507313 and phy9722039 to oakland university . computations were performed at the national center for supercomputing applications ( university of illinois ) .
Abstract:

Answer: numerical studies of the plane symmetric , vacuum gowdy universe on @xmath0 yield strong support for the conjectured asymptotically velocity term dominated ( avtd ) behavior of its evolution toward the singularity except , perhaps , at isolated spatial points . a generic solution is characterized by spiky features and apparent `` discontinuities '' in the wave amplitudes . it is shown that the nonlinear terms in the wave equations drive the system generically to the `` small velocity '' avtd regime and that the spiky features are caused by the absence of these terms at isolated spatial points .    # 1tt@xmath1#1


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: entanglement , an essential feature of the quantum mechanics , has been introduced in many fields of physics . in the field of quantum information , the entanglement has played a key role . the study of entanglement properties in many - body systems have attracted much attention  @xcite-@xcite . the heisenberg chains , widely studied in the condensed matter field , display rich entanglement features and have many useful applications such as in the quantum state transfer  @xcite . most of the systems considered in previous studies are spin - half systems as there exists a good measure of entanglement of two spin - halves , the concurrence  @xcite , which is applicable to an arbitrary state of two spin halves . on the other hand , the entanglements in mixed - spin or higher spin systems are not well - studied due to the lack of good operational entanglement measures . there are several initial studies along this direction  @xcite , however these works are restricted to the case of two particles . for the case of higher spins , a non - entangled state has necessarily a positive partial transpose according to the peres - horodecki criterion  @xcite . in the case of two spin halves , and the case of ( 1/2,1 ) mixed spins , fortunately , a positive partial transpose is also sufficient . thus , the sufficient and necessary condition for entangled state in ( 1/2,1 ) mixed spin systems is that it has a negative partial transpose . this allows us to investigate entanglement features of the mixed spin system . the peres - horodecki criterion give a qualitative way for judging whether the state is entangled or not . the quantitative version of the criterion was developed by vidal and werner  @xcite . they presented a measure of entanglement called negativity that can be computed efficiently , and the negativity does not increase under local manipulations of the system . the negativity of a state @xmath0 is defined as @xmath1 where @xmath2 is the negative eigenvalue of @xmath3 , and @xmath4 denotes + the partial transpose with respect to the first system . the negativity @xmath5 is related to the trace norm of @xmath3 via @xmath6 where the trace norm of @xmath3 is equal to the sum of the absolute values of the eigenvalues of @xmath3 . in this paper , we will use the concept of negativity to study entanglement in ( 1/2,1 ) mixed - spin systems .    as shown in most previous works , models with the nn exchange interactions are considered and it is not easy to have pairwise entanglement between the nnn spins  @xcite . it is true that there exist some quasi - one - dimension compounds offering us systems with nnn interactions . bose and chattopadhyay  @xcite and gu et al . @xcite have investigated entanglement in spin - half heisenberg chain with nnn interactions . in our paper here , we study entanglement properties not only in the ( 1/2,1 ) mixed - spin systems only with nn interactions , but also in the system with nnn interactions . entanglement in a system with a few spins displays general features of entanglement with more spins . for instance , in the anisotropic heisenberg model with a large number of qubits , the pairwise entanglement shows a maximum at the isotropic point  @xcite . this feature was already shown in a small system with four or five qubits  @xcite . so , the study of small systems is meaningful in the study of entanglement as they may reflect general features of larger or macroscopic systems . also , due to the limitation of our computation capability , we only concentrate on small systems such as 4 , 5 and 6-spin models . the paper is organized as follows . in sec . ii , we study the systems with only nn interactions . the analytical results of negativity for the cases of two and three spins are given . the relation between entanglement and the macroscopic thermodynamical function , the internal energy is revealed . also we numerically compute the negativity in more general mixed - spin models up to eight spins , and consider the effects of magnetic fields in this section . in sec . iii , the system with nnn interaction is discussed . for the four - spin case , we analytically calculate the eigenenergy of the system from which we get the analytical results of the negativity of the nn spins . we numerically study negativities versus nnn exchanging coupling , and the case of finite temperature is also considered . for larger system up to eight - spin system , we get some numerical results . the conclusion is given in sec . we study entanglement of states of the system at thermal equilibrium described by the density operator @xmath7 , where @xmath8 , @xmath9 is the boltzmann s constant , which is assumed to be one throughout the paper , and @xmath10 is the partition function . the entanglement in the thermal state is called thermal entanglement . we consider two kinds of spins , spin @xmath11 and @xmath12 , alternating on a ring with antiferromagnetic exchange coupling . the hamiltonian is given by @xmath13 where @xmath14 and @xmath15 are spin-1/2 and spin-1 operators , respectively . the exchange interactions exist only between nearest neighbors , and they are of the same strength which are set to one . we adopt the periodic boundary condition . in fig .  1 , we give the schematic representation of the above hamiltonian . next , we first consider the models with two and three spins , and aim at getting analytical results of entanglement .      for the two - spin case , the hamiltonian  ( [ h01 ] ) reduces to @xmath16 . to have a matrix representation of the hamiltonian , we choose the following basis @xmath17 where @xmath18 is the eigenstate of @xmath19 and @xmath20 with the corresponding eigenvalues given by @xmath21 and @xmath22 , respectively .        in the above basis , the hamiltonian can be written as a block - diagonal form with the dimension of each block being at most @xmath23 . thus , the density matrix @xmath24 for the thermal state is obtained as @xmath25 with the partition function and the matrix elements given by @xmath26 after the partial transpose with respect to the first spin - half subsystem , we can get @xmath27 @xmath28 which is still of the block - diagonal form , and computation of its eigenvalues is straightforward . there are only two eigenvalues which are possibly negative . the negativity is thus given by @xmath29\nonumber\\ + & \frac{1}2\max\big[0,\sqrt{(a_5-a_6)^2 + 4b_1 ^ 2}-a_5-a_6\big].\label{n}\end{aligned}\ ] ]    substituting eqs . ( [ para1])([para2 ] ) leads to the analytical result of negativity @xmath30\nonumber\\ = & 2\max[0,a_2 - 2a_1]\nonumber\\ = & \frac{1}{3}\max\big[0,\frac{e^\beta-4e^{-\frac{1}{2}\beta } } { e^\beta+2e^{-\frac{1}{2}\beta}}\big],\label{n1}\end{aligned}\ ] ] where the second equality follows from eq . ( [ para2 ] ) . we can see that the negativity is a function of the single parameter @xmath31 . in the limit of @xmath32 , the negativity becomes 1/3 and the ground state is entangled . from eq . ( [ n1 ] ) , it is direct to check that the negativity is a monotonically decreasing function when temperature increases . after a certain threshold value of the temperature , the entanglement disappears . this threshold value @xmath33 can be obtained as @xmath34 for a ring of spin - half particles interacting via the hisenberg hamiltonian , it was shown that the pairwise thermal entanglement is determined by the internal energy  @xcite . it is natural to ask if similar relations exist in the present mixed - spin system . the internal energy can be obtained from the partition function as @xmath35 substituting eq . ( [ para1 ] ) into the above equation leads to @xmath36 from eqs . ( [ n1 ] ) and ( [ u ] ) , we obtain a quantitative relation between the negativity and the internal energy @xmath37.\label{u1}\ ] ] the above equation builds a connection between the microscopic entanglement and the macroscopic thermodynamical function , the internal energy . the internal energy completely determine the thermal entanglement . from the equation , we can also read that the thermal state becomes entangled if and only if the internal energy @xmath38 . since @xmath39 , we have @xmath40,\label{u2}\ ] ] which is consistent with the result obtained in ref . @xcite by the group - theoretical technique . we now consider the three - spin case and the schematic representation of the corresponding hamiltonian is given by fig .  1 . in this situation , there are two types of pairwise entanglement , the entanglement between spin 1/2 and spin 1 and the entanglement between two spin halves . the eigenvalue problem can be solved analytically , and after tracing out the third spin - half system the reduced density matrix @xmath24 is still of the same form as in eq . ( [ rho ] ) with matrix elements given by @xmath41 substituting the above equations to eq . ( [ n ] ) leads to the negativity @xmath42/z.\ ] ] it is evident that the negativity becomes 1/3 in the limit of @xmath32 . from the expression of the negativity , the threshold value of temperature after which the entanglement vanishes can be estimated as @xmath43    to examine the entanglement between two spin halves , we trace out the spin 1 system and get the reduced density matrix @xmath44 as follows @xmath45 with the matrix elements given by @xmath46 after taking the partial transpose , we can get @xmath47 @xmath48 then , the negativity is readily obtained as @xmath49.\ ] ] it is straightforward to check that the negativity is always zero . or , from another way , all the eigenvalues of the matrix @xmath47 are obtained as @xmath50 obviously the negativity vanishes here , in other words there is no entanglement between the two spin halves . the ground - state negativity @xmath51=1/3 and @xmath52=0 , here @xmath53 denotes the negativity between the @xmath54 and @xmath55 spin on the chain and @xmath56 denotes the one between the @xmath54 and @xmath57 spin . the equation above can be obtained from the non - degenerate ground state given by : @xmath58 it is interesting to see that the ground - state entanglement between the spin half and spin 1 in the three - spin case is the same as that in the two - spin case .    due to the su(2 ) symmetry in our system , there are following relations between correlation functions and negativities  @xcite @xmath59 where we have removed the max function in the negativity , implying that the negative value of @xmath5 indicates no entanglement . then , we have the relation between the internal energy and the negativities @xmath60 the second equality follows from the exchange symmetry , namely , the hamiltonian is invariant when exchanging two spin halves . so , for the three - spin case , the internal energy is related to two negativities .    to apply the above result , we consider the the ground - state properties ( @xmath61 ) . the hamiltonian can be rewritten as @xmath62.\ ] ] then , by the angular momentum coupling theory , the ground - state energy is obtained as @xmath63 . substituting the ground - state energy and @xmath64 to eq . ( [ relation ] ) , we obtain @xmath65 , indicating that there exists no entanglement between two spin halves . next , we consider more general situations , i.e. , the case of even @xmath66 sites .      except for the su(2 ) symmetry in the system , there exists exchange symmetry for the case of even spins . for instance , in the four - spin model , the hamiltonian is invariant when exchanging two spin halves or two spin ones . thus , for even - spin model , the entanglements between the two nearest - neighbor spins and the correlation functions @xmath67 are independent on index @xmath68 . therefore , the internal energy per spin is equal to the correlation function @xmath67 @xmath69 from eqs . ( [ n12 ] ) and ( [ uuu ] ) , we have @xmath70 this equation indicates that for the case of even spins the entanglement between two nearest neighbors is solely determined by the internal energy per spin . and for the case of zero temperature , the entanglement is determined by the ground - state energy . the less the energy , the more the entanglement . we now apply the above result to the study of ground - state entanglement . we rewrite the four - spin hamiltonian as follows @xmath71.\ ] ] then , by the angular momentum coupling theory , we obtain the ground - state energy per site @xmath72 . thus , from eq . ( [ n1212 ] ) , we have @xmath73 . for @xmath74 , it is hard to get analytical results . we next numerically calculate the entanglement for the case of more spins , and also consider the effects of magnetic fields . having obtained analytical results of entanglement in the heisenberg model with a few spins , we now numerically examine the entanglement behaviors in more general hamiltonian including more spins and magnetic fields . we consider the entanglement versus temperature for different number of spins @xmath66 , and the numerical results are plotted in fig .  2 and fig . it is clear to see that the ground state exhibits maximal entanglement , and with the increase of temperature , the entanglement monotonically decreases until it reaches zero . the decrease of entanglement is due to the mixture of less entangled excited states when increasing the temperature . the existence of the threshold temperature is also evident . for the case of even number of sites ( fig . 2 ) , the threshold temperature decreases with the increase of @xmath66 . in contrast to this behavior , for the case of odd number of spins , as seen from fig .  3 , the threshold temperature increases with the increase of @xmath66 . we now examine the effect of magnetic fields on entanglement . the heisenberg hamiltonian with a magnetic field along @xmath75 direction is given by @xmath76 for @xmath77 , the analytical result of negativity can be obtained ( see appendix a ) . .,scaledwidth=45.0% ]    in fig .  4 , we plot the negativity versus the magnetic field at a low temperature . for the two - spin case , with the increase of the magnetic field , the negativity rapidly reach a platform , and then after a certain magnetic field @xmath78 , it jumps down to zero , indicating no entanglement . when @xmath79 , the ground state is two - fold degenerate with @xmath80 . for a small @xmath81 , the system is no longer degenerate , and the ground state is given by @xmath82 which is of the schmidt form . for a state written as its schmidt form @xmath83 the negativity is obtained as  @xcite @xmath84/2.\ ] ] here , @xmath85 are the schmidt coefficients , and @xmath86 and @xmath87 are bases for subsystems 1 and 2 , respectively . applying this formula to the ground state , we immediately have @xmath88 . when @xmath89 , the ground state is two - fold degenerate , and when @xmath90 the ground state becomes non - degenerate and the corresponding wave function is given by @xmath91 which is obviously of no entanglement . then , the ground - state negativity forms a platform when @xmath92 . the jump of negativity at @xmath93 is due to the level crossing . for @xmath94 , the effects of magnetic fields on entanglement can be also explained by level crossing . for instance , for @xmath95 , there are two level crossing , and the entanglement displays two jumps . we have studied the effects of finite temperature and magnetic fields on entanglement , and now consider the model containing two kinds of spins , spin @xmath11 and @xmath12 , alternating on a ring with antiferromagnetic exchange coupling between both the nn spins and the nnn spins . the hamiltonian can be expressed as @xmath96 where the @xmath14 and @xmath15 are spin-1/2 and spin-1 operators in the @xmath68th cell . @xmath97 characterizes the nn exchange coupling and @xmath98 the nnn coupling . we consider the antiferromagnetic interaction by taking @xmath99 . @xmath66 is the total number of spins and here we choose it be even . also we adopt the periodic boundary condition . the model with four spins is the simplest model with nnn interactions . we first solve the eigenvalue problem of this model . the key step is to write the four - spin hamiltonian in the following form , @xmath100\nonumber\\ & -&2j_2({\bf s}_1 ^ 2+{\bf s}_2 ^ 2+{\bf s}_1 ^ 2+{\bf s}_2 ^ 2 ) \big\},\end{aligned}\ ] ] where @xmath101 denotes the total spin . from the above form , by angular momentum coupling theory , one can readily obtain all eigenvalues of the system as @xmath102 where parameter @xmath103 in expressions @xmath104 and @xmath105 and @xmath106 in expressions @xmath107 and @xmath108 , respectively . then , from eq . ( [ e ] ) , we may find the ground - state energy as @xmath109 clearly , there are two level - crossing points , which will greatly affect behaviors of ground - state entanglement .    from eq . ( [ h4 ] ) , it is obvious that in addition to the the su@xmath110 symmetry , there also exists an exchange symmetry , namely , exchanging two spin halves or two spin ones yields invariant hamiltonian . then , the correlator @xmath111 between any nn spins are the same . therefore , we can get the correlator @xmath112 from the ground - state energy via the hellmann - feynman theorem .    when @xmath113 , after applying the hellmann - feynman theorem to the ground - state energy , we obtain @xmath114 due to the su@xmath110 symmetry in our system , we have the following relation between the negativity and the correlator @xmath115@xcite , @xmath116 where we use @xmath117 to denote the negativity between the nn spins . thus we obtain @xmath118 in other regions , we find that @xmath119 for @xmath120 .    from the above analytical results , we can find that the negativity between the nn spins is not a continuous function of parameter @xmath98 . it jumps from the value @xmath121 down to zero at @xmath122 . hence , we can see that the nnn interaction may deteriorate the entanglement of the nn spins . the critical point of @xmath122 may be considered to be a threshold value , after which the nn entanglement vanishes . having studied ground - state entanglement , now we make a short discussion of entanglement of excited states . we consider the first excited state . in the region @xmath123 , the first excited energy is given by @xmath124 from that we can get @xmath125 and @xmath126 . while in the rest region the negativity of the first excited state is zero .      next , we consider the thermal entanglement . from eq . ( [ e ] ) , the partition function can be obtained as @xmath127 the correlator @xmath128 at finite temperature can be computed from the partition function via the following relation @xmath129    substituting ( [ pf ] ) to eq . ( [ r1 ] ) yields @xmath130.\end{aligned}\ ] ] after substituting the above equation into ( [ cor ] ) , we may get analytical expression of the negativity @xmath117 at finite temperatures . the negativity is a function of @xmath97 , @xmath98 and @xmath131 .    _ low - temperature case _ :     at the temperature @xmath132 in the four - spin system . the ground - state and the first two excited - state energy levels versus @xmath98 are inserted.,scaledwidth=45.0% ]    we now make numerical calculations of entanglement and first consider the low - temperature case . we take @xmath133 in all the following plots . in our system there exist three kinds of negativity , the negativity @xmath117 between spin-1 and spin-1/2 , @xmath134 between two spin-1/2 , and @xmath135 between two spin-1 .    in fig . 5 , we plot the negativity versus @xmath98 in four - spin system at a low temperature of @xmath132 . it is clear to see that @xmath136 keeps a value about @xmath137 when @xmath98 increases from zero until it reaches the critical point , at which the @xmath136 displays a jump to zero . this behavior of entanglement is consistent with that at zero temperature from the analytical results . it is natural to see that increase of nnn exchange interaction will suppress the entanglement of nn spins , and at last completely erase the entanglement .    in comparison with @xmath117 , the negativities @xmath135 and @xmath134 behave distinctly . we see that near the point of @xmath138 , @xmath139 increases quickly to a steady value about @xmath140 , and when @xmath98 reaches the value about @xmath141 , @xmath139 jumps another steady value 1 . these two jumps result from the two level crossing as seen clearly from the figure inserted . the second level crossing also leads to a small dip in the curve of @xmath139 . the negativity @xmath142 displays a jump to a steady value of 1/2 near @xmath143 . it is evident that the entanglement between nnn spins is enhanced by increasing nnn interactions . the competition between nn and nnn interactions leads to rich behaviors of quantum entanglement . another observation is that there is a range of @xmath98 , at which negativities @xmath117 and @xmath142 are zero , and only @xmath139 is not zero . this indicates that the nnn interaction must be strong enough to build up the entanglement of two spin halves .    _ entanglement versus @xmath98 and @xmath131 _ : as temperature increases the entanglement will decrease due to the mixing of less entangled excited states to the thermal state . it is obvious that there exists a threshold temperature after which the negativity is zero . in the frustrated system , there exists the parameter @xmath98 , and with its increase , the negativity @xmath117 will decrease to zero , while @xmath139 and @xmath134 increase from zero to their maxima . so it is clear that there also exists a threshold @xmath144 corresponding to the boundary between zero and nonzero negativities .     versus both the exchange interaction coefficient @xmath98 and the temperature in the four - spin system.,scaledwidth=45.0% ]    in fig .  6 , we show the negativity @xmath117 versus the temperature and @xmath98 . when the temperature approaches zero , the @xmath117 reaches its maximum , and then with the temperature increasing , @xmath117 decreases to zero . on the @xmath145 plane , there is a curve along which the negativity just turns to be zero . it is possible to consider that the curve describes the threshold @xmath146 versus the temperature . obviously , @xmath146 does not behave as a monotonous function of the temperature , and it displays a peak at about @xmath147 , this behavior is in big contrast with the case of non - mixed qubit systems  @xcite . when the temperature rises , the weight of excited states will increase and it may strongly affect the negativity . this behavior of @xmath146 results from both the mixture of excited states to the thermal state and the intrinsic properties of the mixed - spin system . in addition , after crossing the temperature about @xmath148 , @xmath117 will vanish , irrespective of the value of @xmath98 .    from another point of view , we can read the threshold temperature @xmath149 versus different @xmath98 from the curve in the @xmath145 plane . when @xmath98 increases , @xmath149 decreases , and when @xmath98 crosses about @xmath150 , @xmath117 will disappear at any temperature . next , we consider the entanglement between nnn spins . in fig .  7 , we plot the negativity @xmath134 as a function of the temperature and @xmath98 . we can see that , before @xmath98 reaches the value about @xmath151 , @xmath134 keeps being zero at any temperature . and in the region @xmath152 , the @xmath142 can be enhanced by the increasing nnn interaction . this is a result from the competition of two kinds of exchange interactions . the thermal fluctuation all along suppresses the entanglement . so , from the curve lying on the @xmath145 plane which corresponds to the boundary of the nonzero and zero values of @xmath134 , we may find that the higher the temperature is , the larger the threshold @xmath146 will be . from another point of view , the @xmath149 increases as @xmath98 increases . versus both the exchange interaction coefficient @xmath98 and the temperature in four - spin system.,scaledwidth=45.0% ]    in fig .  8 , we plot the negativity @xmath139 versus @xmath131 and @xmath98 . in the region of @xmath153 , @xmath139 is zero at any temperature . when @xmath154 , the increasing nnn exchange interaction @xmath98 enhances the negativity and exhibits two particular flat roofs . with the temperature rises , @xmath135 is suppressed to zero . also we can consider the threshold @xmath149 and @xmath146 from the critical curve on the @xmath145 plane , and the @xmath146 also behaves as an increasing function of the temperature . we can see the nonzero region of @xmath139 is much larger than @xmath142 . but here it should be pointed out that , because @xmath155 only gives a sufficient condition for entangled state , we can not definitely say that the state in the area of zero negativity is not entangled .     versus both the exchange interaction coefficient @xmath98 and the temperature in four - spin system.,scaledwidth=45.0% ]      in this section , we present numerical results of negativity for more spins , and first consider the low - temperature case . at the temperature @xmath156 in six - spin system . the ground and the first two excited energy levels versus the @xmath98 are inserted in the plot.,scaledwidth=45.0% ]     at the temperature @xmath156 in eight - spin system , the ground and the first three excited energy levels versus the @xmath98 are inserted.,scaledwidth=45.0% ]    in fig . 9 , we plot the negativity versus @xmath98 for the case of six spins at a lower temperature . the negativity @xmath136 behaves as a decreasing function of @xmath98 . it decreases to zero at about @xmath157 , which is the special point corresponding to the energy level crossing . on the contrary , around the point @xmath157 , @xmath139 jumps up to a nonzero value , and then increases gradually until approaching the limit about @xmath158 . this behavior is quite different from that in the four - spin model . we also see that @xmath134 is zero all the time . it can be understood as follows . in the six - site system there are three spin halves with the nnn interaction . even for a pure homogeneous three - qubit system , there is no entanglement between two spins , irrespective of the strength of the exchange interactions  @xcite . now , in addition to the interaction among three spin halves , there are also interaction between spin halves and spin ones . so , it is reasonable that the entanglement between two spin halves is zero .    the negativity versus @xmath98 in the eight - site case is shown in fig . after the first sharp jump to a value ( not zero ) at about @xmath122 , @xmath117 goes down to zero gradually . at approximately @xmath159 , the negativity is zero . on the contrary , the negativity @xmath139 jumps up at about @xmath122 , and then goes up gradually and almost linearly until @xmath98 reaches about @xmath160 . then there happens a sharp decrease to nearly zero , and after that it begins to increase gradually . the negativity @xmath134 keeps zero until @xmath98 reaches about @xmath161 , and then it goes up until reaches a steady value .    in fig . 10 , we find that the three kinds of negativity exhibit different properties . from figure inserted , i.e. , the energy levels of the eight - spin system , we can see that in the region from about @xmath159 to @xmath162 the first excited energy is quite close to the ground energy . it is known that the energy level crossing can greatly affect the entanglement . here , the two close energy levels also play an important role in the behavior of entanglement . the approximate degenerate energy levels may remarkably change the probability distribution even at a very low temperature , thereby affect the negativity .       versus both the exchange interaction @xmath98 and the temperature in the six - spin system.,scaledwidth=45.0% ]     versus both the exchange interaction @xmath98 and the temperature in six - spin system.,scaledwidth=45.0% ]    now , we present the entanglement versus @xmath98 and @xmath131 in the six - spin model . the nn negativity @xmath117 is shown in fig . 11 . on the @xmath145 plane , similar to the four - spin case , the @xmath146 does not behave as a monotonous function of @xmath131 and it reaches its maximum at about @xmath163 . from the figure , we may find that the entanglement only exists in the region approximately @xmath164 and @xmath165 . the strong nnn interaction and thermal fluctuation will suppress the nn entanglement to zero .    we do not plot the negativity @xmath134 as it is zero all along for any @xmath98 and @xmath131 . we give the nnn negativity @xmath135 in fig . 12 . in the region of @xmath166 , there is no negativity @xmath139 at any temperature . the threshold @xmath146 is increased by the increasing temperature , similar to the four - spin case . in conclusion , we have studied the entanglement properties of the ( 1/2,1 ) mixed - spin systems described by the heisenberg model . in the systems only with nn exchange interactions , for two - spin and three - spin cases analytical results of the negativity have been obtained , which facilitate our discussions of entanglement . the analytical expression of threshold temperature after which the entanglement vanishes are obtained for the two - site case . for the case of even number of particles , it has been found that the pairwise thermal entanglement is solely determined by the internal energy , and thus builds an interesting relation between the microscopic quantity , entanglement , and the macroscopic thermal dynamical function , the internal energy in the mixed - spin systems . for the odd number of particles such as the three - site case , we also provide a relation between the internal energy and negativities . we have numerically studied the effects of different finite temperature and magnetic fields on entanglement . as a conclusion , the thermal fluctuation suppresses the entanglement , and entanglement may change evidently at some critical points of magnetic field .    in the systems also with nnn interactions , by applying the angular momentum coupling method , we obtained analytical results of all the eigenenergies of the four - spin system , based on which the negativity of the nn spins has been obtained . we also considered the excited - state entanglement . at finite temperature , from the partition function , the analytical results of negativity has been given . we have numerically studied the effects of the nnn interaction on the nn entanglement and nnn entanglement . it is natural to see that the nnn interaction suppresses the nn entanglement , and enhances the nnn entanglement . we found that the negativity between two spin ones is sensitive to @xmath98 and displays some interesting properties . at finite temperature , the thermal fluctuation suppresses both the nn entanglement and the nnn entanglement . the threshold values @xmath146 and @xmath149 are studied in detail . the entanglement displays some peculiar properties , which are quite different from those of the spin - half model . these are due to inherent mixed - spin character of our system . it is more interesting to study entanglements in other mixed systems and explore some universal properties , which are under consideration . the hamiltonian of the two - site heisenberg model with a magnetic field is written explicitly as @xmath167 following the same way as the discussions of subsection ii.a , the density matrix of the thermal state is given by eq . ( [ rho ] ) with the matrix elements @xmath168 and the partition function @xmath169 having obtained the analytical expressions of the matrix , we directly obtain the negativity after substituting the matrix elements to eq . ( [ n ] ) . 99 m. a. nielsen , ph . d thesis , university of mexico , 1998 , quant - ph/0011036 ; m. c. arnesen , s. bose , and v. vedral , phys . lett . * 87 * , 017901 ( 2001 ) . d. gunlycke , v. m. kendon , v. vedral , and s. bose , , 042302 ( 2001 ) . x. wang , phys . a * 64 * , 012313 ( 2001 ) . x. wang , h. fu , and a. i. solomon , j. phys . a : math . gen . * 34 * , 11307(2001 ) ; x. wang and k. mlmer , eur . j. d * 18 * , 385(2002 ) . g. l. kamta and a. f. starace , phys . lett . * 88 * , 107901 ( 2002 ) . k. m. oconnor and w. k. wootters , , 0520302 ( 2001 ) . d. a. meyer and n. r. wallach , j. math . phys . * 43 * , 4273 ( 2002 ) . t. j. osborne and m. a. nielsen , , 032110 ( 2002 ) . a. osterloh , l. amico , g. falci and r. fazio , nature * 416 * , 608 ( 2002 ) . y. sun , y. g. chen , and h. chen , phys . a * 68 * , 044301 ( 2003 ) . l. f. santos , , 062306 ( 2003 ) . y. yeo , , 062312 ( 2002 ) . d. v. khveshchenko , , 193307 ( 2003 ) . l. zhou , h. s. song , y. q. guo , and c. li , , 024301 ( 2003 ) . g. k. brennen , s. s. bullock , , 52303 ( 2004 ) . r. xin , z. song , and c. p. sun , quant - ph/0411177 . g. tth , quant - ph/0406061 . f. verstraete , m. popp , and j. i. cirac , phys . rev . lett . * 92 * , 027901 ( 2004 ) . f. verstraete , m. a. martin - delgado , j. i. cirac , phys . rev . lett . * 92 * , 087201 ( 2004 ) . j. vidal , g. palacios , and r. mosseri , phys . a * 69 * , 022107 ( 2004 ) . n. lambert , c. emary , and t. brandes , phys . lett . * 92 * , 073602 ( 2004 ) . s. j. gu , s. s. deng , y. q. li , and h. q. lin , , 086402 ( 2004 ) . y. chen , p. zanardi , z. d. wang , and f. c. zhang , quant - ph/0407228 . s. b. li and j. b. xu , phys . a * 311 * , 313 ( 2003 ) . k. audenaert , j. eisert , m. b. plenio , and r. f. werner , , 042327 ( 2002 ) . l. a. wu , m. s. sarandy , and d. a. lidar , quant - ph/0407056 .    g. vidal , j. i. latorre , e. rico , and a. kitaev phys lett . * 90 * , 227902 ( 2003 ) . s. bose , , 207901 ( 2003 ) ; v. subrahmanyam , , 034304 ( 2004 ) ; m. christandl , n. datta , a. ekert , and a. j. landahl , , 187902 ( 2004 ) ; y. li , t. shi , z. song , and c. p. sun , quant - ph/0406159 ; m. b. plenio and f. l semio , quant - ph/0407034 . w. k. wootters , , 2245 ( 1998 ) . j. schliemann , , 012309 ( 2003 ) . l. zhou , x. x. yi , h. s. song , and y. q. guo , quant - ph/0310169 . y. q. li and g. q. zhu , submitted . a. peres , , 1413 ( 1996 ) ; m. horodecki , p. horodecki , and r. horodecki , phys . a * 223 * , 1 ( 1996 ) . g. vidal and r. f. werner , , 032314 ( 2002 ) . i. bose and e. chattopadhyay , phys . a * 66 * , 062320 ( 2002 ) . s. j. gu , h. b. li , y. q. li , and h. q. lin , phys . rev . a * 70 * , 052302 ( 2004 ) . x. wang , phys . a * 329 * , 439 ( 2004 ) . x. wang and p. zanardi , phys . a * 301 * , 1 ( 2002 ) .
Abstract:

Answer: by using the concept of negativity , we investigate entanglement in ( 1/2,1 ) mixed - spin heisenberg systems . we obtain the analytical results of entanglement in small isotropic heisenberg clusters with only nearest - neighbor ( nn ) interactions up to four spins and in the four - spin heisenberg model with both nn and next - nearest - neighbor ( nnn ) interactions . for more spins , we numerically study effects of temperature , magnetic fields , and nnn interactions on entanglement . we study in detail the threshold value of the temperature , after which the negativity vanishes .