the study of charged relativistic fluid balls attract the interest of researchers of different areas of physics and astrophysics .
there exists a general consensus that astrophysical objects with large amounts of charge can not exist in nature @xcite , @xcite .
this point of view had been challenged by several researchers @xcite , @xcite , @xcite , @xcite .
it can not be discarded the possibility that during the gravitational collapse or during an accretion process onto a compact object the matter acquires large amounts of electric charge .
this has been considered in @xcite and @xcite .    in this paper
we study the stellar structure and temporal evolution of compact charged fluid spheres , independently on the mechanism by which the matter acquires an electric charge .    from a pure theoretical point of view , the collapse of charged fluid balls and shells are connected with the laws of black hole thermodynamics .
the third law of black hole thermodynamics states that no process can reduce the surface gravity of a black hole to zero in a finite advanced time @xcite , @xcite , @xcite . the surface gravity , @xmath0 , plays the role of a temperature while the area of the event horizon is equivalent to the entropy .
israel @xcite gave a formulation and proof of the third law .
it was demonstrated that the laws of black hole thermodynamics are analogous to the common laws of thermodynamics @xcite .
translated to the physics of fluid collapse , the third law implies the impossibility of forming an extremal black hole , for which @xmath1 .
an extremal black hole has a total charge @xmath2 : @xmath3 , where @xmath4 is the gravitational constant , and @xmath5 is the total mass of the black hole .
the extremal black hole has not an horizon and constitutes a naked singularity , contradicting the cosmic censorship hypothesis .
so , in this paper we explore the possibility of forming an extremal black hole from the collapse of a ( charged ) compact object . we will show that the answer is no .
the relativistic equations for the collapse of a charged fluid ball were obtained by bekenstein @xcite .
so far in the literature , as long as the author knows , it was studied the dynamics of charged shells or scalar fields falling onto an already formed charged black hole ( see @xcite , @xcite , @xcite , and references therein ) , and solutions of the tolman - oppenheimer - volkoff @xcite equations for charged stars ( see @xcite , @xcite , @xcite , @xcite , etc ) . however , in the stellar collapse , the collapsing matter determines the background geometry and the metric of space - time is changing with time as the star collapses .
moreover , in general relativity , the pressure of the fluid contributes to the gravitational field and strengthen it . from the present study it is possible to understand if the coulomb repulsion will prevent or not the total collapse of a charged fluid ball , and which are the stability limits for a charged relativistic star .
we emphasize that we will not be concerned here with the mechanism of electric charge generation , nor with effects that could neutralize or discharge the star .    in the present paper ,
the equations for the evolution of charged fluid spheres in schwarzchild - like coordinates are presented in sec .
[ relatequ ] . in the sec .
[ sec2 ] , we re - derived the equations of hydrostatic equilibrium and the relativistic equations for the temporal evolution of charged spheres in a form closer to that obtained by may and white @xcite , @xcite ( for the case of zero charge ) .
the equations obtained are well suited for performing a finite difference scheme , and allows a direct implementation of very well known numerical techniques @xcite , @xcite , @xcite .
the formalism is also compatible with the bekenstein equations @xcite , and with the misner and sharp equations @xcite for the case of zero charge .
the calculations performed in this study are rather lengthy in general and we tried to be the more self - consistent as possible .
however , we apologize that in most of the paper we only give hints for the computation of the intermediate algebraic steps in order to keep the paper at a reasonable length .
the matching conditions between the exterior and interior solutions of the einstein - maxwell equations are given in the sec .
[ match ] .    for completeness , in sec .
[ seceos ] it is considered an equation of state of a zero temperature neutron gas .
two codes were built and are used in the present study : one to obtain neutron stars models in hydrostatic equilibrium , presented in the sec .
[ sectov ] ; and the other code is to evolve the stars forward in time , integrating the einstein - maxwell equations .
this code is introduced in sec .
[ sechydro ] . in all the simulations ,
the charge distribution is taken proportional to the rest mass distribution .
we found that although the total binding energy grows with the total amount of charge , the binding energy per nucleon tends to zero as the charge tends to the extremal value .
this is related to the impossibility of obtaining bounded solutions of extremely charged fluid configurations .
in addition , it results that the upper limit for the total amount of charge could be slightly lower than the extremal value .
this is discussed in the sec .
[ secresults ] .    in the sec .
[ secresults ] , is discussed the stability and evolution of the stellar models .
we study , as well , the evolution of stars that were suddenly discharged to take into account the case of charged spheres being a metastable state of more complex scenarios .    in a recent work ghezzi and
letelier @xcite presented a numerical study of the collapse of charged fluid spheres with a polytropic equation of state , and with an initial uniform distribution of energy density and charge .
they found that during the evolution a shell of higher density is formed near the surface of the imploding star . from that work
we can not say if a shell will always form when changing the initial conditions and the scenario for the collapse . in the present study
, we checked that the effect is negligible during the collapse of a charged neutron star ( see sec .
[ secresults ] ) .    in the sec .
[ secremarks ] , we end with some final remarks .
considering a spherical symmetric fluid ball , the line element in schwarzchild - like , or standard form , is ( see @xcite , @xcite , @xcite ) : @xmath6 an observer moving radially have a 4-velocity , and the electric 4-current is .
the energy - momentum tensor of the charged fluid is given by : @xmath7\end{aligned}\ ] ] where @xmath8 is the density of mass - energy , @xmath9 is the scalar pressure , and @xmath10 is the electromagnetic tensor .
the electromagnetic field satisfies the maxwell equations : @xmath11 and @xmath12}=0\,.\ ] ] only the radial component @xmath13 is non - zero , and the last equation is satisfied if @xmath14 . the covariant derivative in eq .
( @xmath15 ) can be written as @xmath16_{,\nu}=4 \pi j^{\mu}\,,\ ] ] with @xmath17 .
defining @xcite @xmath18 the eq .
( [ maxwell2 ] ) gives @xcite @xmath19 and @xmath20    integration of eq .
( [ mw1 ] ) gives @xcite @xmath21 where @xmath22 and from eq .
( [ mw2 ] ) @xmath23 from this equation we see that the charge is a constant @xmath24 outside the fluid ball , where @xmath25 . far away from the sphere ,
where @xmath26 ( see below ) , the eq .
( [ campoe ] ) gives the electric field of a classical charged particle .
the einstein s equations to be solved are : @xmath27 here ( and only here ) , @xmath28 denotes the scalar curvature , and @xmath29 is the ricci tensor .    the components of the einstein equations
are ( see @xcite , @xcite ) : @xmath30=   \nonumber\\ & & \hspace{0.8 cm } e^{-\lambda } \,\biggl(\frac{\lambda'}{r}-\frac{1}{r^2}\biggr)+\frac{1}{r^2}\,,\end{aligned}\ ] ]    @xmath31= \nonumber\\ & & \hspace{0.8 cm } -\,e^{-\lambda}\ , \biggl(\frac{\nu'}{r}+\frac{1}{r^2}\biggr)+\frac{1}{r^2}\,,\end{aligned}\ ] ]    @xmath32    multiplying the right hand side of eq .
( [ e1 ] ) by @xmath33 and rearranging , the equation can be cast in the form @xmath34\ , r^2\,.\ ] ] setting @xmath35 , with @xmath36 and @xmath37 : @xmath38 the eq .
( [ re1 ] ) becomes @xmath39\ , r^2\,.\ ] ] integrating this equation , we get the equation for the mass @xcite : @xmath40\ , r^2\ , dr \nonumber\\ & & \hspace{0.8cm}+\frac{1}{2 c^2 } \frac{q^2}{r}+ \frac{1}{2 c^2 } \int^r_0 \frac{q^2}{r^2 } dr+m_0\ , , \end{aligned}\ ] ] where @xmath41 is an integration constant , and can be taken as zero for the purpose of this work .    outside the fluid ball
, the mass do not depend on @xmath42 and then : @xmath43 for @xmath44 , where @xmath45 is the coordinate @xmath42 of the surface of the sphere .
actually the mass do not depend on @xmath46 for @xmath44 ; from eqs .
( [ e3 ] ) and ( [ gamma0 ] ) : @xmath47 as @xmath2 is constant and @xmath48 outside the fluid ball , @xmath49 is independent of @xmath46 .    subtracting eq .
( [ e1 ] ) from eq .
( [ e2 ] ) it is @xmath50 so , @xmath51 is independent of @xmath42 . to get an asymptotic flat solution
it must be @xmath52 as @xmath53 .
thus a solution of eq .
( [ metrica0b ] ) is @xcite @xmath54    using the eq .
( [ gamma0 ] ) , ( [ coefmetrica1 ] ) , and substituting into eq .
( [ metrica0 ] ) , the line element outside the ball is @xmath55 which is the reissner - nordstrm space - time .
the gravitational mass is @xmath56 .
we see that although the ( interior ) mass distribution and the metric depends on time , the external space - time is static .
this is the birkhoff theorem , and implies that a spherical distribution of mass and charge can not emit gravitational waves .
in this section we will specialize the equations for the case of a coordinate system co - moving with the fluid .
we will use a notation closer to that used in the may & white papers @xcite .
the 4-velocity for an observer co - moving with the fluid is @xmath57 , and the measured 4-current is @xmath58 .
the 4-velocity satisfies : @xmath59 .
the line element for co - moving coordinates in  standard " form ( see @xcite , pag .
336 ) is : @xmath60 note that now @xmath61 , @xmath62 , and @xmath28 are functions of @xmath63 and @xmath46 .
we will see in the next subsection that the coordinate @xmath63 can be chosen to be the total rest mass inside a sphere , such that each spherical layer of matter can be labeled by the rest mass it contains .
the stress tensor in co - moving coordinates is obtained using the solution of the maxwell equations ( [ campoe ] ) : @xmath64 from now on we will restore in the equations the gravitational constant @xmath65 and the speed of light @xmath66 .    similarly as was defined above , @xmath67 is the density of energy in @xmath68 $ ] and @xmath9 is the scalar pressure in the same units
. it will be useful to perform the split @xmath69 where @xmath70 is the rest mass density and @xmath71 is the internal energy of the gas .
we want to set the radial coordinate equal to the rest mass enclosed by co - moving spherical layers .
we performed a calculation analogous to that of may and white @xcite , and used the same notation for clarity and comparison .    in this work
is assumed that there is only one species of particles , and there are no particles created nor destroyed .
the number density of particles will be denoted by @xmath72 and the rest mass density by @xmath70 .
if the rest mass of the particles is @xmath73 then : @xmath74 the conservation of the baryon number current is expressed : @xmath75 or equivalently : @xmath76    the line element in spherical coordinates written in
standard " form is : @xmath77 observe that this line element is not identical with that given in eq .
( [ comline ] ) .
the conservation of baryons in this frame is reduced to ( see eq .
[ baryons ] ) : @xmath78 thus the quantity @xmath79 is a function of @xmath28 only @xmath80 choosing @xmath81 : @xmath82 on the other hand , the proper mass is @xmath83 where @xmath84 is the volume element of the space section , and @xmath85 is the volume of integration .
then @xmath86 is the proper mass enclosed by a sphere of circumference @xmath87 , and radial coordinate @xmath88 . with this election of coordinates : @xmath89 so we can perform the coordinate transformation @xmath90 and in this case the metric functions change to    @xmath91    replacing these functions on the line element ( [ comline2 ] ) , dropping the primes , and using the eq .
( [ diffs ] ) we obtain the line element given in the eq .
( [ comline ] ) .
so , the radial lagrangian coordinate was gauged to be the rest mass `` @xmath63 '' of each spherical layer of matter .
the proper time derivative of the charge is ( see eq .
[ dqdt ] ) @xmath92 because it is possible to write the 4-current as a product of a scalar charge density times the 4-velocity , i.e. : @xmath93 .
thus the electric charge is conserved in spherical layers co - moving with the fluid .      the einstein equations
( see the appendix ) for the components @xmath94 are @xmath95 where we use the notation @xmath96 , @xmath97 , etc .
the equation for the component @xmath98 is , must be a plus .
this is immediately evident when obtaining eq .
( 42 ) from eq .
( 35 ) , in that work . ] : @xmath99_{,\mu}\ ,   \nonumber\\ & & \hspace{0.45cm}-\frac{g q q_{,\mu}}{c^2 r}=g m_{,\mu}-\frac{g q q_{,\mu}}{c^2 r}\ , ,   \end{aligned}\ ] ] and the equation for the component @xmath100 is ) to eliminate @xmath101 from the component @xmath100 of the einstein s equations , given in the appendix .
in addition , use that @xmath102 . ] : @xmath103_t\,\nonumber\\ & & \hspace{2.02cm}=-g m_{,t}\ , .
\end{aligned}\ ] ]    we introduced into eqs .
( [ e00 ] ) and ( [ e11 ] ) , the definition of the total mass : @xmath104    we give further definitions that simplify the aspect of the equations : @xmath105 with this definitions an integral of the eq .
( [ e00 ] ) and ( [ e11 ] ) is : @xmath106 where @xmath49 is given by the eq .
( [ totalmass ] ) .
combining eqs ( [ metricgrr ] ) and ( [ r01 ] ) we get the mass conservation equation @xmath107    it can be seen that the quantity @xmath108 ( see eq . [ udefin ] ) is the radial component of the 4-velocity of the fluid in schwarzschild coordinates . using the transformation rule for the contra - variant 4-vector @xmath109 from co - moving to schwarzschild coordinates : @xmath110 , yields @xcite :
@xmath111 @xmath112      from the einstein equations and the contracted bianchi identities follows that : @xmath113 in co - moving coordinates the components are : @xmath114    expanding the eq .
( [ dmetb ] ) and using the eq .
( [ em0 ] ) we get : @xmath115 which is identical to the non - relativistic adiabatic energy conservation equation , and express the first law of thermodynamics
. it could be surprising for the reader that there is not an electromagnetic term on this equation : this is due to the spherical symmetry of the problem .
( [ em1 ] ) can be written : @xmath116_{,\mu}}{a ( \delta+p / c^2)/\rho } = \frac{1}{\bigl[(\delta+p / c^2)/\rho\bigr ] c^2 } \biggl[\epsilon_{,\mu } + \nonumber\\ & & \hspace{3.3cm}p \biggl(\frac{1}{\rho } \biggr)_{,\mu } + \frac{q q_{,\mu}}{4 \pi r^4 \rho } \biggr]\ , , \end{aligned}\ ] ] or using the eq .
( [ split ] ) and the definition of the relativistic specific enthalpy ( see @xcite , @xcite ) @xmath117 : @xmath118 we can write the eq .
( [ gtt1 ] ) in a form closer to that used by may & white @xcite : @xmath119\,.\ ] ]      the equation of motion can be obtained expanding the eq .
( [ e11 ] ) , using the eq .
( [ r01 ] ) to eliminate the factor @xmath120 , the eq .
( [ em1 ] ) to eliminate @xmath101 and the definitions given in the eqs .
( [ udefin ] ) , ( [ gammadef ] ) , ( [ gamma2 ] ) , ( [ enthalpy ] ) . after lengthy calculations we obtain : @xmath121\,.\end{aligned}\ ] ] this equation reduce to the equation of motion obtained by misner and sharp @xcite and may and white @xcite for the case @xmath122 , and is compatible with the equation derived by bekenstein @xcite .
( [ motion ] ) reduce to the newtonian equation of motion of a charged fluid sphere letting @xmath123 , and consequently : @xmath124 ( see eq .
[ gamma2 ] ) ; @xmath125 ( see eq .
[ enthalpy ] ) ; and @xmath126 ( see eq .
[ em1 ] ) .
the equation for charged fluid spheres in hydrostatic equilibrium is a generalization of the tolman - oppenheimer - volkov @xcite equation obtained by bekenstein @xcite .
we can obtain it from the equation of motion ( [ motion ] ) , taking .
in hydrostatic equilibrium , the factor @xmath127 ( see eq .
[ gammadef ] ) is @xmath128    rearranging ) can be simplified using the chain rule : @xmath129 , and using the eq .
( [ gammadef ] ) . ] the eq .
( [ motion ] ) and using the eqs .
( [ enthalpy ] ) and ( [ gammahe ] ) , the tolman - oppenheimer - volkov ( tov ) equation for charged fluid spheres is : @xmath130 this equation gives the well known tolman - oppenheimer - volkov @xcite equation when @xmath122 .
in the section [ relatequ ] we find an exterior solution to the problem which agrees with the reissner - nordstrm solution . in the preceding section we derived an interior solution that we want to integrate with a computer .
however , we must show that the interior and the exterior solutions match smoothly .
the matching conditions between the interior and the exterior solutions are obtained by the continuity of the metric and its derivatives along an arbitrary hypersurface .
first of all we must establish the equality between the exterior and the interior metric at the surface of the star .
this can be done transforming the interior solution given in co - moving coordinates into the schwarzschild frame .
the metric tensors in the two frames are related by the tensor transformation law : @xmath131 the line element in the schwarzschild frame is given by @xmath132 with @xmath133 . using the eq .
( [ tensorlaw ] ) we obtain the non - trivial relations between the metric coefficients :    [ metriceqs ] @xmath134    from eq .
( [ tensor2 ] ) and using eqs .
( [ gammadef ] ) and ( [ udefin ] ) we get : @xmath135 or equivalently , @xmath136 from eqs .
( [ tensor1 ] ) and ( [ tensor3 ] ) : @xmath137 with a little amount of algebra we can transform this equation into : @xmath138 , and using the eq .
( [ b1 ] ) : @xmath139 this last equation can be obtained independently by comparing the induced metric , from interior and exterior solutions , on the hypersurface @xmath140 : @xmath141
. the line element corresponding to the exterior solution is given by : @xmath142 from now on a plus ( minus ) subscript or superscript denotes exterior ( interior ) solutions .
the line element for the interior solution is ( see eq .
[ comline ] ) : @xmath143 due to the spherical symmetry : @xmath144 .
the line element compatible with the induced metric on the hypersurface , from the exterior solution is : @xmath145 and the metric induced on @xmath140 by the interior solution gives : @xmath146 on @xmath140 we have @xmath147 , so we obtain :    @xmath148    and we recovered the eq .
( [ a1 ] ) .
so eq . ( [ a2 ] ) ( or equivalently eq . [ a1 ] )
are conditions for the continuity of the metric along @xmath140 . from now on , we will use indistinctly @xmath149 .
the relation between the exterior and interior metric coefficients constitutes three equations with the four unknowns ( see eq .
[ metriceqs ] ) : @xmath150 .
observe that @xmath151 and @xmath152 are obtained as part of the interior solution .
in addition , the solution exterior to the sphere of matter must satisfy the field equations : @xmath153 , here @xmath154 is the ricci tensor .
so we obtain ( see for example @xcite pag .
180 and pag .
337 ) : @xmath155 the continuity of the metric is almost already established . only for completeness , and after a little algebra
, we can find an equation for @xmath156 : @xmath157 ( although we will not use it ) .
the equation of the hypersurface containing the surface of the star is @xmath158 where @xmath159 is the lagrangian mass at the surface of the star , so the metric coefficient @xmath160 is : @xmath161 with @xmath162 , @xmath163 , and @xmath164 .
the coefficient @xmath165 is @xmath166 , and we recover the reissner - nordstrm solution : @xmath167 it is important to remark that we obtained the reissner - nordstrm solution by performing a tensor transformation of the interior solution . the exterior solution must satisfy the einstein field equations , and so we used them to find the eq .
( [ a3 ] ) .
thus , the exterior and interior metric are consistently matched .
it remains to prove the continuity of the derivatives of the metric along @xmath168 .
this is a more lengthy calculation and we will give only the the main algebraic steps and the final results .
the continuity of the derivatives of the metric is established from the continuity of the second fundamental form , or extrinsic curvature tensor @xmath169 of the hypersurface @xmath168 . by definition
( see @xcite , pag .
59 ) : @xmath170 here @xmath171 is a unit vector normal to @xmath168 ; @xmath172 , are basis vectors on @xmath168 ; the coordinates on @xmath168 are denoted with @xmath173 , while @xmath174 are the space - time coordinates .
equivalently : @xmath175 , here @xmath176 is a lie derivative .
we will choose @xmath177 as coordinates on @xmath168 .
the equation of the hypersurface approaching from the exterior is : @xmath178 the basis vectors on @xmath179 are :    @xmath180    in the coordinate basis @xmath181 , and the normal co - vector is : @xmath182 the equation for @xmath168 as seen from the interior is ( see eq .
[ sigmas ] ) : @xmath183 the basis vectors on @xmath184 are :    @xmath185    in the coordinate basis @xmath186 , and the normal co - vector is : @xmath187 if there are no surface distributions of energy - matter on @xmath168 , it must be : @xmath188 \equiv k^+_{ab}-k^-_{ab } = 0\,.\ ] ] the angular components of this equation are @xmath189 , and @xmath190 . for @xmath191
we have : @xmath192 here @xmath193 is the connection . from now on , in the present section , we are not using the sum rule over repeated indexes . in the last two steps we used the eqs .
( [ a2 ] ) and ( [ a3 ] ) .    for @xmath194 , we find @xmath195 , and using the eqs .
( [ gammadef ] ) and ( [ gamma2 ] ) : @xmath196 , so eq . ( [ keq ] )
gives : @xmath197 or equivalently : @xmath198 here @xmath199 and @xmath200 , with @xmath201 a positive arbitrary small number or zero .
the definition of @xmath108 assures its continuity across @xmath168 , so @xmath202 .
moreover , we already found that @xmath5 and @xmath2 are constants independent of @xmath46 and @xmath63 , so the eq .
( [ gammaspm ] ) implies that @xmath203 .
then , using the eq .
( [ dmdt ] ) in co - movel coordinates ( @xmath204 in both coordinate systems , see eq .
[ dqdt ] and eq .
[ chargeconsintime ] ) , we obtain the boundary condition at the surface : @xmath205 the equation @xmath206=0 $ ] gives an analogous result .
the component @xmath207 is : @xmath208    the external component is : @xmath209 the sum rule is not used here . for clarity
, we will define @xmath210 , so eq . ( [ a2 ] ) becomes @xmath211 and an over - dot means @xmath212 . applying the over - dot operator to this equation
we obtain an equation that will be used later : @xmath213    so , @xmath214 expanding this equation , and using the eq .
( [ ddott ] ) to simplify the result we obtain : @xmath215 further simplifications can be made using : @xmath216 thus we obtain : @xmath217 comparing this equation with the eq .
( [ kminus ] ) we obtain an expresion for @xmath218=0 $ ] : @xmath219 now rearranging and remembering that @xmath220 and @xmath221 , using the eq .
( [ em1 ] ) to replace @xmath222 we get : @xmath223\,.\end{aligned}\ ] ] then , we see that the equation @xmath218=0 $ ] is satisfied by virtue of the equation of motion ( see eq .
[ motion ] ) , and thus the exterior solution matches with the interior solution .
the equations obtained above must be supplemented with an equation of state for the gas ( eos ) .
in particular we choose to represent the gas of neutrons as obeying the quantum statistic of a fermion gas at zero temperature .
this let us use the oppenheimer - volkov numerical results as a test - bed for the code .
however , any other eos can be equally implemented in the two codes presented in this paper .
the energy per unit mass is given by @xcite @xmath224 where @xmath225 is the relativity parameter , and @xmath226 is the momentum of the particles .
the fermi parameter is : @xmath227 , with @xmath228 , where @xmath72 is the number density of neutrons .
it can be written ( we suppress the sub - index f from now on ) : @xmath229 where @xmath230\,.\ ] ] the eq .
( [ eose ] ) can be integrated to give : @xmath231\ , , \end{aligned}\ ] ] the quantity @xmath232 is measured in [ @xmath233 , or equivalently in [ @xmath234 .
this is an expression for the total energy of the gas , so it includes the rest mass energy density .
although this is obvious by definition ( see eq .
[ eose ] ) , it can be useful to check this expanding eq .
( [ eose2 ] ) for @xmath235 ( or @xmath236 ) . in this two cases
an split similar to eq .
( [ split ] ) is obtained .
the pressure is given by @xcite : @xmath237 and performing the integral : @xmath238\ , , \end{aligned}\ ] ] measured in [ @xmath234 .
for charged matter , the equation of state must describe several species of particles , i.e. , neutrons , protons and electrons .
however , the amount of charge in all the models studied here is very small to change the eos . for an extremal charged fluid ball
there are only one unpaired ( not screened ) proton over @xmath239 particles . in this case
, the change on the chemical potential of the neutrons is negligible and charge neutrality can be assumed from a microscopic point of view . ]
relaxing the condition of charge neutrality , it is straightforward to show that for the case of extremely charged matter , the chemical potential of the neutrons is @xmath240 thus , the eos is slightly modified , and charge neutrality can be assumed from a microscopic point of view .
we will come back to this issue in a forthcoming paper . ]
the simulations were extended beyond the densities at which the eos is valid .
however , we want to take into account possible charge and pressure regeneration effects , that occurs at high densities .
the equations that must be integrated to obtain the stars in hydrostatic equilibrium are conveniently written in dimensionless form . from eqs .
( [ metrica0 ] ) and ( [ re1 ] ) it is possible to obtain an equation for the metric component , @xmath241    in addition we use the chain rule on the tov equation ( eq . [ toveq ] ) : @xmath242    the dimensionless equations ( [ toveq2 ] ) and ( [ eqlambda ] ) can be obtained with the replacements : @xmath243 where the bar denotes dimensionless variables , and the subscript @xmath244 indicates the dimensional constants .
the full set of dimensionless equations to be integrated are :    @xmath245    the dimensional constants are given by the set of equations :    @xmath246    it is natural to choose : @xmath247 since @xmath248 is the factor in the expressions for the pressure ( eq . [ press ] ) and energy ( eq . [ eose ] ) . hence :    @xmath249    moreover the set of equations ( [ lambdaadim])-([massadim ] ) are invariant under the transformations : @xmath250 , @xmath251 , and @xmath252 , etc ( the other constants can be obtained from these ) , with @xmath253 an arbitrary number different from zero .
so , it is possible to choose any other set of constants consistent with these transformations . in particular ,
the dimensional constants obtained by oppenheimer and volkov @xcite , are recovered by setting : @xmath254 . in this case : @xmath255 cm ; @xmath256 g , etc , equivalent to the constants obtained in that paper @xcite .    in the eq .
( [ tovadim ] ) the term @xmath257 , is obtained deriving the eq .
( [ press ] ) : @xmath258    it is used a @xmath259 order runge - kutta scheme to integrate simultaneously the set of equations ( [ lambdaadim])-([massadim ] ) .
the charge distribution is chosen proportional to the rest mass distribution : @xmath260 .
for simplicity , in this paper we will concentrate on the case @xmath261 .
the code implemented to integrate the equations above is called he05v1 .
the he05v1 is used to build neutron star models in hydrostatic equilibrium .
this models constitutes an initial data set for another code presented below .
in this section we collect the equations that must be integrated numerically in order to simulate the temporal evolution of the stellar models constructed with the he05v1 code .
the set of equations is :    [ allmotion ] @xmath262\,dt \hspace{1.56cm}[\ref{motion}]\\ % & & r_{,t}=u\,a \\ & & w=1+\frac{\epsilon}{c^2}+\frac{p}{\rho\ , c^2 } \hspace{2.96cm}[\ref{enthalpy } ] \\ & & \gamma^2=1+\frac{u^2}{c^2}-\frac{2 m g}{r c^2 } + \frac{g q^2}{c^4 r^2 } \hspace{1.335cm}[\ref{gammahe } ] \\ & & d\mu=4 \,\pi\ , \rho\ , r^2\ , dr/\gamma\hspace{2.8cm}[\ref{gammadef } ] \\ & & { a w}={a_0 w_0}\,{\rm exp}\,\biggl[\,\int_0^\mu\,\biggl ( d\epsilon + p \,d \biggl(\frac{1}{\rho } \biggr)+ \nonumber\\ & & \hspace{0.8 cm } \frac{q \,d q}{4\ , \pi\ , r^4\ , \rho }   \biggr ) /w\,c^2\biggr ] \hspace{2.45cm}[\ref{gtt2}]\\ & & dm=4\ , \pi\ ,   \rho\ , \biggl(1+\frac{\epsilon}{c^2 } \biggr ) \,r^2\ , r_{,\mu}\ , d\mu + \nonumber \label{hamilton}\\   & & \hspace{0.8 cm } \frac{1}{c^2 } \frac{q q_{,\mu}}{r } d\mu   \hspace{3.25cm}[\ref{totalmass}]\\   & & { \rho\,r^2}={\rho_0\ , r_0 ^ 2}\,\,{\rm exp}\ , \biggr[-\int_0^t a \frac{u_{,\mu}}{r_{,\mu } } \,dt\biggl ] \label{momentum } \hspace{0.57cm}[\ref{massconservation}]\\ & & d\epsilon =- p\ , d\biggl(\frac{1}{\rho}\biggr ) \hspace{3.39cm}[\ref{1stlaw}]\\ & & dq=4\ , \pi\ , \rho_{ch } \,r^2\ , dr/\gamma \label{egamma}\\   & & d{\mu}=4 \,\pi\ , \rho\,{r}^2\,dr/\gamma\,.\hspace{2.7cm}[\ref{gammadef } ] \label{emu}\end{aligned}\ ] ]    where we mean @xmath263 and @xmath264 .
the references to other equations are enclosed in square brackets . comparing the eqs .
( [ allmotion ] ) with the may and white equations @xcite , it can be seen that they share a very similar structure when @xmath265 .
we have chosen a charge distribution proportional to the rest mass distribution so we obtain the eq .
( [ egamma ] ) . the eq .
( [ dynam ] ) is the only equation that have a second time derivative and constitutes the dynamical part of the einstein - maxwell equations .
( [ hamilton ] ) and ( [ momentum ] ) represent the hamiltonian and momentum constraint equations , respectively .
the numerical code to integrate the equations above are called .
the initial conditions at @xmath266 are obtained with the he05v1 for each stellar model , and consist on : the initial mass density distribution @xmath267 ; the initial electric charge density distribution @xmath268 ; the kinetic energy @xmath269 ; the cell spacing @xmath270 and the surface radius @xmath271 . here
, @xmath272 is the mass coordinate of the surface .
then , the output of the he05v1 is used as an initial cauchy data set in the space - time hypersurface @xmath266 , @xmath273 , @xmath274 , @xmath275 , for the .
the boundary conditions are :    @xmath276    the condition given in eq .
( [ bca ] ) , makes the time coordinate synchronized with an observer co - moving with the surface of the star .
the conditions expressed in the eqs .
( [ bcp])-([bcr ] ) , also gives : @xmath277    the finite difference algorithm was implemented using a method similar to that developed by may and white @xcite .
( [ me1])-([emu ] ) were integrated with a leap - frog finite difference method plus a predictor - corrector step .
we call this numerical code as .
the integration is iterated according to a desired error control .
the method is second order accurate in space and time .
each experiment was repeated with different number of particles and with different values of the numerical viscosity parameter , in order to check the convergence of the results . in general , good results are obtained using @xmath278 points along the coordinate @xmath63 .
the compatibility between the two codes is an indirect check of the convergence properties of the , since the he05v1 is @xmath279 order accurate .
for example , at @xmath266 there exists differences in the integrated mass between the two codes of the order of @xmath280 , using @xmath281 points in the .
this difference can be made @xmath282 by taking , instead , @xmath283 points , and so on .
for the integration in the time dimension is needed to take initially small time - steps @xmath284 s , but the code is capable of self - adjust the time - step so as to take small changes in the physical variables .
therefore , with this algorithm it is indifferent to take an initial time - step of @xmath285 s or @xmath286 s , the code always self - adjust the time - step and prove to converge always to the same solution .
the shocks are treated by adding an artificial viscosity , which is non - zero only on discontinuities .
its effects are to smear out the discontinuities over several cells , and to reduce the post - shock oscillations and the numerical errors related to the leap - frog method .
we experienced with distinct functional forms for the viscosity term , probing the method of may and white @xcite to be the best in the present algorithm .
however , there are not strong shocks ( like in ref .
@xcite ) in the set of simulations studied in this paper and with the particular initial conditions chosen .    in the the eqs .
( [ egamma ] ) and ( [ emu ] ) are not integrated in time , and the rest mass and charge are constrained to be constant in the layers of matter ( as it must be ) .
we check that this kind of algorithm reduce the numerical errors .    as was showed in the sec .
( [ relatequ ] ) the total mass - energy at the surface of the star is also a constant @xmath287 .
the numerical errors have an impact on the constancy of the total mass - energy @xmath5 , and this is the physical quantity with which is checked the accuracy of the simulation . in all the runs
@xmath5 is very well conserved , with the exception at the time when a trapped surface forms . in that case , very large gradients forms and serious numerical instabilities appear .
although , it is possible to manage the code to keep the error in mass - energy conservation below @xmath288 at the time the apparent horizon forms .
of course , once the trapped surface forms nothing can change the fate of the matter inside the apparent horizon ( see the sec .
[ subsubsec : ah ] ) .
the simulation must be discontinued soon after the matter cross the gravitational radius because the time - step attains very small values . in general ,
when @xmath289 is roughly @xmath290 s , it is needed a very large number of time - steps to obtain further progress and this consume too much computer time .
the reason is that the large gradients produce huge changes in the physical variables and a very small time - step is required to keep the errors at a low level .
thus , it is not possible to follow the dynamical evolution until a stationary regime is reached ( or the formation of an event horizon , if ever possible ) . when the outer trapped surface forms , some layers of matter external to it can be collapsing or expanding and is not possible in general to say which will be their fate .
however , the matter that falls in the trapped surface can not escape from it and this is enough to say that the star ( or some part of it ) collapsed .
there are two different sources of perturbation introduced in the models : the first source of perturbation is of numerical origin and cames from the mapping of the equilibrium models in the evolution code .
this mapping is not exact , and there is a difference in the integrated rest mass and total mass between the two codes that is taken @xmath291 .
this difference can be made as small as wanted by taking a higher number of points to perform the integration with the ; the second source of perturbation , is introduced by multiplying the kinetic energy by a small constant grater than one , i.e. , @xmath292 . in the present simulations we choose @xmath253 to take the values @xmath293 .
the structure of the space - time of a collapsed charged star is very complex , with the possible existence of time - like singularities and tunnels connecting several disjoint asymptotically flat space - time regions .
this tunnels and the true singularities , pertains to the analytically extended portions of the space - time , and the present code is not prepared to reach that regions of the total manifold .
in fact , with the present code is not possible to reach the cauchy horizon where important physical effects must occur @xcite , @xcite .
however , from an astrophysical point of view is interesting to know what happens outside the apparent horizon , where the astronomers live .
we will not describe the numerical methods in more detail because is not the objective of this paper .
moreover , the numerical methods are very well known and better explained in textbooks and papers ( see for example ref .
the is an extended version of a preliminary code developed by ghezzi ( 2003 ) , earlier test - beds and results obtained with this code were published in collaboration @xcite , @xcite .
.[tab : table1]mass , radius , and central density for neutron star models with zero charge . [ cols=">,>,^,^,^,^,^,^ " , ]      a trapped surface is a space - like orientable two dimensional compact surface , such that inward and outward null geodesics normal to it converge .
the apparent horizon is the outer boundary of all the trapped surfaces ( see @xcite , @xcite , @xcite , @xcite , @xcite , and @xcite , for definitions and related theorems ) .
a black hole in asymptotically flat space - time is the region from where no causal or light signals can reach @xmath294 ( the future null infinity , see for example @xcite and @xcite for definitions ) .
the event horizon @xmath295 is the boundary of the black hole region @xmath296 .
the formation of an apparent horizon is a sufficient , although not necessary , condition for the formation of a black hole ( see ref .
in fact , in the case an apparent horizon forms the theorems of hawking and penrose ( 1970 ) can be applied to know that the outcome will be the formation of a singularity , i.e. : the spacetime contains at least one incomplete timelike or null geodesics ( see the theorems 9.5.3 and 9.5.4 in @xcite , pags
. 239 - 241 , known as `` singularity theorems '' ) .
moreover , the trapped surface ( and the apparent horizon ) is contained within a black hole as a mathematical proposition asserts ( see the propositions 12.2.2 and 12.2.3 in @xcite , pags .
309 - 310 , and proposition 9.2.1 in @xcite , pag .
so , we assumed that if an apparent horizon is formed in the stellar collapse , the result will be the inexorable formation of a black hole .    in a numeric simulation , in general ,
is needed an algorithm to find out the apparent horizon s location . however , in the present study it is easy to guess that the trapped surfaces are spheres . with this in mind , we can describe the apparent horizon formation and evolution .
the vector @xmath297 is tangent to the outgoing null geodesics , and @xmath298 is tangent to the ingoing null geodesics , in co - moving coordinates .
they are also orthogonal to the constant @xmath299 surfaces .
following mashhoon and partovi @xcite , we define the quantities @xmath300 and @xmath301 which are related by a positive multiplicative factor to the expansion factor for radially ingoing and outgoing null geodesics , respectively .
thus , the coordinates of the apparent horizon are found solving the equation : @xmath302 using the eqs . ( [ gammadef ] ) and ( [ udefin ] ) , and squaring , this is equivalent to @xmath303 and using eq .
( [ gamma2 ] ) , this gives : @xmath304 the solution of this equation is : @xmath305 where @xmath306 is the inner boundary , while @xmath307 is the outer boundary of the trapped surfaces .
so , the apparent horizon is the surface ( in schwarszchild coordinates ) : @xmath308 in order to know the space - time coordinates of the apparent horizon in co - moving coordinates , we must solve the eq .
( [ ahrp ] ) for @xmath46 and @xmath63 .
this is done with a numeric subroutine that for each time @xmath46 and for each coordinate @xmath63 asks if @xmath309 .
some of the stars when evolved in time give rise to the formation of an apparent horizon , that is the coordinate @xmath28 contracted to the value @xmath307 for some @xmath46 and @xmath63 .
the table [ tab : table4 ] summarize the results of the simulations : it is shown that some of the stars collapsed , and it is given the time elapsed from the beginning of the simulation and the formation of the apparent horizon .
( [ apparenteq ] ) is also equivalent to : @xmath310 this equation can also be obtained from the equation for light rays @xmath311 . from this equation
we see that outgoing light rays have zero expansion in the co - moving frame when @xmath312 .
however , from the simulation was obtained that the apparent horizon forms before @xmath313 , this means that radial light rays enter the trapped region ( at its boundary @xmath314 ) with @xmath315 .    in the figs .
( [ fig : v08m21]b ) and ( [ fig : v05m21]b ) is shown the evolution of the metric coefficient @xmath61 for model number @xmath316 with @xmath317 , and for model number @xmath316 with @xmath318 ( this model is one of the entries of the table [ tab : table2 ] ) . as we see , @xmath319 for this models , so they collapse .
the event horizon can not be described in the simulations , since it is the result of the complete history of the collapsing matter ( while the simulation lasts a finite amount of time ) .
but we know that the trapped region is contained in the black hole region .
there exists only one situation in which the formation of the event horizon can occur in a finite coordinate time : when all the matter composing the star collapses in a finite time . only in the special case that the star collapse completely , i.e. : if @xmath320 for some finite time @xmath321 , the apparent horizon will coincide with the event horizon of a reissner - nordstrm space - time : @xmath322 , where @xmath323 is the reissner - nordstrm radius of a black hole of total mass - energy @xmath5 and total charge @xmath2 . in all the simulations
performed the apparent horizon is formed for some coordinate @xmath63 with @xmath324 .
so , in the present set of simulations , the apparent horizon never coincides with the event horizon .
we picked up some stable stars and followed their temporal evolution after a sudden discharge of its electric field .
the scenario is a pure hypothetical one , in which the neutron star is leaved in a charged metastable state after its formation , or during an accretion process onto it .
the charged matter could suffer a discharge by recombination of charges of opposite sign , or by effect of the high electrical conductivity .
we emphasize that we are not claiming that this is a possible astrophysical scenario , but we want to explore the different theoretical possibilities for the stellar dynamics .
the simulation is performed with the worst conditions for the stellar stability : the electric field is simply turned off after the first time - steps of temporal evolution of an otherwise stable star .
the energy is maintained constant in the process , corresponding to a conversion of the electromagnetic energy into heat or internal energy .
this is easily performed with the code .
the results of the simulation of the stars suddenly discharged are summarized in the table [ tab : table4 ] .
the models discharged are indicated with primed numbers .
we see that not all the models collapse to form black holes : the model number @xmath325 with an initial charge to mass ratio @xmath326 oscillate , while the same model with @xmath327 collapsed after a sudden discharge ; the model number @xmath328 oscillate if its initial charge to mass ratio is @xmath326 , while the same model collapses if @xmath327 .
the models number @xmath316 always collapse , as they fall on the unstable branch , then it is unnecessary to consider them .      in the reference
@xcite it is described the formation of a shell of higher density formed near the surface of the star . although this effect must happen in newtonian physics ,
its evolution in the strong field regime is highly non - linear and far from obvious .
the weight of the star is supported by the gas pressure and by the coulomb repulsion of the matter .
if the coulomb repulsion is important respect to the gravitational attraction -although not necessarily stronger- a shell of matter can form ( see @xcite for an explanation ) .
the contrast of density between the shell and its interior depends on the energy and charge distribution .
so , it is important to check if the effect take place in the present simulations .
we found that the shell forms only mildly , or do not form , and its late behavior can not be followed clearly .
from all the numerical experiments performed , it is observed that the shell formation effect arise more clearly when the initial density profile is flat ( like the simulations of ref .
@xcite ) . in consequence ,
the shell formation is not an important physical effect in the collapse of charged neutron stars .
however , it could be important for the collapse of the core of super - massive stars @xcite .
@xmath329 \epsfig{width=7cm , height=6 cm , file = prdiimassrad0.eps } & \epsfig{width=7cm , height=6 cm , file = prdiimassrad3.eps }   \\ [ 0.4 cm ] \mbox{\bf ( a ) } & \mbox{\bf ( b ) } \end{array}$ ]    @xmath329 \epsfig{width=7cm , height=6 cm , file = prdiibindingdens.eps } & \epsfig{width=7cm , height=6 cm , file = prdiimassadens3.eps } \\ [ 0.4 cm ] \mbox{\bf ( a ) } & \mbox{\bf ( b ) } \end{array}$ ]    @xmath329 \epsfig{width=7cm , height=6 cm , file = prdiivel08m21.eps } & \epsfig{width=7cm , height=6 cm , file = prdiigtt08m21.eps } \\ [ 0.4 cm ] \mbox{\bf ( a ) } & \mbox{\bf ( b ) } \end{array}$ ]    @xmath329 \epsfig{width=7cm , height=6 cm , file = prdiivel05m21b.eps } & \epsfig{width=7cm , height=6 cm , file = prdiigtt05m21.eps } \\ [ 0.4 cm ] \mbox{\bf ( a ) } & \mbox{\bf ( b ) } \end{array}$ ]    @xmath329 \epsfig{width=7cm , height=6 cm , file = prdiigamma08m21.eps } & \epsfig{width=7cm , height=6 cm , file = prdiigamma05m21b.eps } \\ [ 0.4 cm ] \mbox{\bf ( a ) } & \mbox{\bf ( b ) } \end{array}$ ]    @xmath329 \epsfig{width=7cm , height=6 cm , file = prdiivel05m13.eps } & \epsfig{width=7cm , height=6 cm , file = prdiigtt05m13.eps } \\ [ 0.4 cm ] \mbox{\bf ( a ) } & \mbox{\bf ( b ) } \end{array}$ ]
it is usually assumed in astrophysics that stars havent important internal electric fields . whether an star can have large internal electric fields or
a net total charge is not yet clear .
however , several features could be common to the evolution of rotating collapsing stars , with the angular momentum playing the role of the electric charge .
so , the present study can shed some light on more realistic astrophysical scenarios .    in the models studied in this paper we assumed , for simplicity , that the charge density is proportional to the rest mass density .
we found that in hydrostatic equilibrium the charged stars have a larger mass and radius than the uncharged ones .
this is , as expected , due to the coulomb repulsion .
the mass of the models tends to infinity as the charge approaches the extremal value @xmath330 .
the hydrostatic equilibrium solutions with @xmath331 gives models with an infinite mass and radius ( see also @xcite ) .
this means that in this particular cases the integrated mass and radius diverges within the computer capacity .
all the models with charge less than the extremal ( @xmath332 ) , have a mass limit and there are unstable and stable solutions .
we checked the stability of the solutions integrating forward in time the models , and applying strong perturbations to them .
some of the models collapse directly to form black holes , without ejecting matter .
other models oscillate . for a given model , with fixed central density ,
the frequency of oscillation is lower when the charge is higher .
the frequency of oscillation is weakly dependent on the numerical viscosity .
the models that collapse are solutions with @xmath333 , while the oscillating models have @xmath334 .
so , the stability of the models agree with the predictions of the first order perturbation theory @xcite .
it seems that there is a limit for the charge that a star can have , which is lower than the extremal case .
this limit arise because the binding energy per nucleon of the models with @xmath335 , is lower than the binding energy per nucleon on an atomic nucleus .
then , it is possible that this stars disintegrate to reach a more bounded energy state .
this point deserves further study .    from a pure theoretical point of view , the issue of the collapse of charged fluid spheres , is related to the third law of black hole thermodynamics @xcite , and with the cosmic censorship hypothesis @xcite , @xcite .
the third law of black hole physics states that the temperature of a black hole can not be reduced to zero by a finite number of operations .
the impossibility of transforming a black hole into an extremal one , in a finite number of steps , is related to the impossibility of getting @xmath330 in some particular experiment @xcite , @xcite . in agreement with this law
, we found that it is not possible to form extremal black holes from the collapse of a charged fluid ball .
in fact , any charged ball with a charge to mass ratio greater or equal than one explodes , or its matter spreads .
the black holes with @xmath331 represent naked singularities , and the impossibility of getting black holes with this charge to mass ratio in the present simulations are in agreement with the  cosmic censorship hypothesis " .
it must be remarked that once an apparent horizon forms the formation of an event horizon is inevitable , and as we proved the matching of the interior solution with an exterior reissner - nordstrm solution , we simulated here the formation of a reissner - nordstrm space - time from a gravitational stellar collapse .
although we are not using a realistic equation of state we think that the results are of general validity , at least from a qualitative point of view .    other fields or more exotic physics
must be considered in order to form extremal black holes .
for completeness we reproduce here the einstein equations in co - moving coordinates as derived by landau and lifshitz ( see @xcite , pag .
311 ) . the equations of section [ sec2 ] were derived from these by making the replacements @xmath336 , @xmath337 , and @xmath338 .
the einstein equations are : @xmath339
the author thank professor dr .
jacob d. for clarifications about his equations .
this work was performed with the computational resources of professor dr .
patricio s. , so the author thank him .
the author thank an anonymous referee for useful comments that help to improve this paper .
the author acknowledge the brazilian agency fapesp for support .


