there is increasing evidence that supermassive black holes ( smbhs ) exist at the center of all galaxies , and that they are the sources which power active galactic nuclei and quasars @xcite .
for example , vlbi observations of the keplerian disk around an object in ngc4258 indicate that the central object has a mass @xmath7 .
also , large numbers of observations are provided by the hubble space telescope suggesting that smbhs exist in galaxies such as m31 ( @xmath8 ) , m87 ( @xmath9 ) and our own galaxy ( @xmath10 ) ( see for example , * ? ? ? * for a brief overview ) . although evidence of the existence of smbhs is compelling , the actual formation process of these objects is still uncertain @xcite .
several different scenarios have been proposed , some based on stellar dynamics , others on gas hydrodynamics , and still others which combine the processes . at present , there is no definitive observation as yet which confirms or rules out any one of these scenarios .    here
we discuss the collapse of a supermassive star ( sms ) as one scenario of formation of an smbh .
this subject is also interesting from a viewpoint of general relativity , since the onset of radial instability at the mass shedding limit takes place at @xmath11 ( @xmath1 is the angular momentum , @xmath2 is the total gravitational energy ) in the uniformly rotating sms @xcite , which is also considered as a critical threshold between formation of a stationary black hole ( bh ) and of a naked singularity .
the bh uniqueness theorem says that a complete gravitational collapse of a body always results in a bh rather than a naked singularity , and that the final state of the bh should go into a stationary kerr bh ( e.g. , * ? ? ?
however there is a candidate to except this cosmic censor .
the collapse of a prolate spheroid with large semimajor axis leads to spindle singularities without an apparent horizon @xcite .
therefore , a collapse of a star with the critical value @xmath12 may show us an interesting phenomenon in general relativity .
what happens when the star of @xmath13 collapses ?
what is the final fate of a collapsing star of @xmath13 ?
numerical simulations can clearly show the answers to these questions .
the gravitational collapse of a star of @xmath13 has been investigated most of all in axisymmetric spacetime .
the pioneering study in this field was made by @xcite .
he set up a differentially rotating star , adding radial velocity to induce the collapse .
he found that the criterion of formation of a bh apparent horizon is @xmath14 .
the following study has been done by @xcite .
they set up a uniformly rotating @xmath15 polytropic star and deplete 99 % of the whole pressure to induce the collapse .
they found that the criterion of bh formation is @xmath16 , and when the star exceeds the above value , flattened disk formed ( @xmath17 for the case of 99 % pressure deplete ) .
note that @xmath18 represents the critical kerr parameter that the collapse proceeds to a bh .
finally @xcite performed a collapse of a differentially rotating @xmath15 polytropic star .
he also depleted the pressure of the star to induce the collapse .
he found that when @xmath19 is less than 0.5 , a bh is formed when the rest mass is larger than the maximum mass of the pressure depleted @xmath1 constant sequence .
when @xmath19 is slightly less than @xmath20 , a bh formed when the rest mass is sufficiently larger than the maximum allowed mass on the pressure depleted @xmath1 constant sequence .
he also found , by comparing his results from polytropic evolution and from @xmath21-law evolution of the hydrodynamics , that shock heating prevents the prompt collapse to a bh under the condition of @xmath13 .    at the same time
we have submitted our paper , several groups have also investigated the collapse of a relativistic star of @xmath13 in full general relativity .
@xcite investigated the onset of radial instability of uniformly rotating stars in 2d and find that the criterion of @xmath19 is @xmath22 in a very soft polytropic equation of state ( @xmath23  @xmath24 ; @xmath25 is a polytropic index ) .
@xcite investigated the collapse of a differentially rotating @xmath15 polytropic star in 3d by depleting the pressure and find that the criterion of bh formation is @xmath26 , and when the star exceeds the above value , the collapsing star forms a torus which fragments into nonaxisymmetric clumps .
@xcite studied the collapsing star of @xmath13 in 2d with a polytropic index of @xmath15 and @xmath27 and find that the criterion for bh formation is determined by the central @xmath28 ( @xmath29 ; @xmath30 is the specific angular momentum and @xmath31 is the cylindrical mass ) at the initial state of the star .
the purpose of this paper is the following threefold .
the first is to verify the nature of the gravitational collapse of a differentially rotating radially unstable star from the viewpoint of cosmic censor .
if we collapse a star of @xmath32 , we expect that the star can not directly form a bh of @xmath32 because of cosmic censor .
therefore it is important to find the main cause to prevent bh formation of @xmath33 . from the previous computational results , shock heating @xcite and core bounce @xcite are the dominant phenomena to prevent a star from forming a bh .
however , all of the previous calculations have set up violent initial data sets , adding radial velocity or depleting pressure , to induce the collapse , it may cause abrupt transport of the energy .
we therefore set up a mild , natural situation , that is a collapse of differentially rotating stars from the onset of radial instability , to focus on a graduate transport of the energy . also , we compute the collapsing star in 3d to allow shock propagation or bar formation , if it occurs .
the second is to determine the final outcome of the collapse of differentially rotating smss .
two different types of rotation profile arises during the quasi - static evolution of the rotating star , that depends on the environment .
one is a uniformly rotating star , maintained uniform rotation during the quasi - static evolution by sufficiently strong viscosity or strong magnetic field .
the other is a differentially rotating star , impossible to retain uniform rotation in low viscosity and low magnetic field (
even the star initially takes uniform rotation ) . for the collapse of a uniformly rotating sms from the onset of radial instability
, @xcite studied 3d relativistic hydrodynamic simulation in the post - newtonian gravitation and found that the collapse is coherent and that it is likely to form an smbh with no significant bar nor disk formation .
followup computation has been performed in 2d hydrodynamics in full general relativity and found that approximately 10 % of the total mass can form a disk , while approximately 90 % of that should form a bh @xcite . for the differentially rotating stars ,
the final outcome may strongly depends on the amount and the distribution of initial angular momentum , and the nature of angular momentum distribution of smss is largely unknown .
we also do not know what path does a differentially rotating sms take during the quasi - static evolution .
when the amount of angular momentum of the star is sufficient , the star seems to follow in a quasi - static evolution up to the mass - shedding limit .
@xcite find in the newtonian quasi - static evolution that bar formation is inevitable before reaching mass - shedding limit .
when the amount of angular momentum of the star is not sufficient , the star seems to collapse at the onset of radial instability before reaching secular instability ( @xmath34 for a uniformly rotating incompressible newtonian stars @xcite ) or mass - shedding limit .
the purpose of this paper is to examine the collapse of the sms from the onset of radial instability .
when we consider the collapse of radially unstable differentially rotating stars , the final outcome has a possibility of differing from the collapse of uniform rotation because of the strong centrifugal force at the central core , which may prevent the prompt collapse .
what is the final fate of the collapse ? does the star fragment due to the growth of the degree of differential rotation ? does the disk form during the collapse ?
relativistic simulation can clearly answer to these questions .    finally , it is important to probe whether the collapse of differentially rotating smss could be a promising source of gravitational waves .
direct detection of gravitational waves by ground based and space based interferometers is of great importance in general relativity , in astrophysics , and in cosmology ( e.g. , * ? ? ?
the catastrophic collapse is one of the promising sources of gravitational waves ( e.g. , * ? ? ?
for a gravitational collapse of the star in a dynamical timescale , there are two main reasons that prevent a prompt collapse of the star , which should produce gravitational waves at that moment .
one reason is core bounce and/or shock heating .
suppose gravitational force is balanced to the centrifugal force ( @xmath35 ) in newtonian gravity .
the total mass and the angular momentum ( @xmath36 ) are almost conserved during the collapse , assuming that gravitational radiation takes a little role during the collapse in the absence of nonaxisymmetric deformation prior to form a bh .
note that @xmath37 is angular velocity , @xmath4 is the radius .
we may estimate the bounce radius of the star in a global sense as @xmath38 .
since the horizon radius is roughly the order of @xmath2 , core bounce might take place for the case @xmath32 , which prevents violation of cosmic censor @xcite .
the other is bar formation . from the dimensional analysis
, we can describe the rotational kinetic energy @xmath39 and the gravitational binding energy @xmath40 as @xmath41 and @xmath42 .
we also accept the assumption that total mass and angular momentum are almost conserved during the collapse .
we could estimate the radius of bar formation in terms of the ratio of the rotational kinetic energy to the gravitational binding energy as @xmath43 . since the dynamical instability for a uniformly rotating , incompressible newtonian star sets in at @xmath44 @xcite and for relativistic gravitation as @xmath45 @xcite , bar formation takes place at the radius @xmath46 for the collapse of @xmath13 .
therefore , we may expect that a gravitational collapse of @xmath13 is a promising source of quasi - periodic gravitational waves .
this paper is organized as follows . in   [ sec : cf ] we present basic equations of our conformally flat approximation in general relativity .
we demonstrate our code tests in   [ sec : ctest ] .
we discuss our initial data sets and our numerical results for our catastrophic collapse in
[ sec : jm2 ] . in   [ sec : discussion ] we summarize our findings . throughout this paper
, we use geometrized units ( @xmath47 ) and adopt cartesian coordinates @xmath48 with the coordinate time @xmath49 .
greek and latin indices run over @xmath50 and @xmath51 , respectively .
in this section , we briefly describe the conformally flat spacetime ( e.g. , * ? ? ? * ) .
we solve the fully relativistic equations of hydrodynamics , but neglect nondiagonal spatial metric components @xcite .
we define the spatial projection tensor @xmath52 , where @xmath53 is the spacetime metric , @xmath54 the unit normal to a spatial hypersurface , and where @xmath55 and @xmath56 are the lapse and shift . within a first post - newtonian approximation , the spatial metric @xmath57
may always be chosen to be conformally flat @xmath58 where @xmath59 is the conformal factor ( see * ? ? ?
* ; * ? ? ?
the spacetime line element then reduces to @xmath60 we adopt maximal slicing , for which the trace of the extrinsic curvature @xmath61 vanishes , @xmath62 the gravitational field equations in conformally flat spacetime for the five unknowns @xmath55 , @xmath56 , and @xmath59 can then be derived conveniently from the 3 + 1 formalism ( e.g. , * ? ? ?
* ) .    since the spatial metric is conformally flat , the transverse part of its time derivative vanishes .
the transverse part of the evolution equation of the spatial metric therefore relates the extrinsic curvature to the shift vector , @xmath63 inserting eq .
( [ eqn : evolutionmetric ] ) into the momentum constraint equation , then , yields an equation for the shift @xmath64 @xmath65 where @xmath66 is the flat space laplacian and @xmath67 is the momentum density . in the definition of @xmath68
, @xmath69 is the stress energy tensor .
the conformal factor @xmath59 is determined from the hamiltonian constraint @xmath70 where @xmath71 is the mass - energy density measured by a normal observer .
maximal slicing implies @xmath72 , so that the trace of the evolution equation for the intrinsic curvature yields an equation for the lapse @xmath55 @xmath73 combined with eq .
( [ eqn : hamiltonianc ] ) , where @xmath74 .
therefore , conformally flat gravitational field equations for the five unknowns @xmath59 , @xmath75 , @xmath64 can be derived by eqs .
( [ eqn : momentumc ] ) - ( [ eqn : evolution ] ) .      for a perfect fluid ,
the energy momentum tensor takes the form @xmath76 where @xmath77 is the rest - mass density , @xmath78 the specific internal energy , @xmath79 the pressure , and @xmath80 the four - velocity .
we adopt a @xmath21-law equation of state in the form @xmath81 where @xmath21 is the adiabatic index which we set @xmath82 in this paper .    in the absence of thermal dissipation ,
( [ gammalaw1 ] ) , together with the first law of thermodynamics , implies a polytropic equation of state @xmath83 where @xmath84 is the polytropic index and @xmath85 is a constant .    from @xmath86 together with the equation of state ( eq . [ [ gammalaw1 ] ] ) , we can derive the energy and euler equations according to @xmath87 where @xmath88 and @xmath89 , @xmath90 is the 3-velocity , pressure viscosity , respectively .
note that we treat the matter fully relativistically ; the conformally flat approximation only enters through simplifications in the coupling to the gravitational fields . note also that we include an artificial viscosity , since core bounce might occur in our gravitational collapse and that should produce shocks .
we will explain our form of an artificial viscosity in
[ sec : ctest ] . as a consequence to treat shocks we also need to solve the continuity equation @xmath91 separately .
we solve the matter evolution in second order accurate in space and time with the transfer scheme of @xcite .
there are three key issues to solve gravitational field equations numerically .
the first is to introduce symmetric tensor @xmath92 related to the extrinsic curvature @xmath61 as @xcite @xmath93 therefore , we describe the equations of conformal factor ( eq . [ [ eqn : hamiltonianc ] ] ) and lapse ( eq . [ [ eqn : evolution ] ] ) as @xmath94 note that the dependence value ( such as @xmath59 , @xmath75 ) of the source term is always lower than o(@xmath95 ) , o(@xmath96 ) , which is safe for a convergence during the iteration . , the first term of eq .
( [ eqn : cfactor ] ) behaves as @xmath97 . ]
the second is to decompose the vector type elliptic equation .
the equation for the shift , @xmath98 can be further simplified by introducing a vector @xmath99 and a scalar @xmath100 according to @xcite @xmath101 the shift can then be computed from @xmath102 and will automatically satisfy eq .
( [ eqn : shifteq ] ) .
next , we decompose the symmetric tensor @xmath92 as ( e.g. , * ? ? ?
* ) @xmath103 where @xmath104 therefore , we can rewrite the momentum constraint as @xmath105 we also decompose the vector @xmath106 with the same manner as the shift ; @xmath107 where @xmath108 and @xmath109 satisfies @xmath110 note that @xmath68 only appears in the presence of matter , which means that the source terms of eqs . ( [ eqn : wb ] ) and ( [ eqn : wchi ] ) are compact .
therefore , we can compute these values quite accurately .    to summarize
, we have reduced einstein equations in a conformally flat spacetime to 10 elliptic equations for 10 variables ( @xmath108 , @xmath109 , @xmath59 , @xmath75 , @xmath99 , @xmath100 ) , @xmath111 these poisson - type equations are solved by imposing the following boundary condition at outer boundaries @xmath112 here we briefly explain our order of computing 10 poisson - type equations .
first , we solve @xmath108 and @xmath109 since the source terms do not depend on other unknown spacetime variables .
once we derive the symmetric tensor @xmath113 from @xmath108 and @xmath109 , we solve @xmath59 iteratively till the convergence .
after that we solve @xmath75 iteratively till the convergence .
finally , we solve @xmath99 and @xmath109 .
we use modified incomplete cholesky decomposition conjugate gradient ( miccg ) method @xcite to solve elliptic equations .
finally , we should maintain our grid resolution at the central core of the collapsing star .
the nested grid is one of the most appropriate way to handle this problem ( see , e.g. , * ? ? ?
our method is a mimic version of the nested grid to handle the collapsing star problem , which is based on the rigridding method of @xcite . for the post
homologous collapse of the sms , there are two different timescales in the collapsing star ; one is the collapsing timescale at the center and the other is the one at the envelope , which is sufficiently longer than that at the center .
therefore , @xcite add the grid number at the intermediate stage of the collapse in order to control both regions , center and envelope of the star .
since we perform the collapse of an sms in 3d , we concentrate on the central core of the collapse due to the limitation of computational resource . in order to maintain spatial grid resolution especially at the center ,
we shrink the total grid to one half as the core radius of the star becomes half from the beginning of the present grid resolution .
although this method changes the boundary condition when we `` zoom in '' the computational domain , it is approximately good when the difference of the gravitational mass between the end of the previous computational grid and the beginning is small .
in fact , the difference of the gravitational mass between every two computational domains is less than @xmath114 , and therefore this `` zoom - in '' method should behave as a reasonable approximation for a collapse of differentially rotating smss .
we monitor the gravitational mass @xmath2 and the angular momentum @xmath1 @xmath115 \psi^{5 }   \right .
\nonumber \\ & & \left .
+ \frac{1}{16 \pi } \psi^{5 } k_{ij } k^{ij } \right ] d^{3 } x   , \\ j & = & - \frac{1}{2 \pi } \oint_{\infty } ( x k^{i}_{y } - y k^{i}_{x } ) \psi^{6 } ds_{i } \nonumber \\ & = & \int ( x j_{y } - y j_{x } ) \psi^{6 } d^{3 } x , \end{aligned}\ ] ] during the evolution .
in all cases reported in   [ sec : jm2 ] , the gravitational mass was conserved up to @xmath116% , and the angular momentum up to @xmath117% of their initial values .
we also compute proper mass @xmath118 , rotational kinetic energy @xmath39 , gravitational binding energy @xmath40 in equilibrium as @xmath119 where @xmath37 the angular velocity of the star .
since we use a polytropic equation of state at @xmath120 , it is convenient to rescale all quantities with respect to @xmath85 @xcite . since @xmath121 has dimensions of length , we introduce the following nondimensional variables @xmath122 henceforth , we adopt nondimensional quantities , but omit the bars for convenience ( equivalently , we set @xmath123 ) .
first , we demonstrate 1d relativistic wall shock problem to check whether our code has an ability to treat shock .
we set up the form of artificial pressure viscosity @xcite as @xmath124 where @xmath125 , @xmath126 is the local grid spacing and where @xmath127 is the dimensionless parameter . when evolving the matter equations we limit the stepsize @xmath128 by the following courant condition ( @xmath129 ) , where the second term represents dynamical time of a collapsing star .
we choose the dimensionless parameter @xmath130 .
we have tested the ability of our code to resolve shocks by performing a wall - shock problem , in which two phases of a fluid collide . in fig .
[ fig : wshock ] we compare numerical results with the analytic solution for initial velocities that are similar to those found in our simulations below .
with @xmath131 we find good agreement , and set this value to simulate the catastrophic collapse in  [ sec : jm2 ] .
next , we demonstrate a spherical dust collapse in fig [ fig : oscollapse ] . note that our conformally flat approximation retains all the nonlinear terms to maintain the exact dynamics for a spherical spacetime .
we compare 1d and 3d results to check whether our 3d code has an ability to reproduce a spherical dust collapse .
note that we construct our 1d code following wilson s approach @xcite described in appendix [ sec:1dfield ] .
we choose the grid size of 1d computation as 5,000 , while of 3d one as ( @xmath132 ) .
we find a very good agreement of the central lapse between our 1d result and our 3d conformally flat one within the error of 1% .
we terminate the integration of our 3d code at @xmath133 , since the convergence of the iteration process in our elliptic solver becomes significantly slow .
then , we demonstrate the collapse of a spherical star in fig .
[ fig : spalp ] ( see table [ tbl : sphinitial ] for an equilibrium profile ) . since
a radially unstable spherical star only collapses to form a spherical bh promptly , we should follow the collapse of a spherical star to check whether our code has an ability to follow the collapse and to find a signal of bh formation .
we set up the same grid size as we did for a spherical dust collapse , covering the star with 81 grid points across the diameter .
the central lapse of the star decreases monotonically , which means a collapse of a spherical star directly forms a spherical bh .
although we terminate our integration around @xmath134 , the fact that 1d computation of the spherical star collapse finds an apparent horizon at @xmath135 indicates that a bh might form .
c c c c @xmath136 & @xmath137 & @xmath138 & @xmath139    finally , we check the stability of a uniformly rotating star whether our code has an ability to determine the radial stability of a star . since we determine the radial stability of a differentially rotating star in
[ sec : jm2 ] , we should check the sensitivity of our code to determine the critical onset of radial instability of uniformly rotating stars by comparing the results derived from the turning point method .    to assess the ability of our code to distinguish stable stars from unstable ones with rotation
, we consider an equilibrium sequence of uniformly rotating stars of the fixed angular momentum @xmath1 ( @xmath140 at the turning point ) .
while the turning point criterion strictly identifies the onset of secular instability , the point of onset of dynamical instability nearly coincides with the secular instability point .
we adopt the polytropic index of @xmath141 , which is regarded as radiation pressure dominant , sms sequence , and grid resolution as in the spherical simulations reported above .
we decrease the initial pressure to induce the collapse ( @xmath142 ) .
figure [ fig : roteq ] summarizes our dynamical stability analysis for the rotating sms .
we conclude that with the adopted grid resolution , our code can distinguish stable rotating stars from unstable ones within 0.3% of the maximum gravitational mass .
figure [ fig : rotst ] shows the evolution of the central density for stable and unstable rotating stars .
here we explain our initial data sets for collapsing stars .
since we are interested in a collapsing star of @xmath13 , we have three requirements to construct differentially rotating equilibrium stars .
the first requirement is to construct radially unstable differentially rotating stars .
the critical @xmath21 to the onset of radial instability for slow rotation and weak gravitational field is described analytically as @xcite @xmath143 where @xmath144 is the inertia momenta of the star .
note that relativistic gravitation unstabilizes the star , while rotation , which produces centrifugal force , stabilizes the star .
although we take into account the differential rotation in the equilibrium star , this criterion ( eq . [ [ eqn : critg ] ] ) may indicate the appropriate direction to choose a parameter sets of the initial condition of the rotating collapsing star . from eq .
( [ eqn : critg ] ) , we should at least choose soft equation of state to induce the collapse and set an @xmath145 polytropic star , an sms sequence of the star .
the second requirement is to construct a star which holds @xmath26 . from the critical onset of radial instability in the equilibrium star ,
a uniformly rotating sms takes the maximum of @xmath146 when @xmath147 @xcite . to verify the nature of cosmic censor
, we should at least go beyond @xmath32 .
the main restriction to hold large @xmath148 comes from the mass shedding limit of the star .
therefore , to construct a star of @xmath149 in radially unstable branch , differential rotation of the star is required in sms sequence .
the last requirement is high degree of differential rotation .
it is indeed one path to construct high degree of differential rotation in quasi - static evolution , and reach the onset of radial instability .    to summarize , we need soft equation of state and high degree of differential rotation , and choose an @xmath145 polytropic star ( sms sequences ) and the degree of differential rotation as @xmath150 where we define the rotation profile as @xmath151 in the newtonian limit ( @xmath152 , @xmath153 ) , this rotation law can be written as @xmath154 where @xmath155 is the degree of differential rotation
, @xmath156 is the cylindrical radius of the star .
since @xmath155 has a dimension of length , we normalize it with a proper equatorial radius @xmath157 , ( @xmath158 ) .
hereafter we choose @xmath159 to construct relatively a high degree of differential rotation .
we briefly summarize our method to construct relativistic rotating equilibrium stars in appendix [ sec : rst ] .
we summarize the parameters of our differentially rotating stars at initial in table [ tbl : initial ] .
we slightly perturb our initial equilibrium state according to @xmath160 where @xmath161 .
we install @xmath162 and @xmath163 density perturbation to provide the seed for one - armed spiral and bar formation , if the physical situation should lead to unstable growth .
we adopt a grid size ( @xmath164 ) , so that the star is initially covered by 161 points across the equatorial diameter .
we evolve the rotating sms up to the point at which the conformally flat approximation breaks down .
c c c c c & @xmath165 & @xmath166 & @xmath167 & @xmath168 + @xmath169 & @xmath170 & @xmath170 & @xmath170 & @xmath170 + @xmath171 & @xmath172 & @xmath173 & @xmath174 & @xmath175 + @xmath176 & @xmath177 & @xmath178 & @xmath179 & @xmath180 + @xmath181 & @xmath182 & @xmath183 & @xmath184 & @xmath185 + @xmath2 & @xmath186 & @xmath187 & @xmath188 & @xmath189 + @xmath1 & @xmath190 & @xmath191 & @xmath192 & @xmath193 + @xmath194 & @xmath195 & @xmath196 & @xmath197 & @xmath198 + @xmath19 & @xmath199 & @xmath200 & @xmath201 & @xmath202 + @xmath203 & @xmath204 & @xmath205 & @xmath205 & @xmath206 + stability & unstable & unstable & stable & stable    figure [ fig : rhomax ] shows the results of radial stability in 4 differentially rotating stars . since we do not have a tool to determine radial stability in differentially rotating stars from their equilibrium states ( but for determining the criterion of secular stability in rotating stars , see @xcite ) , the evolution is necessary to determine its stability . the criterion to determine the radial stability is as follows .
when the central density of the star grows exponentially within a few dynamical time , we determine the star radially unstable . on the other hand , when the central density of the star oscillates around its equilibrium state , we determine the star radially stable . from this criterion ,
model i and ii have an exponential growth of the central density that means radially unstable , while model iii and iv have maximum at several dynamical times that means radially stable .
we also show the evolution of central lapse in fig .
[ fig : lapse ] .
the rapid decrease of @xmath207 below @xmath208 indicates that a bh is likely to form .
model ii shows that the central lapse monotonically decreases from @xmath209 to @xmath210 .
this figure shows that we can follow the collapse from the regime of newtonian gravity ( @xmath211 ) to that of relativistic gravity ( @xmath212 ) . on the other hand model iii
shows that the central lapse oscillates around its equilibrium state , and that also means the star is radially stable .
we show the final density snapshots of the stars in the equatorial plane ( fig .
[ fig : qeq ] ) and in the meridional plane ( fig .
[ fig : qmd ] ) . even in the final snapshots of radially unstable stars ,
the collapse is almost axisymmetric .
also , from the snapshots in the meridional plane , the material around the rotational axes collapses faster than the surrounding material due to the strong nonlinear gravitational field .
let us now focus on the collapsing star ( model ii ) to probe the final outcome .
we locally define the cylindrical rest mass @xmath31 , angular momentum @xmath30 , specific angular momentum @xmath213 , keplerian angular velocity @xmath214 as @xmath215 to indicate the transport of angular momentum and the distribution of @xmath216 .
we assume the axisymmetric collapse to investigate the final outcome , since the meridional density snapshots ( fig . [
fig : qmd ] ) behave almost axisymmetry .
figure [ fig : mass ] shows the concentration of mass density profile during the collapse . as the collapse goes on , the increase ratio of the cylindrical mass at the core of the star becomes high .
in fact , @xmath217% of the cylindrical mass is inside the radius of @xmath218 at @xmath219 .
we also show the mean radius , defined as @xmath220 , during the collapse in fig .
[ fig : rm ] .
note that @xmath221 is the rest mass .
we find that the mean radius monotonically decreases , which means that the collapse of the central core is coherent .
figure [ fig : omega ] shows the variation of angular velocity profile during the collapse . since the cylindrical rest mass density shows a coherent collapse , the degree of differential rotation significantly increases at the core during the collapse and approaches to the `` keplerian '' angular velocity defined in eq .
( [ eqn : keplarian ] ) .
we show the specific angular momentum distribution during the collapse in fig .
[ fig : jm ] . in a global sense
, @xmath19 is conserved so that a final bh should violate cosmic censor , if it forms .
however , we find that most of the matter collapses to form a bh , while not a few amount of angular momentum stays at the surface area of the star that prevents forming a bh of @xmath149 .
we show the distribution of @xmath216 in fig .
[ fig : jm2 ] during the collapse to verify cosmic censor .
we find that mass density collapses first in the central part and angular momentum remains in the surface of the star that prevents forming a bh of @xmath149 .
note that the distribution of @xmath222 is only an indicator to interpret the physical cause to prevent bh formation .
although we should define the total gravitational mass locally to discuss the local distribution of @xmath19 , there is no knowing how to define it .
therefore we use the local rest mass instead .
we investigate the collapse of differentially rotating smss by means of hydrodynamic simulations in conformally flat approximation in general relativity .
we start our collapse around the onset of radially instability at @xmath3 to the point where conformally flat approximation breaks down .
we find that cosmic censor even holds for a gravitational collapse of a radially unstable differentially rotating equilibrium sms of @xmath32 .
the main reason to prevent formation of a bh of @xmath32 is that quite a large amount of angular momentum stays at the surface , not core bounce nor shock propagation , bar formation in our model .
note that even a thin disk near the surface of the star can hold relatively a large amount of angular momentum if the radius is large .
the above conclusion is supported by @xcite in full general relativistic simulation that the criterion of bh formation is determined at the central @xmath28 ( @xmath223 ) .
in fact , the central @xmath28 of model ii is @xmath224 , which satisfies their criterion .
the collapse of a differentially rotating , radially unstable sms is coherent and likely leads to formation of an smbh .
this situation is quite similar to the collapse of a uniformly rotating sms ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , since @xmath225 seems still not sufficient to change the path of the collapse from the one of uniform rotation .
however , this final outcome may depend on the equation of state of the star .
@xcite treated the isothermal ( @xmath226 ) collapse of initially homogeneous , uniformly rotating , low entropy clouds via smooth particle hydrodynamics ( sph ) simulations .
they found considerable fragmentation into dense clumps , and disk formation containing @xmath227 of the mass .
they concluded that a seed bh will form at the center and that it likely will grow gradually by accretion in our model .
also @xcite treated the collapse of uniformly rotating polytropic star from the critical onset in the range of @xmath228 and found that the mass of the disk is less than @xmath229 of the gravitational mass .
we can not find any evidence of bar formation nor significant disk formation from the rotating collapse prior to bh formation .
the phenomenon of no bar formation also comes from the fact that mass density collapses first to form a bh . in such case , @xmath194 can not scale in @xmath230 due to the growth of the degree of differential rotation , and as a fact the star of @xmath194 can not reach the dynamical instability point of @xmath231 . since the @xmath162 dynamical instability takes place when the star has a toroidal structure and soft equation of state @xcite , we can not find any evidence of toroidal structure at the time we stop our integration .
the rotating sms collapse is a promising source of burst gravitational waves and of quasi - normal mode ringing waves .
the characteristic frequency for burst ( @xmath232 ) , quasi - normal mode ringing ( @xmath233 ) , and the wave amplitude for burst ( @xmath234 ) , quasi - normal mode ringing ( @xmath235 ) , are @xcite @xmath236 , \\
h_{\rm burst } & \sim &   1 \times 10^{-18 } \left ( \frac{m}{10^{6 } m_{\odot } } \right ) \left ( \frac{1 { \rm g pc}}{d } \right ) \left ( \frac{m}{r } \right ) , \\
f_{\rm qnm } & \sim & 2 \times 10^{-2 } \left ( \frac{10^{6 } m_{\odot}}{m } \right ) { \rm [ hz ] } ,
\\ h_{\rm qnm } & \sim & 6 \times 10^{-19 } \left ( \frac{\delta e_{\rm gw}/m}{10^{-4 } } \right)^{1/2 } \left ( \frac{2 \times 10^{-2 } { \rm [ hz]}}{f_{\rm qnm } } \right)^{1/2 } \nonumber \\ & & \times \left ( \frac{m}{10^{6 } m_{\odot } } \right)^{1/2 } \left ( \frac{1 { \rm gpc}}{d } \right ) , \end{aligned}\ ] ] where @xmath237 is the distance from the observer and @xmath238 is the total radiated energy .
we set @xmath239 , a characteristic mean radius during bh formation . since the main targets of lisa are gravitational radiation sources between @xmath240 and @xmath241 hz , it is possible that lisa can search for the burst and quasi - normal ringing waves accompanying rotating sms collapse and formation of an smbh .
we would like to thank an anonymous referee for his / her critical reading of our manuscript and constructive suggestions .
he also thanks hideki asada , silvano bonazzola , eric gourgoulhon , takashi nakamura , ken - ichi nakao , luciano rezzolla , stu shapiro , nick streigeolous , takahiro tanaka , koji ury , shin yoshida for discussion .
this work has been supported in part by mext grant - in - aid for young scientists ( no . 200200927 ) , astronomical data analysis center , national astronomical observatory of japan , and by mext grant - in - aid for 21coe program at the department of physics , kyoto university .
numerical computations were performed on the nec sx-5 machine in the yukawa institute for theoretical physics , kyoto university , on the vpp-800 machine in the academic center for computing and media studies , kyoto university , and on the vpp-5000 machine in the astronomical data analysis center , national astronomical observatory of japan .
here we briefly explain relativistic hydrodynamics in spherically symmetric spacetime @xcite that we use in    [ sec : cf ] .
we can describe the line element of spherically symmetric spacetime in the isotropic coordinate as @xmath242 where @xmath55 is the lapse , @xmath243 is the shift , @xmath59 is the conformal factor .
we use the maximal slicing condition @xmath244 , which derives the following relation to the components of an extrinsic curvature as @xmath245 as a consequence , we only need to consider @xmath246 to construct an extrinsic curvature .
the momentum constraint derives the equation of the extrinsic curvature as @xmath247 the integration form of eq .
( [ eqn:1dmom ] ) is @xmath248 the restriction of the spatial metric to be conformally flat requires the following equation @xmath249 which gives us an appropriate boundary condition for the extrinsic curvature as @xmath250 .
also the integration form of eq .
( [ eqn:1dbet ] ) is @xmath251 the hamiltonian constraint and the trace of the evolution equation guide @xmath252 =   2 \pi \alpha \psi^{5 } ( \rho_{\rm h } + 2 s )   + \frac{21}{16 } \alpha \psi^{5 } ( k^{r}_{r})^{2 } \label{eqn:1dev } , \end{aligned}\ ] ] with the boundary conditions @xmath253 at the grid edge and @xmath254 at the center ( @xmath255 ) .    therefore , we can determine 4 unknown variables ( @xmath246 , @xmath59 , @xmath256 , @xmath243 ) with 4 equations ( eqs .
[ [ eqn:1dintk ] ] , [ [ eqn:1dintb ] ] , [ [ eqn:1dham ] ] , [ [ eqn:1dev ] ] ) . instead of
@xmath246 from a computational point of view ( see  [ subsec : nt ] . ) . ]
matter equations can be written as @xmath257 we have not taken into account artificial viscosity , since shock does not seem to play an important role in a spherical dust collapse and the collapse of a spherical star . in   [ sec : ctest ] , we adopt our 1d relativistic hydrodynamic code to a spherical dust collapse by setting @xmath258 and @xmath259 @xmath260 equilibrium variables of them in typical polytropic index .
we briefly mention a method to compute an apparent horizon in 1d @xcite .
an apparent horizon is defined as an outermost marginally trapped surface whose future - directed outgoing null geodesics have zero expansion @xmath261 @xcite .
the outgoing null vector @xmath262 can be expressed as @xmath263 , where @xmath264 is unit normal of the spacelike hypersurface defined as @xmath265(@xmath266 , 0 , 0 , 0 ) , @xmath267 is the out - directed spacelike vector orthogonal to the surface .
the projection tensor @xmath268 and the extrinsic curvature @xmath61 are described as @xmath269 the expansion @xmath261 is then written as @xmath270 for spherically symmetric spacetime , we can choose the spacelike vector @xmath271 as @xmath272 , and then the radius @xmath273 of an apparent horizon satisfies the following equation , @xmath274 note that an event horizon always lies outside an apparent horizon , and both horizons coincide with each other when the spacetime is stationary .
also an event horizon always forms before the time that an apparent horizon does @xcite .
here we summarize our method to construct rotating relativistic equilibrium stars , which is based on @xcite ( see also * ? ? ? * for a historical review ) .
first we focus on how to solve the relativistic euler equation .
the equation can be described in axisymmetric spacetime as @xmath275 where @xmath276 is a specific enthalpy .
we also assume specific type of rotation law as @xmath277 where @xmath278 is the central angular velocity , in order to integrate eq .
( [ eqn : reuler ] ) . the bernoulli s equation is driven by integrating the relativistic euler equation ( [ eqn : reuler ] ) as @xmath279 where @xmath280 , \\ k & \equiv & \int \frac{du^{t}}{u^{t } } = \ln u^{t }   = - \frac{1}{2 } \ln   [ \alpha^{2 } - \psi^{4 } [ ( \beta^{x } - \omega y)^{2 } +   ( \beta^{y } + \omega x)^{2 } ] ] , \\ r & \equiv & \int u^{t } u_{\varphi } d\omega   = -\frac{1}{2 } a^{2 } ( \omega_{c } - \omega)^{2 } , \end{aligned}\ ] ] where @xmath281 in this section . ] .
note that we have already assumed polytropic equation of state @xmath282 where @xmath283 , @xmath25 is a polytropic index .
therefore , we can describe the matter distribution equation by using eq .
( [ eqn : bernoulli ] ) and the rotation law ( eq . [
[ eqn : rlaw ] ] ) as @xmath284}{\sqrt{u^{t } } } - 1 \right ] \nonumber \\ & = & \frac{1}{n+1 }   \left [     \frac{c \exp[a^{2}(\omega_{c } - \omega)^{2 } / 2 ] }         { \sqrt{\alpha^{2 } - \psi^{4 } [ ( \beta^{x } - \omega y)^{2 } +            ( \beta^{y } + \omega x)^{2 } ] } } - 1 \right ] \label{eqn : qdis } , \\
a^{2 } ( \omega_{c } - \omega ) & = &   u^{t } u_{\varphi } = \frac{\psi^{4 } [ x ( \beta^{y } + \omega x ) - y ( \beta^{x } - \omega y ) ] } { \alpha^{2 } -     \psi^{4 } [ ( \beta^{x } - \omega y)^{2 } + ( \beta^{y } + \omega x)^{2 } ] } .
\label{eqn : rotlaw}\end{aligned}\ ] ] note that we assume conformally flat spacetime to derive the final part of the equality in eqs .
( [ eqn : qdis ] ) and ( [ eqn : rotlaw ] ) .    from the computational point of view
, we introduce nondimensional quantities rescaled in the equatorial proper radius of the star ( @xmath285 ) as @xmath286 we also rescale the lapse and conformal factor using the nature of scale free in newtonian gravity as @xmath287 to determine the matter distribution and rotation profile of the next iteration step , we have to determine five unknown variables ( @xmath285 , @xmath288 , @xmath278 , @xmath289 , @xmath290 ) from five equations .
therefore , we evaluate eqs .
( [ eqn : qdis ] ) and ( [ eqn : rotlaw ] ) at three locations , the point of the maximum rest mass density , that of polar surface , and that of the equatorial surface .
note that eq .
( [ eqn : rotlaw ] ) becomes an identity at polar surface of the star , and thus we have five equations to be solved .
we use newton - rapson method @xcite to solve these equations .
once we determine 5 unknown variables , we solve eq .
( [ eqn : rotlaw ] ) using newton - rapson method to determine the rotation profile .
after that we can determine the matter distribution ( eq . [ [ eqn : qdis ] ] ) .
next we briefly summarize the gravitational field equations normalized by proper radius as follows .
@xmath291 we choose the same boundary conditions as we choose in our evolution code ( see   [ sec : cf ] ) .
note that we construct the star in 3d using octant symmetry to adopt our initial data to our 3d evolution code smoothly .    to summarize , we first solve the gravitational field equations ( eqs . [
[ eqn : apcfbx ] ]  [ [ eqn : apcfeta ] ] ) .
next , we determine 5 unknown quantities ( @xmath285 , @xmath288 , @xmath278 , @xmath289 , @xmath290 ) from 5 equations ( eqs . [ [ eqn : qdis ] ] and [ [ eqn : rotlaw ] ] ) using newton - rapson method .
finally we determine the matter distribution and rotation profile for the next iteration step .
we continue this iteration cycle till @xmath285 converges ; i.e. @xmath285 goes bellow the error of @xmath292 .
we also check the error rate of physical quantities such as @xmath2 and @xmath1 and find that they are in fact below the error rate of @xmath240 .
99 baumgarte , t. w. & shapiro , s. l. 1999 , , 526 , 941 blanchet , l. , damour , t. & schfer , g. 1989 , , 242 , 289 bonazzola , s. , frieben , j. , & gourgoulhon , e. 1998 , , 331 , 280 bowers , r. l. , & wilson , j. r. 1991 , numerical modeling in applied physics and astrophysics ( boston : jones and bartlett pub . ) , chap . 9 .
centrella , j. m. , new , k. c. b. , lowe , l. l. , & brown , j. d. 2001 , , 550 , l193 chandrasekhar , s. 1965 , , 142 , 1488 . 1969 , _ ellipsoidal figures of equilibrium _ , ( new york : yale univ .
press ) , 61 chandrasekhar , s. & lebovitz , n. r. 1965 , , 142 , 1488 cook , g. b. , shapiro , s. l. , & teukolsky , s. a. 1992 , , 398 , 203 . 1994 , , 422 , 227 duez , m. , shapiro , s. l. , & yo , h. -j .
2004 , , 69 , 104016 hawkings , s. w. , & ellis , g. f. r. 1973 , the large scale structure of space - time , chap .
9.2 hawley , j. f. , smarr , l. l. , & wilson , j. r. 1984 , , 55 , 211 isenberg , j. , & nester , j. 1980 , in general relativity and gravitation vol . 1 : one hundred years after the birth of albert einstein , ed .
a. held ( new york : plenum press ) , 23 komatsu , h. , eriguchi , y. , & hachisu , i. 1989 , , 237 , 355 kormendy , j. 2004 , in carnegie observatories centennial symposium . 1 : coevolution of black holes and galaxies , ed
. l. c. ho ( cambridge : cambridge univ . press ) , in press ( astro - ph/0306353 ) loeb , a. , & rasio , f. a. 1994 , , 432 , 52 leaver , e. w. 1985 , proc .
london a , 402 , 285 murata , k. , natori , r. & karaki , y. 1990 , large - scale numerical simulation ( tokyo : iwanami ) , chap . 5.3 ( in japanese ) nakamura , t. 1981 , prog .
phys . , 65 , 1876 nakamura , t. , oohara , k. , & kojima , y. 1987 , prog .
, 90 , part i , sec .
2 nakamura , t. , shapiro , s. l. , & teukolsky , s. a. 1988 , , 38 , 2972 new , k. c. b. 2003 , living reviews in relativity , 6 , 2 new , k. c. b. , & shapiro , s. l. 2001 , , 548 , 439 oohara , k. , & nakamura , t. 1989 , prog .
phys . , 82 , 535 press , w. h. , flannery , b. p. , teukolsky , s. a. , & vetterling , w. t. 1992 , numerical recipes in fortran ( 2nd ed . ; cambridge : cambridge univ .
press ) , chap .
9.6 rees , m. j. 1998 , in black holes and relativistic stars , ed .
r. m. wald ( chicago : univ .
chicago press ) , 79 . 2001 , in black holes in binaries and galactic nuclei , ed .
l. kaper , e. p.
j. van den heurel , & p. a. woudt ( new york : springer - verlag ) , 351 ruffert , m. 1992 , , 265 , 82 saijo , m. , shibata , m. , baumgarte , t. w. , & shapiro , s. l. 2001 , , 548 , 919 saijo , m. , baumgarte , t. w. , shapiro , s. l. , & shibata , m. 2002 , , 569 , 349 saijo , m. , baumgarte , t. w. , & shapiro , s. l. 2003 , , 595 , 352 sekiguchi , y. -i .
, & shibata , m. 2004 , , submitted ( gr - qc/0403036 ) shapiro , s. l. & teukolsky , s. a. 1980 , , 235 , 199 . 1983 , black holes , white dwarfs and neutron stars , ( new york : john wiley and sons ) , chap . 7 . 1985 , , 298 , 34 . 1991 , , 66 , 994
. 1992 , , 54 , 2006 shibata , m. 1997 , , 55 , 2002 . 2000 , prog .
phys . , 104 , 325 . 2003 , , 595 , 992 . 2004 , , 605 , 350 shibata , m. , baumgarte , t. w. , & shapiro , s. l. 2000 , , 542 , 453 shibata , m & nakamura , t. 1995 , , 52
, 5428 shibata , m & shapiro , s. l. 2002 , , 572 , l39 stark , r. f. & piran , t. 1985 , , 55 , 891 ( erratum 56 , 97 ) . 1986 , in dynamical spacetimes and numerical relativity , ed .
j. m. centrella ( cambridge : cambridge univ .
press ) , 40 stergioulas , n. 2003 , living reviews in relativity , 6 , 3 thorne , k. s. 1998 , in black holes and relativistic stars , ed .
r. m. wald ( chicago : univ .
chicago press ) , 41 van leer , b. 1977 , j. comp .
phys . , 23 , 276 wald , r. m. 1984 , general relativity ( chicago : univ .
chicago press ) , chap .
12.1 wilson , j. r. 1979 , in sources of gravitational waves , ed .
l. l. smarr ( cambridge : cambridge univ .
press ) , 423 wilson , j. r. & mathews , g. j. 1989 , in frontiers in numerical relativity , ed . c. r. evans , l. s. finn , d. w. hobill ( cambridge : cambridge univ . press ) , 306 . 1995 , , 75 , 4161 york , j. w. , jr .
1979 , in sources of gravitational waves , ed .
l. l. smarr ( cambridge : cambridge univ .
press ) , 83