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: various approaches have been used to understand ageing experiments @xcite @xcite @xcite : droplet picture @xcite , mean field model @xcite @xcite , and trap models @xcite @xcite @xcite @xcite . ageing is by definition the property that the two points correlation function @xmath1 depends explicitly on the waiting time @xmath2 . @xmath3 is a macroscopic observable , for instance the magnetization in the zero field cooled experiment @xcite or the thermoremanent magnetization experiment @xcite . on the theoretical level , dynamics can be implemented in several ways . one possibility is langevin dynamics which was used for instance in @xcite where mean field ageing was first obtained , or in @xcite where ageing in the presence of flat directions was studied . another possibility is glauber dynamics @xcite . within the framework of glauber dynamics and given the probability distribution @xmath4 , for the system to be in any microcanonical state @xmath5 at time @xmath6 , the probability distribution at time @xmath7 is given by @xmath8 where @xmath9 is the glauber matrix . by `` microcanonical states '' , we mean a state such as the system is in a given microcanonical configuration with probability one . this article is devoted to the study of the glauber dynamics of various models , phenomenological or microscopic , in connection with ageing . more precisely , for a finite size system of size @xmath10 , we calculate from the glauber matrix a time @xmath0 such as the finite size system equilibrates on a time of the order of @xmath0 , starting from out  of  equilibrium initial conditions at time @xmath6 . in terms of ageing , the correlation function ( [ eqcorr ] ) depends explicitely on the waiting time @xmath2 provided the system is out  of  equilibrium at time @xmath2 , namely @xmath2 is smaller than @xmath0 . in other words , ageing for a finite size system is interupted after @xmath0 [ see @xcite where the denomination `` interrupted ageing '' was proposed ] . in order to understand ageing in the thermodynamic limit ( @xmath11 ) , the variations of @xmath0 with the system size have to be carefully studied . if @xmath0 goes to a constant in the thermodynamic limit , ageing is interrupted , whereas if @xmath0 diverges ageing is not interrupted in the thermodynamic limit : in the thermodynamic limit , the system reaches equilibrium on infinite time scales . we also discuss the possibility of metastability without ageing on experimental time scales , as it is the case for instance in diamond .    in general , @xmath0 depends on the eigenvalues of the glauber matrix , but also on the overlaps between microcanonical states , and the transformation matrix elements between the eigenstates of the glauber matrix and the microcanonical states . however , in some simple cases , with only self  overlaping states , the time @xmath0 can be expressed only in terms of the density of eigenvalues of the glauber matrix . the form of the glauber matrix will be given in each case in the core of the paper . whatever the system , the eigenvalues of @xmath9 are negative . the spectrum is bounded above by zero ( zero is always an eigenvalue , corresponding to the equilibrium boltzmann distribution ) . the eigenvalues of the glauber matrix have the dimension of the inverse of a time , so that they can be viewed as the inverse relaxation times of the system . since we are interested in long time behavior , a special attention will be paid to the smallest eigenvalues . this article is organized as follows . in section [ gm ] , we present a detailed study of the average density of small eigenvalues of the random energy model ( rem ) like trap model proposed and studied in @xcite and inspired by the rem model @xcite .    in the case of the rem  like trap model , two regimes are predicted and checked numerically for the density of small eigenvalues in the glass phase : an arrhenius  like regime and a non arrhenius regime for the smallest eigenvalues . by arrhenius  like regime , we mean a regime in which the density of eigenvalues is the same as the density of eigenvalues calculated from the arrhenius law , assuming a correspondence between the eigenvalues of the glauber matrix and the inverse relaxation times of the system . the eigenvalue density in the arrhenius regime are given in section [ arrh - def ] . we also establish a connection between the behavior of the density of small eigenvalues of rem  like trap models and the spectral dimension of the associated diffusion process .    next , in section [ general ] , we give a general relation for the existence of ageing from a glauber dynamics point of view . this section is the main part of the present article . when applied to the rem  like trap model , our criterium is equivalent to the fact noticed in @xcite that the average trapping time diverges below @xmath12 . the rest of the paper is devoted to examine the case of other models : brownian motion which was shown in @xcite to be an ageing phenomenon , a model introduced by barrat and mzard @xcite ( bm model ) , the one dimensional ising chain ( in this case , our results are consistent with the existence of interrupted ageing ) , and the sherrington  kirkpatrick model ( sk model ) . in all these cases , we find consistent results between the existence of ageing calculated from the behavior of the two times autocorrelation functions and our criterium based on the glauber matrix . we end  up with some final remarks . trap models can be defined without a priori reference to a microscopic model . however , some trap models , such as the rem  like trap model @xcite are related to the parisi solution of spin glasses : even though these models are toy models , some of them retain some feature of the dynamics of microscopic models . the rem  like trap model consists of a trap model in infinite dimension , that is , the system can `` hop '' from one trap to any other trap , according to the rules of the dynamics . other toplogies will be analyzed below , in reduced dimension . the rem  like trap model corresponds to a single replica symmetry breaking step . generalizations of the rem  like trap model including multi replica symmetry breaking steps were studied in @xcite and their relevance to explain experimental data was shown . the models of this section consist of @xmath10 traps of depth @xmath13 , with the @xmath14 chosen among the distribution @xmath15 . the reason why this particular energy trap distribution is chosen is that there exists a finite transition temperature between a high and a low temperature glassy phase @xcite . the dynamics rules are given by @xmath16 for @xmath17 if the traps @xmath18 and @xmath19 are connected and zero otherwise . the diagonal terms are defined by @xmath20 which enforces probability conservation . equation ( [ propcons ] ) is a general constraint for any glauber dynamics . another constraint is the detailed balance condition @xmath21 it is clear that the aforementioned rules for the dynamics of the rem  like trap model satisfy detailed balance . however , as we shall see below , other types of dynamics also satisfy detailed balance and probability conservation . as a consequence of detailed balance , the eigenvectors of the glauber matrix are not physical states ( in the sense that they are not probability distributions ) except for the boltzmann distribution ( the stationary eigenvector corresponding to the @xmath22 eigenvalue ) . notice that if @xmath9 is a glauber matrix @xmath23 , with @xmath24 a positive number , is also a glauber matrix . multiplying @xmath9 by @xmath24 amounts to changing the unit time . the proportionality constant is imposed by physical considerations . for the models considered in this section , this constant is chosen such as the diagonal elements of the glauber matrix do not scale with the number @xmath10 of traps in the large @xmath10 limit . for instance , in the infinite dimensional case , we choose @xmath25 . this dynamics will be studied with different phase space topologies : infinite dimension , where all the traps are connected ( section [ infinite ] ) , one dimensional topology , where a given trap is connected only to its two nearest neighbors ( section [ disor ] ) , and also one dimensional topology with `` ordered '' traps , in the sense that the trap energies have been ordered @xmath26 , and the trap at site @xmath18 has an energy @xmath14 . a useful property that we will use throughout this paper is that glauber matrices can be mapped onto symmetric matrices with the same spectrum . let us call @xmath27 the diagonal matrix @xmath28 . then , because of detailed balance , the matrix @xmath29 is symmetric and has the same spectrum as @xmath9 . using a quantum mechanical analogy , we denote by @xmath30 the microcanonical state corresponding the system in trap @xmath18 with a probability unity . the analogy to quantum mechanics is useful , even though there is no hermitian structure . in particular , any linear combination of microcanonical states is not a probability distribution . we propose ourselves to carry out a detailed study of the density of eigenvalues of the rem  like trap model and related models in reduced dimension . since the arrhenius law is introduced `` by hand '' in the dynamics , it is interesting to question in which regime the density of eigenvalues is given by the arrhenius law , assuming a correspondence @xmath31 between the eigenvalues of the glauber matrix and relaxation times . we will naturally call _ arrhenius regime _ such a regime . as we will see below , the density of eigenvalues is not always in the arrhenius regime .    in this arrhenius regime , and by definition , the off - diagonal terms are neglectable as far as the density of eigenvalues is concerned . the diagonal coefficients of the glauber matrix are proportional to @xmath32 , so that the distribution of eigenvalues is @xmath33 this distribution is the same as the one expected from an arrhenius law with life times @xmath34 . if one assumes an arrhenius picture , the distribution of relaxation times is @xmath35 which is equivalent to ( [ aha ] ) . as noticed in @xcite , the average trapping time diverges if @xmath36 , leading to ageing . distribution of trapping times with no first moment were already studied in @xcite , and shown to lead to anomalous diffusion , namely to a non trivial spectral dimension , even above the upper critical dimension . we will compare our results for the rem  like trap model with the results of @xcite derived of the context of anomalous diffusion on site  disordered hypercubic lattices . the density of eigenvalues is thus described by the arrhenius law provided the off - diagonal terms of the glauber matrix are a small perturbation , so that the density of eigenvalues of the glauber matrix is given by the diagonal terms only . we now examine several cases with different phase space topologies and determine in the different cases on which condition the density of eigenvalues is determined by the arrhenius picture . we also compare our analytical approach to numerical diagonalizations of the symmetrized matrix @xmath37 .      in this section , the system can hop from one trap to any other trap ( infinite dimension ) . the diagonal coefficients are @xmath38 and the off  diagonal coefficients of the symmetrized glauber matrix are @xmath39      we treat the off  diagonal coefficients @xmath40 ( @xmath17 ) as a perturbation . this perturbative treatment is consistent since we are looking for a regime in which the off  diagonal coefficients do not modify the density of eigenvalues with respect to the purely diagonal case . the second order correction to the eigenvalue @xmath41 is @xmath42 the first order corrections being identically zero since the perturbation consists only of off  diagonal elements . we assume that the density of trap energies in the vicinity of @xmath14 is large enough for the discrete summation to be replaced by an integral . if this condition is not satisfied , the second order corrections @xmath43 may have large fluctuations because the divergence for @xmath44 is no more cancelled . the condition that the unperturbed eigenvalue density is large in the vicinity of @xmath14 simply reads @xmath45 . we assume in the rest of this section that this condition is fulfilled , and we can safely replace the discrete sum by an integral : @xmath46 the singularities for @xmath47 and @xmath48 cancel each other , and we can approximate @xmath43 as @xmath49 with @xmath50 . we have used the following approximation : if @xmath51 is small , then @xmath52 and we use @xmath53 as a cut  off . clearly ( [ app ] ) is only an approximation , but we are only interested in the dominant behavior of @xmath54 as a function of @xmath14 . now , we want to compare the approximation ( [ app ] ) of @xmath43 to the typical level spacing between @xmath41 and @xmath55 , where we assume that the traps are ordered such as @xmath56 . to to so , we first calculate the distribution of @xmath57 and we obtain @xmath58 , with @xmath59 notice here that @xmath60 diverges as @xmath61 goes to zero if @xmath36 . this accumulation of small relaxation times is related to the existence a glass transition for this model for @xmath62 . let us now calculate the level spacing statistics @xmath63 , which is the density probability that the spacing between the absolute value of the eigenvalue @xmath9 and the absolute value of the next eigenvalue is @xmath64 . because of the statistical independence of the trap depths , @xmath63 is readily obtained as @xmath65 in this expression , @xmath66 is the number of levels in the interval @xmath67 $ ] ( with @xmath68 ) . the other terms correspond to the condition that there is no eigenvalue in the interval @xmath69 $ ] . assuming that @xmath61 is small compared to unity , and that @xmath10 is large , we get @xmath70 the off - diagonal terms are irrelevant to the density of eigenvalues provided @xmath71 is much smaller than unity . this condition reads @xmath72 , with @xmath73 defined in ( [ deltai ] ) . since we are looking for a necessary condition valid if @xmath36 , we neglect the temperature  dependent prefactor and conclude that the density of states is in the arrhenius regime only for traps such that @xmath74 . notice that the traps such as @xmath75 are deeper than the typical deepest trap ( of the order of @xmath76 ) . the arrhenius picture is thus valid for a _ typical _ trap realization , with a lowest trap not deeper than @xmath76 . however , non typical disorder realizations may be generated . for these realizations , the lowest eigenvalues are not described by an arrhenius picture . the existence or absence of the arrhenius regime even for a very large number of traps is thus sample dependent . this non self averaging behavior is reminiscent on the sample dependence of the quantity @xmath77 with @xmath78 the probability to find the system in trap @xmath18 @xcite . the density of small eigenvalues of the rem  like trap model is pictured on figure [ fig1 ] . for practical purposes , we do not calculate the density @xmath79 since this quantity is too small for small eigenvalues ( see ( [ aha ] ) ) . we rather calculate @xmath80 , where @xmath81 is the boltzmann weight associated to the eigenvalue @xmath82 , which is calculated as the average of the boltzmann weight operator @xmath83 over the eigenstate associated to the eigenvalue @xmath82 . in the arrhenius picture , we have @xmath84 since @xmath36 , @xmath85 increases exponentially as a function of @xmath86 when @xmath87 . we see on figure [ fig1 ] the existence of the arrhenius regime , in good agreement with the condition @xmath74 . we can now fit the density of eigenvalues in the arrhenius regime to the form @xmath88 and see if we recover @xmath89 , as predicted from the arrhenius law . the fit to the arrhenius law is shown on figure [ fig2 ] . the variations of the exponent @xmath90 as a function of the inverse temperature are plotted on figure [ fig3 ] , and we conclude to a good agreement with the arrhenius law . random walks on a finite @xmath91dimensional lattice with a broad distribution of waiting times have been analyzed in the past ( see @xcite and references therein ) . as we shall see , we can relate our results concerning the existence of a arrhenius  like regime for the rem  like trap model to the results of @xcite , where the long time behavior of random walk on such lattices was investigated by means of the real space renormalization method introduced in @xcite . the rem  like trap model falls in the universality class described by machta since the ingredients are the same in the two cases , rather than the universality class of random walks in a random field force @xcite . in order to make contact with the work of machta , we first calculate the probability @xmath92 of being at the origin : @xmath93 where @xmath94 is the probability to find the system in the trap @xmath18 at time @xmath7 , starting from the trap @xmath18 at time @xmath95 , and @xmath96 is the set of eigenvalues of the glauber matrix . we have shown that , in the infinite dimensional case , the typical small eigenvalue density ( relevant to the long time behavior of @xmath97 ) has the arrhenius form @xmath98 . a simple scaling analyzis then shows that @xmath99 in the long time limit . this result was derived for a graph such as any site is connected with any other site , that is we started from the begining from an infinite dimensional situation . in order to extrapolate this result to a finite dimension @xmath91 larger than the upper critical dimension 2 ( as shown in @xcite ) we first notice that , within the framework of the rem  like trap model , the microscopic time @xmath100 scales like @xmath101 , with @xmath102 the lattice spacing , whereas , in the context of random walks , the correct scaling is @xmath103 [ an illustration of the validity of such a scaling will be given in section [ brownian - motion ] ] . this implies that @xmath104 , which gives the correct prescription to go from the rem  like trap model in infinite dimension to the finite @xmath91 dimensional random walk model : one should replace the time @xmath7 in ( [ p0rem ] ) by @xmath105 . this leads to @xmath106 for @xmath91 larger than the upper critical dimension @xmath107 . this is exactly what was found in @xcite using a real space renormalization group approach , with a non trivial spectral dimension @xmath108 . we aim to calculate the structure of the eigenstates of the rem  like trap model . to do so , we decompose a given microcanonical state @xmath30 into the eigenvectors @xmath109 , and calculate the absolute value of the coefficients of this decomposition . the result is plotted on figure [ fig4 ] . we see that a microcanonical state @xmath30 with an eigenvalue @xmath110 is mainly a linear combination of eigenstates with eigenvalues close to @xmath110 . this observation can be understood in the framework of perturbation theory as follows . the first order correction to the eigenvectors are @xmath111 so that , at first order @xmath112 at the first order of perturbation theory , the pure states @xmath30 is thus mainly a linear combination of eigenvectors with an eigenvalue close to the eigenvalue of @xmath113 .    starting from a pure state @xmath30 at time @xmath6 , we can calculate the evolution of the system : @xmath114 \nonumber , \end{aligned}\ ] ] where we have used ( [ per1 ] ) and ( [ per2 ] ) . as expected from the arrhenius law , the relaxation time is of the order of @xmath115 , which diverges in the limit @xmath116 . we will compare below the structure of the eigenvectors of the rem  like trap model to the structure of eigenvectors of another model ( bm model ) . assuming periodic boundary conditions , the diagonal terms of the glauber matrix are given by @xmath117 the perturbation induced by the off - diagonal terms does not affect significatively the density of eigenvalues provided the perturbation @xmath118 is small compared to the level spacing between @xmath41 and @xmath55 . since @xmath119 is uncorrelated from @xmath14 , the matrix element @xmath120 must be averaged over @xmath119 , so that @xmath40 is of the order of @xmath121 , that is @xmath122 the distribution of eigenvalues of the unperturbed system is @xmath123 from what we deduce the level spacing statistics in the long time ( small @xmath61 ) approximation @xmath124 the off - diagonal terms do not affect the density of eigenvalues provided @xmath125 which leads to the condition @xmath126 . the average number of traps generated among a sequence of @xmath10 independent traps chosen among the distribution @xmath127 and fulfilling the aforementioned condition is @xmath128 which is vanishing in the large @xmath10 limit . we thus conclude that the arrhenius regime with this topology does not exist , that is off - diagonal terms will always modify significatively the density of eigenvalues with respect to the arrhenius case . as we will see below , the arrhenius regime is replaced by a renormalized arrhenius regime , where the density of eigenvalues is given by ( [ lawarr ] ) , but with @xmath90 lower than @xmath129 . again , we generate numerically the density of small eigenvalues @xmath85 . the fit is shown on figure [ fig2 ] , and the variations of the exponent in the arrhenius regime are plotted on figure [ fig3 ] . we conclude to the existence of a renormalized arrhenius regime , in the sense that the density of eigenvalues is given by ( [ lawarr ] ) but with an exponent @xmath90 smaller than @xmath129 , indicating a slowing down of the dynamics due to the one dimensional topology compared to the infinite dimensional case : in one dimension , to go from one trap to a deeper trap , the system has to jump to intermediate traps . section [ relation1 ] was devoted to the connection between the spectrum of the glauber matrix in infinite dimension and the long time behavior of random walks on finite @xmath91 dimensional trapped hypercubic lattices above the upper critical dimension . we now analyze the analogy in one dimension . we have shown numerically that the density of eigenvalues behaves like @xmath130 for small @xmath82 , where @xmath131 is different from @xmath9 . in terms of random walks , this result becomes , with the notations of section [ relation1 ] @xmath132 namely the spectral dimension is @xmath133 . we now use the results of @xcite where @xmath134 was calculated : @xmath135 from what we conclude that the coefficient @xmath90 plotted on figure [ fig3 ] is @xmath136 this analytic form is consistent with numerical calculations of the density of eigenvalues . we now assume that @xmath10 energy traps were generated among the distribution @xmath127 , and that these energies have been ordered in such a way that @xmath56 . a site @xmath66 with a trap of energy @xmath137 is connected to the sites @xmath138 ( with an energy @xmath139 ) and @xmath140 ( with an energy @xmath141 , except for the sites @xmath142 and @xmath10 which have no neighbor . intuitively , such a trap ordering should make the dynamics faster with respect to the case analyzed in section [ disor ] since the system can decrease its energy continuously during the relaxation . we have @xmath143 since the density of level spacings @xmath63 ( [ p(s ) ] ) is large for small values of @xmath61 , we expect that @xmath55 is close to @xmath41 at least if @xmath144 is small , so that we can use the following approximation : @xmath145 . the off - diagonal terms will not modify the density of eigenvalues with respect to the arrhenius case provided @xmath146 that is @xmath75 . the number of traps fulfilling this condition in a distribution of @xmath10 traps is on average equal to unity , so that only a small number of eigenvalues contribute to the arrhenius regime . nevertheless , it is still possible to make a statistics over a large ensemble of @xmath10 traps distribution and to test numerically if the arrhenius regime can be observed . we generated numerically the density of small eigenvalues @xmath85 . as shown on figure [ fig2 ] , we still can fit @xmath147 to the form @xmath148 , and the exponent @xmath90 is close to @xmath129 , as in the arrhenius regime . ( see figure [ fig3 ] ) . we conclude that ordering the traps in one dimension restores the arrhenius regime for the lowest eigenvalues . again , the existence of the arrhenius regime is sample dependent . we are now going to give a general relation between ageing and the glauber matrix . this relation is a generalization of a criterium given in @xcite . in the context of the rem  like trap model , the average trapping time @xmath149 diverges if @xmath150 . on the other hand , ageing exists below @xmath12 only , in the limit where the microscopic time scale @xmath100 goes to zero . we are now going to argue that a similar criterium can be found for any statistical mechanical model , provided the master equation is linear . this is not always the case . for instance , the master equation of the backgammon model in a mean field dynamics at zero temperature is not linear @xcite . let us now start our general argument . the idea is to start with a system in a given microcanonical state @xmath30 and write the glauber dynamics equation for the evolution of this pure state . first , we call @xmath151 the transformation matrix @xmath152 where @xmath30 is a microcanonical state and @xmath153 are the eigenstates of the glauber matrix . by microcanonical state , we mean that if the system is in the state @xmath30 , the system is with probability one in the microcanonical state @xmath5 . in order to avoid confusion , all the labels related to the eigenvectors will be denoted by greek letters whereas the labels for pure states will be denoted by latin letters . the evolution of the pure state ( [ pure2 ] ) is readily obtained as @xmath154 in other words , the probability to be in microcanonical state @xmath155 at time @xmath7 starting from the microcanonical state @xmath5 is @xmath156 we now consider the time  dependent correlation function @xmath157 where @xmath158 is the number of microcanonical states . we use @xmath158 instead of @xmath10 to avoid confusion between the number of microstates @xmath158 and the number of sites @xmath10 in a spin model . in a trap model , @xmath159 , whereas in an ising spin system , @xmath160 . in ( [ correlation ] ) , we have carried out an average over all the initial microcanonical states , that is the system is first quenched from an infinite temperature , and the correlation ( [ correlation ] ) is averaged over the realizations of the infinite temperature states . in ( [ correlation ] ) , @xmath161 is the overlap between the pure states @xmath5 and @xmath155 , which has to be defined within each model . as we will see below in section [ brownian - motion ] , even within a given dynamics , the choice of the overlap depends on what type of observables are considered in the two times correlation functions . we now want to calculate from ( [ correlation ] ) a typical time by integrating ( [ correlation ] ) . however , the boltzmann eigenstate leads to a trivial divergence of this integral . we thus need a prescription to eliminate the equilibrium distribution in ( [ correlation ] ) . we first discuss this prescription in the absence of metastability and postpone for later discussion the case of metastability . we define the typical time @xmath162 by @xmath163 where the summation over the eigenstates excludes the boltzmann distribution . the overline denotes an average over the disorder realizations , if necessary . the significance of @xmath0 is clear : for a finite size system with @xmath10 sites , the out  of  equilibrium dynamics occurs before @xmath0 and , for times larger than @xmath0 , the finite size system has equilibriated . in other words , @xmath0 is the cross  over time before which the finite size system exhibits ageing . notice that in this argument , the system is a finite size system . one should thus carefully study finite size effects . two cases may occur : first , @xmath0 tends to a constant in the thermodynamic limit . in this case , ageing is interrupted @xcite in the thermodynamic limit . as we shall see , this is the case for the ising chain at a finite temperature . on the other hand , if @xmath0 diverges in the thermodynamic limit , ageing is not interrupted . in this case , the exitence of a divergent relaxation time with respect to times @xmath95 and @xmath7 implies the existence of a divergent relaxation time with respect to any finite waiting time @xmath2 and @xmath164 , and thus the existence of ageing in the thermodynamic limit . the relation between the existence / absence , the nature of ageing in the thermodynamic limit and the behavior of @xmath0 was derived on heuristic basis . we could not find a rigourous derivation using a general argument . however , in what follows , we investigate the behavior of several models ( rem  like trap model , random walks , bm model , one dimensional ising chain and sk model ) and we find in all these cases that our criterium based on the behavior of @xmath0 is correct . as far as non interrupted ageing in the thermodynamic limit is concerned , we will see below that this is the case for the bm and sk models . this is also the case for the rem  like trap model . we have shown in section [ infinite ] that for all the typical realizations of the trap configurations ( typical is defined in section [ infinite ] ) , the density of eigenvalues is given by the arrhenius law , which means that , in the case of the rem  like trap model , our criterium for the existence of ageing is equivalent to the criterium of bouchaud @xcite that we have recalled at the beginning of this section , provided the overlap between traps is taken to be @xmath165 , where @xmath166 is the edwards  anderson order parameter ( within the one replica symmetry breaking step picture ) . notice that the existence of a non arrhenius regime generated by non typical trap configurations does not modify our conclusion regarding the divergence of @xmath167 since the non arrhenius regime only occurs for traps such as @xmath168 if @xmath11 . the correlation @xmath169 is @xmath170 which is non integrable if @xmath171 . the typical time @xmath172 only depends on the density of eigenvalues @xmath173 for a system of @xmath10 traps , and is independent on the transition matrix , which motivates the study of the density of eigenvalues carried out in section [ gm ] . in general , in any model with only self  overlapping pure states , the existence of ageing is coded only in the behavior of the density of small eigenvalues : there is ageing in the thermodynamic limit if @xmath174 tends to a constant if @xmath82 goes to zero ( an exemple will be provided in section [ bm - model ] ) or if @xmath174 diverges if @xmath82 goes to zero ( this is the case for the rem  like trap model already studied in section [ gm ] ) . let us now discuss the behavior of the correlation function ( [ correlation ] ) in the presence of metastability . we have in mind the case of diamond for instance , which is metastable but in which case there is no ageing . we do not have a glauber matrix description of diamond , but take this case as a generic situation for metastability and draw qualitative conclusions . in this case we expect that for a finite size system , ( [ correlation ] ) will first decay to a finite value over a first time scale @xmath175 . after @xmath175 , the system reaches its metastable configuration ( for instance diamond ) and thermalizes in a portion of its phase space . we thus expect that , for times larger than @xmath175 , the correlation ( [ correlation ] ) reaches a plateau . on larger time scales , the correlation ( [ correlation ] ) should decay to its equilibrium value with a cross  over time @xmath176 , corresponding for instance to the transition from diamond to graphite . the time scale @xmath176 should increase with the number of sites , so that , with a macroscopic number of sites , no transition from diamond to graphite is observable on experimental time scales . on the other hand , @xmath175 may remain small even in the thermodynamic limit , signaling the absence of ageing in the graphite phase . in such a situation where the cross  over time scale @xmath176 is not observable , the only relevant time scale is @xmath175 . in order to calculate @xmath175 one should not use the prescription ( [ prescription ] ) but rather substract to the correlation ( [ correlation ] ) its plateau value for @xmath177 . brownian motion was shown in @xcite to be an ageing phenomenon . this model is simple and thus provides a simple test of the criterium of section [ general ] . brownian motion in an euclidian space of dimension @xmath91 has the following glauber dynamics : @xmath178 and @xmath179 if @xmath18 and @xmath19 are nearest neighbors on a hypercubic lattice . @xmath100 is the microscopic hopping time . in this case , the glauber matrix is symmetric and , because of translational invariance , its eigenstates are plane waves @xmath180 where @xmath181 ( we have assumed cyclic boundary conditions with a lattice spacing @xmath102 ) . the matrix elements of the glauber matrix on the plane wave basis are @xmath182 where @xmath90 runs over the @xmath91 directions and @xmath183 is the component of the wave vector @xmath184 on the @xmath90 direction . the long time behavior of the random walk is readily obtained by expanding the cosine up to the second order , and the probability to find the brownian particle at time @xmath2 on site @xmath185 is , in the long time limit : @xmath186 where the particle is at the origin at time @xmath6 ( out of equilibrium state ) . as noticed first in @xcite , the two times correlation functions of the position are proportional to the waiting time : @xmath187 in the large @xmath2 limit . in order to make contact with @xcite , one should adopt the scaling @xmath188 in the thermodynamic / long time limits , which allows to recover the langevin dynamics results and thus the existence of ageing in any dimension for the position  position correlation function .    in order to make contact between the correlation ( [ c(1 ) ] ) and the general form of the correlation ( [ correlation ] ) , we should choose the overlap @xmath189 . with this choice of the overlap , ( [ correlation ] ) and ( [ c(1 ) ] ) are both related to the same physical quantity : the mean square displacement . using this specific form for the overlap , one can calculate @xmath190 , integrate with respect to @xmath7 and show that @xmath191 , in agreement with the existence of ageing in this system for the position correlations . one could also be interested in the following correlation : @xmath192 in the large @xmath7 , @xmath2 limit , this correlation is found to be @xmath193 in the thernodynamic limit , @xmath194 , @xmath188 , so that @xmath195 is vanishing in the thermodynamic limit , signaling the absence of ageing for the correlation ( [ c(2 ) ] ) . the overlap associated to the correlation ( [ c(2 ) ] ) is the local one : @xmath196 . we can thus calculate @xmath197 using the integration ( [ taurem ] ) over the density of eigenvalues . in order to calculate the behavior of @xmath0 , we consider a finite system of volume @xmath198 with a lattice spacing @xmath102 . the number @xmath10 of sites is @xmath199 . by thermodynamic limit , we mean that @xmath200 is fixed , @xmath194 and @xmath188 . using ( [ taurem ] ) , we have then @xmath201 where the summation is restricted over the wave vectors @xmath184 such as @xmath202 . clearly , the dominant contribution to ( [ tau(2 ) ] ) comes from the smallest eigenvalues only . it is thus legitimate to expand the eigenvalues up to the quadratic order in @xmath184 : @xmath203 one then obtains @xmath204 if @xmath205 , @xmath206 since @xmath188 , @xmath207 in the thermodynamic limit in two dimensions . now if @xmath208 , @xmath209 since @xmath188 , this expression goes to zero in the thermodynamic limit whatever the dimension . we have thus shown that , in the case of random walks , and in the case of the correlations ( [ c(1 ) ] ) and ( [ c(2 ) ] ) , the existence / absence of ageing in the thermodynamic limit is consistent with the general argument of section [ general ] . we assume that the traps have been ordered such as @xmath56 . unlike the rem  like trap model , we do not make any asumption about the trap distribution . the finite temperature dynamics is @xmath210 if @xmath17 and satisfies detailed balance . the zero temperature form of the glauber matrix of this model is very simple since @xmath211 , @xmath212 if @xmath213 and @xmath214 if @xmath215 at zero temperature . the spectrum is @xmath216 . the eigenvector @xmath217 for the eigenvalue @xmath22 is @xmath218 . it is easy to show that the eigenvector @xmath219 associated to the eigenvalue @xmath220 is @xmath221 we can now express the pure states in terms of eigenvectors @xmath222 this means that , at zero temperature , a pure state @xmath223 [ @xmath66 is the @xmath66-ieth traps , with @xmath26 ] overlaps only with eigenstates associated to eigenvalues @xmath225 , and this overlap does not depend on the eigenstate . compared to the result of section [ eigi ] , the structure of the decomposition of a pure state into eigenstates of the glauber dynamics is very different from the rem  like trap model . if one starts from a pure state @xmath226 at time @xmath6 , the state of the system at time @xmath7 is given by @xmath227 we recover the fact that this model can reach its ground state in a time of order @xmath10 at zero temperature , whereas in the limit of zero temperatures , and for the rem  like trap model , the time spent in a trap @xmath18 is proportional to @xmath228 , which goes to infinity if @xmath229 and @xmath10 fixed . the two times autocorrelation functions were calculated in @xcite by means of a laplace transform formalism . the result is @xmath230 if @xmath231 . we rederive this result in the appendix using a treatment different from @xcite , and based only on the diagonalization of the glauber matrix . the bm dynamics thus exhibits ageing . provided only self overlaps are non zero ( @xmath232 ) , the correlation function introduced in section [ general ] is @xmath233 , provided @xmath234 . the typical time @xmath235 below which ageing occurs is thus logarithmically divergent . on the other hand , this can be also shown directly from the spectrum : @xmath236 which diverges logarithmically . as far as the finite temperature behavior of the two times autocorrelation function is concerned , this model exhibits interrupted ageing @xcite since at any finite temperature the correlation length is finite , of the order of @xmath237 , ageing does not exist after a time independent of the system size . it is thus interesting to apply to this case the general argument of section [ general ] . we first begin with the definition of the glauber dynamics . we use the dynamics introduced in @xcite . the master equation of this single spin flip dynamics is @xmath238 as a consequence of detailed balance , @xmath239 notice that , unlike the case of the rem  like trap or bm models , there is no need to multiply the matrix elements by @xmath240 . it is shown in @xcite that in the infinite temperature limit , the eigenvalues of the glauber matrix are negative integers between @xmath241 and @xmath95 , with degeneracies given by the binomial coefficients . the largest relaxation time in the infinite temperature limit is equal to unity whatever the system size , and there is thus no need to rescale the coefficients of the glauber matrix .      in the case of a spin model , the overlap @xmath242 is given by @xmath243 since ageing is interrupted , we should observe a saturation of the time @xmath244 as @xmath10 increases . moreover , @xmath244 should increase as the temperature decreases . in practise , we work with an open chain and we take advantage of the existence of the global @xmath245 symmetry ( invariance of the glauber dynamics under the @xmath246 transformation in a zero magnetic field ) , and the reflection symmetry of the chain . the glauber matrix is diagonalized in each of these four symmetry sectors . the results are plotted on figure [ fig5 ] . the difference @xmath247 decreases as @xmath10 increases , in agreement with the existence of interrupted ageing . this model was proposed in 1975 as an `` exactly solvable '' model of spin glasses @xcite . for reviews on this model and its developments , we refer the reader to the reviews @xcite . the hamiltonian of the sk model is @xmath248 where the sum is carried out over all the pairs of sites ( infinite dimensional model ) and the quenched random interactions are distributed according to a gaussian distribution @xmath249 this model is a spin glass if @xmath250 and a paramagnet otherwise . the two times correlation function were shown in @xcite to exhibit ageing in the glassy phase . we are going to apply our general argument of section [ general ] to the case of the rem  like trap model . to do so , we calculate numerically @xmath251 for a large number of disorder realizations . the result is plotted on figure [ fig6 ] , as well as the differences @xmath252 for @xmath253 sites . we observe that @xmath254 increases at low temperatures and decreases at high temperature . the crossing point on the insert of figure [ fig6 ] should be identified to the glass temperature . this temperature is not equal to unity since the system sizes are very small . notice that even though the sizes are quite small , a large number of disorder configurations is required to have a significant statistics ( of the order of 500000 independent disorder configurations ) . the computer time required was equivalent to 1000 hours of cpu time on a single processor . our results are thus consistent with the existence of ageing in the low temperature phase , and no ageing in the high temperature phase . we have thus carried out a detailed analysis of ageing from a glauber dynamics point of view , with an emphasis on finite size effects . we have shown how to extract information from the glauber matrix about ageing , and how to analyze finite size effects . in some simple cases with only self overlaps , the behavior of the time after which there is no ageing only involves the knowledge of the average density of eigenvalues @xmath173 when @xmath87 . this is the case for the rem  like trap model , in which case we have carried out a careful analysis of the density of small eigenvalues . our analysis was confirmed by numerical diagonalizations of glauber matrices of small clusters . in this case , we were able to related the typical behavior of the density of eigenvalues to the spectral dimension of the corresponding diffusion process .    : r.m . acknowledges interesting discussions with a. barrat , j.p . bouchaud , s. franz and n. wingreen . the authors acknowledge b. douot for his interest to this work and his encouragements . r.m . also acknowledges the hospitality of nec research institute at princeton where part of this work was done . we are going to calculate the two times autocorrelation function @xmath255 using a method based only on the exact diagonalization of the glauber matrix at a zero temperature . starting from an infinite temperature state at time @xmath6 , @xmath256 with @xmath30 a pure state , the natural definition of the two times autocorrelation function within this model is @xmath257 where @xmath258 is the probability to find the system in the trap number @xmath18 at time @xmath2 and @xmath259 is the escape rate from trap @xmath18 , that is the diagonal element of the glauber matrix @xmath260 . this definition of the waiting time dependent autocorrelation function is similar to the one used in @xcite . using the expression ( [ eig ] ) of the eigenvectors and the decomposition ( [ dec1 ] ) and ( [ dec2 ] ) of the pure states into eigenvectors , we can calculate in a straightforward fashion the evolution of the system , starting from the state ( [ initial ] ) : @xmath261 |e_k \rangle .\ ] ] the probabilities @xmath258 are easily recovered : @xmath262 . after straightforward calculations , the dominant behavior of @xmath255 is obtained , in the limit @xmath263 : @xmath264 which is equivalent to the result of @xcite . notice that if @xmath2 is of the order of the number of traps , the system has reached equilibrium : the dynamics is stationary ( the system is stuck in the lowest trap ) and there is no ageing . m. alba , m. ocio , j. hamman , europhys . lett . , 45 ( 1986 ) ; j. phys . ( 1985 ) l-1101 ; m. alba , j. hamman , m. ocio , p. refregier , j. appl . phys . , 3683 ( 1987 ) ; e. vincent , j. hamman , m. ocio , _ recent progress in random magnets _ , d.h . ( world scientific , singapore , 1992 ) . m. mzard , g. parisi and m. virasoro , _ spin glass theory and beyond _ , world scientific ed . ( 1993 ) ; k.h . fisher and j.a . hertz , _ spin glasses _ , cambridge studies in magnetism ( 1991 ) ; v.s . dotsenko , m.v . feigelman and l.b . ioffe , _ spin glasses and related problem _ , physics reviews , hardwood academic publishers ( 1990 ) . density of eigenvalues in the arrhenius regime for the rem  like trap model and fit to the exponential form ( @xmath265,@xmath266 traps ) . the insert shows the same quantity , but with also the smallest eigenvalues . for these small eigenvalues , the arrhenius picture is not valid . density of small eigenvalues of the models of section [ gm ] , for @xmath267 . ( a ) : infinite dimensional case ( rem  like trap model ) ; ( b ) : one dimensional case , with `` disordered '' traps ; ( c ) : one dimensional case with `` ordered '' traps . the lines in cases ( b ) and ( c ) are parallel , and in good agreement with the arrhenius law , whereas we observe a significant deviation to the arrhenius law in ( b ) case . exponent @xmath90 deduced from the fit ( [ lawarr ] ) of @xmath85 . the circles correspond to the infinite dimensional case ( rem  like trap model ) , the squares to the one dimensional `` disordered '' case , and the diamonds to the `` ordered '' one dimensional case . the solid line corresponds the the arrhenius regime ( in other words , to site  disordered random walks above the upper critical dimension 2 ) , and the dashed line corresponds to the one dimensional renormalized arrhenius regime ( in other words , to site  disordered random walks in one dimension ) . structure of the eigenstates of the rem  like trap model for 20 traps . the label @xmath18 is an eigenvalue label , with eigenvalues in growing order and @xmath19 labels the traps ( the depth increases if @xmath19 increases ) . ( a ) represents @xmath268 for the @xmath269 pure states . ( b ) represents the depth of the traps . ( c ) represents the 20 eigenvalues between @xmath270 and @xmath95 for 4 different temperatures ( circles : @xmath271 , squares : @xmath267 , diamonds : @xmath272 , triangles : @xmath273 ) .
Abstract:

Answer: the glauber dynamics of various models ( rem  like trap models , brownian motion , bm model , ising chain and sk model ) is analyzed in relation with the existence of ageing . from a finite size glauber matrix , we calculate a time @xmath0 after which the system has relaxed to the equilibrium state . the case of metastability is also discussed . if the only non zero overlaps between pure states are only self  overlaps ( rem  like trap models , bm model ) , the existence or absence of ageing depends only on the behavior of the density of eigenvalues for small eigenvalues . we have carried out a detailed numerical and analytical analysis of the density of eigenvalues of the rem  like trap models . in this case , we show that the behavior of the density of eigenvalues for typical trap realizations is related to the spectral dimension of the equivalent random walk model .    16.8cm-.4cm-1 cm    plus 0.2pt minus 0.1pt


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: in recent years the inner structure of dark matter haloes has been the topic of some discussion . observations suggest that the dark matter distribution in disk galaxies has a roughly constant density core , with a typical size of order a few kpc . a recent analysis suggests a shallow power - law mass density distribution @xmath1 with @xmath2 @xcite . on the other hand , numerical simulations based on the ( @xmath0)cold dark matter ( cdm ) paradigm suggest a very steep inner mass density distribution , a so - called `` cusp '' @xcite [ nfw ] . the most recent simulations suggest slopes with @xmath3 @xcite . there is strong consensus that this slope can not be as shallow as the observations suggest . the rotation curves of low surface brightness ( lsb ) galaxies are considered particularly clean tests of the cdm paradigm at galaxy scales . lsb galaxies are dominated by dark matter @xcite , and their dynamics should thus give a more or less direct map of the dark matter distribution . lsb galaxy rotation curves , together with the missing dwarf or substructure problem ( e.g.  @xcite ) , currently form the most serious challenges for the cdm model . there has been much debate on this topic in recent years . this is however not the place to review this debate , and the reader is referred to the summaries by e.g.  @xcite and @xcite and references therein .    on the theoretical side the `` universal mass density profile '' ( with its characteristic @xmath4 cusp ) as proposed by @xcite , has long been a corner - stone of cdm ( but see @xcite ) . recent simulations presented in @xcite and @xcite [ hereafter h04 ] show that simulated cdm haloes may actually contain a range of slopes in their inner regions ( though the slopes remain steep ) . the shallowest inner slope measured is @xmath5 , with an average over the simulated haloes of @xmath6 ( at the innermost reliably resolved radius for galaxy mass haloes ) . with this scatter in mind , h04 argue that simple fitting functions ( such as `` the '' nfw halo ) do not capture the full variety and diversity in shapes and slopes , and suggest that this may partly explain the discrepancy between observations and simulations . they claim that the inner slope is difficult to constrain , both observationally as well as theoretically and point out that rotation curve constraints are strongest where numerical simulations are least reliable . this , as well as the potential for systematic effects in the observations , has , at least according to h04 , led to an `` unwarranted emphasis '' on the value of the inner slope , at the cost of evaluation of the data and simulations over their full radial extent . h04 present a novel method to fit rotation curves derived from their simulations as well as lsb galaxy rotation curves from the literature . their method is able to describe the variety present in data and simulations and from the distributions of respective fitting parameters h04 conclude that the observed curves _ are _ consistent with the cdm paradigm ; a conclusion that contradicts a large body of observational work .    for their analysis , h04 make use of the rotation curves presented in @xcite , @xcite and @xcite , and derive circular velocity curves for haloes from their simulations . to quantify both sets of observed and simulated rotation curves they use a three - parameter fitting formula ( see e.g.  @xcite ) :    @xmath7    here @xmath8 is the asymptotic velocity of the flat part of the rotation curve , @xmath9 is a scale radius ( the transition radius between the rising and flat part of the rotation curve ) and the parameter @xmath10 describes the abruptness of the turn - over between the rising and flat parts of the rotation curve . a higher value of @xmath10 results in an abrupt transition , a low value in a more gradual turn - over . note that this function was designed to fit observed rotation curves : it specifically assumes a solid - body rise in the inner parts and a flat rotation curve in the outer parts . h04 fit the two sets of rotation curves with this function , constraining the fit parameters slightly by insisting that    1 . @xmath11 , 2 . @xmath12 , 3 . @xmath13 ,    with @xmath14 the maximum observed velocity in the rotation curve . fits with @xmath15 [ rejected by condition ( 1 ) ] correspond with very abrupt transitions between the rising and flat parts of the curve . condition ( 3 ) prevents some of the solid - body rotation curves where the asymptotic velocity is not well - constrained , from producing fits with extremely large @xmath8 . conditions ( 1)(3 ) thus help removing unrealistic values for the fit parameters . figure [ fig : hayashidata ] ( after h04 , their fig . 9 ) shows the results presented in h04 . the lsb rotation curves show a broad distribution in @xmath10 with two pronounced peaks : one around @xmath16 and a narrower one at @xmath17 . some 70 per cent of the fits have @xmath18 . the @xmath10 distribution of the simulated cdm haloes is markedly different . it is narrow and centered around @xmath19 with a dispersion of @xmath20 .    at first sight the two distributions seem quite incompatible , and one might conclude that the lsb rotation curves are not consistent with the cdm simulations . h04 argue however that for most of the lsb galaxy rotation curves the reduced @xmath21 distribution is quite broad and shallow , so that by constraining the fit parameters slightly more , a fit can be found that is consistent with the simulations , at the cost of only a small increase in @xmath22 . they define these cdm - compatible constraints as :    1 . @xmath23 2 . @xmath12 3 . @xmath13 4 . @xmath24 . here @xmath25 is the average density within the radius where the velocity reaches half its maximum value ( alam et al 2002 ) . see h04 for a justification of these additional constraints , here it is sufficient to note that for the lsb galaxies @xmath25 is derived from the fit of eq . ( 1 ) to the rotation curves , whereas @xmath26 is the predicted value for a galaxy with @xmath27 in the @xmath0cdm model presented in @xcite and @xcite ( see also fig .  2 in @xcite and fig . 11 in h04 ) . the constraint on @xmath25 ensures that only fits are considered with average densities close to those predicted in a standard @xmath0cdm universe . h04 find that for the majority ( @xmath28 70 per cent ) of the galaxies fits can be found for which both the best - fit @xmath29 [ using constraints ( 1)-(3 ) ] and the constrained @xmath30 [ using constraints ( i)-(iii ) ] . this group of galaxies ( group a using the terminology in h04 ) is thus considered to be consistent with cdm . for the remaining galaxies no cdm - compatible fit could be found ( group b in h04 ) , or the curves were too irregular to be well fitted by eq .  ( 1 ) ( group c in h04 ) . the main claim in h04 is thus that most of the lsb galaxies can still be fitted reasonably well with these extra constraints and without a large increase in @xmath22 . they conclude , based on the agreements of the constrained , cdm - compatible @xmath10 and @xmath25 distributions with those of the simulations , that `` this sample of lsb rotation curves is not manifestly inconsistent with the predictions of @xmath0cdm cosmological models . '' this conclusion thus contradicts most observational work on lsb rotation curves which argues the opposite .    in this paper i discuss these conclusions and investigate the method with which they are derived . i find that the method presented in h04 can not distinguish between shallow and steep mass - density slopes , and illustrate this by showing that with this method even core - dominated , pseudo - isothermal halos would be inferred to be consistent with cdm . i introduce an additional criterion that does take the inner mass - density slope into account , and find that only a quarter of the lsb galaxies investigated are possibly consistent with cdm . however , most of these galaxies fit parameters are so weakly constrained that this is not a very strong conclusion . only 3 out of the 51 galaxies investigated have tightly constrained solutions that are not significantly inconsistent with cdm . however , two of them are high surface brightness galaxies that are likely dominated by their stellar population . the net result is that in the total sample there is only one galaxy that is likely to be dark matter dominated and not significantly inconsistent with cdm . lsb galaxies that are consistent with cdm seem to be rare indeed .    in sect . 2 i derive an analytical expression for the slope of a rotation curve , based on eq . 3 compares the slopes in both data and simulations at identical physical radii . in sect . 4 i show that the constraints imposed by h04 are not sufficient to prove or disprove agreement with cdm . in sect .  5 a new constraint is introduced that does take the slope into account . the results are summarised in sect . equation  ( 1 ) is a flexible , analytical function that can describe a wide range of rotation curves ( see @xcite for more examples ) . one can , analogous to h04 , use the fitting function to investigate the properties of the underlying data . assuming a spherical halo dominated by dark matter we can use the inversion used in e.g.  @xcite to derive the corresponding mass - density distribution :    @xmath31    in combination with eq . ( 1 ) yields    @xmath32.\ ] ]    the logarithmic slope of the mass density distribution can now be derived from eq . ( 3 ) :    @xmath33    the logarithmic , or power - law slope @xmath34 thus _ depends on both @xmath10 and @xmath9_. fig . [ fig : courslope ] shows @xmath34 for a number of values of @xmath10 and @xmath9 , where the values have been chosen to cover the range shown by the lsb rotation curve fits from h04 . changing @xmath10 has the effect of changing the shape of the slope - radius relations ; changing @xmath9 shifts the relations horizontally without affecting the shape for any particular @xmath10 . most of the plotted slopes are remarkably _ shallow _ at small radii . an inner asymptotic value of @xmath35 is predominant . the cdm slope @xmath5 can only be reproduced with eqs . ( 1 ) and ( 4 ) for small values of @xmath9 and/or extremely small values of @xmath10 . similarly , at large radii the slope converges to @xmath36 . these are not naturally occurring values for @xmath0cdm mass profiles , with asymptotic slopes @xmath37 in the inner parts and @xmath38 in the outer parts . equation ( 1 ) is really built to mimic the same striking features of rotation curves that led to the introduction of the pseudo - isothermal ( iso ) halo , namely a solid - body rise in the inner parts , and a constant ( `` flat '' ) rotation velocity in the outer parts . one needs to make very specific choices of fitting parameters to make this function resemble an nfw profile . the h04 simulations predict steep slopes , even at the smallest reliably resolved radii , defined by the so - called convergence radius @xmath39 @xcite . for the h04 dwarf galaxy models @xmath40 kpc or @xmath20 kpc for @xmath41 . the slope at @xmath39 in the simulations varies between @xmath5 and @xmath42 . equation  ( 4 ) can be used to derive the slopes implied by the fits to the observed lsb rotation curves . for this one needs to choose a radius at which the slopes will be evaluated . in the rest of this paper two choices will be used . firstly , the radius of the innermost measured point of the rotation curve as given in @xcite and @xcite will be considered . here a value of @xmath43 is used ( converted using the appropriate distance ) , which is the typical spacing between data points for the curves presented there . ] . one might argue that at @xmath44 a different radius is probed for each rotation curve , so , secondly , a constant radius @xmath45 kpc @xmath46 kpc for all curves will be used as well . one should not expect dramatically different results from both choices though : the average value of @xmath47 is @xmath48 kpc , very much comparable with @xmath39 . it should be kept in mind that even though we evaluate the slope at these particular choices for the radius , the @xmath10 and @xmath9 parameters ( and therefore the slope ) are constrained by the _ entirety _ of the curve . figure [ fig : hayslope3 ] shows the distribution of inferred mass density slopes derived using the lsb rotation curves and the best - fit @xmath10 and @xmath9 parameters listed in table 2 of h04 . the slopes for both @xmath49 and @xmath50 kpc are shown . it is clear that the distributions are heavily biased towards shallow slopes and very different from the values found in the simulations at these radii @xmath51 . the small values of @xmath52 derived in h04 indicate that eq . ( 1 ) fits the data well , and the problem is thus not with the quality of the fits . in fig . [ fig : hayslope3 ] is also indicated the distribution of slopes for the galaxies that h04 claim are consistent with cdm ( their group a with @xmath53 and @xmath54 ) . this distribution is not markedly different from that of the entire sample . from fig . [ fig : hayslope3 ] one could thus conclude that the majority of slopes are inconsistent with the cdm prediction . how can this be reconciled with the h04 conclusions ?    as discussed above , h04 add an extra constraint to the @xmath10-parameter , and argue that the resulting agreement with the range of @xmath10-values shown by the simulations implies consistency with cdm . this point will be explored in more detail below , but one can already see here that this does not solve the problem of the discrepancy between observed and simulated slopes . consider the distribution of slopes of the subset of group a for which @xmath55 ( and both @xmath56 ) , i.e. galaxies for which the best fit is ( apparently ) already consistent with cdm . the distribution of slopes of this sub - group is also plotted in fig . [ fig : hayslope3 ] . the selection criteria for this sub - group favour galaxies with steeper slopes , and it should thus come as no surprise that the dominant peak at @xmath57 has disappeared . nevertheless it is remarkable that the distributions are still dominated by slopes @xmath58 . the distribution using @xmath59 peaks at @xmath60 . the distribution using @xmath61 kpc peaks at @xmath62 . as we are probing the slope at @xmath63 kpc and not in the centre , we should not expect all iso haloes to have a flat slope at this radius . the distribution of slopes found is consistent with that expected at @xmath45 kpc for an ensemble of pseudo - isothermal ( iso ) halos with core - radii between @xmath64 and @xmath65 kpc . the decrease in mass - density at this radius is thus much less steep than expected for a cdm halo . the steepest slope expected for a realistic iso halo ( i.e.  consistent with the smallest value @xmath66 kpc measured in real lsb galaxies ; see results in @xcite ) is @xmath67 , with the large majority of iso haloes having slopes less steep than that . the slopes measured in the h04 simulations at @xmath45 kpc vary between @xmath68 and @xmath42 , firmly inconsistent with the observed distributions . the reason for the discrepancy can best be appreciated by realizing that according to eq . ( 4 ) the slope is a function of both @xmath10 and @xmath9 . this is shown graphically in fig . [ fig : gammartfits ] which shows iso - slope contours in the @xmath69 plane , as evaluated at @xmath63 kpc . it is immediately obvious that @xmath70 is not a unique function of @xmath10 , but depends equally strongly on @xmath9 . most combinations of @xmath10 and @xmath9 with @xmath71 and @xmath72 yield slopes shallower than @xmath73 . only a very small part of parameter space results in steep , cdm - compatible slopes ( essentially only the area with @xmath74 and/or @xmath75 ) . constraining @xmath76 does not actually constrain the slope to steep , cdm - like values , but still allows values up to @xmath77 . over - plotted in fig . [ fig : gammartfits ] are the best - fit @xmath10 and @xmath9 values of the lsb rotation curves as determined by h04 . a distinction is made between the entire sample ( groups a , b and c ) , the `` cdm - compatible '' group a , and the sub - group discussed above where the best fit is apparently already cdm - compatible . the distribution of the points is not very different from group to group , though the @xmath78 subgroup tends to have slightly steeper slopes , as already shown in fig . [ fig : hayslope3 ] . regardless of how the total sample is divided up , the very large majority of the galaxies have slopes inconsistent with the values @xmath37 derived in the h04 simulations . note again that the slopes in simulations and observations are compared at _ identical radii _ , and that the @xmath10 and @xmath9 values were derived using the _ entirety of the rotation curves_. the extra constraints introduced by h04 are thus insufficient conditions for agreement with cdm , as constraining to @xmath76 does not automatically imply a mass - density slope consistent with cdm . unless one is prepared to allow for shallow slopes in cdm at these radii ( making it inconsistent with its own simulations ) an extra constraint on the slope is needed to ensure agreement with cdm adding this extra constraint to the four introduced by h04 is only necessary if these prior constraints are insufficient to make a unique distinction between cdm and non - cdm models . if , for example , rotation curves known to be inconsistent with cdm are not rejected as such , the prior constraints are clearly insufficient to address cdm - related questions . one can test the strength of these constraints by `` observing '' a sample of simulated iso haloes . this will show that the h04 constraints are not sufficient as they imply that the majority of these iso haloes , all with easily detectable shallow slopes , are consistent with cdm . iso haloes form of course the one model known to be certainly inconsistent with cdm ( if they were not , there would be no cusp / core controversy ) . the modeling procedure as described in @xcite [ hereafter db03 ] is used . the reader is referred to that paper for an extensive description of these models . the aim there was to estimate the impact of systematic effects on the observations . for that purpose many different rotation curves were simulated and `` observed '' after the addition of observational uncertainties due to e.g.  resolution effects and random motions . db03 used the same datasets as the basis of their models as h04 use for their analysis , and the db03 procedure can be directly applied here to test the fitting method and make comparisons with the h04 results . the iso halo parameters are chosen in such a way that the resulting curves resemble minimum disk iso fits to the observed rotation curves . fig .  10 from @xcite shows that mass models of lsb galaxies have asymptotic velocities @xmath79 ) < 2.2 $ ] ( between @xmath80 and @xmath81 km s@xmath82 ) . central dark matter core densities vary between @xmath83 ) < 0 $ ] ( between @xmath84 and @xmath85 @xmath86 pc@xmath87 $ ] ) . the core radius @xmath88 follows from these two parameters , and has a range @xmath89 ) < 0.83 $ ] ( between @xmath90 and @xmath91 kpc ) . for the models a random @xmath92 and @xmath93 are chosen from the respective logarithmic intervals . note that these intervals encompass the entire observed parameter range . no account is taken of any correlations that might exist between these parameters , nor of any other uncertainties that might be present in the determination of halo parameters . this is reflected in the range of core - radii , which is somewhat larger than that found from observations ( typically @xmath94 kpc ) . the choice of parameter ranges is therefore liberal , but sufficient for our purposes : this is not an attempt to model individual galaxies but a _ reductio ad absurdum _ proof of principle .    for each model halo an inclination , resolution and sampling interval are chosen randomly from the observed distributions described in db03 . error - bars are assigned to each sampled point , again modeled on the observed distribution of uncertainties , with a minimum value of 4 km s@xmath82 . the curve is then corrected for inclination with a corresponding increase in the size of the error - bars . to each data point a random velocity component between @xmath95 and @xmath96 km s@xmath82 is added . this is slightly simpler than the method used in db03 , but results in similarly sized random motions . lastly , as real rotation curves do not all extend out to the same radius , for each galaxy a random outer radius between 3 and 20 kpc is defined , out to which the galaxy is `` observed '' . figure [ fig : panel ] shows a small selection of the artificial iso curves . over - plotted is the eq . ( 1 ) fit using the best - fit constraints ( 1)-(3 ) , as given in sec .  1 . the reduced @xmath21 and the fit parameters are also given in the figure . figure  [ fig : modelgamma ] shows the distribution of the best - fit @xmath10 parameters for 100 different realisations of the model . some @xmath97 per cent of the iso galaxies have @xmath18 . this number should be compared with h04 where it was found that @xmath98 per cent of the lsb rotation curves had @xmath18 . the @xmath10-distribution of lsb galaxies is at least in this respect consistent with that of iso haloes . [ fig : modelslopes ] shows a histogram of the slopes at @xmath63 kpc , both for the input iso models and the `` observed '' output models , the latter derived using the best - fit @xmath10 and @xmath9 parameters . both the input and output slopes show a clear peak at @xmath35 with a tail towards steeper slopes where the data are probing the edge of the dark matter core kpc and not in the center , we do not expect to find flat slopes for all iso haloes . the steepest iso slopes expected are nevertheless significantly less steep than the predicted cdm values . cf .  the remarks in sect .  3 the output distribution has a less pronounced peak at @xmath35 , and a longer tail towards steeper slopes , but as the discussion in @xcite shows , this is fully consistent with the effects of the limited resolution used to `` observe '' the input distribution . the range in slopes found is slightly larger than that exhibited by real lsb galaxies shown in fig . [ fig : hayslope3 ] . this is entirely due to the larger range in core - radii used . fits of eq . ( 1 ) to a set of model iso rotation curves `` observed '' under the same conditions as the samples of lsb galaxies under investigation here are thus able to retain and retrieve the presence of shallow inner slopes .    in order for the h04 criteria to make a valid discrimination between cdm and non - cdm models they thus need to be able to reject the majority of the iso models presented here . applying these criteria to the iso models and again fitting using eq . ( 1 ) , one finds that reducing the upper limit of @xmath99 to @xmath100 has only a modest effect on @xmath22 and almost the entire iso sample can still be fitted reasonably well with eq . ( 1 ) and constraints ( i)-(iii ) , as can be seen from the representative examples in fig . [ fig : panel ] . cdm constraint ( iv ) turns out to be the most restrictive one . it rejects about @xmath101 per cent of the simulated iso haloes with @xmath102 ( note that h04 show that a similar fraction of the lsb galaxies have incompatible densities as well ) . the bottom row of fig . [ fig : panel ] shows a few of the rejected models . these do not look any different from the cdm - compatible iso curves : condition ( iv ) merely removes some fraction of the iso haloes that do not have a cdm - compatible density , without any particular selection against specific values of @xmath10 , @xmath9 , @xmath8 or other relevant parameters . the most surprising outcome though is that of this sample of random iso haloes almost 60 per cent match the h04 cdm criteria . a comparison of the distribution of slopes in this cdm - compatible sample of haloes with the total distribution shows that there has been no selection against shallow slopes ( fig . [ fig : slopescomp3 ] ) . of the total sample of 100 observed iso models , @xmath103 per cent has an observed slope @xmath104 . once condition ( iv ) is applied , 57 haloes remain , of which 82 per cent have best - fit slopes @xmath105 . even when the slopes are derived using the constrained `` cdm compatible '' fit parameters , one still finds that 70 per cent has a shallow slope @xmath104 . as the range in core - radii used for the models is slightly larger than found observationally , the fractions mentioned here should be considered lower limits . regardless of how the slope is defined , shallow slopes are still profusely present ( and detectable ! ) in the models . a comparison with the input iso halo parameters shows that the rejected haloes are not abnormal in any way . in fig . [ fig : rho0 ] @xmath25 is plotted against the input iso halo central density @xmath93 . there is a clear relation , as might be expected : the average density of a halo is related to its central density . galaxies are plotted with different symbols depending on whether they meet the @xmath25 constraint or not . it is clear that the halo models that do not meet the cdm criteria are only guilty of having average densities that are too high or too low . they are not degenerate , nor do they suffer from e.g.  larger error bars . a more finely - tuned choice of halo parameters could bring the rejection rate down significantly , but as the input ranges of @xmath93 and @xmath106 are based on observed distributions , it is more likely that real galaxies have a wider range of densities than predicted by cdm , and condition ( iv ) merely selects from a much broader spectrum of galaxy properties those galaxies that happen to lie in the cdm - compatible range . figure [ fig : chi2 ] compares the values of @xmath22 for the best - fit parameters and the cdm - compatible parameters . haloes that were rejected by the @xmath25 criterion are indicated in the right - hand side of the plot . following h04 one can divide the haloes in three categories : _ ( group a ) : _ rotation curves that are well fit by eq .  ( 1 ) and for which a good fit with apparently cdm - compatible parameters can also be found , i.e. @xmath107 and @xmath108 . _ ( group b ) : _ rotation curves that are well fit by eq . ( 1 ) but for which no good fit with cdm - compatible parameters can be found , usually because of a discrepant value of @xmath25 ; these curves have @xmath109 and @xmath110 or undefined ( i.e. , with @xmath25 out of range ) . _ ( group c ) : _ rotation curves that are poorly fit by eq . ( 1 ) , i.e. @xmath111 and @xmath110 . the latter category only contains a few galaxies , as there are no large - scale asymmetries in the model rotation curves , and most rise smoothly and monotonically with radius . note the numerical agreement between the observational @xmath112 values presented in h04 and those of the iso models , showing that the additional uncertainties introduced in those models are a good representation of the true observational uncertainties . models in group a are , according to h04 , `` consistent with @xmath0cdm in terms of their inferred central densities and the shape of their rotation curves . '' one now has to add several caveats to this statement . firstly , even with @xmath113 and @xmath114 a steep slope is not guaranteed . the constraints used to define group a are not sufficient to filter out the shallow slopes that are still present and detectable in the iso models , even after the addition of observational uncertainties . secondly , as the inner slopes are well - defined and not necessarily steep , the shapes of the rotation curves are also not necessarily consistent with @xmath0cdm . constraints ( i)-(iii ) thus seem to have had hardly any effect in improving the agreement with cdm . the only restrictive constraint is condition ( iv ) , the @xmath25 criterion , but even this seems to have only selected models with a certain range in average density from amongst a much larger range of models . this is illustrated in fig . [ fig : bullock ] where a comparison is made between average @xmath25 densities of the iso models from group a and b with those expected in a @xmath0cdm universe . galaxies in groups a and b are indicated separately . in the h04 interpretation the iso halos from group a are consistent with cdm in terms of shape and average density , while group b is consistent in terms of shape but not density , even though groups a and b both consist of pure iso haloes , with easily detectable shallow slopes . the cdm constraints ( i)-(iv ) fail to take this into account and are thus insufficient to distinguish between cdm and non - cdm models . they can not exclude the possibility that lsb galaxies have shallow inner mass - density slopes , even at @xmath115 . additional constraints are needed to make a unique identification of galaxies compatible with cdm . now that it is established that constraints ( i)-(iv ) are not sufficient to uniquely identify cdm - compatible galaxies , the next step is to identify additional conditions that can . figures  [ fig : courslope ] and [ fig : gammartfits ] show that the inner slope and therefore the shape of the rotation curve do not depend on @xmath10 alone , but that there is an equally strong dependence on @xmath9 . these figures also explain why h04 came to the conclusion that their four constraints were sufficient : implicit in their analysis is the assumption that the shape of the rotation curve depends on @xmath10 alone . while it is true that for @xmath116 kpc a smaller @xmath10 implies a steeper slope , this does not ensure a steep slope in an absolute sense . for @xmath76 , depending on the value of @xmath9 , slopes as shallow as @xmath117 are still possible . furthermore , it is shown above that the h04 method does not take these shallow slopes into account , and thus allows models with easily detectable and prominent shallow slopes purpose - built in to pass the cdm test with flying colours .    in order to identify galaxies that are really compatible with cdm , an additional constraint is needed , in such a way that only @xmath118 combinations are allowed that yield a steep slope [ as well as pass conditions ( i)-(iv ) ] . in practice this means that @xmath69 parameter space needs to be searched for the minimum @xmath21 value that still results in a steep slope , one needs to take that parameter into account as well . this is done by stepping through values of @xmath8 at fixed @xmath10 and @xmath9 , and choosing the value that yields the lowest @xmath21 . therefore @xmath119 $ ] . ] . in order to compare directly with the simulations the slopes are again evaluated at @xmath63 kpc . at this radius the h04 simulations show slopes @xmath120 , and strictly speaking one ought to restrict the search of @xmath121 space to slopes @xmath122 . however , to take into account uncertainties in data and simulations , and give the cdm models as much leeway as possible , a more liberal range of @xmath123 will be used . [ fig : slopeonly ] shows that this is the steepest slope we can expect at this radius from an ultra - compact iso halo with @xmath124 kpc , and this value thus cleanly separates the iso and cdm domains .    at this point a short digression on the significance of @xmath112 is useful . in testing the goodness - of - fit of a model , the reduced @xmath21 , that is , the total @xmath21 divided by the degrees of freedom of the data , is often used . a value of @xmath125 is then taken as an indication that the model is a good description of the data . the reduced @xmath21 is an appropriate statistic when comparing different models to a single set of data , or data with a similar number of degrees of freedom . an example would be the testing of an nfw model versus an iso model for a single rotation curve , or a set of curves with a similar number of data points . h04 determine for each lsb rotation curve the reduced @xmath21 for the best - fitting @xmath118 model , as well as the minimum @xmath112 value using their cdm constraints . they argue that for most of the galaxies the cdm constraints only cause a small increase in the _ reduced _ @xmath21 , and that therefore the cdm - constrained models are not much worse than the best - fit models . this is not necessarily so . remember that the reduced @xmath21 is the @xmath21 divided by the degrees of freedom of the data ( the number of data points minus the number of parameters in the model ) . it is easy to see that in order to cause the same change in _ reduced _ @xmath21 the change in _ absolute _ @xmath21 must be much larger for a rotation curve with a large number of data points than for a rotation curve with a small number of data points . in other words , a small change in @xmath22 can be very significant ( that is , unlikely to be consistent with the uncertainties in the data points ) for a curve with a large number of data points , while it can be insignificant ( entirely consistent with the uncertainties ) for a curve with only a small number of data points . the observed sample we are dealing with here contains a large number of rotation curves , with a wide range in the number of data points , varying from @xmath126 to over 350 . the relevant parameter when judging whether a change in fit parameters yields a model that is still consistent with the data is thus _ not _ the reduced @xmath21 , but the _ change in absolute @xmath21_. the confidence level of any parameter combination is given by the corresponding change in @xmath21 with respect to the minimum value . for a model with three free parameters @xmath127 as used here , the @xmath128 confidence intervals for any fit are given by @xmath129 , respectively . using these criteria , a rotation curve can be defined to be consistent with cdm if the best - fit values of @xmath10 and @xmath9 lie within @xmath130 of the closest @xmath118 combination with a @xmath131 slope . in order for the fit to give any meaningful constraints , one should furthermore require that the area enclosed by the @xmath130 confidence level is a small fraction of the total area of parameter space investigated is not crucial . we have tested the analysis with a @xmath132 level as well . though the number of badly constrained fits increases , the relative proportions of cdm consistent _ vs _ inconsistent fits does not change appreciably . ] . if the latter is not the case , then any conclusions regarding agreement or disagreement with cdm are only very weak at best .    here the group a rotation curves , defined in h04 to be consistent with their definition of cdm , are investigated again using the new constraints . groups b and c are already inconsistent with cdm using the more relaxed h04 constraints , and will therefore not be investigated further . for each rotation curve from group a the minimum @xmath21 was found using eq . ( 1 ) , and the change in @xmath21 necessary to be consistent with a slope @xmath123 was determined . the @xmath133-confidence interval corresponding to this @xmath134 was also calculated , as well as the ratio @xmath135 of the area enclosed by the @xmath130 contour and the total area of @xmath121 parameter space searched @xmath136 kpc ) . as described above , the ratio @xmath135 is used to determine whether a solution is well - constrained or not . small values of @xmath135 indicate that the @xmath130 area is only a small fraction of total parameter space , and therefore yield significant solutions . large values of @xmath135 indicate that the fit does not constrain models in a significant way . the results are given in table [ table : fit ] , where the galaxies are sorted in order of the number of data - points in each individual rotation curve . [ fig : chi2surface ] shows the @xmath128 confidence contours for a number of representative galaxies as described below . from studying the plots of the @xmath21 contours for all galaxies , it was found that a value @xmath137 is a reasonable discriminator between well - constrained and badly constrained fits , but the precise value of @xmath135 is not critical ( see fig . [ fig : chi2lim ] for some examples of @xmath138 fits ) . the galaxies in table [ table : fit ] fall in four distinct categories :    1 .   _ well - constrained solutions , consistent with cdm : _ objects with best fits within @xmath130 of a steep slope @xmath123 parameter combination , and with a well - constrained @xmath139 solution . an example is ngc3274 , as shown in fig . [ fig : chi2surface ] . this category consists of 3 out of the 36 galaxies in group a. 2 . _ well - constrained solutions , inconsistent with cdm : _ objects with best fits more than @xmath130 away from a steep slope @xmath123 parameter combination , and with a well - constrained @xmath139 solution . an example is ugc 11819 , as shown in fig . [ fig : chi2surface ] . this category contains 17 out of 36 galaxies . badly constrained solutions , consistent with cdm : _ objects with best fits within @xmath130 of a steep slope @xmath123 parameter combination , but with a badly constrained @xmath140 solution . an example is f563 - 1 , as shown in fig . [ fig : chi2surface ] . this category contains 12 out of 36 galaxies . badly constrained solutions , inconsistent with cdm : _ objects with best fits more than @xmath130 away from a steep slope @xmath123 parameter combination , and with a badly constrained @xmath140 solution . an example is ddo 185 , as shown in fig . [ fig : chi2surface ] . this category contains 4 out of 36 galaxies . it is clear that only galaxies from the first two categories , i.e.  with @xmath139 , can help constrain halo models . a glance at table [ table : fit ] shows that most of the solutions that are consistent with cdm are actually badly constrained . this is investigated further in fig . [ fig : resolve ] . here @xmath135 is plotted against the number of data points in each rotation curve . a distinction is made between rotation curves that are consistent with cdm and those that are inconsistent with cdm using the constraints introduced here . table [ table : fit ] and fig . [ fig : resolve ] lead to several conclusions . firstly , of the 15 galaxies that are consistent with cdm , 13 have less than 30 data points . furthermore , of these 13 galaxies , only 1 has @xmath139 ( i.e.  is well - constrained ) . looking at the 21 galaxies that are inconsistent with cdm , only 10 have less than 30 data points . of these 10 galaxies , all except one have @xmath139 . rotation curves with more than 30 data points tend to be better constrained on average , as well as have a tendency to be inconsistent with cdm . it is also interesting to note that even amongst galaxies with only a limited number of data points , the most constrained solutions ( small @xmath135 ) tend to be those inconsistent with cdm . the majority of solutions @xmath141 per cent ) that are consistent with cdm have only a small number of data points _ and _ are ill - constrained . these ill - constrained solutions therefore do not decide between any model one way or the other , and statements that these galaxies are consistent with cdm are thus not very strong . of the well - constrained solutions the large majority is inconsistent with cdm . furthermore , of the 13 galaxies with more than 30 data points only 2 are consistent with cdm . these two galaxies ( ngc 3274 and ngc 4455 ) also happen to have the smallest optical scale - lengths and highest surface brightnesses from all galaxies with available photometry in the @xcite , @xcite and @xcite samples : both galaxies have @xmath142 mag arcsec@xmath143 . these galaxies are therefore unlikely to be dominated by dark matter in their inner parts and it is very likely that the steep slopes in these two galaxies simply reflect the stellar mass distribution ( see also the discussion on ngc 3274 in sec .  9.2.2 of @xcite ) .    amongst the rotation curves with less than 30 points there is only one galaxy with @xmath144 that seems consistent with cdm . this is ugc 731 , a bona - fide lsb dwarf , and is the best and only candidate in the entire observed sample for a cdm - consistent , dark - matter - dominated galaxy with a well - constrained solution . the level of agreement is not perfect though , ugc 731 is consistent with a @xmath145 slope at the @xmath146 confidence level , with the confidence level decreasing to @xmath147 for an @xmath4 slope , which is still only the upper limit of the range of slopes expected from the simulations . once the steep slope constraint is added to the h04 constraints there are thus only 3 galaxies left that have a well - constrained fit consistent with a steep slope and @xmath0cdm . of these 3 galaxies , 2 are of very high surface brightness and are likely to be dominated by stars in the inner parts , leaving only one possible candidate with a well - constrained cdm - compatible fit . the large majority of the well - constrained fits strongly prefer shallow slopes . the 2@xmath133 limit used here already indicates that for many galaxies it is difficult to adjust the fit parameters to get a cdm - compatible fit . restrictive as this criterion is , it really only gives an upper limit to the true probability . while one has the freedom to adjust the parameters of an individual fit to test the likelihood of an hypothesis for one particular rotation curve , this becomes problematic when dealing with many fits . if every fit is individually adjusted to fit a particular model , does that still constitute a proper test ? an analogue can perhaps be found in some of the first supernova results that indicated the existence of a cosmological constant ( e.g.@xcite ) . the majority of the supernovae in that early work were each individually consistent within @xmath130 with the _ _ non-__existence of @xmath0 . even though the data showed a systematic trend , each measurement could have been corrected individually to fit the no-@xmath0 model . however , the joint result of these small inconsistencies has formed the basis for the discovery of the effects of @xmath0 on the expansion of the universe . while it would be presumptuous to compare those results with the current analysis , it does show that it is important to pay attention to the global conclusion that is forced upon us by the data . one thing that is clear is that galaxies that are significantly consistent with @xmath0cdm in terms of their inferred central densities and the shapes of their rotation curves seem to be rare indeed . i have investigated the claim by h04 that the majority of rotation curves of lsb galaxies are not inconsistent with cdm . i have shown that the method used is unable to distinguish between shallow and steep slopes , even when these are obviously present and easily detectable . i illustrate this by showing that with the h04 method one would infer even pseudo - isothermal , core - dominated haloes to be consistent with cdm . an extra constraint based on the mass - density slope is thus needed to make more definitive statements on any ( dis)agreement with cdm . using this extra criterion i show that only a quarter of the lsb galaxies can be said to be consistent with cdm , as opposed to the three quarters found by h04 . the majority of these cdm - consistent lsb galaxies do however have fit parameters that are so ill - constrained , that they are consistent with virtually anything . restricting the analysis to galaxies with well - constrained solutions , i find that only 3 out of the 20 well - constrained galaxies are possibly consistent with cdm ; two of these are high surface brightness dwarf galaxies that are likely to be dominated by stars .    in summary , a comparison of the circular velocity profiles of cdm haloes with rotation curves of lsb galaxies indicates that the shapes and inferred central densities of most lsb galaxies are mostly _ inconsistent _ with those of simulated haloes within the limitations imposed by observational error .      alam , s.  m.  k. , bullock , j.  s. , weinberg , d.  h.  2002 , , 572 , 34 binney , j.  2004 , iau symposium , 220 , 3 bullock , j.  s. , kolatt , t.  s. , sigad , y. , somerville , r.  s. , kravtsov , a.  v. , klypin , a.  a. , primack , j.  r. , & dekel , a.  2001 , , 321 , 559 courteau , s.  1997 , , 114 , 2402 de blok , w.  j.  g.  2004 , iau symposium , 220 , 69 de blok , w.  j.  g. , bosma , a.  2002 , , 385 , 816 de blok , w.  j.  g. , mcgaugh , s.  s.  1997 , , 290 , 533 de blok , w.j.g . , mcgaugh , s.s . , rubin , v.c . 2001 , , 122 , 2396 de blok , w.  j.  g. , mcgaugh , s.  s. , bosma , a. , & rubin , v.  c.  2001 , , 552 , l23 de blok , w.  j.  g. , bosma , a. , & mcgaugh , s.  2003 , , 340 , 657 fukushige , t. , kawai , a. , & makino , j.  2004 , , 606 , 625 ghigna , s. , moore , b. , governato , f. , lake , g. , quinn , t. , & stadel , j.  2000 , , 544 , 616 hayashi , e. , navarro , j.f . , power , c. , jenkins , a. , frenk , c.s . , white , s.d.m . , springel , v. , stadel , j. , quinn , t.r . 2004 , , 355 , 794 jing , y.  p. , mo , h.  j. , borner , g. , & fang , l.  z.  1995 , , 276 , 417 klypin , a. , kravtsov , a.  v. , bullock , j.  s. , & primack , j.  r.  2001 , , 554 , 903 mcgaugh , s.  s. , rubin , v.  c. , & de blok , w.  j.  g.  2001 , , 122 , 2381 moore , b. , governato , f. , quinn , t. , stadel , j. , & lake , g.  1998 , , 499 , l5 moore , b. , ghigna , s. , governato , f. , lake , g. , quinn , t. , stadel , j. , & tozzi , p.  1999 , , 524 , l19 navarro , j.f . , frenk , c.s . , white , s.d.m . 1996 , , 462 , 563 navarro , j.f . , frenk , c.s . , white , s.d.m . 1997 , , 490 , 493 navarro , j.  f. , et al . 2004 , , 349 , 1039 power , c. , navarro , j.  f. , jenkins , a. , frenk , c.  s. , white , s.  d.  m. , springel , v. , stadel , j. , & quinn , t.  2003 , , 338 , 14 riess , a.  g. , et al . 1998 , , 116 , 1009 swaters , r.  a. , madore , b.  f. , van den bosch , f.  c. , & balcells , m.  2003 , , 583 , 732 wechsler , r.  h. , bullock , j.  s. , primack , j.  r. , kravtsov , a.  v. , & dekel , a.  2002 , , 568 , 52      f563 - 1 & 10 & 1.92 & 0.67 & 0.54 & m01 + eso 4250180 & 11 & 24.94 & 0.25 & 4.31 & m01 + eso 4880490 & 12 & 1.07 & 0.27 & 0.28 & m01 + eso 3020120 & 12 & 1.72 & 0.62 & 0.48 & m01 + ddo 189 & 14 & 0.75 & 0.79 & 0.18 & db02 + eso 1870510 & 15 & 0.02 & 0.45 & 0.0011 & m01 + ugc 4115 & 15 & 0.36 & 0.87 & 0.064 & m01 + f568 - 3 & 16 & 43.27 & 0.13 & 5.99 & m01 + f583 - 1 & 17 & 15.12 & 0.43 & 3.14 & m01 + f571 - 8 & 19 & 93.60 & 0.10 & @xmath148 & m01 + f563 - 1 & 19 & 0.02 & 0.56 & 0.0011 & db02 + ugc 11583 & 19 & 4.31 & 0.57 & 1.20 & m01 + eso 1200211 & 19 & 2.09 & 0.45 & 0.59 & m01 + f583 - 4 & 19 & 19.67 & 0.13 & 3.72 & m01 + eso 2060140 & 20 & 9.57 & 0.06 & 2.28 & m01 + ugc 11819 & 23 & 181.61 & 0.07 & @xmath148 & m01 + f730-v1 & 24 & 18.66 & 0.04 & 3.60 & m01 + ugc 4173 & 24 & 0.00 & 0.62 & 0.00 & db02 + ugc 11616 & 26 & 105.32 & 0.05 & @xmath148 & m01 + ddo 47 & 26 & 0.05 & 0.91 & 0.0040 & db02 + ugc 731 & 28 & 4.95 & 0.11 & 1.35 & db02 + ugc 11454 & 28 & 218.96 & 0.05 & @xmath148 & m01 + ugc 5005 & 28 & 1.47 & 0.87 & 0.40 & db02 + ugc 1230 & 34 & 20.27 & 0.31 & 3.79 & db02 + ugc 628 & 37 & 8.92 & 0.05 & 2.17 & db02 + ddo 185 & 42 & 19.23 & 0.44 & 3.67 & db02 + ugc 10310 & 49 & 32.78 & 0.20 & 5.09 & db02 + ugc 3371 & 57 & 11.68 & 0.57 & 2.63 & db02 + ngc 3274 & 66 & 0.00 & 0.01 & 0.00 & db02 + ugc 3137 & 72 & 24.17 & 0.15 & 4.23 & db02 + ngc 4455 & 80 & 1.63 & 0.10 & 0.45 & db02 + ugc 4325 & 84 & 76.41 & 0.24 & @xmath148 & db02 + ddo 52 & 86 & 41.8 & 0.02 & 5.87 & db02 + ngc 100 & 92 & 21.16 & 0.03 & 3.90 & db02 + ngc 5023 & 149 & 95.96 & 0.01 & @xmath148 & db02 + ngc 4395 & 355 & 23.51 & 0.01 & 4.16 &
Abstract:

Answer: a recent study has claimed that the rotation curve shapes and mass densities of low surface brightness ( lsb ) galaxies are largely consistent with @xmath0cdm predictions , in contrast to a large body of observational work . i demonstrate that the method used to derive this conclusion is incapable of distinguishing the characteristic steep cdm mass - density distribution from the core - dominated mass - density distributions found observationally : even core - dominated pseudo - isothermal haloes would be inferred to be consistent with cdm . this method can therefore make no definitive statements on the ( dis)agreement between the data and cdm simulations . after introducing an additional criterion that does take the slope of the mass - distribution into account i find that only about a quarter of the lsb galaxies investigated are possibly consistent with cdm . however , for most of these the fit parameters are so weakly constrained that this is not a strong conclusion . only 3 out of 52 galaxies have tightly constrained solutions consistent with @xmath0cdm . two of these galaxies are likely dominated by stars , leaving only one possible dark matter dominated , cdm - consistent candidate , forming a mere 2 per cent of the total sample . these conclusions are based on comparison of data and simulations at identical radii and fits to the entire rotation curves . lsb galaxies that are consistent with cdm simulations , if they exist , seem to be rare indeed .