the isoperimetric problem is a classical subject in differential geometry with its origin in ancient greece .
it consists in finding on a riemannian manifold @xmath2 the regions that minimize the perimeter among sets enclosing the same volume .
the solutions are called isoperimetric regions and their boundaries isoperimetric hypersurfaces .
the euclidean plane is historically the first space where the problem started to be investigated rigorously .
it is now a well known fact that the round circles are the optimal curves for the problem .
this geometric fact is often seen through the following classical inequality : @xmath3 where @xmath4 and @xmath5 stand for the length and the enclosed area of a simple closed curve @xmath6 respectively .    the framework of geometric measure theory and its tools work successfully well in tackling the aspects of existence and regularity of this variational problem .
when @xmath7 is closed or homogeneous , then isoperimetric hypersurfaces do exist and are smooth up to a closed set of hausdorff dimension @xmath8 .
the regular part is a stable hypersurface of constant mean curvature .
this major contribution to the isoperimetric problem was achieved thanks to the efforts of many people , including f. almgren , r. schoen , l. simon , f. morgan , and others ( see @xcite for a comprehensive list ) . despite the long history of the problem
, it remains largely open with few 3-manifolds where the problem is completely understood .
the simply connected space forms , @xmath9 , @xmath10 and @xmath11 , are the most appealing spaces to begin the study of the isoperimetric problem .
it turns out that their symmetries are enough to characterize the geodesic spheres as the isoperimetric hypersurfaces .
a complete solution on @xmath12 with the standard product metric can be found in @xcite .
for other homogeneous manifolds with certain product structure , such as @xmath13 , @xmath14 , and @xmath15 , see @xcite , @xcite , and @xcite respectively . the case @xmath16 is also treated in @xcite and
they show that when @xmath17 unduloids are minimizers rather than cylinders for certain volumes .
one key idea exploited in the results listed above is the use of symmetry to reduce the problem to an ode analysis .
the case @xmath18 where @xmath19 is a flat torus is not solved in full ; great progress can be found in @xcite and in @xcite .
more generally , it is known that boundaries of small isoperimetric regions in closed manifolds are nearly round spheres , see @xcite .
to finish this brief and not exhaustive account of results on the isoperimetric problem we mention that bray and morgan ( @xcite and @xcite ) classified the horizon homologous isoperimetric surfaces in the schwarzschild manifold .
the work in @xcite highlighted the interesting relationship between isoperimetric surfaces and the concept of mass in general relativity .
we will be interested in spherical space forms in this paper .
a significant result in this direction was given in @xcite where a. ros and m. ritor solved the isoperimetric problem in the projective space @xmath20 .
they show that the solutions are geodesic spheres or flat tori .
later , a. ros @xcite used the above result to give a proof of the willmore conjecture in @xmath21 for the special case of surfaces that are invariant by the antipodal map .
the real projective space is a special case of an important family of riemannian manifolds , namely the lens spaces @xmath22 .
these are elliptic space forms obtained as a quotient of @xmath21 by a finite group of isometries that are isomorphic to @xmath23 but which also depend on @xmath24 .
they are , along with @xmath25 , characterized by having heegaard genus one .
we give a complete solution for the isoperimetric problem in the lens spaces with large fundamental group :    [ theorem 1 ] _ there exists a positive integer @xmath26 such that for every @xmath27 and every @xmath28 the isoperimetric surfaces in @xmath22 are either geodesic spheres or tori of revolution about geodesics . _    the isoperimetric problem in dimension three was previously solved for only a finite number of non - diffeomorphic 3-manifolds .
we also add to the literature the observation that the proof of the willmore conjecture by marques and neves @xcite can be applied to extend the work of ros and ritor @xcite , on the classification of stable cmc surfaces in @xmath20 , to @xmath0 and @xmath1 : and @xmath1 are diffeomorphic but not isometric . ]    [ theorem 2 ] _ the immersed stable cmc surfaces in @xmath0 and @xmath1 are either geodesic spheres or flat tori .
moreover , the minimal clifford torus is , up to ambient isometries , the only index one minimal surface in @xmath0 and @xmath1 .
_    the idea of the proof of theorem [ theorem 1 ] is as follows .
stability implies that every isoperimetric surface is connected and its genus is @xmath29 or @xmath30 .
if follows from a classical result of hopf that if the genus is 0 , then it is a geodesic sphere .
from @xcite we know that if the genus is 1 , then it is flat , and this forces the surface to be of revolution about a geodesic .
we are left to rule out other topological types .
to do so we argue by contradiction .
we assume that there exists a sequence of lens spaces with increasing fundamental group containing isoperimetric surfaces of genus @xmath31 .
after a suitable rescaling on the metrics we use compactness results to obtain a limit for the sequence of lens spaces which will be a flat three manifold of rank one . in the same way , the sequence of isoperimetric surfaces will converge to a flat surface in the respective ambient manifold .
the topology of the surfaces will force the limit to be an union of planes . on the other hand , the minimization property rules this configuration out .
the arguments in the proof of theorem [ theorem 1 ] generalize naturally for the berger spheres @xmath32 .
this is a well known one parameter family of homogeneous metrics on the 3-sphere ; the case @xmath33 corresponds to the round metric .
all the hopf fibers in @xmath32 have the same length @xmath34 .
when @xmath35 is not too small , spheres are the only solutions of the isoperimetric problem , see @xcite . when @xmath35 is very small , it is pointed out in @xcite that some tori are better candidates to solve the isoperimetric problem rather than spheres for certain volumes .
_ there exists @xmath36 such that for every @xmath37 the isoperimetric surfaces in the berger spheres @xmath32 are either rotationally invariant spheres or tori . _
it follows from the work in @xcite and @xcite that the minimal clifford torus in @xmath38 ( resp .
@xmath32 ) is isoperimetric for every @xmath39 ( resp .
@xmath40 ) , see section 3 .
_ acknowledgements .
_ i am grateful to my advisor andr neves for his constant encouragement and support during the course of this work .
i also thank jason lotay for many helpful comments on an earlier version of this manuscript .
let @xmath41 be an orientable riemannian manifold of dimension @xmath42 .
the @xmath42-dimensional hausdorff measure of a region @xmath43 is denoted by @xmath44 .
similarly , we denote the @xmath45-dimensional hausdorff measure of the hypersurface @xmath46 by @xmath47 .
the class of regions considered here are those of finite perimeter , see @xcite .
a region @xmath43 is called an _ isoperimetric region _
if @xmath48 in this case , the hypersurface @xmath49 is called an _ isoperimetric hypersurface_.    in general , the existence of isoperimetric surfaces is handled by a compactness theorem from geometric measure theory . for non - compact manifolds one needs to be careful since a minimizing sequence of regions of fixed volume may drift off to infinity .
we recommend @xcite for a recent reference on the regularity of isoperimetric hypersurfaces :    _ let @xmath41 be a closed riemannian manifold .
for every @xmath50 there exists an isoperimetric region @xmath51 satisfying @xmath52 .
moreover , @xmath49 is smooth up to a closed set of hausdorff dimension @xmath8 .
_      isoperimetric hypersurfaces are stable critical points of the area functional for variations that preserve the enclosed volume ; thus , the regular part of isoperimetric hypersurfaces has constant mean curvature .
more generally , we say an isometric immersion @xmath53 has _ constant mean curvature _
( cmc ) if it is a critical point of the area functional for volume preserving variations .
critical points are called _ stable cmc _ if the second derivative of the area is non - negative for such variations .
equivalently , @xmath54 is stable if for every @xmath55 with compact support such that @xmath56 , we have @xmath57 @xmath58 is the unit normal vector of @xmath59 and @xmath5 is the _ second fundamental form _ of the immersion @xmath54
mean curvature _ of @xmath59 , denoted by @xmath60 , is defined by @xmath61 .
the study of immersed stable cmc hypersurfaces started in @xcite and @xcite with a new characterization of the geodesic spheres in the simply connected space forms @xmath10 , @xmath9 and @xmath11 . the classification of stable cmc surfaces is often a way to approach the isoperimetric problem in reasonable spaces . with this purpose in mind ,
a. ros and m. ritor in @xcite used the hersch - yau trick to study orientable stable cmc surfaces on three manifolds of positive ricci curvature .
[ up bound w energy , g=2 ] _ let @xmath62 be a closed three manifold with positive ricci curvature . if @xmath63 is a stable cmc immersion , then @xmath64 .
moreover , if @xmath65 or @xmath30 , then @xmath66 _    [ wilmore energy genus 2 ]    the first statement follows from the previous theorem since @xmath67 .
let @xmath68 be the induced map in fundamental groups .
as @xmath69 has finite index there exists a finite covering @xmath70 such that @xmath71 and @xmath72 .
this means there exists a lifting of this map into @xmath21 and we denote it by @xmath73 .
it follows that @xmath74 .
if @xmath54 is not an embedding , then @xmath75 is not embedding either . by the work of li and yau @xcite
the willmore energy of the immersed surface @xmath76 , i.e. @xmath77 , is strictly greater is discussed in @xcite . ] than @xmath78 .
therefore , if @xmath79 , we obtain a contradiction and @xmath75 is an embedding .
moreover , for closed surfaces with genus greater than or equal to @xmath80 in @xmath21 the willmore conjecture , recently proved in @xcite , states that @xmath81 .
let @xmath82 be the pre - image of @xmath59 by the universal covering map , then @xmath83 therefore , if @xmath84 , then @xmath85 and @xmath86 .
in particular , if @xmath87 or @xmath30 , then there exist no stable cmc surface @xmath59 with @xmath88 .    for each @xmath89
we define the clifford torus @xmath90 as : @xmath91    every flat tori with constant mean curvature in @xmath21 is congruent to a clifford torus @xmath90 .
this follows from the rigidity theorem , pg .
49 in @xcite , for they have the same second fundamental form .    _
the stable cmc surfaces in @xmath0 and @xmath1 are totally umbilical spheres or flat tori .
in addition , the index one minimal surfaces in @xmath0 and @xmath1 are congruent to the the projection of minimal clifford torus .
_    let @xmath92 , @xmath93 , be in the conditions of the corollary . by proposition [ wilmore energy genus 2 ] , @xmath94 or @xmath80 .
if @xmath94 , then it follows from the hopf holomorphic quadratic differential that @xmath59 is totally umbilical .
if @xmath95 , then it is proved in @xcite that @xmath59 is flat .
the isoperimetric properties of @xmath2 can be encapsulated in a single function called the _
isoperimetric profile_. this is the function @xmath96 \rightarrow [ 0,+ \infty)$ ] defined by @xmath97 we finish the section with some well known facts on the analytic nature of @xmath98 .
these will be used later in section 3 .
let @xmath51 be an isoperimetric region in @xmath2 such that @xmath99 for some @xmath100 . the function @xmath98 has left and right derivative @xmath101 and @xmath102 . in addition , if @xmath60 is the mean curvature of @xmath49 in the direction of the inward unit vector , then @xmath103 the second derivative also exists but weakly in the sense of comparison functions .
more precisely , we say @xmath104 weakly at @xmath105 if there exists a smooth function @xmath106 such that @xmath107 , @xmath108 , and @xmath109 . in this sense
we have @xmath110 the equations ( [ first derivative ] ) and ( [ second derivative ] ) are first presented on @xcite .
we sketch the proof of ( [ first derivative ] ) and ( [ second derivative ] ) .
let @xmath111 be the variation @xmath112 of @xmath59 reparametrized in terms of the enclosed volume @xmath113 .
in addition , let @xmath114 ( resp .
@xmath115 ) be the area of @xmath116 ( resp .
@xmath117 ) . by the first variation formula for the area and volume we have @xmath118 and @xmath119 respectively . since @xmath120 , we conclude that @xmath121 and also that @xmath122 . on the other hand ,
the second derivative of area for general variations implies the following : @xmath123 hence , in the sense of comparison functions , ( [ second derivative ] ) follows from : @xmath124
in order to define the lens spaces , we first recall the round three sphere as : @xmath125 fix @xmath126 integers with the following property @xmath127 and @xmath128 . let @xmath23 be the group @xmath129 acting on @xmath21 as follows : @xmath130 the group @xmath23 acts freely and properly discontinuous on @xmath21 .
the orbit space @xmath131 is a closed three manifold called the _ lens space _ , it is denoted by @xmath22 .
the hopf fibration , which is the riemannian submersion @xmath132 defined by @xmath133 , can be extended naturally to @xmath22 .
indeed , the group @xmath23 acts on the set of _ hopf fibers _ through the following cyclic action on @xmath134 : @xmath135 the hopf fibration for @xmath22 is then defined as @xmath136 .
the set @xmath137 is a two dimensional orbifold with conical singularities at the north and south pole .
the preimage of horizontal slices of @xmath137 via @xmath138 correspond to the clifford torus described in ( [ clifford torus ] ) .
they are natural candidates to solve the isoperimetric problem in @xmath22 .
steiner and schwarz symmetrization theorems were proved in @xcite for certain fiber bundles such as the lens spaces . to explain this symmetrization procedure we restrict to the case @xmath38 where the hopf fibration @xmath139 is a smooth riemannian submersion .    the symmetrization consist of associating to each set of finite perimeter @xmath140 the set @xmath141 in the product manifold @xmath142 defined by replacing the slice of @xmath143 in each fiber with a ball of the same volume about the respective fiber in the product .
the coarea formula for riemannian submersions implies that @xmath141 encloses the same amount of volume as @xmath143 .
it is proved in @xcite that @xmath141 has no greater perimeter than @xmath143 .
one immediate consequence is that @xmath144 . applying the classification of the isoperimetric problem on @xmath142 , @xcite ,
we conclude that @xmath145 in a interval around @xmath146 . in particular , the minimal clifford torus is isoperimetric in @xmath38 for every @xmath39 .
the isoperimetric profiles , however , do not coincide as the profile of geodesic spheres on the respective spaces are different . therefore , this technique is not enough to completely solve the isoperimetric problem .
d_{\mathbb{s}^3}[(\theta^{kq } z ,
\theta^k w),(z , w)]&=&2\arcsin(\frac{a}{2})\geq 2\arcsin\big({\frac{2}{2}\sin{\frac{2\pi}{2p}}}\big)=\frac{2\pi}{p}.\end{aligned}\ ] ] however , is not true in general that @xmath153 as @xmath154 .
let s consider @xmath155 , @xmath156 .
we show that the injectivity radius at points far away from the critical fibres are @xmath157 .
if the round metric is scaled by the factor @xmath158 , then we have the riemannian submersion : @xmath159 the fibers have constant length @xmath160 except the critical fibres which have length @xmath161 .
the right hand side will converge , as @xmath162 , to @xmath163 .
it follows from the coarea formula for riemannian submersions that the volume of the geodesic ball @xmath164 in @xmath165 is bounded from below .
therefore , by cheeger s inequality , lemma 51 in @xcite , the injectivity radius of the sequence @xmath165 is bounded from below .
this sequence converges to a flat @xmath166 .
let @xmath173 and recall that @xmath174 .
without loss of generality , let s assume that @xmath175 . hence , under the rescaled metric @xmath176 , the extrinsic diameter of @xmath170 is greater than or equal to @xmath177 .
let @xmath178 be a geodesic segment realizing the intrinsic diameter of @xmath170 .
thus , we can find disjoint balls @xmath179 , with @xmath180 , @xmath181 , and @xmath182 + 1 $ ] .
hence , @xmath183 + 1}\mathcal{h}^2(b_r(x_{i_0})\cap t_{\frac{\pi}{4}}/\mathbb{z}_p ) \leq |t_{\frac{\pi}{4}}/\mathbb{z}_p|= \lambda_p^2\,\frac{2\pi^2}{p}.\ ] ] therefore , there exists @xmath184 + 1\}$ ] such that : @xmath185 the first inequality follows from the monotonicity formula , proposition [ monotonicity ] in the appendix , applied to @xmath186 .
this finishes the proof of the lemma .
let s use the notation @xmath187 to represent the group @xmath23 acting on @xmath21 and its dependence on the parameter @xmath24 . by using the toroidal
coordinate system for @xmath21 , @xmath188 the action of @xmath189 on @xmath90 corresponds to the following action on @xmath190 : @xmath191 in these coordinates , the @xmath189 orbit at the point @xmath192 is given by : @xmath193      to prove the lemma it is enough to consider @xmath194 .
if there is a subsequence for which the diameter of @xmath195 is going to zero as @xmath154 , then @xmath196 is clearly becoming dense on @xmath197 and item 1 is proved .
let s consider now the case where the diameter of @xmath197 in @xmath22 is bounded away from zero . from the equivalence between extrinsic and intrinsic distance and by lemma [ injective radius ]
we conclude that the euclidean injectivity radius satisfies @xmath198 .
in particular , there exist @xmath199 such that @xmath200 therefore , there exists @xmath201 such that @xmath202 infinitely often and @xmath203 is the euclidean injectivity radius of @xmath170 at @xmath204 for this subsequence .
hence , the sub - orbit generated by the translation @xmath205 is contained in the line @xmath206 .
it follows that the @xmath207 is contained in a union of equidistant lines parallel to the one described above by homogeneity .
modulo @xmath208 the number of such lines is finite , let s denote it by @xmath209 .
modulo @xmath208 there are @xmath210 points of @xmath207 in each of these lines .
hence , @xmath211 therefore , @xmath209 divides @xmath212 and is independent of @xmath213 .
in other words , @xmath196 is contained in @xmath214 integral curves of @xmath215 .
[ main theorem ] _ there exists a positive integer @xmath26 such that for every @xmath27 and every @xmath28 the isoperimetric surfaces in @xmath22 are either geodesic spheres or tori of revolution about geodesics .
_    [ proof of theorem 1.1 ] arguing by contradiction , let s assume that there exists an infinite sequence of lens spaces @xmath22 containing an isoperimetric surface @xmath216 of genus @xmath217 or @xmath30 for each @xmath213 .        1 .
there exists an exhaustion of @xmath2 by compacts @xmath220 , i.e. @xmath221 and @xmath222 .
2 .   there exists a family of diffeomorphism onto their images , @xmath223 , such that @xmath224 and @xmath225 in the @xmath226 topology .
we consider pointed manifolds @xmath227 with base points @xmath228 belonging to @xmath229 . by the cheeger - gromov compactness theorem
we have that @xmath230 , where @xmath231 is a flat three manifold . the inclusion of @xmath216 into @xmath2 through the diffeomorphism @xmath232 is still denoted by @xmath216 .
let @xmath235 be such that @xmath236 and define @xmath237 . arguing by contradiction , let s assume that @xmath238 . in local coordinates around @xmath239 we consider the surface @xmath240 on the euclidean ball @xmath241 endowed with the rescaled metric @xmath242 .
therefore , @xmath243 converges to @xmath244 as @xmath154 .
the surface @xmath245 now has the property that @xmath246 .    by the strong compactness for a sequence of isoperimetric surfaces with bounded second fundamental form ,
see corollary [ compactness iso ] in the appendix , there exists a subsequence converging to a properly embedded surface @xmath247 , the convergence is in the sense of graphs and with multiplicity one .
moreover , @xmath248 is also stable , i.e. : @xmath249 if @xmath248 is compact , then it has to be a round sphere by alexandrov s theorem , which is a contradiction since strong convergence preserves topology .
if @xmath248 is non - compact , then it has infinite area by the monotonicity formula : indeed , by proposition [ monotonicity ] in the appendix there exists a positive constant @xmath250 such that @xmath251 in particular , @xmath252 . as @xmath248 is properly embedded , it has infinite extrinsic diameter and the claim follows .
therefore , @xmath248 is totally geodesic by da silveira s theorem [ silveira ] in the appendix , which is a contradiction since @xmath253 .
below we denote @xmath254 by @xmath90 .
let @xmath255 be the clifford torus through @xmath228 enclosing a region @xmath256 . under the scaling by @xmath257 we have that @xmath258 .
if @xmath259 , then @xmath260 . moreover , the second fundamental form @xmath261 of @xmath255 satisfies @xmath262 .
the critical fiber @xmath263 is distant from @xmath228 by @xmath264 since @xmath256 is converging to a compact region in @xmath2 . instead of using base points @xmath228
we choose new base points @xmath265 ; it follows that @xmath266 and @xmath267 .
we claim that rank of @xmath58 is at most one : @xmath268 let s assume now that @xmath269 , consequently @xmath270 and @xmath271 .
recall the function @xmath272 , the distance from the clifford torus through @xmath149 to the critical fiber @xmath273 with respect to the round metric .
the unit vector field @xmath274 is orthogonal to @xmath90 for every @xmath275 and it is well defined on @xmath276 .
let @xmath277 be the geodesic whose velocity is @xmath274 and such that @xmath278 .
consider @xmath279 . by the triangle inequality @xmath280 under the metric @xmath281 .
applying the coarea formula for @xmath282 , @xmath283 , we obtain : @xmath284 where @xmath285 $ ] is from the mean value theorem for integrals .
the last inequality is justified as follows .
either the extrinsic diameter of @xmath286 is going to infinity and @xmath286 is converging with multiplicity to a flat surface or the extrinsic diameter of @xmath286 is bounded .
the former implies that @xmath287 .
the latter implies that @xmath288 , which is a contradiction since @xmath289 .
we have conclude that @xmath290 and rank of @xmath2 is at most one .
[ isoperimetric small volume ] _ for @xmath213 large enough there exist @xmath292 and @xmath293 such that @xmath291 is given by the profile of spheres on @xmath294 $ ] and by the profile of flat tori on @xmath295 .
moreover , if @xmath216 is an isoperimetric surface such that @xmath296 , then @xmath297 .
_    for each @xmath213
we consider the first volume , @xmath292 , for which there is transition on topology of isoperimetric surfaces from spheres to something else . if @xmath298 is the volume for which the profile of geodesic spheres intersect the profile of flat tori , then @xmath299 . the value of @xmath298 is computed by solving the following system of equations : @xmath300 the left hand sides ( right hand sides ) correspond to the enclosed volume and area of the clifford torus @xmath90 ( geodesic spheres @xmath301 of radius @xmath302 ) , respectively .
it follows that @xmath303 ; another way to see this is by recalling that the injectivity radius of @xmath38 is @xmath304 at every point
. therefore , @xmath305 .
let @xmath216 be an isoperimetric surface with the property that @xmath296 . by lemma
[ bound second fundmental form lemma ] the sequence @xmath306 has bounded second fundamental form in @xmath227 ; thus , it strongly converges to a properly embedded surface @xmath59 of finite area in some orientable flat three manifold @xmath231 of rank at most one by lemma [ rank one ] . by the monotonicity formula , proposition [ monotonicity ] in the appendix , @xmath59 is a closed surface .
it follows that the pre - image @xmath307 of @xmath59 in @xmath308 is contained in a solid cylinder .
hence , @xmath307 is an union of round spheres by alexandrov s theorem or is a surface of revolution about the axis of the cylinder by theorem [ kks ] in the appendix .
therefore , @xmath94 or @xmath80 . from the strong compactness for isoperimetric surfaces ,
see corollary [ compactness iso ] , we have that @xmath297 , @xmath309 and the existence of the desired @xmath293 .      by the strong compactness for isoperimetric surfaces , corollary [ compactness iso ] , there exists @xmath310 such that if @xmath59 is an isoperimetric surface enclosing volume @xmath311 $ ] , then @xmath59 is a flat torus .
it follows from lemma [ isoperimetric small volume ] and for large @xmath213 that the isoperimetric profile @xmath291 is given by the area of geodesic spheres for volumes in @xmath294 $ ] and by the area of flat tori for volumes in @xmath312\cup [ \widehat{v},\frac{\pi^2}{p}]$ ] . in other words , if @xmath313 is the function defined by @xmath314 , where @xmath315 is the clifford torus enclosing a volume equal to @xmath316 , then @xmath317 on @xmath312\cup [ \widehat{v},\frac{\pi^2}{p}]$ ] .
it follows that @xmath318 has a local maximum point at @xmath319 $ ] .
the claim will follow by exploring the weak differential equation for @xmath291 . from ( [ second derivative area ] )
we have @xmath320 let @xmath59 be an isoperimetric surface such that @xmath321 . if @xmath117 is the unit normal variation of @xmath59 parametrized by the enclosed volume , then we define @xmath322 . since @xmath323 , we have that @xmath324 has also a local maximum point at @xmath325 .
hence , @xmath326 , i.e. , @xmath327 , and @xmath328 .
applying equation ( [ second derivative area ] ) for @xmath138 together with the gauss equation and the gauss - bonnet theorem we obtain : @xmath329 therefore , @xmath95 and @xmath330 in @xmath331 $ ] .
we claim that the injectivity radius of @xmath22 at every point is @xmath264 . indeed , let @xmath207 be the orbit of @xmath334 with respect to @xmath189 . as before , @xmath335 points of @xmath207 lie on the curve @xmath336 . here
, @xmath337 is the one parameter family of diffeomorphisms associated to @xmath338 : @xmath339 when ordered according the orientation of @xmath340 , those points determine a piecewise closed geodesic @xmath178 with @xmath341 .
as @xmath154 , @xmath178 converges to @xmath340 .
the claim now follows from : @xmath342    by the cheeger - gromov compactness theorem we obtain , up to subsequence , the following convergence : @xmath343 where @xmath344 is a flat three manifold of rank one by lemma [ rank one ] . by ( [ injective radius ] )
we conclude that @xmath345 .
the curves @xmath346 represent integral curves of @xmath338 through @xmath148 ; they have bounded geodesic curvature and @xmath347 .
hence , the integral curves of @xmath338 converge to closed geodesics in @xmath2 under ( [ convergence ] ) . as sets ,
they coincide with the standard fibers of @xmath348 .    by the poincare - hopf index theorem there
exists a zero for the vector field @xmath349 since @xmath350 .
hence , we can choose the base points @xmath228 to satisfy @xmath351 , here @xmath58 is the unit normal vector of @xmath216 .      by lemma [ bound second fundmental form lemma ] the sequence of isoperimetric surfaces
@xmath356 has uniformly bounded second fundamental form . applying the strong compactness theorem for isoperimetric surfaces , corollary [ compactness iso ] ,
we conclude that @xmath216 converges smoothly and with multiplicity one to a properly embedded stable cmc surface @xmath357 .
if @xmath358 , then the monotonicity formula , proposition [ monotonicity ] , implies that the extrinsic diameter of @xmath216 and @xmath354 are bounded .
this is impossible since the sequence @xmath216 separates @xmath22 in two regions of the same volume that goes to infinity as @xmath359 .
therefore , @xmath354 is a complete properly embedded stable cmc surface in @xmath355 with infinite area . applying da silveira s theorem [ silveira ]
once more , we obtain that @xmath354 is totally geodesic .
as @xmath351 , we conclude that @xmath354 is orthogonal to the standard fibers of @xmath355 .    we claim that @xmath354 separates @xmath355 .
if it does not separate , then there exists a loop @xmath360 intersecting @xmath354 at a single point .
as @xmath361 with multiplicity one , the same conclusion holds for @xmath216 , which contradicts the fact that @xmath216 separates @xmath22 .
therefore , there is @xmath362 such that @xmath363 , where @xmath364 is a flat plane for each @xmath365 .
we regard @xmath355 as a slab in @xmath308 with height @xmath160 .
now , we construct a deformation of @xmath354 which decreases its area as follows .
first , we cut off the @xmath366 solid cylinders obtained from the intersection of @xmath367 with a vertical solid torus of radius @xmath143 . to balance the enclosed volume
we add a vertical solid torus of radius @xmath275 , see figure 1 for the case @xmath368 .
if @xmath369 is the distance between @xmath370 and @xmath371 , then the radius @xmath275 is given by : @xmath372 the boundary of this new region is denoted by @xmath373 and @xmath374      if @xmath143 is large enough , then @xmath373 has less area than @xmath354 .
this is impossible since the strong multiplicity one convergence allow us to carry out this deformation of @xmath354 to @xmath216 which contradicts the fact that @xmath216 is an isoperimetric surface .
it is proved in @xcite that if @xmath49 has positive mean curvature , then @xmath377 applying this formula to the sequence of isoperimetric surfaces , we conclude that the mean curvature of @xmath216 satisfy @xmath378 .
we have used that @xmath216 separates @xmath22 into regions of the same volume .
as @xmath379 has area and mean curvature bounded , we apply _
allard s compactness theorem _ ,
theorem 42.7 and remark 42.8 in @xcite , to obtain an integral varifold @xmath380 that , up to subsequence , is @xmath381 .
recall that @xmath382 acts on @xmath21 via @xmath383 .
we claim that @xmath384 is @xmath382 invariant : indeed , if @xmath385 and @xmath386 , then for each @xmath213 there is @xmath387 such that @xmath388 . on the other hand , as @xmath389 , there is @xmath390 such that @xmath391 and @xmath392 .
it follows that @xmath393 and @xmath394 .
in particular , @xmath395 where @xmath396 is a clifford torus .
the monotonicity formula implies that the convergence @xmath397 is also in hausdorff distance ; hence , we have that @xmath398 and @xmath399 . since @xmath245 is @xmath332 invariant there exists @xmath400 for which @xmath401 . by taking the limit
we obtain @xmath402 for some @xmath403 .
as this is impossible , we conclude that @xmath368 and all @xmath245 are connected for @xmath213 large .
now we consider @xmath404 , the space of 2-dimensional integral _ currents _ on @xmath21 .
each @xmath405 and @xmath406 .
as @xmath407 is a region of finite perimeter ( @xmath408 is bv function with uniform bounded variation ) , then @xmath409 and @xmath410 as currents , @xmath411 is an open set of finite perimeter , see theorem 6.3 and proof of theorem 37.2 in @xcite . since @xmath412 we conclude that @xmath413 . applying the _ constancy theorem _ ,
theorem 26.27 in @xcite , we conclude that @xmath411 is the handlebody bounded by the clifford torus @xmath414 , and consequently @xmath415 .
we proved that @xmath416 for some positive integer @xmath417 .
since @xmath216 is isoperimetric , it follows that @xmath418 .
indeed , @xmath419 as @xmath197 is smooth , we have for @xmath420 sufficiently small that the density @xmath421 , where @xmath422 is from theorem [ allard regularity ] in the appendix . on the other hand , as @xmath245 is converging to @xmath197 with multiplicity one , then @xmath423 for @xmath213 large enough .
now we invoke the smooth version of _
allard s regularity theorem _ ,
theorem [ allard regularity ] , to concluded that the convergence @xmath424 is strong , i.e. , graphical with multiplicity one . as strong
convergence preserves topology , we conclude that @xmath425 .
this completes the proof of theorem [ theorem 1 ] .
the _ berger metrics _ are riemannian metrics @xmath430 on @xmath21 defined as : @xmath431 the riemannian manifolds @xmath432 are called the _ berger spheres _ , they are denoted by @xmath32 .
geometrically , the metric @xmath430 shrinks the hopf fibers to have length @xmath34 .
the berger metrics are also homogeneous and their group of isometries has dimension four .
it follows from the work of abresch and rosemberg @xcite that every constant mean curvature surface in @xmath32 admits a holomorphic quadratic differential .
in particular , every cmc sphere in @xmath32 is rotationally invariant .    a precise study of closed orientable surfaces with constant mean curvature on the berger spheres is given in @xcite .
it is proved there the existence of @xmath433 with the following property : if @xmath434 $ ] , then every stable constant mean curvature surface in @xmath32 has genus zero or one . moreover ,
if @xmath435 $ ] , then these stable cmc surfaces are totally umbilical spheres or the minimal clifford torus , the latter only occurring when @xmath436 .
in particular , rotationally invariant spheres are the only solutions of the isoperimetric problem in @xmath32 for @xmath435 $ ] .
we rescale the metric @xmath430 of @xmath32 by the factor @xmath441 .
the hopf fibers have constant length equal to @xmath160 under the new metric @xmath442 .
it follows that the injective radius of @xmath32 at a point @xmath213 is equal to @xmath443 for every @xmath444 .
since @xmath138 is a local trivial fibration , we have that for each @xmath445 there exist a neighbourhood @xmath446 of @xmath447 and a diffeomorphism @xmath448 such that @xmath449 , where @xmath450 given by @xmath451 .
moreover , @xmath452 in the @xmath226 topology .
therefore , in the sense of cheeger - gromov we have : @xmath453    we pick the points @xmath454 with the property that @xmath455 , this means @xmath427 and @xmath456 are parallel at @xmath457 .
these points exist by the poincar - hopf index theorem . by lemma [ bound second fundmental form lemma ]
the inclusion of @xmath439 in @xmath458 has the following property : @xmath459 by the strong compactness theorem for isoperimetric surfaces , corollary [ compactness iso ] in the appendix , we can extract a subsequence , @xmath460 , which converges with multiplicity one to a properly embedded surface @xmath461    if @xmath462 , then the monotonicity formula , proposition [ monotonicity ] , implies that @xmath354 is compact .
we apply theorem [ kks ] to conclude that @xmath354 is either a round sphere or torus .
this is impossible since we have strong convergence and @xmath463 .
therefore , @xmath354 is a complete non - compact surface with infinite area .
moreover , @xmath354 is also a stable cmc surface in @xmath355 : @xmath464 it follows from theorem [ silveira ] that @xmath354 is totally geodesic . by the choice of @xmath457
we conclude that @xmath354 is orthogonal to the @xmath465 fibers of @xmath355 .
since @xmath354 separates @xmath355 , we also conclude that @xmath354 is an union of at least two totally geodesic planes .
as shown in the proof of theorem 1.1 , this configuration can not be a limit of isoperimetric surfaces .
[ monotonicity ] _ let @xmath466 be a three manifold with bounded curvature , @xmath467 , and with positive lower bound on the injectivity radius @xmath468 . if @xmath469 is a smooth surface with mean curvature which satisfies @xmath470 , then there exists a positive constant @xmath471 such that @xmath472 for every @xmath473 and @xmath474 . _
[ allard regularity ] _ let @xmath466 be a closed manifold and @xmath476 .
there exist @xmath477 and @xmath478 with the following property : if @xmath479 is a smooth embedded cmc surface satisfying @xmath480 for every @xmath481 and @xmath482 , then its second fundamental form is uniformly bounded , i.e. , @xmath483 . _
let @xmath484 be a sequence of surfaces in a manifold @xmath2 .
we say that @xmath485 converge to @xmath59 in the sense of graphs if near any point @xmath486 and for large @xmath45 the surface @xmath485 is locally a graph over an open set of @xmath487 and these graphs converge smoothly to the graph of @xmath59 .
in addition , we say that @xmath488 satisfy _ local area bounds _ if there exist @xmath420 and @xmath169 such that @xmath489 for every @xmath490 .
[ compactness ] _ let @xmath491 be a sequence of embedded surfaces with constant mean curvature satisfying local area bounds and such that @xmath492 .
let s assume that @xmath493 converges to a metric @xmath494 in the @xmath226 topology . if @xmath495 has an accumulation point , then we can extract a subsequence that converges to a properly weakly embedded cmc surface @xmath59 in @xmath231 . _
let s first recall the constant mean curvature equation for graphs .
if @xmath248 is a surface with constant mean curvature @xmath60 in @xmath62 , then @xmath248 can be written locally as a graph over a neighbourhood @xmath496 : @xmath497 in coordinates @xmath498 where @xmath499 is the coordinate base associated to @xmath500 .
let @xmath501 be the coordinate base for @xmath485 , i.e. , @xmath502 .
the induced metric @xmath138 is expressed by @xmath503 .
a simple computation gives : @xmath504 we also have @xmath505 where @xmath506 are the christoffel symbols of @xmath106 .
therefore , @xmath507 finally , the mean curvature equation is written as : @xmath508 since @xmath509 and @xmath510 , then @xmath511 .
the equation ( [ mean curvature equation ] ) is uniformly elliptic as long as @xmath512 .
let @xmath513 be an accumulation point for the sequence @xmath488 .
by the upper bound on the second fundamental form , @xmath514 , there exists @xmath515 such that for every @xmath516 we have that @xmath517 is locally a graph @xmath518 over a neighbourhood @xmath519
. moreover , there exists @xmath520 for which @xmath521 . as @xmath522 ,
the schauder estimates for solutions of elliptic equations , see @xcite , imply that @xmath518 have @xmath523 estimates on @xmath524 , i.e. , @xmath525 .
therefore , @xmath526 are uniformly bounded and equicontinuous .    as @xmath484
satisfy local area bounds , then @xmath527 . on the other hand , the monotonicity formula , proposition [ monotonicity ] ,
gives that @xmath528 , where @xmath529 is a connected component of @xmath530 .
it follows that the number of components of @xmath531 is finite and independent of @xmath45 . by the arzel - ascoli theorem
we can extract a subsequence for which @xmath532 converges to @xmath533 for every @xmath534 .
moreover , @xmath533 also satisfies the constant mean curvature equation ( [ mean curvature equation ] ) .
as the set of accumulation points of @xmath488 is compact in @xmath535 we can cover this set by finite balls @xmath536 with @xmath537 . repeating the arguments in each of these balls and applying a diagonal argument
we obtain a properly immersed surface @xmath59 on @xmath538 with constant mean curvature @xmath60 .
since the surfaces @xmath485 are embedded , it follows that @xmath59 does not cross itself though it may have tangential self - intersections .
therefore , @xmath59 is properly weakly embedded in @xmath2 .
[ compactness iso ] _ let @xmath539 be a sequence of isoperimetric surfaces with @xmath540 . assume that @xmath541 converges , in the sense of cheeger - gromov , to a three manifold @xmath542 .
there exists a properly embedded surface @xmath543 such that @xmath544 in the sense of graphs and the convergence is with multiplicity one . _
first we remark that @xmath488 satisfy local area bounds .
indeed , @xmath545 . by proposition [ compactness ]
we only need to rule out possibly multiplicities for the convergence @xmath544 and points where @xmath59 fails to be embedded .
if the multiplicity of the limit is bigger than two , then @xmath546 has several components getting arbitrarily close .
this allow us to do a local cut and past deformation , as shown in figure 2 , that preserves the enclosed volume . if @xmath494 is the euclidean metric , then @xmath547 and @xmath548 .
thus , if @xmath549 , then @xmath550        the construction to deal with points where @xmath59 has tangential self intersections is similar to the multiplicity two case .
the corollary now follows since these constructions contradict the fact that @xmath485 is isoperimetric for every @xmath45 .
[ silveira ] _ let @xmath554 be a complete orientable surface conformally equivalent to a compact riemann surface punctured at a finite number of points .
let @xmath555 be an operator on @xmath59 satisfying @xmath556 and @xmath557 .
if @xmath59 has infinite area , then there exists a piecewise smooth function @xmath558 with compact support such that @xmath559 _
u. abresch and h. rosenberg .
_ generalized hopf differentials_. mat .
contemp . 28
( 2005 ) , 1 - 28 .
h. bray , _ the penrose inequality in general relativity and volume comparison theorems involving scalar curvature _ ( thesis ) .
arxiv:0902.3241v1 , ( 1998 ) .
d. gilbarg and n.s .
trudinger . _ elliptic partial differential equations of second order_. second edition .
fundamental principles of mathematical sciences , 224 .
springer - verlag , berlin , 1983 .
xiii+513 pp .
p. petersen .
_ riemannian geometry_. second edition . graduate texts in mathematics , 171 .
springer , new york , 2006 .
xvi+401 pp .
r. pedrosa and m. ritor .
_ isoperimetric domains in the riemannian product of a circle with a simply connected space form and applications to free boundary problems_. indiana univ .
j. 48(4 ) , ( 1999 ) , 1357 - 1394 .