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: the phonon spectrum of defect - free crystals is well understood , since the advent of modern solid - state physics in the mid-20th century  @xcite . at low frequency and long wavelength , the linear dispersion relation between frequency and momentum results from the breaking of translational symmetry due to the periodic lattice , a manifestation of the goldstone theorem , and gives rise to the @xmath1 debye law in the density of states ( dos ) , in 3d . at higher frequencies , phonon propagation through brillouin - zone boundaries may appear as sharp peaks in @xmath2 , known as van hove singularities  @xcite . in the presence of structural disorder , the spectrum of vibrational modes presents very different features which remain poorly understood . the most striking anomaly in glasses is the deviation from the debye law which manifests itself as the well - documented excess of low - frequency modes visible as a peak in the the normalized dos @xmath3 . this effect is widely known as the boson peak anomaly , and is a universal feature in glasses  @xcite , although it has often been observed in crystals as well  @xcite . the ioffe - regel crossover  @xcite defines the frequency @xmath4 at which the phonon mean - free path becomes equal to its wavelength . very close to this frequency , is the crossover frequency @xmath5 from ballistic phonons in the linear regime @xmath6 to quasi - localized modes characterised by diffusive propagation . this crossover is supposed to play an important role for the boson peak in glasses , where local disorder gives rise to scattering at sufficiently small wave - vector @xmath7 , as well as in defective crystals where vacancies and interstitials act as local scattering centres  @xcite . this effect has been attributed , among other mechanisms , to the lowest van hove singularity in the spectrum of the reference crystalline system , shifted to lower frequencies by disorder - induced level - repulsion effects  @xcite . however , no clear or unifying understanding of the role of local structure has emerged for the boson peak in glasses and defective crystals using standard tools such as e.g. the bond - orientational order  @xcite , to characterize the effect of structural disorder . intriguingly , a recent experimental study has shown that the low-@xmath8 peak in the dos is very similar for the silica glass and the @xmath9-quartz crystal with matched density  @xcite . it was also found that the low-@xmath8 peak in the reduced dos of @xmath9-quartz could be accurately reproduced by dft calculations with just standard phonon physics and interpreted as the lowest van hove singularity , as shown also in ref . @xcite .    here we present numerical results for both model glasses and defective crystals with randomly - cut bonds , and a new conceptual framework to explain the boson peak based on a unifying local symmetry principle . we identify the micro - structural key - factor , which controls the boson peak in both these models of glasses and crystals , with the _ local _ center - inversion symmetry measured not with respect to the center of the unit cell , but , crucially , with respect to any atom taken as a local center of symmetry . given the small number of physical parameters in our model , we discuss each of them carefully in terms of whether they correlate or not with the boson peak . the local breaking of inversion symmetry appears to be the only microscopic structural signature that correlates with the emergence of the boson peak in both the glass and the defective crystal , without leading to contradictions or paradoxes . we show that model glasses and defective crystals having the same average atomic connectivity @xmath10 , as well as the same density and interatomic interaction , display the same boson peak in spite of having very different values of bond - orientational order . the proposed framework thus naturally provides a unifying framework to explain the boson peak in glasses and defective crystals . over the last decades , the vibrational density of states of glasses has been intensely studied from the point of view of theory and simulations , and it is impossible to give a full account of all models . here we limit ourselves to some of the models that have been applied to experimental data and which connect to our main interest . among the early models , thorpe  @xcite used a scalar elastic potential for central - force components of interatomic interaction supplemented with bond - bending terms to achieve a description of the dos of amorphous silicon . subsequently , a model was proposed which is known as the soft - potential model  @xcite , and is largely based on anharmonicity . within this framework , the anharmonic part of the soft potential causes a redistribution of local oscillator frequencies with a @xmath11 scaling in the dos which then crosses over into a linear @xmath12 scaling . the crossover gives rise to a boson peak in the dos . another model developed and extensively studied since the 90 s is based on the concept of shifted and smeared van hove singularity . in crystals , van hove singularities are discontinuities in the slope of the dos which occur when the dispersion curve is flat , i.e. when @xmath13 at the brillouin - zone boundaries . in ref . @xcite , a scalar model of a crystal with disorder in the spring constants was studied and the disorder was gradually increased to the point that the van hove singularity appeared shifted to much lower frequency and much smeared . since scalar models can not distinguish between longitudinal and transverse modes , and bear little resemblance to real materials , this model was refined for a vector model by taraskin et al . @xcite . here it became clear that it is the transverse van hove singularity which gets shifted to low frequency thus producing a close resemblance with the boson peak observed experimentally in glasses . subsequent calculations by zorn  @xcite suggest that the eigenvalues corresponding to the lowered transverse van hove singularity are presumably those which contribute the most to the boson peak . in these models , it is clear that the local atomic packing and structure play a major role in transforming the eigenvalue statistics and distribution of a perfect crystal into the one typical of glasses and defective crystals . one of the most popular models of the boson peak was proposed by schirmacher  @xcite , and is based on the concept of heterogeneous elasticity . the starting point is the assumption that the shear modulus is a spatially varying quantity due to structural disorder . the microscopic details of atomic arrangements are not specified and the model is entirely a macroscopic one . yet , the model is very powerful because the assumption of gaussian quenched disorder allows one to employ the standard tools of statistical field theory to arrive at elegant expressions for both longitudinal and transverse elastic green functions that are used to calculate the dos . the assumption of gaussian disorder has been later generalized to non - gaussian disorder within the scheme of coherent potential approximation  @xcite . overall , this model provides a picture of the boson peak as a consequence of the crossover from phonon physics ( at very large wavelengths ) into random - matrix physics at length - scales where the effect of disorder becomes important . more recently , important studies , both experimental and computational , have focused on the comparison between crystal polymorphs and glassy phases with the same atomic composition , such as e.g. vitreous silica and @xmath9-quartz , at matched densities  @xcite . it was found that the two systems have a very similar reduced dos in the boson peak region , and that the boson peak of the glass appears as a smoothened version of the van hove peak in the crystal . the latter peak in the crystal was shown to be as a lowered and smeared van hove singularity , as clarified by first - principles numerical calculations  @xcite . also the specific heat for the two systems was found to be identical . while most theoretical models are limited in their ability to describe the microscopic disorder , numerical simulations have been used in an attempt to identify the microscopic signature of structural disorder which is directly responsible for the boson peak . the most common parameter used to this aim is the bond - orientational order parameter  @xcite , which quantifies the spread in the angular orientations of the bonds connecting the atoms in the lattice . our main contribution below is to show that in model systems , which share many features of real systems in terms of dos and elastic properties , the bond - orientational order parameter does not correlate with the boson peak . instead , the breaking of local inversion symmetry displays a much stronger correlation and proves to be a good candidate for a microscopic structural signature of disorder which could be used as a universal order parameter to link the low - frequency non - debye peak in the dos to the underlying atomic structure . this microscopic interpretation is proved here for random network models of glasses and for defective crystals . it remains to see in future studies if this concept can prove useful also for defect - free non - centrosymmetric crystals such as @xmath9-quartz  @xcite . finally , the boson peak in glasses has been found to be dominated by transverse modes and to correlate strongly with softening of the shear modulus  @xcite . we will show below that both the boson peak and the softnening of the shear modulus can be linked to the phenomenon of nonaffine displacements which is caused by the breaking of local inversion symmetry . in our simulations , we use a random network created by first randomly placing @xmath14 soft spheres in a box and letting them interact via a truncated lennard - jones ( lj ) potential @xmath15 . the system is brought to a metastable lower energy state by a monte carlo energy - relaxation algorithm  @xcite . bonds are formed only between nearest neighbours and the bond length is distributed around the mean value @xmath16 . the volume of the box is chosen such to create a dense network with an average coordination number @xmath17 which is almost delta - distributed . the fact that the coordination is the same for all atoms implies the absence of regions which are locally under - coordinated or over - coordinated ( with respect to the average @xmath10 ) , hence the local rigidity is uniform throughout the sample  @xcite . to simulate systems with lower @xmath10 , we randomly cut bonds from the initial configuration , while keeping a narrow distribution of @xmath10 . we studied systems with coordination numbers from @xmath18 down to @xmath19 . the density is kept at a constant value @xmath20 and we implemented periodic boundary conditions to avoid surface effects . to reduce noise we calculate our results for ten independent realizations , over which we then take averages .    . the dots represent atoms which are connected by harmonic springs . local inversion - symmetry is evidently broken by the randomness and lack of correlation in bond orientations . ( b ) schematic 2d picture of a regular lattice where locally the removal of a bond breaks the inversion - symmetry on atom @xmath21 . the main consequence of inversion - symmetry breaking induced by the cutting of the bond is the imbalance of nn forces ( arrows ) acting on atom @xmath21 when it reaches its affine position under strain . the net force acting on atom @xmath21 in its affine position has to be released through an additional _ nonaffine _ displacement . the atom @xmath21 also acts as a scattering and quasi - localization center  @xcite for incoming vibrational excitations . ]    in the final step to create our model glass , we then used the so obtained configurations to generate harmonic random - network ( rn ) model glasses , where the lennard - jones interactions between nearest neighbours are all replaced by harmonic springs with pair potential @xmath22 , with spring constant @xmath23 . a sample realization of the model rn glass is shown in fig.1 for the marginally stable ( isostatic ) network @xmath24 . we also generated harmonic fcc crystals with the same density and spring constant as the rn glass , and randomly cut bonds to vary @xmath10 and to induce the breaking of inversion - symmetry . this procedure allows us to use all the tools of lattice dynamics and nonaffine linear response theory to analyse the dos and the shear modulus . the lattice dynamics is governed by the equation of motion for the displacement @xmath25 of atom @xmath21 , @xmath26 , with oscillating solutions @xmath27 , leading to @xmath28 . here @xmath29 denotes the unit vector pointing from atom @xmath21 to atom @xmath30 . using the latter relation , the time - independent part of the displacement is related to the dynamical ( hessian ) matrix  @xcite @xmath31 , where @xmath32 and its eigenvalues @xmath33 , via @xmath34 . hence , @xmath35 , upon recalling that the atomic mass is @xmath36 . in this way , the phonon density of states @xmath2 is obtained from the diagonalization of the hessian matrix , from which one obtains the set of eigenvalues @xmath33 . different eigenvalues are obtained from different realizations of the same sample , and averaging over the realizations leads to the distribution @xmath37 . the dos is thus calculated for different values of connectivity @xmath10 , for both the model rn glasses and the fcc crystals with randomly - cut bonds , which allows us to vary @xmath10 by keeping the density constant . also , bonds are severed to always keep a very narrow distribution of @xmath10 in all the samples , which ensures that spatial fluctuations of the affine part of the elastic constants is negligible . the hessian matrix is also a key quantity to evaluate the _ nonaffine _ elastic response of disordered solids . the latter is closely connected with the local inversion symmetry of the lattice  @xcite . in glasses , when applying shear stress to the solid , the atoms tend to reach a new position ( _ affine _ position ) proportional to the applied shear strain @xmath38 . in the affine position , the forces transmitted to any atom @xmath21 by its nearest neighbours ( nn ) cancel each other out only if atom @xmath21 is a local center of symmetry . if the atom is not a center of symmetry for the nn bonds , as schematically depicted in fig.1(b ) , the nn forces ( arrows in fig.1b ) can not cancel each other out and a net force acting on atom @xmath21 in the affine position has to be released via an additional _ nonaffine _ displacement . this is always true for glasses ( fig.1a ) , but also for crystal lattices with defects or with randomly - cut bonds ( fig.1b ) , and also for intrinsically non - centrosymmetric crystals like e.g. quartz  @xcite . in the latter systems , however , the covalent character of interatomic bonding , with its non - central component of interaction , makes the applicability of analytical theories of nonaffine lattice dynamics not yet established due to the difficulty of analytically evaluating hessian matrices with non - central interactions . in the harmonic approximation , the total nn force acting on @xmath21 under a strain @xmath38 can be expressed as @xmath39 , where   @xmath40 ( see ref . @xcite ) . the vector @xmath41 plays a very important role because it encodes the local inversion - symmetry of the lattice . as one can easily verify , @xmath42 if atom @xmath21 is a local center of symmetry , while @xmath43 if the lattice does not have local inversion - symmetry . hence , in non - centrosymmetric and disordered lattices , a net total force @xmath44 acts on any atom @xmath21 in its affine position . under the action of this force , the atoms have to perform an additional nonaffine displacement into their final _ nonaffine _ equilibrium positions , which is an internal work contributing negatively to the free energy of deformation , @xmath45 . the first term , @xmath46 , is the contribution from the affine displacements , which is the sum of all the bond - stretching energies , and can be calculated based on the born - huang lattice dynamics . the second term , @xmath47 , contains the reduction of the elastic free energy due to the nonaffine relaxation of the system caused by the lack of local inversion - symmetry . recalling that the shear modulus is given by @xmath48 , the local inversion - symmetry breaking thus causes the shear modulus of disordered solids to be lower compared to defect - free centrosymmetric crystals , and as shown in several previous works  @xcite : @xmath49    here the second , nonaffine ( negative ) contribution to the shear modulus @xmath50 is identically zero only for defect - free centrosymmetric crystal lattices . next we shall use this formalism to evaluate the shear modulus for our model glasses and randomly - cut fcc crystals as a function of the atomic connectivity @xmath10 . the results are shown in fig.2 . both our model systems follow the well known @xmath51 scaling with respect to the isostatic point @xmath24 found in many previous works . while it is well established  @xcite that @xmath52 , we note here , importantly , that the nonaffine contribution @xmath53 also depends on @xmath10 , and , in particular , it decreases with increasing @xmath10 . in fact , while the affine contribution @xmath54 for the glass is in nearly exact quantitative agreement with analytical mean - field predictions for random isotropic networks  @xcite , the nonaffine contribution decreases linearly upon increasing @xmath10 , which deviates from the mean - field theory  @xcite . from fig.2 we find that the following law is obeyed : @xmath55 for both the rn glass and the defective fcc crystal , where @xmath56 and @xmath57 are numerical coefficients . for the fcc crystal the nonaffinity vanishes in the limit of the perfect crystal with @xmath58 , and thus @xmath59 . overall , the nonaffinity decreases with increasing @xmath10 in qualitatively the same way for both rn lattice and defective crystal , which suggests a common microscopic structural origin for this behaviour , as discussed below . we shall now consider the density of states of both rn glass and defective fcc crystal , for the same conditions investigated for the shear modulus above . the results are shown in fig.3 . at large @xmath10-values we observe that , at the lowest frequencies , the parabolic debye law @xmath1 is visible , for both glass and crystal . the only difference in the spectrum is at higher frequencies where two peaks emerge in the fcc spectrum which are reminiscent of the typical peaks in the phonon spectrum of perfect fcc crystals  @xcite ; the latter spectrum is eventually recovered at @xmath58 , which we checked . at lower @xmath10 , where breaking of local inversion symmetry becomes important , the debye regime shrinks and the boson peak becomes more prominent . both spectra are quite similar to those of harmonic , stress - free random packings  @xcite . we also verified that the boson peak frequency scales with connectivity as @xmath60 . the latter scaling can be explained in terms of the crossover between the elastic - continuum regime ( ballistic phonon propagation ) and the quasi - localized random - matrix - dominated regime ( diffusive - like propagation ) , as suggested in ref.@xcite . calculated for the rn glass ( solid line ) and for the randomly - cut fcc crystal ( dotted line ) , for @xmath61 different values of atomic connectivity @xmath62 . the solid arrow indicates the approximate position of the boson peak frequency , @xmath63 , while the dashed arrow indicates the position of the lowest van hove peak , @xmath64 . for @xmath24 , @xmath65 . the low - energy part of the spectrum , including the boson peak , appears practically identical for the rn glass and for the randomly - cut fcc crystal . ]    , for the same data of fig.3 . ] ( which marks the start of the random - matrix regime linear in @xmath8 in the dos ) , and the boson peak frequency @xmath63 , both plotted as a function of the average connectivity @xmath10 . the data points refer to both random network and defective fcc systems , since both systems have exactly the same values of @xmath5 and @xmath63 . the dashed lines are guides to the eye . ]    the most striking fact , in fig.3 and fig.4 , is that the low - frequency part of the spectrum , including the boson peak , is practically identical for the rn glass and for the randomly - cut fcc crystal . this is an important observation which calls for a mechanistic explanation . the structural origin of the boson peak in our system can not be traced back to spatial fluctuations of the local elastic modulus , because the @xmath10 distribution is very narrow by construction , hence the connectivity and the shear modulus are spatially homogeneous . furthermore , anharmonic effects , also invoked in the past to explain the boson peak  @xcite , obviously play no role because atoms , in our simulations , are connected by strictly harmonic springs . one is then tempted to look for an explanation based on microstructure . what is really puzzling , however , is that the glass and the randomly - cut fcc crystal display the same boson peak and low - frequency spectrum , in spite of having widely different microscopic structure and disorder . in the rn lattice , nn bonds can have any orientation and the nn unit vector orientation is distributed nearly at random ( isotropically ) in the solid angle ( apart from some weak correlations due to the self - organization of the network ) . in the randomly - cut fcc crystal , instead , the nn bonds , basically with no exceptions , can have very few orientations only , which are dictated by the crystallographic structure . this important microstructural difference between the glass and the randomly - cut crystal becomes evident upon quantifying the bond - orientational order in the two systems . to this aim , we employ the standard bond - orientational order parameter @xmath66 , which has been used many times on glasses and defective crystals  @xcite . for each pair of nn atoms @xmath21 and @xmath30 , one first defines the correlator of nn orientations , @xmath67 , where @xmath68 is the usual definition of the local bond - orientational order parameter in terms of spherical harmonics  @xcite . one then defines the local bond - orientational order parameter as @xmath69 $ ] , where @xmath70 is a threshold equal to @xmath71 , as discussed in  @xcite , while @xmath72 is the connectivity of atom @xmath21 and @xmath73 the heaviside function . we finally average @xmath74 over all atoms in the system to obtain @xmath75 . the latter parameter measures the degree of correlation among bond orientations , or in simple words , how many bonds are aligned along the same directions . hence , @xmath75 has its largest value and is equal to 1 for crystal lattices where all bonds are aligned along the crystallographic orientations . we thus find @xmath76 for our randomly - cut fcc crystal under all conditions , as shown in fig.4 . this was expected from the fact that practically all surviving ( non - severed ) bonds in our randomly - cut crystal are perfectly aligned with the crystallographic directions , and thus have a very high degree of correlation reflected in the @xmath66 being close to 1 . this is different from other defective crystals like those studied in  @xcite , where bond - orientational disorder is important because e.g. interstitial atoms introduce bond - orientations which differ from those prescribed by the crystal lattice . the fact that some bonds are not oriented along the crystallographic axes leads , in that case , to @xmath66 values significantly below 1 . we also calculated @xmath66 for our model rn glass , and in this case we find a much smaller value , about @xmath77 , consistent with the large degree of bond - orientational disorder in our rn glass . we thus face the question of why such widely different degrees of bond - orientational order , for glass and crystal , can coexist with the same boson peak .    in effect , it appears that the microstructural mechanisms proposed in the past to explain the boson peak , can not be responsible for the boson peak in our randomly - cut fcc crystal . we have already showed above that the key mechanism which controls the softening of the shear modulus in disordered solids is the local inversion - symmetry breaking which is active in both the glass and the randomly - cut crystal , see fig.1 . within the nonaffine response formalism used in our analysis , the nonaffine part of the modulus @xmath53 is also closely related to the density of states @xmath2 , and hence to the boson peak , via eq.(34 ) of ref . this fact strongly supports the concept we propose here that the local inversion - symmetry breaking is directly related to the emergence of the boson peak . in order to confirm this hypothesis , we shall now quantify the degree of inversion - symmetry breaking in the two systems . bond - orientational disorder is known to play a role in the glassy behaviour of orientationally disordered organic crystals  @xcite . it is important to notice , however , that , in those systems , orientational disorder is also coupled to anisotropy of the organic molecules . the combination of local orientational disorder and molecular anisotropy leads to breaking the local inversion symmetry , such that the same mechanism of force - imbalance due to local inversion - symmetry in the affine positions ( fig.1 ) and the ensuing nonaffine softening must be important in those systems as well and should be analysed quantitatively in future work . we thus believe that , also in orientationally disordered crystals , inversion - symmetry breaking , resulting from bond - orientational disorder coupled to molecular anisotropy , can be identified as the source of soft modes and boson peak behaviours .    ) and the standard bond - orientational order parameter ( @xmath66 ) for the rn glasses and for the randomly - cut fcc crystals as a function of connectivity @xmath10 . there is a perfect collapse of @xmath78 for the two systems onto a master curve as a function of @xmath10 . ] to this aim we propose a new order parameter which , unlike the standard @xmath66 , is sensitive to the degree of local inversion - symmetry breaking of the lattice and we shall test how it correlates with both the shear modulus and the boson peak . a good starting point is the absolute value of the sum of all nearest - neighbour force vectors ( squared ) in the affine configuration ( affine force vectors ) @xmath79 , which is identically zero for perfect centrosymmetric crystal lattices and has its largest values for lattices where the local inversion symmetry is completely absent . to measure the degree of symmetry breaking independent of the direction of deformation , we additionally sum over all possible carthesian coordinate pairs @xmath80 . the order parameter for local inversion - symmetry is thus defined as @xmath81 , where @xmath82 indicates the limit in which inversion symmetry is completely broken and there can not be any correlations whatsoever between bond orientations . for the latter case , we found @xmath83 , a result derived in appendix c. assuming that each lattice site has the same coordination number @xmath10 , we can simplify the denominator to @xmath84 . hence , @xmath85 for any _ perfect _ centrosymmetric lattice , while @xmath86 for the limiting configuration at which the local breaking of inversion - symmetry is maximum . the new @xmath78 order parameter has a strong correlation with the nonaffine part of the shear modulus , reflected in the empirical relation @xmath87 , with @xmath88 for the fcc case , which we obtain from the simulations . importantly , the values of the is order parameter for both fcc and random network appear to be basically the same in fig.5 . this crucial observation lends further support to the conclusion that the boson peak is controlled by inversion symmetry and this is the only possible explanation to the fact that the boson peak is exactly the same for fcc and rn lattices . the order parameter for local inversion symmetry , @xmath78 , is plotted in fig.6 , in comparison with the standard bond - orientational order parameter @xmath66 . it is seen that @xmath78 displays the linear trend with @xmath10 which correlates well with both the @xmath10-dependence of boson peak frequency , and with the nonaffine shear softening , also linear in @xmath10 . further , @xmath78 displays very similar values , for both the glass and the crystal , at any given @xmath10 , which also appears consistent with the boson peaks being the same for both systems in fig.3 and fig.4 . no such correlation is displayed by @xmath66 , which remains always constant with @xmath10 , and has widely different values for the rn glass and the defective crystal , in fig.6 . in our model system there are very few physical parameters at play . it is therefore possible to check for each of them separately whether they correlate with the boson peak or not . this parameter can not control the boson peak . the density , defined as @xmath89 , i.e. the total number of atoms per unit volume , is constant in all our simulations , and is constant with @xmath10 . the connectivity @xmath10 is decreased not because we decrease the number of atoms per unit volume , which remains always the same , but because we cut bonds connecting atoms . this fact proves that the density plays no role and does not affect the boson peak . the value of force constant @xmath90 for our harmonic springs is never varied in our simulations . hence , this parameter can not control the boson peak . this parameter does not correlate with the boson peak frequency either . we have added a dashed arrow in fig.3 to mark the position of the lowest van hove peak in the dos . at the highest @xmath10 considered , both van hove peaks in our defective fcc lattice are clearly visible , and they occur at the same frequencies expected for the perfect fcc lattice ( as one can easily check ) . upon decreasing @xmath10 , the boson peak develops , and its frequency shifts to lower values . the van hove peak , instead , remains at the same frequency @xmath64 , which is a much higher frequency ( about three times higher ) than the boson peak frequency . as is shown in the inset of fig . 5 , @xmath64 remains constant upon decreasing @xmath10 in our defective fcc crystal , hence it does not correlate with the boson peak . another important observation which confirms that our boson peak is not due to the van hove singularity comes from the scaling of the dos at frequencies just above the boson peak . the scaling is @xmath91 in this regime , which means that the eigenvalue distribution scales as @xmath92 , upon recalling the definition @xmath93 . below the boson peak , and for @xmath94 , the behaviour is instead @xmath95 , i.e. fully consistent with debye law . upon approaching @xmath24 , the coefficient @xmath96 vanishes and we reproduce the random - matrix scaling @xmath97 found analytically in the marcenko - pastur distribution of random matrix theory . for @xmath98 the debye regime arises and alters this scaling . this important observation shows that the boson peak coincides with the crossover from phonon wave propagation to a regime dominated by scattering and by random - matrix eigenvalue statistics . a physical explanation for this may be found in the fact that atoms which are no longer local centers of inversion symmetry act as scattering centers for incoming waves . the same type of crossover from debye regime to random - matrix behaviour across the boson peak has been found by zamponi , parisi and coworkers for a mean - field model of glasses  @xcite . for the crossover from ballistic phonons to quasi - localized modes see also ref . @xcite . a further indication that the boson peak in our systems is of a different nature from the van hove singularity , comes from the analysis of localization . upon considering the plot of the participation ratio of vibrational modes @xmath99 , reported in in appendix a below , it is evident that both the random network and the defective fcc have localized or quasi - localized modes at about the same frequency at which the boson peak is observed . in particular , for @xmath100 , both the random network and the defective fcc have values of participation ratio approaching zero near the boson peak . it is clear that such a strong localization of vibrational modes at the boson peak frequency is not compatible with an explanation based on standard phonon physics and the van hove singularity . there is an obvious correlation between @xmath10 and the boson peak ( both bp frequency and amplitude ) , as there is also an evident correlation between z and the inversion - symmetry order parameter @xmath78 . however , we believe that the boson peak depends on @xmath10 mainly because the boson peak depends on the local inversion - symmetry , which in turn is controlled by @xmath10 due to the bond - cutting method . the reason why we do not believe that @xmath10 can be the ultimate cause of the boson peak is that , if this was true , then we should observe a very strong boson peak also in the simple cubic ( sc ) lattice which has @xmath24 , exactly . this is plainly impossible . in our systems ( random network and defective fcc ) at @xmath24 there is no trace left of debye behaviour ( see our fig . 3 above ) , whilst there is a very strong boson peak at vanishing frequency , @xmath101 . instead , in the nearest - neighbour sc lattice with @xmath24 there is obviously perfect debye behaviour and no trace of boson peak ; the dispersion relation for the sc lattice can be calculated analytically which gives @xmath102 , where @xmath7 is the wave - vector . one should therefore conclude that connectivity alone , or even isostaticity  @xcite , can not explain the boson peak because there is no boson peak in isostatic structures such as the sc lattice with @xmath24 .      as discussed with reference to fig.6 , bond - orientational disorder does not correlate with the boson peak , as in fact it remains constant while the boson peak and its frequency change significantly . this is the only physical parameter which correlates strongly with the boson peak and does not imply any paradox or contradiction . we have studied two numerical models of disordered solids : a disordered glass with bond - orientational disorder ( @xmath103 ) , and an fcc crystal with randomly - cut bonds and perfect bond - orientational order ( @xmath76 ) . in spite of the widely different bond - orientational disorder , the two systems exhibit exactly the same boson peak and almost the same nonaffine softening of the shear modulus . in particular , we showed that in both cases the boson peak frequency and the shear modulus display the same scaling with connectivity , which correlates strongly with the degree of local inversion symmetry . since this observation in our system can not be explained based on other mechanisms invoked in previous models , we arrived at the conclusion that the most likely microscopic origin of both boson peak and nonaffinity resides in the local inversion - symmetry breaking in the lattice , which is very important for both the glass and the randomly - cut crystal . this conclusion is supported by a new order parameter for centrosymmetry which displays a strong correlation with both the boson peak and the nonaffine modulus , for both the glass and the crystal . within this new framework , the boson peak is caused by the scattering of vibrational modes on atoms which are not local centers of symmetry ; such scattering and quasi - localization  @xcite events become important at nanometric length - scales ( frequencies ) comparable to the first coordination shells , as shown in previous simulation studies  @xcite . finally , our analysis identifies new local structural signatures of soft modes , which can not be traced back to purely structural quantities ( e.g. the structure factor or bond - orientational order parameters ) , yet they bear an important relation to dynamical heterogeneities because local configurations lacking inversion symmetry would be inherently unstable already under thermal stresses , e.g. in supercooled liquids  @xcite . these local mechanical instabilities could correlate with local regions of dynamical activity . in future work it will be important to ascertain the existence of quantitative correlations between our proposed inversion - symmetry order parameter and the dynamical activity order parameter used to quantify and map dynamical heterogeneity and elastic heterogeneity  @xcite at the glass transition in relation with soft modes  @xcite . this proposed shift in paradigm , and the proposed new order parameter , can be used in future studies , with the aid of new theoretical concepts  @xcite , to arrive at a unified understanding of amorphous materials . for the random network ( rn ) glass and the randomly - cut defective fcc crystals studied in this work , we also calculated the participation ratio of vibrational modes , in order to determine to which extent the modes are localized or delocalized , as a function of frequency . the participation ratio is defined as follows  @xcite : @xmath104^{-1}.\ ] ] in this expression , @xmath105 is the projection of the normalized eigenvector with frequency @xmath8 , onto atom @xmath21 , or in other words , the displacement on atom @xmath21 belonging to the collective vibrational mode @xmath8 . by construction , @xmath106 for ballistic phonons , while it is equal to zero for completely localized modes . the participation ratio is plotted below for the different values of connectivity @xmath10 , for both rn glass and fcc crystal . the qualitative behaviour is very similar to the one observed in simulations of random packings  @xcite and harmonic lattices with spring - constant disorder  @xcite . in the low frequency regime , where the linear dispersion relation and the debye law are valid , the participation ratio is always large and very close to 1 , as expected for phonons . then the participation ratio goes through a minimum corresponding approximately to the boson peak frequency , which is also close to the ioffe - regel frequency at which the physics changes from phonons to random - matrix transport . this frequency also corresponds to the wavevector or length scale at which scattering of collective vibrational modes due to local inversion - symmetry breaking becomes important .    . ]    these scattering events cause the modes to become quasi - localized  @xcite , which is reflected in much lower values of @xmath99 . the part of the spectrum just above the boson peak is dominated by randomness , and by an eigenvalue distribution characterized by the scaling @xmath107 , typical of random - matrix models  @xcite . finally at the highest frequencies of the spectrum , close to the debye frequency , the participation ratio approaches zero for anderson - localized modes . we present here a derivation of the analytical expression for the inversion - symmetry order parameter @xmath78 for the case of defective fcc crystals with randomly - depleted bonds . we start from a generic defective fcc system with a distribution of bond angles @xmath108 and @xmath109 , which define the orientation of a bond unit vector @xmath110 between two atoms @xmath21 and @xmath30 . in the framework of the affine force field , for every bond vector @xmath29 there exists a vector @xmath111 with the same probability @xmath112 in the solid angle . we now write the general expression of the total affine force field @xmath113 , as @xmath114 where @xmath115 are the cartesian coordinates . we can carry out those sums and regroup the terms to get    @xmath116    now we implement the difference between the most asymmetric configuration where inversion symmetry is completely broken , which we call the isb , and any other configuration that we want to calculate the order parameter for , such as e.g. a defective fcc crystal . + if there are no constraints whatsoever on the angular correlations between bonds connecting to the same atom @xmath21 the center of a unit cell , the second term in is zero . we can explain this by the fact that , as mentioned above , the probability to have any bond vector according to a given angular distribution is equal to the probability to have the negative of this vector ( same orientation , opposite direction ) . in the framework of the scalar product , this means that , for the probability of the quantity in the second right - most sum in eq . , the following equality must hold    @xmath117    where the averaging denotes the isotropic angular averaging @xmath118 . in a hard sphere system , one has the constraint that @xmath119 , since two bonds both connected to the same atom @xmath21 can not have an angle smaller that @xmath120 ( ultimately due to excluded volume ) . this constraint shifts the average in from zero to a negative value and lowers the final value of @xmath113 . this is so because the excluded volume correlations raise the average degree of inversion symmetry in the system with respect to a system where the excluded volume constraint on the angles the bonds is absent . in a system where no correlations exist between bond orientations such that the breaking of inversion symmetry is maximum and the local bond orientations are completely asymmetric , the only term which remains in the expression of @xmath113 is @xmath121 therefore our order parameter becomes @xmath122 this expression can be easily evaluated numerically for different lattices and provides the correctly normalized limit used in the main article to plot @xmath78 for both the fcc crystal and the rn lattice as a function of @xmath10 . we will now derive the analytical value for the affine force field of the depleted fcc lattice in order to get an analytical expression for @xmath78 . to this aim , we have to calculate @xmath123 for the two cases @xmath124 and @xmath125 . we start with the general definition of the affine force field on a generic atom @xmath21 : @xmath126 where @xmath127 are carthesian directions . since no carthesian direction or plane is in any way special , we can pick one example for each of the two cases . so we explicitly calculate @xmath128 and @xmath129 . in the first case the x component of the affine force field is :      here @xmath131 is the number of particles in the system . in the x - component we have 8 allowed bond orientations that can contribute to the affine force field @xmath132 . out of a given value of @xmath10 bonds in the unit cell , only @xmath21 contribute in the x direction . @xmath30 out of those @xmath21-contributing bonds give a positive contribution in the sum of , thus @xmath133 give a negative contribution to the sum . the absolute value of the sum is then @xmath134 times the value that each bond contributes , which is @xmath135 . since we want to calculate the absolute square of the affine force field , we have to consider the square of this value , which gives @xmath136 . for the y- and z - component we get similar expressions with the difference that now only 4 bonds contribute , in each of these two directions . for example , for the @xmath137 component we get :      and we get exactly the same for the @xmath139 component . now we just sum up the x-,y- and z - component to get @xmath140 we can use these results to easily calculate @xmath129 . the x- and y - components are equal to the y- and z - component of , as they correspond to a sum in which two of the indexes @xmath141 in are equal while one is different . this means that we have 4 contributing bonds and can apply eq .. the z - component is 0 , since any product of the three different components of the each unit bond vector vanishes in this system . so we get :          if we insert this into eq.(6 ) of this appendix upon evaluating the denominator in mean - field approximation , the expression @xmath144 leads to the following simple analytical relation : @xmath145 this situation , where @xmath86 and @xmath146 could be achieved for example in a liquid where most nearest - neighbours are short - lived and highly fluctuating , and only one mechanical bond , on average per atom , is active . u. buchenau , y. m. galperin , v. l. gurevich , and h. r. schober , phys . rev . b 43 , 5039 ( 1991 ) ; v. l. gurevich and h. r. schober , phys . b 57(11 ) , 295 ( 1998 ) ; v. l. gurevich , d. a. parshin , and h. r. schober , phys . rev . b 67 , 094203 ( 2003 ) .
Abstract:

