the study of stars is one of the most fascinating topics in astrophysics because it involves many areas of physics : gravitation , thermodynamics , hydrodynamics , statistical mechanics , quantum mechanics , relativity , ... and it has furthermore a very interesting history . in his classical monograph _ the internal constitution of the stars _ , eddington ( 1926 ) @xcite lays down the foundations of the subject and describes in detail the basic processes that govern the structure of ordinary stars . at that stage , only elements of `` classical '' physics are used and difficulties with such theories to account for the structure of high density stars such as the companion of sirius are pointed out .
soon after the discovery of the quantum statistics by fermi ( 1926 ) @xcite and dirac ( 1926 ) @xcite , fowler ( 1926 ) @xcite uses this `` new thermodynamics '' to explain the puzzling nature of white dwarf stars .
he understands that low mass white dwarf stars owe their stability to the quantum pressure of the degenerate electron gas . ) .
the case of a partially degenerate self - gravitating fermi gas at arbitrary temperature has been discussed more recently by hertel & thirring ( 1971 ) @xcite and chavanis ( 2002 ) @xcite in the context of statistical mechanics .
they describe the phase transition , below a critical temperature @xmath9 , from a gaseous configuration to a condensed state ( with a degenerate core surrounded by a halo ) .
this provides a simple physical mechanism showing how the system can reach highly degenerate configurations as a result of gravitational collapse ; see chavanis ( 2006 ) @xcite for a review on phase transitions in self - gravitating systems . ]
the resulting structure is equivalent to a polytrope of index @xmath10 so that the mass - radius relation of classical white dwarf stars behaves like @xmath11 ( chandrasekhar 1931a @xcite ) .
the next step was made by chandrasekhar , aged only nineteen , who was accepted by the university of cambridge to work with fowler . in the boat that took him from madras to southampton @xcite
, chandrasekhar understands that relativistic effects are important in massive white dwarf stars and that einstein kinematic must be introduced in the problem . in his first treatment ( chandrasekhar 1931b @xcite )
, he considers the ultra - relativistic limit and shows that the resulting structure is equivalent to a polytrope of index @xmath12 .
applying the theory of polytropic gas spheres ( emden 1907 @xcite ) , this leads to a unique value of the mass @xmath13 that he interprets as a limiting mass , nowdays called the chandrasekhar limit .
quantum mechanics can not prevent the system from collapsing to a point .
however , he did not take this collapse very seriously and noted that , in reality , there exists stars with mass @xmath14 that do not show this `` ridiculous tendency '' , so that they must possess regions in which the laws of quantum mechanics are violated . ] the complete mass - radius relation of relativistic white dwarf stars was given later ( chandrasekhar 1935 @xcite ) and departs from fowler s sequence as we approach the limiting mass . at this critical mass , the radius of the configuration vanishes .
above the critical mass , the equation of state of the relativistic degenerate fermi gas of electrons is not able to balance gravitational forces and , when considering the final evolution of such a star , chandrasekhar ( 1934 ) @xcite `` left speculating on other possibility '' .
chandrasekhar s result was severly criticized by eddington ( 1935 ) @xcite who viewed this result as _ a reductio ad absurdum _ of the relativistic formula and considered the combination of special relativity and non - relativistic quantum theory as an `` unholly alliance '' .
because of the querelle with eddington , it took some time to realize the physical implication of chandrasekhar s results . however , progressively , his early investigations on white dwarf stars were extended in general relativity to the case of neutron stars ( oppenheimer & volkoff 1939 @xcite ) and finally led to the concept of `` black holes '' , a term coined by wheeler in 1967 . without the dispute with eddington , this ultimate stage of matter resulting from gravitational collapse could have been predicted much earlier from the discovery of chandrasekhar @xcite .
although these historic elements are well - known , it is less well - known that the concept of a maximum mass for relativistic white dwarf stars had been introduced earlier by anderson ( 1929 ) @xcite and stoner ( 1930 ) @xcite .
these studies are mentioned in the early works of chandrasekhar but they have been progressively forgotten and are rarely quoted in classical textbooks of astrophysics .
these authors investigated the equation of state of a relativistic degenerate fermi gas and predicted an upper limit for the mass of white dwarf stars .
stoner ( 1930 ) @xcite uses a uniform mass density to model the star while chandrasekhar ( 1931 ) @xcite considers a more realistic @xmath12 polytrope .
however , as noted in chandrasekhar ( 1931 ) @xcite , the value of the limiting mass found by stoner with his simplified model is relatively close to that obtained with the improved treatment .
interestingly , nauenberg ( 1972 ) @xcite introduced long after a simplified treatment of relativistic white dwarf stars in order to obtain an analytical approximation of the mass - radius relation .
it turns out that this model , which gives a very good agreement with chandrasekhar s numerical results ( up to some normalization factors ) , is similar to that introduced by stoner .
leaving aside these interesting historical remarks , the result of chandrasekhar @xcite concerning the existence of a limiting mass is very profound because this mass can be expressed in terms of fundamental constants , similarly to the bohr radius of the hydrogen atom @xcite .
hence , the mass of stars is determined typically by the following combination @xmath15 where @xmath16 is the constant of gravity , @xmath17 the planck constant , @xmath18 the velocity of light and @xmath19 the mass of the hydrogen atom ( @xmath20 is the solar mass ) .
this formula results from the combination of quantum mechanics ( @xmath17 ) , special relativity ( @xmath18 ) and gravity @xmath21 . since dimensional
analysis plays a fundamental role in physics , it is of interest to investigate how the preceding results depend on the dimension of space @xmath1 of the universe . in our previous investigations @xcite , we considered the case of classical white dwarf stars in @xmath1 dimensions and found that they become unstable in a space of dimension @xmath8 . in that case , the star either evaporates or collapses
. therefore , quantum mechanics can not stabilize matter against gravitational collapse in @xmath8 contrary to what happens in @xmath4 @xcite .
interestingly , this is similar to a result found by ehrenfest @xcite at the atomic level for coulomb forces in bohr s model and for the planetary motion ( kepler problem ) .
the object of this paper is to extend these results to the case of relativistic white dwarf stars and exhibit particular dimensions of space which play a special role in the problem when we combine newtonian gravity , quantum mechanics and special relativity .
we shall see that the problem is very rich and interesting in its own right .
it shows to which extent the dimension @xmath4 of our universe is particular , with possible implications regarding the anthropic principle @xcite .
we note that a similar problem has been considered in @xcite on the basis of dimensional analysis .
our approach is more precise since we generalize the _ exact _ mathematical treatment of chandrasekhar @xcite to a space of dimension @xmath1 .
a connection with other works investigating the role played by the dimension of space on the laws of physics is made in the conclusion . in appendix
[ sec_stab ] , we propose different derivations of the stability limits of polytropic spheres and consider applications of these results to classical and relativistic white dwarf stars .
following chandrasekhar @xcite , we model a white dwarf star as a degenerate gas sphere in hydrostatic equilibrium .
the pressure is due to the quantum properties of the electrons and the density to the protons . in the completely degenerate limit ,
the electrons have momenta less than a threshold value @xmath22 ( fermi momentum ) and their distribution function is @xmath23 where @xmath17 is the planck constant .
there can only be two electrons in a phase space element of size @xmath24 on account of the pauli exclusion principle .
therefore , the number of electrons per unit volume is @xmath25 where @xmath26 is the surface of a unit sphere in @xmath1-dimensions .
the mean kinetic energy per electron is given by @xmath27 where @xmath28 is the energy of an electron with impulse @xmath29 .
in relativistic mechanics , @xmath30 the pressure of the electrons is @xmath31 using eq .
( [ es3 ] ) , the pressure can be rewritten @xmath32 finally , the mass density of the star is @xmath33 where @xmath19 is the mass of the proton and @xmath34 the molecular weight . if we consider a pure gas of fermions ( like , e.g. , massive neutrinos in dark matter models ) , we just have to replace @xmath35 by their mass @xmath36 .    introducing the notation @xmath37 , we can write the density of state parametrically as follows @xmath38 where @xmath39 @xmath40 the function @xmath41 has the asymptotic behaviors @xmath42 @xmath43 the classical limit corresponds to @xmath44 and the ultra - relativistic limit to @xmath45 .
explicit expressions of the function @xmath41 are given in appendix [ sec_exp ] for different dimensions of space .
we provide here the explicit expression of the function @xmath41 defined by eq .
( [ es9 ] ) for different dimensions of space @xmath261 and @xmath262 , respectively : @xmath271 @xmath272 @xmath273 @xmath274
for a spherically symmetric distribution of matter , the equations of hydrostatic equilibrum are @xmath46 @xmath47 they can be combined to give @xmath48 expressing @xmath49 and @xmath50 in terms of @xmath51 and setting @xmath52 , we obtain @xmath53 we denote by @xmath54 and @xmath55 the values of @xmath51 and @xmath56 at the center .
furthermore , we define @xmath57 @xmath58 note that the scale of length @xmath59 is independent on @xmath60 for @xmath2 . substituting these transformations in eq .
( [ c4 ] ) , we obtain the @xmath1-dimensional generalization of chandrasekhar s differential equation @xmath61 with the boundary conditions @xmath62 the radius @xmath63 of the star is such that @xmath64 .
this yields @xmath65 the density can be expressed as @xmath66 where the central density @xmath67 finally , we find that the mass is related to @xmath55 by @xmath68 note that @xmath60 does not
_ explicitly _ enter in this expression for @xmath4 but it is of course present implicitly .
in the classical case @xmath44 , we find that the equation of state takes the form @xmath69 with @xmath70 therefore a classical white dwarf star is equivalent to a polytrope of index @xcite : corresponds to the value of the polytropic index of a classical white dwarf star in @xmath4 .
the index @xmath71 refers to its corresponding value in @xmath1 dimensions .
the same convention is adopted for the other indices @xmath72 , @xmath73 and @xmath74 . ] @xmath75 polytropic stars are described by the lane - emden equation @xcite . this can be recovered as a limit of the chandrasekhar equation . for @xmath44 , we have @xmath76 . we define @xmath77 . to leading order , @xmath78 . setting @xmath79 and combining the foregoing results , we find that eq .
( [ c7 ] ) reduces to the lane - emden equation with index @xmath80 : @xmath81 @xmath82 note that the condition at the origin is @xmath83 instead of @xmath84 as in the ordinary lane - emden equation .
however , using the homology theorem for polytropic spheres @xcite , we can easily relate @xmath85 to @xmath86 , the solution of the lane - emden equation with index @xmath87 and condition at the origin @xmath84 .     in an appropriate system of coordinates ( see @xcite for details ) .
complete polytropes correspond to the terminal point in the series of equilibria .
the series becomes dynamically unstable with respect to the euler - poisson system ( saddle point of the energy functional ) after the turning point of mass @xmath88 which appears for @xmath6 .
thus , classical white dwarf stars are stable for @xmath89 and unstable for @xmath8 . ,
width=302 ]    the structure and the stability of polytropic spheres in various dimensions of space has been studied by chavanis & sire @xcite . for @xmath90 ,
this study exhibits two important indices : @xmath91 according to this study , a classical white dwarf star is self - confined ( complete ) if @xmath92 , i.e. @xmath93 . in that case
, the density vanishes at a finite radius @xmath63 identified as the radius of the star . on the other hand , it is nonlinearly dynamically stable with respect to the euler - poisson system if @xmath94 , and linearly unstable otherwise . therefore , classical white dwarf stars are stable for @xmath89 and unstable for @xmath8 ( see fig .
[ alphac ] and appendix [ sec_stab ] ) . for @xmath95 ,
classical white dwarf stars are always self - confined and stable .
using the results of @xcite , the mass - radius relation for complete polytropes with index @xmath87 in @xmath1 dimensions is @xmath96 where we have defined @xmath97 where @xmath98 is such that @xmath99 . using eq .
( [ class2 ] ) , we find that the mass - radius relation for classical white dwarf stars in @xmath1 dimensions is @xmath100 for @xmath101 , the mass @xmath102 decreases as the radius @xmath63 increases while for @xmath103 and for @xmath104 it increases with the radius ( see fig .
[ classique ] ) .     and the mass is independent on radius in @xmath6 . the dimensional factors @xmath105 and @xmath106 are defined in sec .
[ sec_gen ] .
, width=302 ]    for @xmath6 the mass is independent on the radius and given in terms of fundamental constants by @xmath107 we recall that the value of the gravitational constant @xmath16 depends on the dimension of space so that we can not give an explicit value to this limiting mass .
the central density is related to the radius by @xmath108 for @xmath3 , the radius is independent on mass and given in terms of fundamental constants by @xmath109 furthermore , for @xmath3 , the lane - emden equation ( [ class4 ] ) with index @xmath110 can be solved analytically , yielding @xmath111 and @xmath112 where @xmath113 is the first zero of bessel function @xmath114 .
we obtain a density profile @xmath115 where the central density is related to the total mass by @xmath116 finally , for @xmath4 , using the mass - radius relation ( [ class9 ] ) , we find that the average density is related to the total mass by @xmath117 ( for @xmath118 ) . historically , this result was obtained by chandrasekhar ( 1931 ) @xcite who first applied the theory of polytropic gas spheres with index @xmath10 to classical white dwarf stars .
it improves an earlier result @xmath119 obtained by stoner ( 1929 ) @xcite on the basis of his model of stars with uniform density ( see sec .
[ sec_ana ] ) .
in the ultra - relativistic limit @xmath45 , we find that the equation of state takes the form @xmath120 with @xmath121 therefore , an ultra - relativistic white dwarf star is equivalent to a polytrope of index @xmath122 this also directly results from the chandrasekhar equation . for @xmath45
, it reduces to @xmath123 @xmath124 where we have set @xmath125 and @xmath126 .
this is the lane - emden equation with index @xmath127 .     in an appropriate system of coordinates ( see @xcite for details ) .
complete polytropes correspond to the terminal point in the series of equilibria .
the series becomes dynamically unstable with respect to the euler - poisson system ( saddle point of the energy functional ) after the turning point of mass @xmath88 which appears for @xmath4 .
thus , ultra - relativistic white dwarf stars are stable for @xmath128 and unstable for @xmath129 . ,
width=302 ]    using the results of @xcite for @xmath90 , we deduce that an ultra - relativistic white dwarf star is self - confined if @xmath130 , i.e. @xmath131 .
in addition , it is nonlinearly dynamically stable with respect to the euler - poisson system if @xmath132 and linearly unstable otherwise . therefore , ultra - relativistic white dwarf stars are stable for @xmath133 and unstable for @xmath129 ( see fig .
[ alphar ] and appendix [ sec_stab ] ) . for @xmath95
, ultra - relativistic white dwarf stars are self - confined and stable .
on the other hand , using the results of @xcite , the mass - radius relation for complete polytropes with index @xmath127 in @xmath1 dimensions is @xmath134 where we have defined @xmath135 using eq .
( [ ur2 ] ) , we find that the mass - radius relation for ultra - relativistic white dwarf stars in @xmath1 dimensions is @xmath136 for @xmath137 , the mass @xmath102 increases as the radius @xmath63 increases while for @xmath138 it decreases with the radius ( see fig .
[ relativiste ] ) .     and the mass is independent on radius in @xmath4 ( chandrasekhar s mass ) . , width=302 ]    for @xmath2 ,
the radius is independent on mass and given in terms of fundamental constants by @xmath139 furthermore , for @xmath2 , the lane - emden equation ( [ ur4 ] ) with index @xmath140 can be solved analytically , yielding @xmath141 and @xmath142 .
we obtain a density profile @xmath143 where the central density is related to the total mass by @xmath144 for @xmath4 , the mass is independent on radius and given in terms of fundamental constants by @xmath145 this is the chandrasekhar mass @xmath146 coming back to eq .
( [ c7 ] ) , we can show that for this limiting value , the radius @xmath63 of the configuration tends to zero ( see sec .
[ sec_gen ] ) .
historically , the existence of a maximum mass for relativistic white dwarf stars was first published by anderson ( 1929 ) @xcite who considered a relativistic extension of the model of stoner ( 1929 ) @xcite for classical white dwarf stars .
he obtained a limiting mass @xmath147 ( for @xmath118 ) .
the relativistic treatment of anderson was criticized and corrected by stoner ( 1930 ) @xcite who obtained a value of the limiting mass @xmath148 .
the uniform density model of stoner was in turn criticized and corrected by chandrasekhar ( 1931 ) @xcite who applied the theory of polytropic gas spheres with index @xmath12 to relativistic white dwarf stars and obtained the value ( [ ur13 ] ) of the limiting mass @xmath149 .
it seems that these historical details are not well - known because the references to the works of anderson and stoner progressively disappeared from the literature .
collecting together the results of sec .
[ sec_c ] , the mass - radius relation for relativistic white dwarf stars in @xmath1 dimensions can be written in the general case under the parametric form @xmath150 where we have defined @xmath151 and @xmath152 the mass @xmath105 and the radius @xmath106 can be expressed in terms of fundamental constants as @xmath153 @xmath154 we can now obtain the mass - radius curve @xmath155 by the following procedure .
we fix a value of the parameter @xmath60 and solve the differential equation ( [ c7 ] ) with initial condition ( [ c8 ] ) until the point @xmath156 , determined by eq .
( [ c9 ] ) , at which the density vanishes . the radius and the mass of the corresponding configuration
are then given by eq .
( [ gen1 ] ) . by varying @xmath55 ,
we can obtain the full curve @xmath157 parameterized by the value of the central density @xmath158 given by eq .
( [ c11 ] ) . to solve the differential equation ( [ c7 ] )
, we need the behavior of @xmath159 at the origin . expanding @xmath160 in taylor series and substituting this expansion in eq .
( [ c7 ] ) we obtain for @xmath161 : @xmath162 where @xmath163 we note in particular that @xmath164 .    .
the radius increases with the mass .
the configurations are always stable and there exists a maximum radius achieved in the ultra - relativistic limit for @xmath165 . ,
width=302 ]    .
the configurations are always stable and there exists a maximum radius for partially relativistic distributions . , width=302 ]    . the radius decreases with the mass .
the configurations are always stable and there exists a maximum radius achieved in the classical limit for @xmath166 . , width=302 ]    .
the radius decreases with the mass .
the configurations are always stable . , width=302 ]    .
the radius decreases with the mass .
the configurations are always stable and there exists a limiting mass ( chandrasekhar s mass ) achieved in the ultra - relativistic limit for @xmath167 . for @xmath168 ,
quantum mechanics can not arrest gravitational collapse . , width=302 ]    finally , we give the asymptotic expressions of the mass - radius relation . in the classical limit , using eq .
( [ class9 ] ) , we obtain @xmath169 in the ultra - relativistic limit , using eq .
( [ ur8 ] ) , we get @xmath170 in figs .
[ chandra1]-[chandra4 ] , we plot the mass - radius relation of relativistic white dwarf stars ( full line ) for different dimensions of space .
the asymptotic relations ( [ gen8 ] ) and ( [ gen9 ] ) valid in the classical ( c ) and ultra - relativistic ( r ) limits are also shown for comparison ( dashed line ) together with the analytical approximation ( dotted line ) derived in sec .
[ sec_ana ] .    .
there exists a limiting mass for partially relativistic distributions .
classical configurations are stable and the series of relativistic equilibria becomes unstable after the turning point of mass .
thus , highly relativistic configurations in @xmath129 can not be in hydrostatic equilibrium.,width=302 ]    . for @xmath171
the ultra - relativistic configurations ( equivalent to polytropes of index @xmath172 ) are not self - confined anymore and a spiral develops in the mass - radius relation .
this is somehow similar to the classical spiral occurring in the @xmath173 plane in the thermodynamics of isothermal self - gravitating systems @xcite and to the spiral occurring in the @xmath174 plane in the general relativistic treatment of neutron stars @xcite .
the series of equilibria becomes unstable at the first turning point of mass and new modes of instability occur at the secondary turning points @xcite.,width=302 ]    .
there exists a limiting mass achieved in the classical limit for @xmath175 .
in fact , for @xmath8 , all the configurations ( classical and relativistic ) are unstable . quantum mechanics can not balance gravitational attraction even in the classical limit .
, width=302 ]
we present here an analytical approximation of the mass - radius relation in various dimensions of space based on the treatment by nauenberg @xcite in @xmath4 .
this treatment amounts to considering that the density is uniform in the star and the mass - radius relation is obtained by minimizing the energy functional ( see appendix [ sec_ef ] ) with respect to the radius ( or to the density ) at fixed mass . as mentioned in the introduction , this is similar to the simplified model of relativistic white dwarf stars made by stoner @xcite before chandrasekhar s treatment @xcite .
following nauenberg @xcite , we approximate the kinetic energy @xmath176 by the form @xmath177 where @xmath178 is the number of electrons and @xmath29 is an average over the star of the momentum of the electrons . ) of the kinetic energy in his treatment .
this leads to more complicated expressions with , however , qualitatively similar conclusions .
since the models are based on the ( unrealistic ) approximation that the stellar density is homogeneous , the results can not be expected to be more than qualitatively correct .
therefore , the approximation ( [ ana1 ] ) made by nauenberg for the kinetic energy is sufficient for the purposes of this simplified approach .
] we assume that it is determined by an appropriate average value of the density by the relation @xmath179 based on the pauli exclusion principle .
now , for the density , we write @xmath180 where @xmath181 is a dimensionless parameter .
we also write the potential energy in the form @xmath182 where @xmath183 is another dimensionless parameter . by writing eq .
( [ ana4 ] ) , we have assumed that @xmath184 but we shall see that the following results pass to the limit for @xmath185 .
we introduce two dimensionless variables @xmath186 and @xmath187 and two fixed constants @xmath188 and @xmath189 such that @xmath190 we determine @xmath191 and @xmath189 by the requirement that the relativity parameter @xmath192 have the form @xmath193 and that the potential energy can be written @xmath194 this yields @xmath195 @xmath196 comparing with eqs .
( [ gen4 ] ) and ( [ gen5 ] ) , we find that @xmath197 now , the energy @xmath198 of the star can be written @xmath199 in the classical limit @xmath44 , we get @xmath200 and in the ultra - relativistic limit @xmath45 , we obtain @xmath201 we shall consider the energy as a function of the radius @xmath63 , with the mass @xmath102 fixed .
thus , the mass - radius relation will be obtained by minimizing the energy @xmath202 versus @xmath51 with fixed @xmath203 .
writing @xmath204 , we obtain the equations @xmath205 @xmath206 defining the mass - radius relation in parametric form in the framework of the simplified model . note that @xmath207 represents the value of the density in units of @xmath208 , i.e. @xmath209 .
therefore , eq . ( [ ana12 ] ) can be viewed as the relation between the mass and the density . eliminating the relativity parameter @xmath51 ( or density ) between eqs .
( [ ana12 ] ) and ( [ ana13 ] ) , we explicitly obtain @xmath210 @xmath211 @xmath212 @xmath213 in @xmath2 , there exists a maximum radius @xmath214 achieved for @xmath215 , and we have the scaling @xmath216 in @xmath3 , there exists a maximum radius @xmath214 achieved for @xmath217 , and we have the scaling @xmath218 in @xmath4 , there exists a maximum mass @xmath219 ( chandrasekhar s mass ) achieved for @xmath220 , and we have the scaling @xmath221 in @xmath6 , there exists a maximum mass @xmath219 achieved for @xmath222 , and we have the scaling @xmath223 in order to compare eqs .
( [ ana14])-([ana17 ] ) with the exact results , we need to estimate the values of the constants @xmath181 and @xmath183 . this can be done by considering the limiting forms of eqs .
( [ ana12 ] ) and ( [ ana13 ] ) . in the classical limit @xmath44
, we obtain @xmath224 and in the ultra - relativistic limit @xmath45 , we get @xmath225 comparing with eqs .
( [ gen8 ] ) and ( [ gen9 ] ) , we find that @xmath226 where the quantities @xmath227 and @xmath228 can be deduced from the numerical study of the lane - emden equation .
the analytical approximations of the mass - radius relation ( [ ana14])-([ana17 ] ) are plotted in dotted lines in figs .
[ chandra1 ] , [ chandra2 ] , [ chandra3 ] , and they give a fair agreement with the exact results ( full line ) .
of course , they can not reproduce the spiral in @xmath6 which requires the resolution of the full differential equation ( [ c7 ] ) .
let us address the stability of the configurations within this simplified analytical model .
if we view a white dwarf star as a gas of electrons at statistical equilibrium at temperature @xmath229 , stable configurations are those that minimize the free energy @xmath230 at fixed mass , where @xmath202 is the total energy ( kinetic @xmath231 potential ) and @xmath232 is the fermi - dirac entropy ( see appendix [ sec_ef ] ) . for a completely degenerate gas at @xmath0 ,
stable configurations are those that minimize the energy @xmath202 at fixed mass . as we have seen previously , the first order variations @xmath233 determine the mass - radius relation .
then , the configurations is stable if the energy is minimum , i.e. @xmath234 .
now , we have @xmath235 for @xmath133 , all the solutions are stable while for @xmath8 , all the solutions are unstable . for @xmath236 , the solutions with @xmath237 are stable and the solutions with @xmath238 are unstable where @xmath239 we now show that the onset of instability precisely corresponds to the turning point of mass , i.e. to the point where the mass is maximum in the series of equilibria @xmath240 . in terms of reduced variables , this corresponds to @xmath241 , or equivalently @xmath242 . taking the logarithmic derivative of eq .
( [ ana12 ] ) , we have @xmath243 so that the condition @xmath242 yields @xmath244 . thus , instability sets in precisely at the maximum mass as could have been directly inferred from the turning point criterion @xcite .
it is also useful to discuss the classical and the ultra - relativistic limits specifically . in the classical limit , using eq .
( [ ana11class ] ) , we get @xmath245 so that a classical white dwarf star is stable for @xmath89 and unstable for @xmath8 . in the ultra - relativistic limit , using eq .
( [ ana11relat ] ) , we get @xmath246 so that an ultra - relativistic white dwarf star is stable for @xmath133 and unstable for @xmath129 .
this returns the results obtained in secs .
[ sec_class ] and [ sec_ur ] .    finally , considering the function @xmath247 given by eq .
( [ ana11 ] ) and assuming that the system evolves so as to minimize its energy ( this requires some source of dissipation ) , we have the following results : because they are elementary , but drawing them is helpful to visualize the results described in the text . ] for @xmath248 , there exists a stable equilibrium state ( global minimum of @xmath202 ) with radius @xmath63 for all mass @xmath102 . for @xmath4
, there exists a maximum mass @xmath13 so that : ( i ) for @xmath249 , there exists a stable equilibrium state ( global minimum of @xmath202 ) with radius @xmath250 ( ii ) for @xmath251 , the system collapses to a point ( @xmath167 ) but its energy remains finite ( lower bound ) ( iii ) for @xmath14 , the system collapses to a point @xmath252 and @xmath253 . for @xmath254 , there exists a maximum mass @xmath13 so that : ( i ) for @xmath249 , there exists a metastable equilibrium state ( local minimum of @xmath202 ) with @xmath237 and an unstable equilibrium state ( global maximum of @xmath202 ) with @xmath238 .
the system can either reach the metastable state or collapse to a point ( @xmath252 , @xmath253 ) ; the choice probably depends on a notion of basin of attraction ( ii ) for @xmath255 , the system collapses to a point ( @xmath252 , @xmath253 ) . for @xmath256 , there exists a critical mass @xmath13 so that : ( i ) for @xmath249 , there exists an unstable equilibrium state ( global maximum of @xmath202 ) so the system either collapses ( @xmath252 and @xmath253 ) or evaporates ( @xmath175 and @xmath257 ) ( ii ) for @xmath14 , the system collapses to a point ( @xmath252 and @xmath253 ) . for @xmath258 , there exists an unstable equilibrium state ( global maximum of @xmath202 ) for all mass @xmath102 so the system either collapses ( @xmath252 and @xmath253 ) or evaporates ( @xmath175 and @xmath257 ) .
_ `` why is our universe three dimensional ? does the dimension @xmath4 play a special role among other space dimensions ? '' _    several scientists have examined the role played by the dimension of space in determining the form of the laws of physics .
this question goes back to ptolemy who argues in his treatise _ on dimensionality _ that no more than three spatial dimensions are possible in nature . in the @xmath259 century
, kant realizes the deep connection between the inverse square law of gravitation and the existence of three spatial dimensions .
interestingly , he regards the three spatial dimensions as a consequence of newton s inverse square law rather than the converse . in the twentieth century
, ehrenfest @xcite , in a paper called `` in what way does it become manifest in the fundamental laws of physics that space has three dimensions ? ''
argues that planetary orbits , atoms and molecules would be unstable in a space of dimension @xmath8 .
this idea has been followed more recently by gurevich & mostepanenko @xcite who argue that if the universe is made of metagalaxies with various number of dimensions , atomic matter and life are possible only in @xmath260-dimensional space .
other investigations on dimensionality are reviewed in the paper of barrow @xcite .
we have found that the relativistic self - gravitating fermi gas at @xmath0 ( a white dwarf star ) possesses a rich structure as a function of the dimension of space .
we have exhibited several characteristic dimensions @xmath2 , @xmath3 , @xmath4 , @xmath5 , @xmath6 and @xmath7 . for @xmath248 , there exists stable configurations for any value of the mass . for @xmath4 ,
the sequence of equilibrium configurations is stable but there exists a maximum mass ( chandrasekhar s limit ) above which there is no equilibrium state . for @xmath254 ,
the sequence of equilibrium configurations is stable for classical white dwarf stars but it becomes unstable for relativistic white dwarf stars with high density after the turning point of mass .
therefore , the dimension @xmath4 is special because it is the largest dimension at which the sequence of equilibrium configurations is stable all the way long ; for @xmath129 , a turning point of mass appears so that ultra - relativistic white dwarf stars become unstable .
therefore , the dimension @xmath4 is _ marginal _ in that respect .
finally , for @xmath8 , the whole sequence of equilibrium configurations is unstable .
therefore , as already noted in our previous papers @xcite , the dimension @xmath6 is critical because at that dimension quantum mechanics can not stabilize matter against gravitational collapse , even in the classical regime , contrary to the situation in @xmath4 @xcite .
interestingly , this result is similar to that of ehrenfest @xcite although it applies to white dwarf stars instead of atoms .
our exact description of @xmath1-dimensional white dwarf stars based on chandrasekhar s seminal paper @xcite shows that relativistic white dwarf stars become unstable in @xmath129 and that classical white dwarf stars become unstable in @xmath8 .
therefore , for @xmath8 , a self - gravitating fermi gas forms a black hole or evaporates ( see appendix [ sec_stab ] ) .
these conclusions have also been reached by bechhoefer & chabrier @xcite on the basis of simple dimensional analysis .
this suggests that a @xmath8 universe is not viable ( see also appendix [ sec_fourd ] ) and gives insight why our universe in apparently three dimensional .
we note that extra - dimensions can appear at the micro - scale , an idea originating from kaluza - klein theory .
this idea took a renaissance in modern theories of grand unification which are formulated in higher - dimensional spaces .
cm ) and it is a theoretical challenge to explain why only three dimensions are expanded while the others are compact .
we note that a universe with @xmath129 dimensions would be very different from ours since compact objects like white dwarfs or neutron stars would be unstable and replaced by black holes .
in fact , _ all the matter would collapse to a single point in a finite time _
( see appendix [ sec_fourd ] ) .
this observation may be connected to the fact that only three dimensions are extended . ]
our approach shows that already at a simple level , the coupling between newton s equations ( gravitation ) , fermi - dirac statistics ( quantum mechanics ) and special relativity reveals a rich structure as a function of @xmath1 . in this respect , it is interesting to note that the critical masses ( [ class10 ] ) ( [ ur12 ] ) and radii ( [ class11 ] ) ( [ ur9 ] ) that we have found occur for simple _ integer _
dimensions @xmath261 and @xmath262 , which was not granted a priori .
it is interesting to develop a parallel between the mass - radius relation @xmath240 of white dwarf stars and the caloric curve @xmath263 giving the temperature as a function of the energy in the thermodynamics of self - gravitating systems @xcite . in this analogy ,
the chandrasekhar mass in @xmath4 is the counterpart of the critical temperature for isothermal systems in @xmath3 ( in both cases , the equilibrium density profile is a dirac peak containing all the mass ) .
for relativistic white dwarf stars in @xmath4 ( corresponding to a polytrope @xmath264 ) , the critical mass @xmath265 of bacterial populations described by the keller - segel model in @xmath3 ( corresponding to @xmath266 ) and the critical mass @xmath267 or critical temperature @xmath268 of self - gravitating isothermal systems in @xmath3 ( corresponding to @xmath266 ) @xcite . at these critical values
@xmath251 or @xmath269 , the system forms a dirac peak .
this analogy , sketched in sec .
8 of astro - ph/0604012v1 , is addressed specifically in @xcite .
] on the other hand , for @xmath129 , the mass - radius relation for white dwarf stars exhibits turning points , and even a spiraling behavior for @xmath270 , which is similar to the spiraling behavior of the caloric curve for isothermal systems in @xmath4 . in this analogy , the maximum mass , corresponding to a critical value of the central density which parameterizes the series of equilibria , is the counterpart of the antonov energy ( in the microcanonical ensemble ) @xcite or of the emden temperature ( in the canonical ensemble ) @xcite .
the series of equilibria becomes unstable after this turning point .
in addition , there is no equilibrium state above this maximum mass , or below the minimum energy or minimum temperature in the thermodynamical problem . in that case , the system is expected to undergo gravitational collapse .    as a last comment ( notified by the referee )
, it should be emphasized that the conclusions reached in this paper concerning dimensionality implicitly assume that the laws of physics that we know remain the same in a universe of arbitrary dimension @xmath1 .
this is of course not granted at all .
there may be a new cosmological theory , a new theory of star formation and stellar evolution in higher dimensions .
we also emphasize that our approach does not take into account general relativistic effects that can sensibly modify the results @xcite .    _
acknowledgements : _ i am grateful to e. blackman and g. chabrier for pointing out their works after a first version of this paper was placed on astro - ph/0603250 .
special thanks are due to n. barros for interesting discussions .
we generalize in @xmath1 dimensions the usual stability criteria for polytropic gaseous spheres , and apply them to classical and ultra - relativistic white dwarf stars .
let us consider the euler - poisson system @xcite : @xmath275 @xmath276 @xmath277 describing the dynamical evolution of a barotropic gas with an equation of state @xmath278 .
the euler - poisson system conserves the mass @xmath102 and the energy @xcite : @xmath279 in the following , we shall essentially consider the case of a polytropic equation of state @xmath280 where @xmath281 . for @xmath282 ,
the energy functional @xmath283 is given by @xmath284 it can be rewritten @xmath285 for an isothermal equation of state @xmath286 , corresponding to @xmath287 or @xmath215 , we have @xmath288 for a general polytropic equation of state @xmath289 with arbitrary @xmath290 , it is convenient to write the energy of the gas in the form @xmath291 we have added a constant term @xmath292 proportional to the total mass ( which is a conserved quantity ) so as to recover the energy ( [ qw3 ] ) of an isothermal gas in the limit @xmath287 @xcite .
we first consider the linear dynamical stability of a polytropic star described by the euler - poisson system and generalize the approach developed in chavanis @xcite .
we consider a steady solution of the euler - poisson system satisfying @xmath293 and the condition of hydrostatic balance @xmath294 . for a polytropic equation of state @xmath295 , the equilibrium density profile is solution of the @xmath1-dimensional lane - emden equation @xcite .
we consider _ complete polytropes _ such that the density vanishes at a finite radius @xmath63 . for @xmath296
, there exists a unique steady state for any mass @xmath102 .
the mass - radius relation is @xcite : @xmath297 where @xmath298 is a constant given in terms of the solution @xmath299 of the lane - emden equation of index @xmath186 in @xmath1 dimensions . for the critical index @xmath300
, steady state solutions exist only for a unique value of the mass @xmath301 these solutions have the same mass @xmath13 but an arbitrary radius @xmath63 .
we already anticipate that the index @xmath300 will correspond to a case of marginal stability separating stable and unstable solutions .
linearizing the equations of motion ( [ stab1])-([stab3 ] ) around a stationary solution in hydrostatic balance and writing the perturbation in the form @xmath302 , we obtain after some calculations @xcite the eigenvalue equation @xmath303 where we have restricted ourselves to spherically symmetric perturbations and defined @xmath304 . for a polytropic gas with an equation of state @xmath289 ,
the foregoing equation becomes @xcite : @xmath305 the polytrope is stable if @xmath306 ( yielding oscillatory modes with pulsation @xmath307 ) and unstable if @xmath308 ( yielding exponentially growing modes with growth rate @xmath309 ) . considering the point of marginal stability ( @xmath310 ) and introducing the emden variables @xcite , eq .
( [ stab5 ] ) reduces to @xmath311 this equation has the exact solution @xcite : @xmath312 the point of marginal stability is then determined by the boundary conditions .
one can show @xcite that the velocity perturbation is given by @xmath313 .
therefore , if the density of the configuration vanishes at @xmath314 , one must have @xmath315 to avoid unphysical divergences .
thus , if @xmath316 denotes the value of the normalized radius @xmath63 of the star @xcite such that @xmath317 , the natural boundary condition for the eigenvalue equation ( [ stab5 ] ) is @xmath318 . substituting this condition in eq .
( [ stab7 ] ) , we obtain the critical index @xmath319 corresponding to a marginally stable gaseous polytrope ( @xmath310 ) .
our method gives the form ( [ stab7 ] ) of the neutral perturbation @xmath320 at the critical index @xmath300 .
we conclude that the infinite family of steady state solutions with equal mass @xmath13 and arbitrary radius @xmath63 that exists for @xmath300 is marginally stable .
we now consider the nonlinear dynamical stability problem .
since the euler - poisson system ( [ stab1])-([stab3 ] ) conserves the mass @xmath102 and the energy @xmath283 , a maximum or a minimum of the energy functional @xmath321 $ ] at fixed mass @xmath322=m$ ] determines a steady state of the euler - poisson system that is formally nonlinearly dynamically stable @xcite . because of the presence of the kinetic term @xmath323=(1/2)\int \rho { \bf u}^{2}d{\bf r}$ ] , the functional @xmath321 $ ] has no absolute maximum .
thus , we need to investigate the possible existence of a minimum of @xmath324 $ ] at fixed mass @xmath322=m$ ] . a barotropic star that minimizes the energy functional @xmath283 at fixed mass @xmath102 is nonlinearly dynamically stable with respect to the euler - poisson system .
the cancellation of the first order variations @xmath325 , where @xmath326 is a lagrange multiplier , yields @xmath293 and the condition of hydrostatic balance @xmath294 .
then , the condition of nonlinear dynamical stability is @xmath327 for all perturbations that conserve mass , i.e. @xmath328 .
after some calculations @xcite , this can be put in the quadratic form @xmath329 we are led therefore to consider the eigenvalue problem @xmath330 if all the eigenvalues @xmath331 are negative , then @xmath332 and the configuration is a minimum of @xmath283 at fixed mass .
this implies that it is nonlinearly dynamically stable .
if at least one eigenvalue @xmath331 is positive , the configuration is a saddle point of @xmath283 and the star is dynamically unstable .
the marginal case is when the largest eigenvalue @xmath331 is equal to zero .
now , for @xmath310 , eqs .
( [ fdf2 ] ) and ( [ stab4 ] ) coincide .
this implies that the conditions of linear and nonlinear dynamical stability are the same . in the case of polytropic stars ,
the case of marginal stability corresponds to the critical index @xmath300 . at that index ,
the equilibrium configurations with mass @xmath13 and radius @xmath63 _ all _ have the same value of energy @xmath333 ( see eq .
( [ stab22 ] ) later ) .
this is therefore a very degenerate situation .
the nonlinear dynamical stability of gaseous polytropes can also be investigated by plotting the series of equilibria @xmath334 ( mass vs central density ) of box - confined configurations and using the turning point argument of poincar @xcite .
we can thus determine whether the last point on the series of equilibria , which corresponds to a _
complete polytrope _ whose density vanishes precisely at the box radius , is stable or not . for @xmath95
, there is no turning point of mass ( we restrict ourselves to @xmath335 ) , implying that the gaseous polytropes are always stable . for @xmath90 , a turning point of mass @xmath334
appears precisely for @xmath336 ( see fig . 5 of @xcite ) .
this method shows that , for @xmath90 , complete polytropes with @xmath337 are nonlinearly dynamically stable ( they are minima of @xmath283 at fixed mass @xmath102 ) while complete polytropes with @xmath338 are dynamically unstable ( they are saddle points of @xmath283 at fixed mass @xmath102 ) .
complete polytropes with @xmath336 and @xmath251 are marginally stable .
we can also investigate the linear dynamical stability of gaseous polytropic spheres by using the method introduced by eddington @xcite and ledoux @xcite .
if we introduce the radial displacement @xmath339 we can rewrite the eigenvalue equation ( [ stab5 ] ) in the form @xcite : @xmath340 this is the eddington equation of pulsations , which has been written here in the form of a sturm - liouville problem .
it must be supplemented by the boundary conditions @xmath341 @xmath342 where @xmath343 is the lagrangian derivative of the pressure . since @xmath344 at the surface of the star , it is sufficient to demand that @xmath345 and @xmath346 be finite in @xmath314 .
multiplying eq .
( [ stab10 ] ) by @xmath345 and integrating between @xmath347 and @xmath63 , we obtain @xmath348 the system is linearly dynamically stable if @xmath349 and unstable otherwise .
since @xmath350 , a sufficient condition of stability is @xmath351 , i.e. @xmath352 it can be shown furthermore ( see below ) that the system is unstable if @xmath353 so that the criterion ( [ stab12 ] ) is a necessary and sufficient condition of dynamical stability ( the case @xmath354 is marginal ) . in terms of the index @xmath186 , a complete polytrope is stable with respect to the euler - poisson system in @xmath2 for @xmath335 and for @xmath355 , in @xmath3 for @xmath356 , in @xmath4 for @xmath357 and in @xmath90 for @xmath358 .
from the theory of sturm - liouville problems , it is known that expression ( [ stab11 ] ) , which can be written @xmath359 , forms the basis of a variational principle . the function @xmath360 which maximizes the functional @xmath361 is the fundamental eigenfunction and the maximum value of this functional gives the fundamental eigenvalue @xmath362 .
furthermore , any trial function under - estimates the value of @xmath362 so this variational principle may prove the existence of instability but can only give approximate information concerning stability . as shown by ledoux & pekeris @xcite
, we can get a good approximation of the fundamental eigenvalue by taking @xmath360 to be a constant ( note that @xmath363 , i.e. @xmath364 , is the _ exact _ solution of the sturm - liouville equation ( [ stab10 ] ) at the point of marginal stability @xmath310 for a polytropic equation of state ) . for the trial function @xmath365 , expression ( [ stab11 ] )
gives @xmath366 according to the condition of hydrostatic balance , @xmath367 we have @xmath368 where we recognize the virial @xmath369 inserting @xmath370 in eq .
( [ stab16 ] ) , interchanging the dummy variables @xmath371 and @xmath372 and taking the half - sum of the resulting expressions , we get @xmath373 therefore , the virial can be written @xmath374 @xmath375 where @xmath376 is the potential energy and @xmath377 is the gravitational potential .
therefore , we can rewrite eq .
( [ stab13 ] ) in the form @xmath378 where @xmath379 is the moment of inertia . since @xmath380
, we conclude that the system is unstable if @xmath381 , which completes the proof above . on the other hand , eq .
( [ stab17 ] ) provides an estimate of the pulsation period @xmath307 when the system is stable .
this is the @xmath1-dimensional generalization of the ledoux stability criterion .
note finally that , for a spherically symmetric system , using the gauss theorem @xmath382 we get @xmath383 this expression may be useful to calculate the potential energy .      the virial theorem associated with the barotropic euler - poisson system ( [ stab1])-([stab3 ] )
is @xcite : @xmath384 where @xmath385 is the kinetic energy of the macroscopic motion and @xmath386 . for a polytropic equation of state , the energy functional ( [ wkre ] ) , which is a conserved quantity ,
can be written @xmath387 where @xmath388 is the internal energy .
we note that @xmath389 for a polytrope . at equilibrium , @xmath390
, we get @xmath391 and @xmath392 . for @xmath3
, this implies that @xmath393 . for @xmath394
, this implies that @xmath395 where we recall that @xmath396 for @xmath90 and @xmath335 , the system is dynamically stable if @xmath397 and unstable otherwise ( poincar argument @xcite ) . since @xmath398 , we find from eq .
( [ stab22 ] ) that the polytropic star is stable if and only if @xmath399 . at the point of marginal stability @xmath300 , where several steady configurations exist with the same mass @xmath13 given by eq .
( [ dim13b ] ) and an arbitrary radius @xmath63 , the energy of these configurations is @xmath333 for all @xmath63 .    using @xmath400 and eliminating the pressure term in eq .
( [ vir1 ] ) , the virial theorem for @xmath184 can be put in the form @xmath401 alternatively , eliminating the kinetic energy in eq .
( [ vir1 ] ) , the virial theorem for @xmath184 can be written @xmath402 for @xmath403 , corresponding to classical white dwarf stars , the virial theorem becomes @xmath404 for @xmath405 , corresponding to relativistic white dwarf stars , we have @xmath406 finally , considering the polytropic index @xmath300 of marginal stability , the virial theorem ( [ vir2 ] ) takes the form @xmath407 if we consider the dimension @xmath6 , it reduces to @xmath408 this equation describes the case of classical white dwarf stars at the critical dimension @xmath6 where @xmath409 .
it can be integrated into @xmath410 for @xmath411 , we find that @xmath412 for @xmath413 so that the system evaporates . alternatively , for @xmath414
, we find that @xmath415 in a finite time , so that the system collapses and forms a dirac peak in a finite time .
if we consider the dimension @xmath4 , eq . ( [ vir3 ] ) reduces to @xmath416 this equation describes the case of relativistic white dwarf stars at the critical dimension @xmath4 where @xmath417 . for @xmath411 ,
we find that @xmath412 for @xmath413 indicating that the system evaporates .
finally , we show that the instability criterion for polytropic stars can be obtained from simple dimensional analysis .
we shall approximate the internal energy ( [ vir1c ] ) and the potential energy ( [ stab19b ] ) by @xmath418 @xmath419 where @xmath102 is the total mass of the configuration and @xmath63 its radius ( @xmath183 and @xmath181 are dimensionless parameters ) .
for homogeneous spheres , the values of @xmath181 and @xmath183 are given by eq .
( [ dim1c ] ) . with these expressions , the energy functional ( [ qw1 ] ) becomes @xmath420 we now need to minimize this functional with respect to @xmath63 for a given mass @xmath102 ( a minimum necessarily requires @xmath421 ) .
we first look for the existence of critical points ( extrema ) .
the cancellation of the first order variations @xmath422 yields the mass - radius relation @xmath423 this relation determines the radius @xmath63 of the star as a function of its mass @xmath102 .
this expression is consistent with the exact mass - radius relation ( [ dim13 ] ) deduced from the lane - emden equation ( of course , our simple approach can only model compact density profiles corresponding to @xmath424 for @xmath90 ) @xcite . for @xmath296 ,
there is one , and only one , extremum of @xmath425 for each mass @xmath102 . in order to have a true minimum
, we need to impose @xmath426 therefore , the system is stable if @xmath427 and unstable otherwise .
this simple dimensional analysis returns the exact stability criterion ( [ stab12 ] ) .
we can be a little more precise . for @xmath428 ,
the functional @xmath425 has a global minimum reached for a finite , and non - zero , value of @xmath63 .
this solution is stable . for @xmath429 ,
the functional @xmath425 has an unstable global maximum : for @xmath103 , @xmath430 when @xmath175 ( evaporation ) and @xmath431 when @xmath252 ( collapse ) ; for @xmath3 , @xmath430 when @xmath252 and @xmath175 ; for @xmath90 and @xmath432 , @xmath430 when @xmath252 and @xmath433 ; for @xmath90 and @xmath434 , @xmath430 when @xmath252 and @xmath431 when @xmath433 . since we have constructed a particular configuration ( homogeneous sphere ) which makes the energy @xmath283
diverge to @xmath435 , the above arguments prove that the exact functional @xmath436 $ ] given by eq .
( [ wkre ] ) has no absolute minimum at fixed mass for @xmath429 .
since the functional @xmath324 $ ] has only one critical point ( cancelling the first variations ) at fixed mass , we conclude that , for @xmath429 , this critical point is a saddle point .
therefore , there is no stable steady state of polytropic spheres for @xmath429 . for the critical index @xmath300 , the relation ( [ dim3 ] )
shows the existence of a critical mass @xmath437 the functional ( [ dim1 ] ) can be rewritten @xmath438 for @xmath251 , @xmath439 for all @xmath63 ( at equilibrium @xmath421 ) .
this determines an infinite family of solutions with the same mass @xmath13 and different radii .
these solutions have the same energy and are marginally stable .
this returns the result of the exact model where the configurations are solution of the lane - emden equation ( see sec .
[ sec_eigen ] ) . for @xmath249 ,
the function @xmath425 is monotonically decreasing with @xmath63 ( so the system tends to evaporate ) : for @xmath103 , @xmath425 goes from @xmath347 to @xmath435 and for @xmath90 , @xmath425 goes from @xmath440 to @xmath347 . for @xmath14 ,
the function @xmath425 is monotonically increasing with @xmath63 ( so the system tends to collapse ) : for @xmath103 , @xmath425 goes from @xmath347 to @xmath440 and for @xmath90 , @xmath425 goes from @xmath435 to @xmath347 .
we can obtain a simple dynamical model by using the virial theorem ( [ vir1 ] ) . at equilibrium ( @xmath390 )
, we have @xmath441 inserting the expressions ( [ dim1a ] ) and ( [ dim1b ] ) in eq .
( [ dim6 ] ) , we recover the mass - radius relation ( [ dim3 ] ) . for @xmath296 and for any given mass @xmath102 ,
there is only one steady state .
its radius @xmath106 given by eq .
( [ dim3 ] ) and its energy @xmath442 is given by eq .
( [ stab22 ] ) .
it corresponds to the extremum value of @xmath425 . estimating the moment of inertia by @xmath443 and inserting the expressions ( [ dim1a ] ) and ( [ dim1b ] ) in eq .
( [ vir2b ] ) , we obtain @xmath444 this equation determines the evolution of the radius of the star for a fixed mass @xmath102 and a fixed energy @xmath283 .
the evolution of the kinetic energy is then given by @xmath445 .
the solution of eq .
( [ dim8 ] ) depends on two control parameters @xmath102 and @xmath283 and on the initial condition @xmath446 and @xmath447 .
we shall consider the case where @xmath283 is equal to the value @xmath442 corresponding to the steady state , such that the r.h.s . of the above equation is equal to zero at equilibrium .
then , we can rewrite eq .
( [ dim8 ] ) as @xmath448 where @xmath106 is the radius of the star at equilibrium . since the dynamics is non - dissipative , in the euler equation ( [ stab2 ] ) @xcite . in that case @xmath449 , so that the system evolves so as to minimize the functional @xmath321 $ ] at fixed mass . ]
the system does not evolve towards the minimum of @xmath425 ( unless @xmath450 initially ) .
the system can either oscillate around the minimum ( stable ) or evolve away from it ( unstable ) .
considering small perturbations around equilibrium , writing @xmath451 with @xmath452 , and linearizing the foregoing equation , we obtain @xmath453 where we have used the mass - radius relation ( [ dim3 ] ) to simplify the expression .
this is the equation for a harmonic oscillator with pulsation @xmath454 the star is stable if @xmath455 and unstable otherwise .
this returns the exact stability criterion ( [ stab12 ] ) .
furthermore , using eqs .
( [ dim1b ] ) and ( [ dim7 ] ) , the pulsation can be rewritten in the form @xmath456 which exactly coincides with the ledoux formula ( [ stab17 ] ) .
therefore , our simple dimensional model allows to obtain a lot of interesting results .
the case of arbitrary perturbations around equilibrium will be considered in a future work .
finally , we note that for @xmath300 , eq . ( [ dim8 ] ) becomes @xmath457 and eq .
( [ dim9 ] ) corresponds to @xmath333 and @xmath251 yielding @xmath458 .
hence , @xmath459 corresponding to a marginal evolution .    in our dimensional analysis , the constants @xmath181 , @xmath183 and @xmath326 are dimensionless parameters that could be chosen to fit the exact results at best . alternatively
, we can try to obtain quantitative predictions by calculating their values for homogeneous spheres .
this yields @xmath460 note that the potential energy of a homogeneous sphere in @xmath1 dimensions can be easily calculated from eq .
( [ stab19e ] ) . on the other hand , in @xmath3 ,
a direct calculation gives @xmath461 where @xmath462 is a reference radius where @xmath463 .
all the results given above pass to the limit @xmath185 provided that we take @xmath464 . using eq .
( [ dim1c ] ) , we find that the approximate value of the pulsation ( [ dim11 ] ) becomes @xmath465 on the other hand , if we compare the approximate mass - radius relation ( [ dim3 ] ) with the exact mass - radius relation ( [ dim13 ] ) , we obtain an estimate of the constant @xmath466 in the form @xmath467 this has to be compared with the exact value ( [ dim15 ] ) given in terms of the solution @xmath299 of the emden equation of index @xmath186 in @xmath1 dimensions @xcite .
let us consider the case @xmath4 . for @xmath10 ,
corresponding to classical white dwarf stars , we find @xmath468 instead of the exact value @xmath469 .
for @xmath12 , corresponding to relativistic white dwarf stars , we find @xmath470 instead of the exact value @xmath471 this suggests that the homogeneous star model will provide a fair description of relativisitic white dwarf stars and a poorer description of classical white dwarf stars .
we shall come back to these different issues ( static and dynamics ) in a future work .
in this appendix , we show that the condition of thermodynamical stability in the canonical ensemble is equivalent to the condition of nonlinear dynamical stability with respect to the barotropic euler - poisson system .
we apply this result to white dwarf stars .
we consider a barotropic gas with an equation of state @xmath278 described by the euler - poisson system @xcite .
we introduce the energy functional @xmath472 the first term @xmath473 is the work @xmath474 done in compressing the system from infinite dilution , the second term @xmath475 is the gravitational energy and the third term @xmath476 is the kinetic energy associated with the mean motion . for a gas in local thermodynamic equilibrium , the equation of state
is @xmath477 or @xmath478 and the first law of thermodynamics can be written @xmath479 where @xmath480 is the entropy density and @xmath481 is the density of internal energy . for a gas without interaction ( apart from the long - range gravitational attraction ) ,
the internal energy is equal to the kinetic energy .
there are two important cases where the gas is barotropic .
the first case is when @xmath482 .
this concerns either adiabatic ( or isentropic ) fluids so that @xmath483 or fluids at zero temperature so that @xmath0 . when @xmath482 , the first law of thermodynamics reduces to @xmath484 .
it can be integrated into @xmath485d\rho'$ ] .
then , the work @xmath473 done by the pressure force ( first term in eq .
( [ ef1 ] ) ) coincides with the internal energy @xmath486 of the gas . in that case , we get @xmath487 and the energy functional ( [ ef1 ] ) can be written @xmath488 thus , at @xmath0 or for an adiabatic evolution , the total energy of the gas @xmath202 is conserved by the euler - poisson system ( since @xmath283 is conserved ) . alternatively , for an isothermal gas @xmath489 , the first law of thermodynamics @xmath479 can be written @xmath490 .
it can be integrated into @xmath491d\rho'$ ] .
therefore , the work @xmath473 done by the pressure force ( first term in eq .
( [ ef1 ] ) ) coincides with the free energy @xmath492 of the gas . in that case , we get @xmath493 and the energy functional ( [ ef1 ] ) can be written @xmath494 thus , for an isothermal evolution , the free energy of the system @xmath495 is conserved by the euler - poisson system ( since @xmath283 is conserved ) . at @xmath0
, we recover the conservation of the energy @xmath202 .
let us apply these results to white dwarf stars .
we can view a white dwarf star as a barotropic gas described by an equation of state @xmath278 . according to the discussion of appendix [ sec_eigen ]
( see also @xcite ) , it is nonlinearly dynamically stable with respect to the euler - poisson system if it is a minimum of the energy functional @xmath283 at fixed mass . if a minimum exists , it is necessary that @xmath496 . as a result
, a steady state of the euler - poisson system is nonlinearly dynamically stable if , and only if , it is a minimum of @xmath497={\cal w}[\rho,{\bf u}]-\theta[{\bf u}]$ ] at fixed mass @xmath102 . in conclusion , the condition of nonlinear dynamical stability can be written @xmath498\quad |\quad m[\rho]=m\rbrace . \label{ef4}\ ] ] for a white dwarf star at zero temperature ( @xmath0 ) , according to eqs .
( [ ef2 ] ) and ( [ ef3 ] ) , the functional @xmath499 reduces to the energy @xmath500 , where @xmath486 is the kinetic energy .
we can also view a white dwarf star as a gas of self - gravitating relativistic fermions at statistical equilibrium in the canonical ensemble .
it is thermodynamically stable if , and only if , it is a minimum of free energy @xmath495 at fixed mass @xmath102 .
the free energy is given by @xmath501 where @xmath502 is the maximum value of the distribution function fixed by the pauli exclusion principle .
for a white dwarf star at zero temperature ( @xmath0 ) , the free energy @xmath495 reduces to the energy @xmath500 .
the condition of thermodynamical stability in the canonical ensemble can be written @xmath503\quad |\quad m[f]=m\rbrace . \label{ef6}\ ] ] the critical points of free energy at fixed mass , determined by the variational principle @xmath504 , correspond to the relativistic mean - field fermi - dirac distribution @xmath505 where @xmath506 is a lagrange multiplier ( related to @xmath326 ) determined by the mass @xmath102 and @xmath28 is given by eq .
( [ es3 ] ) . using eqs .
( [ es1 ] ) , ( [ es4 ] ) and ( [ ef7 ] ) , the density and the pressure are of the form @xmath507 and @xmath508 . eliminating @xmath509 between these two expressions , we find that @xmath278 so that the gas is barotropic .
the equation of state is parameterized by @xmath229 and is fully determined by the entropic functional in eq .
( [ ef5 ] ) , which is here the fermi - dirac entropy . at @xmath0
we obtain the explicit relations of sec .
[ sec_es ] .
furthermore , the condition that @xmath510 is a function of the energy @xmath511 implies that the corresponding barotropic gas is at hydrostatic equilibrium . indeed , using eq .
( [ es5 ] ) , @xmath512 so that @xmath513    the relativistic mean - field fermi - dirac distribution ( [ ef7 ] ) is just a critical point of free energy at fixed mass @xmath102 . to determine whether it corresponds to a true minimum of free energy
, we can proceed in two steps .
we first minimize @xmath514 $ ] for a fixed density profile @xmath515 .
since the potential energy @xmath516 $ ] and the mass @xmath322 $ ] are entirely determined by the density profile , this is equivalent to minimizing @xmath517-ts[f]$ ] at fixed @xmath515 , where @xmath486 is the kinetic energy .
this gives a distribution @xmath518 where @xmath519 is a local lagrange multiplier determined by the density @xmath515 , using @xmath520 .
since @xmath521 , the distribution ( [ ef8 ] ) is a true minimum of @xmath514 $ ] at fixed @xmath515 . substituting the optimal distribution function ( [ ef8 ] ) in eq .
( [ ef5 ] ) , we obtain a functional @xmath522\equiv f[\tilde{f}]$ ] of the density @xmath515 . using eqs .
( [ es1 ] ) , ( [ es4 ] ) and ( [ ef8 ] ) , we note that the density and the pressure are of the form @xmath523 and @xmath524 . eliminating @xmath519 between these expressions
, we find that @xmath278 where the equation of state is the same as the one determined from the fermi - dirac distribution ( [ ef7 ] ) at equilibrium . now
, we can show that @xmath525=\tilde{\cal w}[\rho ] , \label{ef9}\end{aligned}\ ] ] where @xmath526 $ ] is the functional defined in sec .
[ sec_efa ] . the relation ( [ ef9 ] ) can be obtained by an explicit calculation which extends the proof given in appendix b of @xcite for classical particles ( this will be shown in a future work ; see also particular cases in sec . [ sec_efc ] ) but it also results from a straightforward argument .
we note that the distribution function ( [ ef8 ] ) corresponds to a condition of local thermodynamical equilibrium with uniform temperature and zero average velocity .
thus , it locally satisfies the first law of thermodynamic @xmath479 with @xmath489 . integrating this relation like in sec .
[ sec_efa ] , we find that @xmath527-ts[\rho]={\cal w}_{1}[\rho]$ ] where @xmath527=u[\tilde{f}]$ ] and @xmath528=s[\tilde{f}]$ ] .
this directly yields the identity ( [ ef9 ] ) .
we can now conclude that @xmath514 $ ] has a minimum @xmath529 at fixed mass @xmath102 if , and only if , @xmath522=\tilde { \cal w}[\rho]$ ] has a minimum @xmath530 at fixed mass @xmath102 .
in that case , @xmath529 is given by eq .
( [ ef8 ] ) where @xmath531 is determined by @xmath530 , writing @xmath532 .
therefore , a system is thermodynamically stable in the canonical ensemble if , and only if , the corresponding barotropic gas with the same equilibrium distribution is nonlinearly dynamically stable with respect to the barotropic euler - poisson system . said differently , a system that minimizes the functional ( [ ef9 ] ) is ( i ) thermodynamically stable in the canonical ensemble and ( ii ) nonlinearly dynamically stable with respect to the barotropic euler - poisson system .
this result applies to white dwarf stars at arbitrary temperature .
although the above results are general , it may be useful to explicitly compute the functionals ( [ ef1 ] ) and ( [ ef5 ] ) for white dwarf stars at @xmath0 and check the relation ( [ ef9 ] ) .
if we view a white dwarf star as a barotropic gas described by an equation of state @xmath533 , its energy @xmath499 can be written @xmath534 with @xmath535 a white dwarf star is nonlinearly dynamically stable with respect to the euler - poisson system if it is a minimum of @xmath499 at fixed mass @xmath102 . at @xmath0 , the equation of state
is given by eq .
( [ es7 ] ) . in the classical limit ( polytrope @xmath87 ) we have @xmath536 and in the ultra - relativisitic limit ( polytrope @xmath127 ) we have @xmath537 in the general case ,
( [ es7 ] ) , the function ( [ ef11 ] ) can be written @xmath538 integrating by parts and using eqs .
( [ es8 ] ) and ( [ es9 ] ) , we find after straightforward calculations that @xmath539    alternatively , we can view a white dwarf star as a gas of self - gravitating relativistic fermions at statistical equilibrium in the canonical ensemble . at @xmath0 ,
its free energy @xmath230 coincides with its energy @xmath202 . using eq .
( [ ef5 ] ) , it is given by @xmath540 with @xmath541 a white dwarf star at @xmath0 is thermodynamically stable if it is a minimum of the energy @xmath202 at fixed mass @xmath102 . in the classical limit , @xmath542 and @xmath543 so @xmath544
therefore , the free energy can be written @xmath545 using eq .
( [ class1 ] ) , this is equivalent to eq .
( [ ef12 ] ) . in the ultra - relativistic limit , @xmath546 and @xmath547 so @xmath548 .
therefore , the free energy can be written @xmath549 using eq .
( [ ur1 ] ) , this is equivalent to eq .
( [ ef13 ] ) . in the general case
, we have @xmath550 using eq .
( [ es3 ] ) , it can be put in the form @xmath551 now , it is straightforward to check that the two expressions in the r.h.s .
( [ ef15 ] ) and ( [ ef21 ] ) are equal so that @xmath552 implying the relation ( [ ef9 ] ) .    finally , if we consider a classical isothermal self - gravitating gas at temperature @xmath229 with an equation of state @xmath286 , its free energy ( [ ef5 ] ) can be written @xmath553 comparing with eq .
( [ qw3 ] ) , we find that relation ( [ ef9 ] ) is indeed satisfied .
therefore , a classical isothermal self - gravitating gas that minimizes the functional ( [ qw3 ] ) or ( [ ef22 ] ) at fixed mass is ( i ) thermodynamically stable in the canonical ensemble and ( ii ) nonlinearly dynamically stable with respect to the barotropic euler - poisson system . we had already made this observation in @xcite . in @xcite ,
this result was extended to an arbitrary form of entropic functional , including the fermi - dirac entropy , for non - relativistic systems .
the present paper shows that , due to relation ( [ ef9 ] ) , the equivalence between nonlinear dynamical stability with respect to the barotropic euler - poisson system and thermodynamical stability in the canonical ensemble is general .
if we consider , in a @xmath1-dimensional universe , a hamiltonian system of self - gravitating classical point masses whose dynamics is described by the newton equations @xmath554 the scalar virial theorem reads @xcite : @xmath555 where @xmath176 is the kinetic energy and @xmath556 the trace of the potential energy tensor for the @xmath178-body system @xcite . for @xmath184 , using eq .
( [ stab18 ] ) and introducing the total energy @xmath198 , which is a constant of the motion for an isolated system , we have @xcite : @xmath557 we note that the dimension @xmath6 is critical . in that case , @xmath558 which yields after integration @xmath559 . for @xmath560 , @xmath561 indicating that the system evaporates . for @xmath562
, @xmath563 goes to zero in a finite time , indicating that the system forms a dirac peak ( `` black hole '' ) in a finite time . more generally , for @xmath8 , since @xmath564 , we have @xmath565 so that the system forms a dirac peak in a finite time if @xmath562
. therefore , self - gravitating systems with @xmath562 are not stable in a space of dimension @xmath8 .
the study of the present paper indicates that this observation remains true if quantum ( pauli exclusion principle for fermions ) and relativistic effects are taken into account . in this sense ,
a universe with @xmath8 is not viable .