Answer: it is well known that amorphous solids display a phonon spectrum where the debye @xmath0 law at low frequency melds into an anomalous excess - mode peak ( the boson peak ) before entering a quasi - localized regime at higher frequencies dominated by scattering . the microscopic origin of the boson peak has remained elusive despite various attempts to put it in a clear connection with structural disorder at the atomic / molecular level . using numerical calculations on model systems , we show that the microscopic origin of the boson peak is directly controlled by the local breaking of center - inversion symmetry . in particular , we find that both the boson peak and the nonaffine softening of the material display a strong positive correlation with a new order parameter describing the local inversion symmetry of the lattice . the standard bond - orientational order parameter , instead , is shown to be a poor correlator and can not explain the boson peak in randomly - cut crystals with perfect bond - orientational order . our results bring a unifying understanding of the boson peak anomaly for model glasses and defective crystals in terms of a universal local symmetry - breaking principle of the lattice .


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: it is well known that anti - de sitter ( ads ) gravity coupled to a scalar field with mass at or slightly above the breitenlohner - freedman bound @xcite admits a large class of boundary conditions , defined by an essentially arbitrary real function @xmath0 . for all ( regular ) @xmath0 , the conserved charges are well defined and finite @xcite , despite the fact that the scalar field falls off slower such that @xmath5 corresponds to turning off the independent subleading term in the asymptotic expansion of the scalar field . the usual boundary conditions instead keep this term and turn off the leading term . ] than usual . theories of this type have been called designer gravity theories @xcite , because their dynamical properties depend significantly on the choice of @xmath0 ( see e.g.  @xcite ) .    in supergravity theories with a dual conformal field theory ( cft ) description , the ads / cft duality @xcite relates @xmath0 to a potential term @xmath6 in the cft action , where @xmath7 is the field theory operator that is dual to the bulk scalar for @xmath5 boundary conditions @xcite . this led @xcite to conjecture that ( a ) there is a lower bound on the gravitational energy in those designer gravity theories where @xmath0 is bounded from below , and that ( b ) the solutions locally minimizing the energy are given by the spherically symmetric , static soliton configurations found in @xcite . more recently , the stability of designer gravity theories has been studied using purely gravitational arguments . in particular , a lower bound on the conserved energy in terms of the global minimum of @xmath0 was rigorously proven within a specific ads gravity theory by relating the hamiltonian charges to spinor charges  @xcite . arguments were given in  @xcite suggesting that these bounds hold more generally . however , it was subsequently discovered  @xcite that solutions with arbitrarily negative energy can be constructed numerically for certain theories with @xmath8 . this raised a puzzle , which we resolve in this paper . our resolution focuses on an auxiliary construct , the  superpotential " @xmath9 built from the bulk scalar potential @xmath10 . this superpotential is an important ingredient in constructing the spinor charge . in @xcite , it was shown that any @xmath10 satisfying the breitenlohner - freedman bound _ perturbatively _ admits an appropriate superpotential @xmath9 . it was of course recognized that the global existence of @xmath9 was required for the proof of an energy bound , and that this global existence may impose constraints on @xmath10 . what is interesting about the counter - examples of @xcite is that they _ do _ admit a globally defined superpotential , but nevertheless violate the proposed bound . the issue turns out to be that superpotentials associated with a given @xmath10 can be of two types , which we call @xmath4 and @xmath3 . a particular @xmath10 may admit ( distinct ) superpotentials of both types , or it may admit only the @xmath4 type . the proof requires existence of a @xmath3-type superpotential , while the examples of @xcite admit only the @xmath4-type can admit both @xmath4 and @xmath3 superpotentials was used to study asymptotically ads domain wall solutions in @xcite . ] . we verify this claim below and also confirm various other details of the arguments of @xcite . after briefly listing our conventions in section [ prelim ] , we review the proof of the lower bound in section [ two ] , illustrating why a @xmath3-type superpotential is essential ( and why a @xmath4-type superpotential is not sufficient ) . section [ solitons ] confirms that the counter - examples of @xcite admit _ only _ @xmath4-type superpotentials and numerically explores the energy bound in a number of examples . we find evidence that the global existence of @xmath3 may be necessary ( as well as sufficient ) for any lower bound to hold . section [ spinors ] investigates other aspects of the proof from @xcite . arguments for positivity of the spinor charge require witten spinors , whose global existence can be difficult to demonstrate . questions about this global existence were raised in @xcite . however , we show that the argument for global existence given in @xcite for a particular @xmath3 extends to the general case . as a check , we also explicitly demonstrate the existence of witten spinors in the context of spherical symmetry by analyzing an associated ordinary differential equation . we evaluate these spinors numerically for a particular example and use the results to check the relation between the spinor and hamiltonian charges derived in @xcite . we close with some discussion in section [ disc ] . our conventions and our definition of asymptotically anti - de sitter spacetimes follow those of @xcite . in particular , we consider gravity theories minimally coupled to a scalar field with lagrangian density given by @xmath11   \ , , \ ] ] where we have set @xmath12 . here the scalar potential @xmath10 is of the form @xmath13 near @xmath14 . for simplicity we assume that @xmath15 is even . the constant @xmath16 is the cosmological constant , given by @xmath17 with @xmath18 a positive length that we may set to one by rescaling the metric and scalar field . the spacetime dimension is denoted @xmath19 , and we assume that @xmath20 . we will furthermore assume that the scalar field is tachyonic ( @xmath21 ) , with mass in the breitenlohner - freedman range  @xcite @xmath22 where @xmath23 . the metric of exact ads space of unit radius ( and @xmath14 ) , given by @xmath24 is an exact solution of the theory , where @xmath25 is the unit - radius round metric on @xmath26 . as part of our boundary conditions , we assume that the metric of a general solution asymptotically approaches ( [ pureads ] ) in the manner described in @xcite , and that the scalar field is asymptotically of the form @xmath27 where @xmath28    to obtain a well - defined dynamics for the linearized theory , it is necessary to impose a boundary condition at @xmath29 on the scalar field , i.e. we must impose a relation between @xmath30 and @xmath31 in  . for example , one can impose @xmath32 , leaving @xmath31 totally unspecified . we refer to this option as  fast fall - off boundary conditions . " alternatively , one may set @xmath33 , leaving @xmath30 unspecified . one may also impose more general boundary conditions of the form @xmath34 where @xmath35 is an arbitrary smooth function . under the ads / cft duality , this function @xmath0 appears as an additional potential term in the action for the dual field theory @xcite . an elegant way to prove energy bounds is witten s spinor method  @xcite , which proceeds by constructing a manifestly positive `` spinor charge , '' and then comparing it to the energy of the gravitational solution . witten s argument was originally given in the context of asymptotically flat spacetimes , but it can be generalized to the asymptotically ads situation . when the matter satisfies the dominant energy condition apart from the negative cosmological constant term ( regarded as a contribution @xmath36 to @xmath37 ) , witten s argument may simply be modified by the addition of a term @xmath38 to the covariant derivative which acts on the spinor , where @xmath1 is a constant proportional to @xmath39 . this term is needed to deal with the negative energy associated with the cosmological constant . our interest here is in tachyonic scalars @xmath40 , whose potential energy may in fact be unbounded below and does not satisfy the dominant energy condition . as shown by @xcite for @xmath41 ( based on @xcite , and extended to higher dimensions by @xcite ) , many such settings may be addressed by generalizing the constant @xmath1 to a real `` superpotential '' @xmath9 satisfying @xmath42 where @xmath10 is the scalar potential .    taking @xmath43 to be the ads vacuum , we are interested in potentials for which @xmath44 . the value @xmath45 then determines the cosmological constant . with fast fall - off boundary conditions ( @xmath32 ) , the proof @xcite requires only that there be a solution to ( [ vtop ] ) with @xmath46 and @xmath47 . perturbatively , i.e. , in the sense of formal power series , a solution of this form always exists when @xmath48 sets the scalar field mass to satisfy the breitenlohner - freedman bound . in fact , there are two such perturbative solutions : @xmath49 where @xmath50 are given by ( [ roots ] ) . however , the proof of @xcite requires that @xmath9 be well - defined ( and real ) for all @xmath40 , and this imposes further restrictions on @xmath10 .    in both the fast @xcite and slow fall - off cases @xcite , the proof proceeds by using @xmath51 to construct a spinor charge @xmath52 . in the fast fall - off case one may show ( e.g. , following the basic method outlined by @xcite ) that @xmath53 , where @xmath54 is the conserved energy . the proof of @xcite is identical no matter which superpotential is used , and it is sufficient that only @xmath4 existboth @xmath4 and @xmath3 exist near @xmath14 . however , a real solution ceases to exist if ( [ vtop ] ) forces @xmath55 to become imaginary ; i.e. , if @xmath56 . since @xmath57 near @xmath14 , and since this in turn implies @xmath58 , one finds @xmath59 for all @xmath60 and global existence of @xmath3 implies global existence of @xmath4 , but not vice versa . ] . however , in order to derive an energy bound with slower fall - off conditions , @xcite assumed the existence of @xmath3 . we verify in subsection [ critical ] below that this choice is critical in this context , and that @xmath4 alone does not lead to an energy bound . this turns out to resolve the issue raised in @xcite . a simple example is discussed in section [ simple ] . we now quickly repeat the derivation of the energy bound from @xcite , examining both the original argument ( using @xmath61 ) and an analogous argument based on @xmath62 . for either superpotential , the spinor charge is defined as @xmath63 where the integrand is the hodge dual of a suitably defined nester two - form  @xcite @xmath64 } \widehat{\nabla}^e \psi + \textrm{h.c.}\ , , \ ] ] and @xmath65 is a surface at spatial infinity that bounds a spacelike surface @xmath66 . in ( [ nester ] ) , @xmath67 is a dirac spinor and the covariant derivative is defined in terms of @xmath1 ( @xmath68 ) as @xmath69 we require that the spinor field @xmath67 approaches a covariantly constant spinor ( i.e. a killing spinor ) of pure ads , @xmath70 , at infinity . we furthermore assume that asymptotically @xmath71 . using gauss s theorem we can rewrite the spinor charge @xmath72 ( @xmath73 ) as @xmath74 where @xmath75 is the integration element on @xmath66 , and @xmath76 the unit normal . letting @xmath77 denote directions in the tangent space of the surface @xmath66 , one can then show @xcite that the integrand of ( [ q2 ] ) is @xmath78 \ , , \ ] ] where @xmath79 the first and third terms in  ( [ divb ] ) are manifestly nonnegative as written . a negative contribution from the second term can be avoided by imposing the witten condition @xcite @xmath80 which is essentially the spatial dirac equation . in section  [ spinors ] , we recall the argument  @xcite that globally smooth spinors satisfying   with the above boundary conditions exist in designer gravity . this establishes @xmath81 . however , an energy bound can be derived only once we relate @xmath72 to the physical energy @xmath54 . while these coincide for fast fall - off boundary conditions , they differ in the slow fall - off case . one consequence of this is that @xmath72 is not in general conserved . another is that @xmath72 may now depend on whether the spinor charge is defined using @xmath4 or @xmath3 . that is , @xmath82 .    in  @xcite , the covariant phase space method of  @xcite , was used to show ( following @xcite ) that energy in designer gravity takes the form ) describes the generic case . additional terms are present for special cases . see @xcite . ] : [ generale ] e = - _ c e_ab ^a ds^b - ( _ + - _ - ) _ c _ b ds^b , where @xmath83 is the time translation conjugate to @xmath54 , @xmath84 is the integration element on the cut @xmath85 of the ads conformal boundary @xmath86 , and @xmath87 is the suitably rescaled electric part of the weyl tensor , which is smooth at @xmath86 as a consequence of our boundary conditions . the numerical value of @xmath54 is independent of the particular choice of the cut . one can derive the relation between @xmath52 and the energy @xmath54 @xcite by expanding the metric and the spinors in an asymptotic series in @xmath88 , using the witten equation , and using einstein s equation . the result is : @xmath89 d\omega + \frac{1}{2 } \lim_{r \to \infty } ( \lambda_{\pm}-\lambda_- ) \ , r^{d-1 - 2\lambda_- } \int_{c_r } \alpha^2\ , d\omega \ ,   , \ ] ] where for simplicity we have chosen @xmath90 to be a cut of constant @xmath91 so that @xmath92 is the integration element of the unit sphere @xmath26 . in the final term , @xmath93 is a large sphere of radius @xmath94 in @xmath66 . choosing @xmath3 and using @xmath95 yields the bound @xmath96 where @xmath97 is the volume of the unit - radius @xmath26 . on the other hand , for @xmath98 choosing @xmath4 causes the final term in ( [ qew ] ) to diverge . since @xmath54 and @xmath0 are manifestly finite , it follows that @xmath99 is infinite . it is also infinitely larger than @xmath54 , and its positivity yields no lower bound on @xmath54 . thus , the energy bound  ( [ bound ] ) has been established only if the theory admits a @xmath3 superpotential ; the existence of @xmath4 alone is not sufficient . let us consider a simple example shown in @xcite to have an energy @xmath54 which is unbounded below . we take @xmath41 and the superpotential @xmath100 this choice corresponds to @xmath101 ( that is , @xmath102 ) , and since the coefficient of @xmath103 is @xmath104 , we see from ( [ series ] ) that ( [ simpleex ] ) is of the @xmath4 type . to see that the corresponding @xmath3 does not exist , let us write  ( [ vtop ] ) as @xmath105     for the example presented in the text . the dashed line corresponds to the @xmath4 solution , which is stationary at the global minimum of the potential , @xmath106 . the solid line corresponds to the @xmath3 solution , whose derivative vanishes at @xmath107 . hence , @xmath3 does not meet the global existence criterion . , width=336 ]    one may solve ( [ dp ] ) by integrating out from @xmath108 and matching to the expansion ( [ series ] ) for small @xmath40 . such solutions exist until the quantity @xmath109 becomes negative . now , note that  ( [ series ] ) implies @xmath110 for all nonzero @xmath40 and similarly for their derivatives . in our example , @xmath4 has a maximum at the global minimum of the potential occurring at @xmath111 , so @xmath112 vanishes there . but this means that @xmath113 must vanish at some @xmath114 . in fact , it must vanish linearly at @xmath115 , since @xmath116 would force @xmath117 by ( [ vtop ] ) . thus a real @xmath3 can not exist for @xmath118 . this is illustrated in figure  [ pprime ] using numerical solutions of  ( [ dp ] ) .    since only @xmath4 exists for all @xmath40 , one does not expect ( [ bound ] ) to hold for solutions which explore large values of @xmath40 . indeed , @xcite found solutions with @xmath119 and arbitrarily negative energy for this potential . in order to confirm the above resolution of the puzzle raised in @xcite , we now numerically investigate energy bounds in a simple class of examples . for designer gravity boundary conditions that preserve the full ads symmetry group , the existence of solitons has proven to be a reliable indicator whether or not a theory satisfies a positive energy theorem @xcite . when the theory has a static spherical soliton solution , no such theorem can hold because ads - invariant gravitational solitons can always be rescaled to obtain solutions with arbitrarily negative energy that obey the same boundary conditions .    in this section we consider a one - parameter class of scalar potentials and construct spherical solitons numerically . we find that if the boundary conditions are specified by a @xmath0 which is i ) ads - invariant and ii ) bounded below , then solitons exist _ if and only if _ @xmath3 does not exist . this provides a strong test of the above stability proof , and it suggests that the existence of @xmath3 is both a sufficient and a necessary condition for the energy to be bounded from below in designer gravity . we consider the following class of potentials in @xmath41 dimensions : [ pot ] v()=-3- ^2 -2 3 ^6 + a ^8 , where @xmath120 is a free parameter . these yield scalar potentials with a negative maximum at @xmath14 , and with two global minima at @xmath121 . small fluctuations around @xmath14 have @xmath122 , which is within the range ( [ range ] ) . hence asymptotically the scalar generically decays as [ genfall2 ] = r + r^2 +  and the asymptotic behavior of the @xmath123 metric component reads [ asmetric2 ] g_rr = 1 r^2 - ( 1+^2/2 ) r^4 + o ( r^-5 ) . ( 0,0 ) ( 90,201)@xmath124 ( 388,22)@xmath125     designer gravity boundary conditions . the dashed line shows a potential that violates the positive energy theorem for designer gravity boundary conditions , but satisfies such a theorem for fast fall - off boundary conditions ( @xmath126).,width=336 ]    to construct the corresponding superpotentials @xmath2 we solve ( [ vtop ] ) , starting with @xmath127 and with @xmath128 . a solution to ( [ vtop ] ) exists unless @xmath109 becomes negative . as we integrate out from @xmath14 , @xmath1 is initially increasing and @xmath109 remains positive because the scalar satisfies the breitenlohner - freedman bound . for sufficiently large values of @xmath129 the global minima at @xmath130 will not be very much lower than the local maximum at @xmath14 , so one expects global solutions for @xmath2 to exist , with @xmath131 . by contrast , if the global minima are too deep , the quantity under the square root may become negative before the global minimum is reached , and consequently a real superpotential @xmath1 will not exist . clearly there are two distinct critical potentials ( with parameters @xmath132 ) , corresponding to those where @xmath133 just vanishes as the global minimum is reached . numerically solving ( [ vtop ] ) for a range of values for @xmath129 , one finds @xmath134 , and @xmath135 . we plot the latter critical potential in figure 2 ( full curve ) . the dashed line in figure 2 shows a potential with @xmath136 . we now look for static spherical soliton solutions of the theory ( [ pot ] ) satisfying ads - invariant boundary conditions with @xmath0 bounded below . since @xmath137 and @xmath138 , ads - invariant boundary conditions are given @xcite by @xmath139 for real @xmath140 . we see that @xmath0 is bounded below only for @xmath141 , so that only solitons with @xmath33 lead to violations of ( [ bound ] ) . the metric of a general such solution takes the form [ metric ] ds^2=-h(r)e^-2(r)dt^2+h^-1(r)dr^2+r^2d_2 ^ 2 , and the field equations read [ hairy14d ] h_,rr+(+_,r^2h+h_,r ) _ , r = v _ , , [ hairy24d ] 1-h - rh_,r-_,r^2h = r^2v ( ) , [ hairy34d ] _ , r = -1 2r _ , r^2 . regularity at the origin requires @xmath142 and @xmath143 at @xmath144 . rescaling @xmath91 shifts @xmath145 by a constant , so its value at the origin is arbitrary . thus soliton solutions can be labeled by the value of @xmath40 at the origin and the set of all soliton solutions of a particular potential with a negative maximum is found by integrating the field equations ( [ hairy14d])-([hairy34d ] ) for different values of @xmath40 at the origin . for @xmath146 the scalar asymptotically behaves as ( [ genfall2 ] ) with constant @xmath147 . the soliton therefore defines a point in the @xmath148 plane for each such @xmath149 . repeating for all @xmath150 yields a curve @xmath151 . given a choice of boundary condition @xmath152 , the allowed solitons are simply given by the points where the soliton curve intersects the boundary condition curve : @xmath153 . here we are interested in the existence of @xmath154 solitons , for potentials of the form ( [ pot ] ) . figure 3 shows how the value of @xmath125 for the @xmath154 soliton changes when one increases the potential parameter @xmath129 , from @xmath155 to its critical value @xmath156 . this corresponds to deforming @xmath15 from the potential given by the dashed line in figure 2 , to the critical potential at which @xmath3 begins to exist . one sees that @xmath157 precisely when @xmath158 . furthermore , for @xmath159 the @xmath151 curve intersects the @xmath154 axis only at the origin , so no nontrivial soliton solution exists in this parameter regime is further decreased to @xmath160 , an intersection point at finite @xmath124 appears between @xmath161 and the @xmath126 axis . this supports the claim of @xcite that the theory admits a positive energy theorem for @xmath126 scalar boundary conditions only when @xmath15 can be derived from a superpotential ( @xmath162 ) . ] . it appears , therefore , that regular spherical @xmath5 soliton solutions cease to exist precisely when @xmath158 . ( 0,0 ) ( 90,201)@xmath129 ( 388,22)@xmath125     for the @xmath5 soliton in a range of theories with different values for @xmath129 , the coefficient of the @xmath163 term in ( [ pot ] ) . one sees that regular spherical @xmath5 soliton solutions cease to exist precisely when @xmath164 . , width=336 ]    the existence of scalar solitons with ads - invariant boundary conditions for @xmath165 implies there are negative energy solutions in these theories . this can be seen as follows @xcite . starting with a static soliton @xmath166 , consider the one parameter family of configurations @xmath167 . because we chose conformally invariant boundary conditions , these are again satisfied by the rescaled data . it then follows from the constraint equations that the total energy of the rescaled configurations takes the form e_= ^-3 e_1 + ^-1 e_2 where @xmath168 is independent of the potential and is manifestly positive , and both @xmath169 are finite and independent of @xmath170 . furthermore , because the static soliton extremizes the energy @xcite , one has 0=d e_d|_=1 = -3 e_1 -e_2 and hence @xmath171 . therefore the contribution to the energy that scales as the volume , which includes the potential and scalar terms , is negative . this means that rescaled configurations are initial data for time dependent ( but time - symmetric ) solutions . for sufficiently small @xmath170 one has a large central region where @xmath40 is essentially constant and away from an extremum of the potential . hence one expects the field to evolve to a spacelike ` big crunch ' singularity @xcite . ] @xmath172 with @xmath173 must have negative total energy , and hence violate the energy bound ( [ bound ] ) . thus , the energy is unbounded below for potentials with @xmath174 and @xmath5 , indicating that existence of @xmath3 is necessary for energy in designer gravity to be bounded below . ( note that we have already proven that it is sufficient . ) as explained in section  [ two ] , the derivation of the lower bound for the energy @xmath54 uses the positivity of the spinor charge . this requires the global existence of smooth spinors satisfying the witten condition l = 0 on a constant time surface @xmath66 with suitable asymptotic behavior , where @xmath175 is the operator in the witten equation , [ witteneq ] l = ^i _ i . here , @xmath176 is the first order covariant differential operator defined in eq .  , taken in directions @xmath177 tangent to the spacelike surface @xmath66 . global existence of such solutions is non - trivial to establish , and indeed it was suggested in @xcite that some subtle failure of such spinors to exist might be responsible for the lack of energy bounds in the systems studied there . in contrast , we have now proposed that the failure of those examples to satisfy ( [ bound ] ) was due solely to the fact that such potentials do not admit a @xmath3-type superpotential . in order to dispel any remaining concerns about the global existence of spinors , we reexamine this issue below . section [ again ] essentially follows @xcite and proves global existence of the witten spinors in general , but we deviate slightly from  @xcite to correct an error in this proof . in order to correct this proof , we need to work on a maximal slice @xmath66 ( i.e. , one for which the trace of the extrinsic curvature satisfies @xmath178 ) in the spacetime that intersects scri transversely . this proof is in fact quite abstract and draws on certain very non - trivial results in the mathematics literature . we therefore restrict attention in section [ sphere ] to the particular case of maximal ( e.g. , time - symmetric ) , spherically symmetric hypersurfaces , which is enough to address any concerns raised by @xcite . in this case we are able to give a more explicit construction of the witten spinors by reducing the problem to ordinary differential equations . existence of solutions is then straightforward to show . finally , we numerically solve for the witten spinors in section [ numerics ] and evaluate the corresponding spinor charges , verifying the relationship between @xmath72 and @xmath54 derived in section [ two ] .      if @xmath66 is a maximal slice in the spacetime , then the witten operator @xmath175 in eq . can be written as [ witteneq1 ] l = ^i ( d_i + p ( ) _ i ) , where @xmath179 is now the spin - connection intrinsic to @xmath66 . on slices with @xmath180 , there would be another term proportional to @xmath181 in the above expression for @xmath175 . for the proof to work , we must assume that such a slice exists . below , we will briefly outline how one might construct such a @xmath66 . we would like to find a smooth spinor field @xmath67 on @xmath66 with the property that @xmath182 , and that @xmath67 tends to a suitable exact killing spinor of ads - space near infinity . following @xcite , the first step in the proof is to find a formal power series solution @xmath183 to the witten spinor equation of the form [ asympt ] ~_0 + r^-1 _ 1 + r^-2 _ 2 +  + r^-n _ n +  , where @xmath184 , and where @xmath70 is a killing spinor in exact ads - space ( i.e. , @xmath185 in exact ads with vanishing @xmath40 ) such that - ( ) ^a = _ 0 ^a _ 0^ . the spinor fields @xmath186 are determined recursively by the condition that @xmath183 be a formal power series solution to the witten equation . it can be proven that @xmath187 are uniquely determined for @xmath188 by making a split of the witten equation into a part containing an @xmath94-derivative and a part containing derivatives tangent to a sphere of constant @xmath94 . the explicit form of the first coefficients @xmath187 in terms of the spacetime curvature near infinity is given in  @xcite . the formal power series solution satisfies the equation l= j , where @xmath189 is a smooth source vanishing faster than any inverse power of @xmath94 near @xmath86 . we may terminate the expansion   at some large finite @xmath190 , thus obtaining a @xmath189 vanishing faster than @xmath191 . the next step is to obtain from the formal power series solution @xmath183 a global solution @xmath67 to the witten equation . this step requires some global analysis . the idea is to write the global solution @xmath67 that we seek as = + , where @xmath192 is a smooth spinor field to be determined . thus , @xmath192 should satisfy @xmath193 , and @xmath192 should vanish at least as fast as @xmath194 in order that @xmath67 have the same asymptotics as @xmath183 to leading order . the global existence of such a @xmath192 can be established using the following key inequality   which holds whenever @xmath195 is a solution to the einstein equations satisfying the asymptotic conditions given in @xcite . to state the inequality , consider any smooth spinor field @xmath196 on @xmath66 that vanishes outside some compact region in @xmath66 , and introduce the following norm on such @xmath196 : [ normu ] u ^2 = _ ds . then there exists a positive constant independent of @xmath196 such that [ ineq ] const.^-1 u ^2 _ |l u|^2 ds u ^2 . the proof of the inequality uses a similar technique as in the `` hardy - inequality , '' combined with the key identity  . it may be found in  @xcite . we now use the inequality to establish a solution to the equation @xmath193 . we will do this by first constructing a distributional solution @xmath192 , and then showing that it is in fact smooth . let @xmath197 be the linear functional on smooth compactly supported spinors @xmath196 defined by f(u ) = -_(nj)^u ds , where @xmath198 , and @xmath76 is the unit timelike normal to @xmath66 . note that @xmath199 , and that [ hermitian ] n l^= l n due to the fact that @xmath178 . by the cauchy - schwartz inequality and the first half of the inequality  , we get |f(u)| & & ( _ |u|^2 ( 1+r^2)^-1 ds ) ^1/2 ( _ |j|^2 ( 1+r^2 ) ds ) ^1/2 + & & const . ( _ |lu|^2 ds ) ^1/2 . the constant is finite because by construction @xmath189 drops off sufficiently rapidly at infinity . we interpret this inequality as saying that @xmath197 is a bounded functional with respect to the positive definite scalar product given by v | u = ( lv)^lu ds . let @xmath200 be the hilbert space defined by this inner product ( which , by the second half of inequality   is identical to the hilbert space obtained from the norm @xmath201 ) . by the riesz representation theorem , there is hence an element @xmath202 such that @xmath203 for all @xmath204 . again by the inequality  , every element in @xmath200 is locally square integrable , so in particular a distribution on @xmath66 . thus , @xmath205 is a solution in the distributional sense of the equation @xmath206 . next , we apply @xmath207 to this equation , use eq .  , and put @xmath208 . it follows that this @xmath192 is the desired distributional solution to @xmath193 , and @xmath209 is a global solution to the witten equation . it remains to prove that @xmath67 is smooth , and that it satisfies the desired boundary conditions . this will follow if we can show that @xmath192 is smooth and vanishes sufficiently fast at infinity . it follows from our construction so far only that @xmath210 . but since @xmath189 is smooth and vanishes quickly , one can now use the mapping properties of the parametrix of the operator @xmath175 established in  @xcite to prove that @xmath192 is indeed smooth and vanishes sufficiently fast . the details of this argument are given in  @xcite . finally , we outline how one might construct the desired @xmath178 slice @xmath66 in the type of spacetimes that we consider . first , let @xmath211 be any ( weakly ) spacelike hypersurface in spacetime , @xmath212 its domain of dependence . if @xmath211 is such that the closure of @xmath212 is contained in a globally hyperbolic subregion of spacetime ( note that the full spacetime need _ not _ be globally hyperbolic ) , then a theorem of bartnik  @xcite guarantees that there exists a sufficiently regular maximal slice @xmath66 with prescribed extrinsic curvature and boundary @xmath213 . unfortunately , it is not immediately clear how to apply this theorem when the surface @xmath211 touches the conformal boundary . however , by considering a sequence @xmath214 eventually touching the conformal boundary transversely , we obtain a corresponding sequence @xmath215 of maximal slices . presumably it can be shown that these have a subsequence that is convergent in a sufficiently strong sense . we now explicitly describe the solution to the witten equation for a slice @xmath66 which is both spherically symmetric and maximal ( i.e , @xmath178 ) though not necessarily static . such a slice automatically exists when the initial data are time - symmetric . since this case is sufficient to address the concerns of @xcite , we explain the analysis in some detail . our conventions will be as follows . for simplicity , we restrict attention in this subsection to the case @xmath41 . we denote by @xmath216 indices on a flat internal space ( while @xmath217 are spacetime indices ) . the gamma matrices satisfy @xmath218 and we choose the explicit matrix representation @xmath219 where @xmath220 is the @xmath221 identity matrix , and where @xmath222 are the standard pauli matrices . we wish to address spherical spacetimes . let @xmath66 be a maximal ( @xmath178 ) , spherically symmetric spacelike slice on which we set @xmath223 , e.g. , a surface of time - symmetry . though such spacetimes are not in general static , on @xmath66 they may nevertheless be written in the form ( [ metric ] ) , which may be described by the orthonormal frame @xmath224 because @xmath178 , the witten equation   only depends on the geometry intrinsic to @xmath66 , i.e. , @xmath225 . we further assume that the scalar field on @xmath66 is a function of @xmath94 alone , and that @xmath226 , @xmath227 are known functions . the intrinsic covariant derivative acting on spinor fields is in general given by @xmath228 where the intrinsic spin connection is determined by the relation @xmath229 . for the metric  ( [ metric ] ) , we find @xmath230    with the above ingredients in hand , the witten condition becomes the explicit equation @xmath231\psi   + \frac{3}{2 } p(\phi ) \psi \,.\ ] ] writing @xmath67 in terms of two - component spinors as @xmath232 ( [ witcon ] ) implies @xmath233\psi_2   + \frac{3}{2 } p(\phi ) \psi_1 , \\ \label{witcon22 } 0 & = & i \left[\sigma_{3 } \sqrt{h } \left(\partial_r + \frac{1}{r } \right ) + \frac{1}{r } \ , \rlap/\!\nabla_{\ ! \!s^2 } \right]\psi_1   - \frac{3}{2 } p(\phi ) \psi_2 \,.\end{aligned}\ ] ] here the operator @xmath234 is the dirac operator on the 2-dimensional sphere @xmath235 . to decouple the differential equations  ( [ witcon21 ] ) and  ( [ witcon22 ] ) , we define @xmath236 then the @xmath237 satisfy @xmath238 \psi_+ & = & 0 \\ \label{witcon32 } \left [ \sigma_{3 } \left(\partial_r + \frac{1}{r } \right ) + \frac{1}{r \sqrt{h } } \ , \rlap/\!\nabla_{\ ! \!s^2 } - \frac{3}{2 \sqrt{h } } p(\phi ) \right ] \psi_- & = & 0 \ , .\end{aligned}\ ] ]    next we wish to expand @xmath237 in spinor spherical harmonics , which are eigenfunctions of @xmath234 . these are given , for example in  @xcite , by @xmath239 @xmath240 where @xmath241 are integers such that @xmath242 and @xmath243 . the @xmath244 are jacobi polynomials and the spherical harmonics satisfy @xmath245 separating variables as @xmath246 @xmath247 we find from  ( [ witcon31 ] ) and  ( [ witcon32 ] ) that the radial functions satisfy the following differential equations : @xmath248 @xmath249    to pick the desired solution to the witten equation , we now impose our boundary conditions . we demand that our solution be regular in the interior , and that it asymptotically approaches a suitable killing spinor @xmath70 of pure ads . as we will see shortly , this @xmath70 has only @xmath250 components . since the remaining modes are decoupled , they may be consistently set to zero . thus we consider only @xmath251 below .    near the origin @xmath144 , the solutions behave as @xmath252 where @xmath253 are constants . regularity at the origin then requires boundary conditions of the form @xmath254 the analysis for the @xmath255 solutions near the origin is exactly the same . as @xmath256 , we require that the solution to the witten condition approaches a killing spinor in exact ads space . for large @xmath94 , solutions to  ( [ radialde ] ) and  ( [ tradialde ] ) behave as @xmath257 the solution which grows as @xmath258 matches the known asymptotic behavior of killing spinors in exact ads space . in particular , an explicit solution to the killing spinor equation is given by @xcite @xmath259 where @xmath260 is a constant spinor and @xmath261 the spinor @xmath260 is chosen so that @xmath70 satisfies the additional requirement @xmath262 , or equivalently @xmath263 . one can achieve this normalization by taking @xmath264 with the choice  ( [ norm ] ) , one can rewrite  ( [ ks ] ) in terms of the spinor spherical harmonics as @xmath265 hence , the solution for the witten spinor is @xmath266 \end{array } \right ) \,.\ ] ] furthermore , comparing to  ( [ ks2 ] ) we see that the asymptotic conditions on the radial functions are @xmath267    comparing ( [ wcasym ] ) and ( [ wcbc ] ) with ( [ near0 ] ) and ( [ largerr ] ) demonstrates that the desired witten spinors exist . one first notes that @xmath268 decouple from @xmath269 . to construct the desired @xmath268 , one simply enforces the boundary condition ( [ wcbc ] ) at @xmath144 for some @xmath270 and integrates ( [ radialde ] ) outward . from ( [ largerr ] ) , any solution will grow asymptotically as @xmath258 , and in fact will satisfy ( [ wcasym ] ) up to an overall scale . rescaling @xmath270 yields the desired solutions @xmath268 . the analysis for @xmath271 is similar . thus , witten spinors exist for any spherically symmetric maximal hypersurface @xmath66 satisfying our boundary conditions . a final crucial aspect of the argument of section [ two ] is the detailed relation ( [ qew ] ) between the spinor charges @xmath52 and the energy @xmath54 . since ( [ qew ] ) was derived via a tedious calculation , it is useful to verify the relation numerically . let us consider the theory with potential @xmath272 which is just ( [ pot ] ) for @xmath273 . this theory admits both a @xmath4 and a @xmath3 superpotential . we take as initial data a soliton with @xmath274 and boundary condition @xmath275 for some constant @xmath276 ( so that @xmath0 is not conformally invariant ) . we find that a soliton exists with @xmath277 and @xmath278 . the energy of this soliton is @xmath279 . numerical results for @xmath3-type witten spinors are shown in figures  [ rrm ] and  [ irm ] . ) with initial data given by the rescaled configurations @xmath280 discussed in section  [ solitons ] . ]    ) confirm the existence of witten spinors that are regular at the origin and have the correct behavior at large @xmath94 . the solid line is @xmath281 , which by the arguments in the text asymptotically approach @xmath282 , shown here as the dashed line.,width=336 ]    , which by the arguments in the text asymptotically approach @xmath283 , shown here as the dashed line . , width=336 ]    using these solutions , we can also calculate the spinor charge and check  ( [ qew ] ) . we obtain @xmath284 , thus confirming the relation @xmath285 . similarly , for the @xmath4 case we find @xmath286 , in agreement with ( [ qew ] ) . we have resolved the puzzle raised in @xcite concerning energy bounds in designer gravity . while the arguments of @xcite are correct as written , global existence of an appropriate superpotential is a subtle requirement . in particular , two types of superpotentials may ( or may not ) exist for a given scalar potential @xmath10 . the proof of @xcite requires the global existence of a @xmath3-type superpotential . if one attempts to follow the same argument with a @xmath4-type superpotential , one finds that the difference between the conserved energy @xmath54 and the associated spinor charge @xmath99 diverges , and in particular that @xmath99 diverges . thus , positivity of @xmath99 does not yield a lower bound for @xmath54 , and the existence of a @xmath4-type superpotential alone is not sufficient to yield a positive energy theorem . numerical explorations support this resolution .    a specific question raised in @xcite concerned the global existence of witten spinors . we have demonstrated ( section [ spinors ] ) that no such difficulties arise , even for the models considered in @xcite . the existence theorem of @xcite was shown to hold in general and , for the special case of spherical symmetry , existence of witten spinors on maximal ( e.g. , time - symmetric ) hypersurfaces was again demonstrated using simple arguments based on ordinary differential equations . the spherical time - symmetric context is sufficient to address the concerns of @xcite . this reinforces our claim that , when @xmath0 has a global minimum , the existence of a @xmath3-type superpotential is sufficient for the energy to be bounded below . interestingly , our numerical explorations of particular examples ( section [ solitons ] ) suggest that global existence of @xmath3 may also be _ necessary _ for the energy to be bounded below . such a property would be analogous to the claim of @xcite that the existence of a superpotential is necessary for an energy bound to hold in the case of fast fall - off boundary conditions @xmath287 . however , in the fast fall - off case the relevant superpotential is @xmath4 , since either superpotential is sufficient and existence of @xmath3 implies existence of @xmath4 ( footnote [ exist ] ) . indeed , for @xmath32 the spinor charges @xmath52 satisfy @xmath288 and , in particular , @xmath289 when both potentials exist . thus we find that each type of superpotential provides a stability criterion for ads gravity - scalar theories , with @xmath4 controlling the @xmath32 case and @xmath3 controlling cases with slower fall - off conditions on the scalar . finally , a key tool in our investigation of particular examples was the argument from @xcite showing that the existence of designer gravity solitons satisfying ads - invariant boundary conditions implies that the energy is unbounded below . we may thus interpret the lower bounds of @xcite in terms of the soliton content of such theories . we conclude that designer gravity theories where @xmath3 exists can have no solitons when @xmath0 is both ads - invariant and bounded below . the authors would like to thank gary horowitz for useful discussions . s.h . would also like to thank piotr chrusciel and bob wald for discussions and the physics department at ucsb for its hospitality . d.m . also thanks jim isenberg and robert bartnik for discussions . a.a . and d.m were supported in part by nsf grant phy0354978 , and by funds from the university of california . p.  breitenlohner and d.  z.  freedman , `` stability in gauged extended supergravity , '' annals phys . * 144 * ( 1982 ) 249 ; `` positive energy in anti - de sitter backgrounds and gauged extended supergravity , '' phys . b * 115 * ( 1982 ) 197 . m.  henneaux , c.  martinez , r.  troncoso and j.  zanelli , `` black holes and asymptotics of 2 + 1 gravity coupled to a scalar field , '' phys .  rev . d * 65 * , 104007 ( 2002 ) [ arxiv : hep - th/0201170 ] . m.  henneaux , c.  martinez , r.  troncoso and j.  zanelli , `` asymptotically anti - de sitter spacetimes and scalar fields with a logarithmic branch , '' phys . d * 70 * , 044034 ( 2004 ) [ arxiv : hep - th/0404236 ] . t.  hertog and k.  maeda , `` stability and thermodynamics of ads black holes with scalar hair , '' phys . d * 71 * , 024001 ( 2005 ) [ arxiv : hep - th/0409314 ] . j. m. maldacena , `` the large n limit of superconformal field theories and supergravity , '' adv . * 2 * ( 1998 ) 231 [ arxiv : hep - th/9711200 ] . d.  z.  freedman , c.  nunez , m.  schnabl and k.  skenderis , `` fake supergravity and domain wall stability , '' phys . d * 69 * , 104027 ( 2004 ) [ arxiv : hep - th/0312055 ] . a.  celi , a.  ceresole , g.  dallagata , a.  van proeyen and m.  zagermann , `` on the fakeness of fake supergravity , '' phys . d * 71 * , 045009 ( 2005 ) [ arxiv : hep - th/0410126 ] . m.  zagermann , `` n = 4 fake supergravity , '' phys . d * 71 * , 125007 ( 2005 ) [ arxiv : hep - th/0412081 ] . k.  skenderis and p.  k.  townsend , `` hidden supersymmetry of domain walls and cosmologies , '' phys . * 96 * , 191301 ( 2006 ) [ arxiv : hep - th/0602260 ] . i.  papadimitriou ,  non - supersymmetric membrane flows from fake supergravity and multi - trace deformations , " [ arxiv : hep - th/0606038 ] . k.  skenderis and p.  k.  townsend , `` pseudo - supersymmetry and the domain - wall / cosmology correspondence , '' arxiv : hep - th/0610253 . s.  hollands , a.  ishibashi and d.  marolf , `` comparison between various notions of conserved charges in asymptotically ads - spacetimes , '' class .  quant . * 22 * , 2881 ( 2005 ) [ arxiv : hep - th/0503045 ] . m. heusler , n. straumann , `` scaling arguments for the existence of static , spherically symmetric solutions of self - gravitating system , '' class . * 9 * ( 1992 ) 2177 ; m. heusler , `` a no hair theorem for selfgravitating nonlinear sigma models , '' j. math * 33 * ( 1992 ) 3497 . r. mazzeo : `` elliptic theory of differential edge operators i , '' commun . part . diff . equations * 16 * ( 10 ) ( 1991 ) , 1615 ; r. mazzeo : `` the hodge cohomology of a conformally compact metric , '' j. diff * 28 * ( 1988 ) 309
Abstract:

Answer: we investigate the stability of asymptotically anti - de sitter gravity coupled to tachyonic scalar fields with mass at or slightly above the breitenlohner - freedman bound . the boundary conditions in these `` designer gravity '' theories are defined in terms of an arbitrary function @xmath0 . previous work had suggested that the energy in designer gravity is bounded below if i ) @xmath0 has a global minimum and ii ) the scalar potential admits a superpotential @xmath1 . more recently , however , certain solutions were found ( numerically ) to violate the proposed energy bound . we resolve the discrepancy by observing that a given scalar potential can admit two possible branches of the corresponding superpotential , @xmath2 . when there is a @xmath3 branch , we rigorously prove a lower bound on the energy ; the @xmath4 branch alone is not sufficient . our numerical investigations i ) confirm this picture , ii ) confirm other critical aspects of the ( complicated ) proofs , and iii ) suggest that the existence of @xmath3 may in fact be necessary ( as well as sufficient ) for the energy of a designer gravity theory to be bounded below .