Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: in general relativity , a black hole may be formed during the evolution from asymptotically flat initial data where none is present . consider a one - parameter family of regular asymptotically flat initial data . it is not difficult to find such families which form a black hole for some range of the parameter ( strong data ) but disperse for another range ( weak data ) . the boundary between the two regimes is the black hole threshold . in what is called type ii critical collapse , the black hole mass can be made arbitrarily small by adjusting the parameter @xmath0 of the initial data to its critical value @xmath1 . near the threshold , the final black hole mass @xmath2 then scales as @xmath3 where @xmath4 is a constant . @xmath4 depends on the family , but the transcendental number @xmath5 is universal  it depends on the type of matter but not on the family of initial data . type ii critical collapse was originally discovered by choptuik in the spherically symmetric massless scalar field @xcite , but has since been found in many simple matter systems in spherical symmetry , and also in axisymmetric gravitational waves @xcite . a review of the field is @xcite . type ii critical phenomena can be described in dynamical systems terms : the phase space of the system has ( at least ) two attracting fixed points , namely black holes and dispersion . the boundary between the two basins of attraction , the critical surface , contains a critical point : it is an attractor within the boundary surface , and a repeller only in the one remaining direction . this means that it must have precisely one unstable linear perturbation , with the property that adding a bit of that perturbation with one sign leads to collapse , while adding it with the opposite sign leads to dispersion . in type ii critical collapse the critical point is either a discretely self - similar ( dss ) or a continuously self - similar ( css ) spacetime . type ii critical collapse is interesting , among other things , because the maximum value of the curvature in a subcritical evolution , and the maximal value of the curvature outside the black hole in a supercritical evolution , both diverge as @xmath6 as the fine - tuning is improved . ( this is basically the same result as the black hole mass scaling , and similar results hold for any curvature invariant . ) from the dynamical systems picture it is clear that the end point of type ii critical collapse in the limit of perfect tuning of @xmath0 to its critical value @xmath1 is not a `` zero mass black hole '' but the critical solution itself . this solution has a naked singularity . it is therefore interesting to examine the global spacetime structure of the critical solution , and in particular its cauchy horizon . here we do this for the spherically symmetric massless scalar field , where the critical solution is dss . we focus on this system because css can be viewed as a limiting case of dss , and because the critical solution in the most interesting system in which type ii critical phenomena have been found , axisymmetric pure gravity , is also dss @xcite . in section [ section : kinematics ] we discuss the global structure of choptuik s critical solution kinematically . section [ section : scalarfield ] sets out the field equations for the real massless scalar field in spherical symmetry , in coordinates adapted to our problem , and describes the mathematical structure of the solution at the cauchy horizon of the singularity . sections [ section : numerics ] , [ section : continuation ] and [ section : globalimages ] show the results of our numerical integration of the critical solution , and section [ section : conclusions ] contains our conclusions . some technical details have been removed to appendices . a spacetime is discretely self - similar ( dss ) if there is a conformal isometry @xmath7 of the spacetime such that @xmath8 the value of the dimensionless `` logarithmic scale period '' @xmath9 is a geometric property of the spacetime , independent of coordinates . it is often useful to work in coordinates adapted to the symmetry at hand . a generic self - similar and spherically symmetric metric can be written as @xmath10 where @xmath11 is the metric on the unit 2-sphere , and where @xmath12 , @xmath13 , @xmath4 and @xmath14 are functions of @xmath15 and @xmath16 . this metric is dss if and only if they are periodic in @xmath15 with period @xmath9 . it is continuously self - similar ( css ) if these functions are completely independent of @xmath15 . we assume that the signature is @xmath17 , and that the metric is non - degenerate . this leads to the inequality @xmath18 . we also assume @xmath19 for definiteness . any four - dimensional spacetime splits into a product of a two - dimensional spacetime ( the reduced manifold ) and a round two - sphere of area @xmath20 . the area radius @xmath21 is a scalar in the reduced manifold . here the coordinates on the reduced manifold are @xmath15 and @xmath16 , and the area radius is given by @xmath22 . geodesics in the reduced spacetime are radial ( constant @xmath23 and @xmath24 geodesics in the full spacetime ) . the hawking mass @xmath25 is defined by @xmath26 . it is a scalar on both the full and the reduced spacetime . from @xmath25 we define the two dimensionless scalars @xmath27 and @xmath28 . it is easy to show that a spherical surface where @xmath29 is a closed trapped surface , and one where @xmath30 is an apparent horizon . in a dss spherical spacetime , @xmath31 and @xmath32 are periodic in @xmath15 .    radial null geodesics which are invariant under the symmetry ( [ dssgeometric ] ) are called self - similarity horizons ( ssh ) . they are the key to determining the causal structure . all coordinate systems of the form ( [ spherical_metric ] ) are related by coordinate transformations of the form @xmath33 where @xmath24 and @xmath34 are periodic in @xmath15 with period @xmath9 . we use this coordinate freedom to make all lines in the reduced manifold where @xmath35 into lines of constant @xmath16 . ( these can be either regular centers or central singularities . ) we also make all sshs into lines of constant @xmath16 where @xmath36 .    in order to discuss the global structure of the choptuik spacetime and its possible continuations we briefly review the kinematical results of @xcite . in a spherically symmetric dss spacetime , two kinds of singularities can be distinguished . from dimensional analysis it can be seen that the kretschmann scalar scales as @xmath37 for constant @xmath16 , and therefore the set @xmath38 is a central ( because @xmath21 scales as @xmath39 ) curvature singularity . we call this the kinematical singularity . geometrically , this singularity is either a point or a null line in the reduced spacetime . there are two types of self - similarity horizons that in @xcite we have called fans and splashes . the kinematical singularity is null if there is at least one splash . additional central singularities can arise where @xmath35 for all @xmath15 . ( we call these dynamical . ) because @xmath15 takes all values up to @xmath40 they are connected to the kinematical singularity . there are at most two of them , connected to the kinematical singularity at its ends . topologically , they are lines in the reduced manifold . all known type ii critical solutions in spherical gravitational collapse can be defined by the properties of self - similarity ( css or dss ) , analyticity at the past lightcone , and the requirement that they have a single unstable perturbation mode . a generic spherically symmetric dss scalar field solution is singular either at the center or the past lightcone . imposing analyticity at both the center and the past lightcone defines a nonlinear pde boundary value - problem which admits at most discrete solutions @xcite . only one such solution has been found , and it empirically turns out to have only one unstable mode , and to agree perfectly with the critical solution found previously in collapse simulations by choptuik @xcite .    the global structure of the choptuik solution up to the future light cone of the kinematical singularity ( which is a cauchy horizon ) is sketched in fig . [ figure : prech ] , together with the @xmath16 and @xmath15-lines of the three coordinate patches that we shall use in the numerical calculations . this structure is the same as for all other known type ii critical solutions in spherical symmetry . these solutions have a regular center @xmath41 in the past . as @xmath16 increases , the @xmath16 lines are at first timelike , so that @xmath42 . they become null at the fan @xmath43 where @xmath36 , @xmath44 and @xmath45 . as @xmath16 increases further , they are spacelike , so that @xmath46 spacelike . somewhere in the spacelike region @xmath13 changes sign . the @xmath16-lines become null again at the second fan @xmath47 , where @xmath36 , @xmath48 and @xmath49 . approaching the kinematical singularity @xmath38 from the range of `` angles '' @xmath50 it is a single point . @xmath51 and @xmath52 are its past and future lightcones .     are shown dashed ( we have assumed @xmath53 while in reality @xmath54 ) . note that the numerical domain is @xmath55 ( shaded ) , with periodic boundary conditions . we have illustrated the three coordinate patches we use for numerical work : in the past patch between @xmath56 and @xmath51 @xmath15-lines are spacelike . in the outer patch between @xmath51 and @xmath52 @xmath15-lines are timelike . in the future patch beyond @xmath52 they are null . the three patches together cover the entire spacetime without overlapping , and the coordinates @xmath16 and @xmath15 are continuous at the interfaces.,width=264 ]    we demonstrate below _ numerically _ that the curvature at the cauchy horizon @xmath47 is finite in choptuik s scalar field critical solution . furthermore all geodesics cross it in finite affine parameter . the spacetime can therefore be extended beyond , but this continuation is not unique . mathematically speaking we shall see that the solution is not analytic at the cauchy horizon in the limit coming from the past , and so there is no preferred analytic continuation to the future . ( the curvature is only @xmath57 from the past . ) the family of continuations in which the curvature is @xmath57 across the cauchy horizon and which are dss is parameterized by one free periodic function @xmath58 . physically speaking this function can be interpreted as data on the naked singularity which determine the continuation , in addition to the null data on the cauchy horizon . we now discuss the possible continuations that are allowed kinematically .     and @xmath52.,height=302 ]    in the simplest case , the @xmath16-lines to the future of the cauchy horizon are all timelike ( @xmath42 ) until a timelike singularity @xmath35 is reached . in our classification this is a dynamical singularity , while the kinematical singularity is a single point . this conformal structure is shown in fig .  [ figure : timelikesingular ] . if @xmath59 , or equivalently @xmath60 as @xmath61 , the conformal structure is the same , but with the singularity replaced by a regular center . nolan @xcite has drawn spacetime diagrams for spherically symmetric css spacetimes in which the point singularity at the origin of the cauchy horizon is only the starting point of an ingoing central null singularity . in our classification this is an extended kinematical singularity , which requires the existence of at least one splash , that is a line where @xmath36 again . [ figure : spacelikesingular ] shows the simplest generic possibility , in which the splash is followed by a spacelike dynamical singularity , which covers part of the naked null singularity . this spacetime structure can actually be realized in spherical css scalar field solutions @xcite . the same , and more exotic structures , can also be realized in spherical css perfect fluid solutions @xcite . contrary to our expectations , we have found that none of these exotic possibilities are realized as spherical dss continuations of the choptuik solution . there is a unique dss continuation with a regular center , and all other dss continuations have a timelike singularity . in hindsight , the reason appears to be that the null data on the cauchy horizon are extremely weak . with stronger data ( not associated with the choptuik spacetime ) we find different kinds of continuations .     and @xmath52 and a splash at @xmath62 . note that this extension is _ not _ realized as a dss continuation of the choptuik solution . , height=302 ] in spherical symmetry , there are four algebraically independent einstein equations , which can be taken to be @xmath63 , @xmath64 , @xmath65 and @xmath66 . the fourth of these is a combination of derivatives of the first three and can therefore be disregarded . the first three equations contain first and second derivatives of @xmath14 , but only first derivatives of @xmath12 , @xmath13 , and @xmath4 . in the following we assume that @xmath14 is a given function of @xmath16 and @xmath15 . the einstein equations then determine three independent linear combinations of the six first derivatives of @xmath12 , @xmath13 and @xmath4 . a non - degenerate basis of this three - dimensional space is @xmath67 , @xmath68 and @xmath69 . ( with @xmath70 given , this last term contains only first derivatives of @xmath12 , @xmath13 and @xmath4 . )    here we investigate massless real scalar field matter . the scalar field @xmath71 obeys the wave equation @xmath72 and the einstein equations can be written as @xmath73 the wave equation , when written as a second order pde , will in general contain both first derivatives of all four metric coefficients . however , by writing the wave equation in a geometrically defined first - order form , all metric derivatives except one can be eliminated . define two null derivative operators @xmath74 and @xmath75 with the usual convention that both point towards the future and that @xmath74 points inwards while @xmath75 points outwards . normalize them by imposing that @xmath76 and @xmath77 . we define the null derivatives of the scalar field as @xmath78 and @xmath79 . the massless wave equation in spherical symmetry can then be written as @xmath80 with the scalar @xmath81 . using the einstein equations , the curvature can be given in terms of @xmath82 and @xmath83 : @xmath84 ( other components of the ricci tensor vanish ) and @xmath85    the most general form of the scalar field compatible with dss is @xmath86 where @xmath87 and @xmath88 is a global constant . in the choptuik solution , @xmath89 empirically , so that @xmath71 is periodic in @xmath15 . this means that @xmath82 and @xmath83 have zero average in @xmath15 . moreover , @xmath82 and @xmath83 in the choptuik solution obey @xmath90 and so for @xmath83 . ( this of course implies zero average . ) as a consequence @xmath91 and so for other suitable metric fields . we shall assume these extra symmetries in our numerical work , but all our analytical expressions remain valid if these symmetries are dropped . we now describe three coordinate patches that cover the critical solution . we demand that both the past and the future lightcones of the singularity occur at lines of constant @xmath16 . this makes it easier to impose regularity at the center and past lightcone , and to investigate the behavior at the future lightcone ( cauchy horizon ) . it also allows us to match the coordinate patches without overlap . subject to these requirements , we have tried to make the field equations in each patch as simple as possible . based on the three patches , it is straightforward to construct a single smooth coordinate system covering the whole spacetime ( see section [ section : globalimages ] ) , but using it from the beginning would unnecessarily complicate our numerical work . we summarize a number of coordinate systems for spherically symmetric css and dss spacetimes , and their advantages and disadvantages , in the appendix .      on the past patch , which extends from the regular center to the past lightcone , we write the metric in terms of two free functions @xmath92 and @xmath93 as @xmath94 for @xmath95 . here @xmath96 is the scalar we defined above . these coordinates can be derived from the schwarzschild - type metric @xmath97 through the coordinate change @xmath98 and @xmath99 , and defining @xmath100 . the remaining gauge freedom of time relabelling is fixed by imposing that the past lightcone is on a constant @xmath43 line . note that we do not follow the convention of @xcite that @xmath101 at the center , which in @xcite forced us to introduce an additional free function of time in the definition of @xmath16 . the einstein equations are @xmath102a^{-2}-1 , \end{aligned}\ ] ] and the matter equations are @xmath103-{xu_{,\tau}}\over x(f+x ) } , \\ v_{,x}&=&{f[(1-a^2)v+u]+{xv_{,\tau}}\over x(f - x)},\end{aligned}\ ] ]    at the regular center @xmath104 we impose elementary flatness , that is the absence of a conical singularity . in order to do this , we define @xmath105 and @xmath106 and impose that both are regular even functions of @xmath16 at @xmath104 . at the past lightcone we have @xmath107 , which by our gauge choice means that @xmath108 there . we also impose the physical regularity condition @xmath109 on the past lightcone . the conditions of dss , regularity at the center and regularity at the past lightcone select a solution . the equations on the past patch are form - invariant under the linear coordinate transformation @xmath110 , @xmath111 , and @xmath32 , @xmath82 and @xmath83 unchanged . in the numerical results presented here we have set @xmath112 on the past patch . note that the regularity condition ( [ past_patch_regularity ] ) is coordinate - independent , as @xmath82 , @xmath83 and @xmath32 and @xmath21 are all scalars and @xmath113 on the lightcone .      on the outer patch , which extends from the past to the future lightcone , we write the metric in terms of @xmath93 , @xmath114 , and the auxiliary function @xmath115 as @xmath116 where @xmath117 is a function of @xmath15 only . @xmath118 is the scalar defined above . we fix the remaining gauge freedom by imposing that the past lightcone @xmath119 occurs at @xmath43 , and the future lightcone @xmath120 at @xmath47 . putting both lightcones on a line of constant @xmath16 requires @xmath115 to be non - constant . there is no outer coordinate patch that does not have at least one function like @xmath115 . as @xmath113 everywhere on the outer patch , @xmath15 is continuous between the past and outer patches . the metric equations are @xmath121\ , a^{-2}-1,\end{aligned}\ ] ] and the matter equations are @xmath122 note that the @xmath123 constraint equation is again linear in @xmath124 . the equations on the outer patch are form - invariant under the linear coordinate transformation @xmath125 , @xmath126 , and @xmath32 , @xmath127 , @xmath82 and @xmath83 unchanged . in the numerical results presented here we have set @xmath128 and @xmath129 on the outer patch . we shall see that at the cauchy horizon the solution is mildly singular . naive finite differencing breaks down there . instead we expand the generic solution around the cauchy horizon in terms of two free periodic functions @xmath130 and @xmath131 , and match this expansion to the numerical evolution at a small finite distance to the past of the cauchy horizon . before describing the full expansion , we focus on the origin of the singular behavior . equation ( [ outer_u ] ) becomes singular at the future lightcone because the denominator of the right hand side vanishes there . in contrast to the past lightcone , we do not have any freedom left to enforce the vanishing of the numerator as well . therefore we expect the solution to have some kind of singularity at @xmath52 . we rewrite ( [ outer_u ] ) @xmath132 where we have defined the metric function @xmath133 @xmath134 is positive on the outer patch and vanishes on the future light cone . the characteristics @xmath135 of eq . ( [ u_singular_eq ] ) are given by @xmath136 recall that we impose the gauge condition that the ch is at @xmath47 , or @xmath137 . we define the shorthand @xmath138 . we now make one fundamental assumption , namely that the spacetime admits regular null data on the cauchy horizon @xmath139 . this assumption is clearly necessary if we want to continue the spacetime through the ch , but here we make it simply because we have not been able to find a more general ansatz . we shall see later that it is sufficiently general to be matched to the critical solution that we have obtained numerically on the past patch . we therefore assume that @xmath83 is continuous , or @xmath140 by substituting this into eq . ( [ outer_a ] ) in the limit @xmath120 we find that @xmath32 is also continuous and @xmath141 obeys @xmath142 because we impose periodic boundary conditions in @xmath15 this ode has a unique solution . the physical significance of this is that the null data @xmath130 determine the geometry of the ch . similarly , from ( [ u_singular_eq ] ) in the limit @xmath143 we find that @xmath82 must be continuous , and @xmath144 is the unique periodic solution of @xmath145 this condition follows from the assumption of dss . finally , from ( [ ddef ] ) , ( [ outer_first ] ) and ( [ outer_second ] ) we find that @xmath134 is once differentiable , so that @xmath146 , and @xmath147 is given by @xmath148    at this point we introduce more shorthand notation . if @xmath149 is any periodic function ( with period @xmath9 ) , let @xmath150 be its average value , and let @xmath151 be its oscillatory part . let @xmath152 be the definite integral @xmath153 where @xmath154 is chosen so that @xmath152 has vanishing average . we can now integrate ( [ ucharacteristics ] ) for the @xmath82-characteristics to leading order and obtain @xmath155 we see that on a characteristic @xmath156 as @xmath157 . because @xmath158 is periodic in @xmath15 with period @xmath9 , an infinite number of oscillations in @xmath159 at constant @xmath15 pile up at the cauchy horizon @xmath139 . we can solve ( [ u_singular_eq ] ) to leading order by the method of characteristics . the general solution is @xmath160 where @xmath161 is an arbitrary periodic function with period @xmath9 and @xmath162    rewriting ( [ a0eqn ] ) as @xmath163 we see that @xmath164 , and so we can express @xmath165 and @xmath166 as @xmath167 our initial assumption that @xmath82 and @xmath83 are continuous at the light cone is therefore justified either if @xmath168 , so that @xmath169 , or if @xmath170 . we shall show numerically that @xmath165 is small but positive on the ch . in this case @xmath82 and @xmath83 are just @xmath57 , and the scalar field is therefore @xmath171 . in spherical symmetry the riemann tensor is determined completely by the ricci tensor , which in turn is quadratic in the partial derivatives of @xmath71 , see ( [ ricci ] ) . the curvature is therefore quadratic in @xmath82 and @xmath83 and so is @xmath57 . a similar analysis , with @xmath82 and @xmath83 interchanged , applies to the past lightcone . in the notation we have introduced here we can describe the past lightcone by saying that @xmath172 there ( because the null data @xmath173 on the past lightcone are large , @xmath174 ) but that the free coefficient @xmath175 vanishes identically ( because we have imposed analyticity as a boundary condition ) . we can not apply fuchsian techniques to our system of equations because they require the _ simultaneous _ vanishing on @xmath139 of the coefficients of @xmath176 and @xmath177 in equation ( [ u_singular_eq ] ) . however , the form ( [ uepsilon ] ) of the leading terms in @xmath82 suggests that the full non - linear solution can be written as a regular part , containing only integer powers of @xmath159 , plus a singular part which contains powers of @xmath178 . we can in fact construct a formal solution near the cauchy horizon in the form of an asymptotic double series : @xmath179 where @xmath180 stands for @xmath82 , @xmath83 , @xmath32 and @xmath127 , and @xmath181 here @xmath165 and @xmath182 are the quantities defined above in terms of @xmath130 . the expansion depends on the two free periodic functions @xmath183 and @xmath161 . by function counting we can therefore match any dss solution to this expansion . the first non - integer term appears in each variable at the following orders : @xmath184 in the previous section we obtained the coefficients of expansions ( [ expansion0][expansion3 ] ) up to @xmath185 . stopping there , the first order we neglect is @xmath186 . this truncation already depends on both free functions @xmath130 and @xmath131 and shows the singular behavior . it is also a sensible truncation numerically because @xmath165 turns out to be very small in the choptuik solution . going further , for the same reason there would be no point in including @xmath186 terms without also including all @xmath187 terms . it turns out that we need to go to @xmath188 . we have used the expansion to that order to check convergence . the expressions are given in appendix [ singular_expansion ] . our analysis of the possible continuations of the critical solution in section  [ section : kinematics ] has shown that we can cover the entire future of the cauchy horizon in a single patch if we make @xmath15 an ingoing null coordinate . this means setting @xmath189 and @xmath49 . in order to put the center @xmath190 at a known coordinate location , we also set @xmath191 . we choose @xmath192 here so that @xmath16 increases as we extend the spacetime away from the cauchy horizon . we parameterize this metric in terms of the scalar @xmath32 and a coefficient @xmath180 ( not the same as @xmath180 in the past patch ) : @xmath193 regularity of the metric requires @xmath118 and @xmath194 . the field equations are @xmath195a^{-2}-1 , \\ u_{,x}&=&{f[(1-a^2)u+v]-xu_{,\tau}\over x(f+x ) } , \\ v_{,x}&=&{u+(1-a^2)v\over x}. \label{futureeq5}\end{aligned}\ ] ] at the cauchy horizon @xmath196 we impose the coordinate condition @xmath197 . the equations are invariant under @xmath110 , @xmath111 . we use this to set @xmath198 on this patch . as a consequence @xmath113 on the lightcone , and so @xmath15 is continuous between the outer and future patches . we find an asymptotic expansion around the cauchy horizon in terms of two free functions @xmath130 and @xmath131 . @xmath130 is given by the null data on the cauchy horizon , and so is the same as on the outer patch . @xmath165 is therefore the same on both sides . @xmath173 obeys the same ode , ( [ u0eqn ] ) , as on the outer patch , and so is the same function . there is no need , however , to make @xmath131 the same on both sides , as we would not gain any differentiability by doing so . instead we consider @xmath131 as free `` data on the naked singularity '' , and we shall find experimentally how the global structure of the spacetime is influenced by this choice . we have carried out a brand new numerical computation of the choptuik spacetime , improving the precision of our previous calculations @xcite by roughly four orders of magnitude . this was mainly needed to assert without doubt that @xmath165 is different from zero , even though it is extremely small . essentially , our new scheme uses shooting methods on the @xmath16 axis , instead of relaxation methods . this allows us to improve the treatment of the regular singular points of the equations ( the center and the lightcones ) by using taylor expansions . we still work with pseudo - spectral fourier techniques in @xmath15 because the solution is periodic . we shall see however that the particular structure of the fourier transform of the choptuik spacetime poses an unexpected problem when combining taylor and fourier expansions .    in this section we explain the numerical scheme and present the results for the first two patches , with special emphasis on convergence properties and error analysis .        for definiteness , let us suppose we work on the past patch . following @xcite we discretize our @xmath9-periodic fields @xmath199 , where @xmath200 stands for any of the set @xmath201 , using @xmath202 equidistant points in one period : @xmath203 for @xmath204 . in this way we transform our 1 + 1 pde problem for @xmath199 into an ode problem for the modes @xmath205 . the essential idea of pseudo - spectral methods is to carry out algebraic operations pointwise on the @xmath206 and @xmath15-differentiation / integration on the @xmath207 , switching from one to the other with a fast fourier transform algorithm . the main drawback of the method is the aliasing problem : pointwise products of fields ( nonlinearities ) generate high frequency modes which can not be sampled with only @xmath202 points . we partially solve that problem by doubling the @xmath207 ( padding with zeros ) before going to the @xmath206 , carrying out the necessary algebraic operations on the doubled @xmath206 and going back to the @xmath207 , then halving the @xmath207 and thus throwing away high frequency noise . we have tried other possibilities , such as padding with @xmath208 or @xmath209 instead of @xmath202 zero fourier components , or extrapolating the fourier coefficients using the observed fact that high frequency modes have a simple exponential dependence on frequency ( see below ) , but the results are not improved . aliasing can only be reduced by going to higher @xmath202 . from a numerical point of view , we are only safe from aliasing when the amplitude of the modes we are cutting off is below machine precision .    because all our fields are real @xmath210 . furthermore , the metric fields @xmath32 and @xmath180 are even [ in the sense @xmath211 and therefore their fourier transform only contains even @xmath212 modes , while the matter fields @xmath82 and @xmath83 are odd [ in the sense @xmath90 ] and therefore their fourier transform only contains odd @xmath212 modes . taking this symmetry into account an even or odd field @xmath200 sampled with @xmath202 points per period is encoded by @xmath214 independent non - zero complex modes . as we have 4 independent variables @xmath215 , the ode system we solve comprises @xmath202 complex or @xmath216 real variables . in our calculations we have used @xmath217 and @xmath218 . previous investigations used @xmath219 . see appendix c of ref . @xcite for a complete discussion of our fourier pseudo - spectral method . we can not cross the lightcones during the integration of the ode in the @xmath16-axis because they are regular singular points of the equations . therefore we perform consecutive shooting calculations on the past , outer and future patches , in this order because we need information from the first to shoot the second and from the second to shoot the third . the issue of error propagation becomes very important . again , we describe the past patch for definiteness . given the equations ( [ pastmetriceqs][pastmattereqs ] ) and free data @xmath220 , @xmath221 at the center , we calculate the solution at @xmath222 slightly larger than @xmath56 using a second - order power expansion [ leaving errors of order @xmath223 . from these data we integrate the ode system forward in @xmath16 using finite differencing , up to @xmath224 . in the same way , given free data @xmath225 and the gauge condition @xmath226 , we calculate the solution at @xmath227 slightly smaller than @xmath51 [ this time with errors of order @xmath228 and integrate backward in @xmath16 up to the same @xmath229 . finally we use newton s method to look for the free data which brings the mismatch between both integrated solutions at @xmath229 down to a minimum , typically of order @xmath230 . ( this is the machine precision of @xmath231 reduced by a factor @xmath232 due to the calculations at each step and a factor @xmath233 from the ode - integration along @xmath234 points . ) we use a grid that becomes logarithmic near the regular singular points , with maximum stepsize @xmath235 at some intermediate position . that grid is constructed using the transformation @xmath236 from a grid of equidistant points in @xmath237 between the values @xmath238 and @xmath239 corresponding to @xmath222 and @xmath240 respectively . near the two endpoints we have @xmath241 and @xmath242 . we integrate on a fixed grid in @xmath16 rather than using a variable stepsize method in order to check for convergence with @xmath235 . this gives us a good estimate of the underlying discretization error . concerning the ode - integrator , we have tried several runge - kutta methods , both explicit and implicit ( gauss - legendre ) , with different convergence orders @xcite . in general , implicit methods are better suited to our problem than explicit methods , particularly for high @xmath202 , because high frequencies make the problem stiffer . we choose an irk2 method ( implicit , 2 stages , order 4 ) , which is a compromise between the accuracy of a high order method and the clarity of convergence of a low order method . our implicit schemes are implemented by iteration until @xmath243-differences between successive iterations converge below @xmath244 . we can not get closer to actual machine precision ( @xmath231 ) due to the accumulation of roundoff error . an implicit step typically takes 5 to 15 iterations to converge . the natural choice for the coordinate of the center is @xmath245 , but there is no preferred value for the past lightcone coordinate and we choose @xmath112 . with the set of parameters @xmath246 our newton s method converges to the free data plotted in fig . [ freedata ] ( see also table [ fouriertable ] ) , with a value @xmath247 which improves the precision of our previous result @xmath248 more than four orders of magnitude . the metric and matter fields integrated from those free data are given in fig . [ 3dinnerplots ] .     best free data on the singular points : @xmath249 ( continuous line ) , @xmath250 ( dotted line ) and @xmath251 ( dash - dot line ) . ]    [ cols="^,<,<,<,<,<,<",options="header " , ]     finally , it is interesting to see that going to @xmath188 is not essential for obtaining an accurate result . table [ eps_conv ] compares the results for @xmath165 using a zeroth - order expansion ( that is , we only include the terms of orders @xmath252 and @xmath178 ) with those from the first order expansion ( which also includes orders @xmath159 and @xmath253 ) . it is clear that they both converge to the same number , even though with very different rates of convergence . we can not continue the solution as flat empty minkowski spacetime after the cauchy horizon because we have a small amount of ingoing radiation there , given by @xmath254 . the simplest continuation to look for is one with a regular timelike center , so that the conformal diagram is the same as for minkowski spacetime . in this case both @xmath82 and @xmath83 are small on the entire future patch , and we can obtain an approximate solution in perturbation theory around minkowski space , using the magnitude @xmath255 of @xmath256 as the small parameter . to leading order in @xmath165 we obtain the dalembert solution on flat spacetime : @xmath257 + o(\epsilon^{3/2 } ) ,   \label{upert } \\ v & = & \epsilon^{1/2 } \left[-g'(\tau ) + ( 1+x)\frac{f(\hat\tau)+g(\tau)}{x } \right ] + o(\epsilon^{3/2 } ) , \label{vpert}\end{aligned}\ ] ] with @xmath258 . in order to match the null data on the cauchy horizon @xmath259 , we need @xmath260 . recall that @xmath254 is given numerically by ( [ v_f_numerical ] ) . if we want to have a regular solution at the center we need @xmath261 . the solution is then completely determined , and it is clear that a nonlinear solution exists of which this is the leading order , and which can be found numerically .    because the null data @xmath256 on the cauchy horizon are dss , it appears highly unlikely that there is another continuation with a regular timelike center that is not dss . we have obtained the ( probably ) unique regular continuation by shooting from expansions around the cauchy horizon and a regular center . the free data for the shooting algorithm ( given by @xmath262 at the ch and @xmath180 , @xmath176 at the center ) were obtained from the flat spacetime approximation . in this case we use an irk1 integrator and @xmath263 with very good convergence , that we do not show again . the fields @xmath32 , @xmath180 , @xmath82 and @xmath83 are shown in fig . [ choptuon_future ] . we now consider the other possible continuations , in particular those of figs . [ figure : timelikesingular ] and [ figure : spacelikesingular ] . before performing a numerical search of the possibilities in the @xmath131 space , we study the equations of motion in the future patch . we proceed in four steps :    [ [ a - necessary - condition - for - another - ssh - is - a22 . ] ] a necessary condition for another ssh is @xmath264 . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    in order to obtain fig . [ figure : spacelikesingular ] , or any even more exotic continuation , we must have a self - similarity horizon before the central singularity occurs at @xmath265 . a self - similarity horizon is a dss line @xmath266 ( periodic ) where @xmath36 . the only factor in @xmath267 that can vanish ( while the metric is regular ) is @xmath268 . we have @xmath269 and so @xmath270 at the cauchy horizon @xmath197 and @xmath271 because @xmath272 is small , and so @xmath273 and @xmath42 at least to the immediate future of the cauchy horizon . if there are more self - similarity horizons to the future of the cauchy horizon , we must have @xmath197 again there , and therefore @xmath274 in some intermediate region . @xmath275 must increase from @xmath276 to @xmath264 in order to achieve this . [ [ once - a21-for - any - tau - a2to-0-for - that - tau . ] ] once @xmath277 for any @xmath15 , @xmath278 for that @xmath15 . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    @xmath32 is given by the constraint ( [ a_constraint ] ) which reduces at the cauchy horizon to @xmath279 and this means that @xmath280 there . now @xmath32 obeys the evolution equation @xmath281 recall that in the future patch @xmath192 and increases towards the future . therefore @xmath282 is a repeller in the absence of matter , but any outgoing radiation @xmath82 drives @xmath32 to smaller values . hence , once @xmath32 becomes smaller than 1 for some value of @xmath15 , it will afterwards decrease to 0 for that @xmath15 whatever happens . because @xmath283 , this means that we are in a negative mass regime and remain there until we reach the central singularity @xmath104 , which is therefore timelike and has negative infinite mass at least for this @xmath15 . the singularity can not be reached with a finite value of @xmath284 . [ [ is - a21-at - the - singularity - possible - at - least - for - some - values - of - tau ] ] is @xmath285 at the singularity possible at least for some values of @xmath15 ? + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    if we choose @xmath131 much larger than @xmath286 , the @xmath287 term will drive @xmath32 to @xmath288 almost immediately . however , we know that the regular solution does not have negative mass and therefore it seems plausible that functions @xmath289 close to the regular case could generate a singularity with positive mass . we shall now argue that this is not true . @xmath32 can only increase from its ch value if @xmath290 is very small . in order to explore this regime we return to the approximation of perturbing around minkowski spacetime , but now dropping the assumption that @xmath261 . the result can be summarized as : @xmath291 ^ 2   \\ & + &   \epsilon \ ; o(x^{-1 } ) + o(\epsilon^{3/2 } ) . \nonumber\end{aligned}\ ] ] the function @xmath292 is positive and always smaller than 2 , and is independent of @xmath14 . all singular terms vanish for @xmath293 . for every other @xmath14 , the function @xmath32 becomes smaller than 1 , which corresponds to negative mass , as the center @xmath104 is approached . this is due to the divergent terms @xmath294 in ( [ upert ] ) and ( [ vpert ] ) . the crucial point is that @xmath275 is integrated from @xmath295 and therefore the term in the perturbation expansion that makes the mass nonzero has coefficient @xmath296 , and so the mass can not be positive . the negative mass regime is reached while perturbation theory still applies , and we have seen that afterwards the mass must decrease indefinitely .    with these arguments we have ruled out the possibility of having positive mass at the singularity for very large @xmath131 and for very small ( in particular , close to the regular case ) @xmath131 . however we could be missing some intermediate regime where perturbation theory does not apply . for these intermediate cases we observe that the minimum of @xmath32 is far from 1 ( and therefore we can not apply perturbation theory ) , but its maximum is positive and close to 1 , and only becomes negative in the vicinity of the singularity . there could be solutions where @xmath32 is positive at the singularity for some values of @xmath15 . numerically , however , we do not find that . we have performed a large number of numerical evolutions in this intermediate regime starting from @xmath297 , for @xmath298 and @xmath299 . even though it is difficult to evolve the system near the singularity , we always observe a final decay of @xmath32 to 0 . see fig . [ negativemass ] . therefore we believe that , starting from small null data @xmath256 at the ch it is not possible to form another ssh . therefore , either we have a regular center or a timelike singularity , and this singularity always has negative mass . in the approach to the center with @xmath300 we always have @xmath301 with both @xmath302 finite . it should be stressed that this scenario is a consequence of the small amplitude of the null data @xmath256 on the cauchy horizon in the choptuik solution . we have performed other evolutions starting from large null data on the cauchy horizon ( not @xmath256 of the choptuik solution ) where @xmath275 increases beyond 2 . maximum and minimum of @xmath32 for different initial data functions @xmath131 from numerical evolutions of the nonlinear equations . note that the @xmath16 axis becomes logarithmic when approaching the boundaries , while the @xmath32 axis is linear between the dotted lines at @xmath303 , @xmath304 , and logarithmic elsewhere . the initial @xmath131 have been obtained multiplying the regular data ( black ) by different constants : 0 ( red ) , -1 ( blue ) and 5 ( green ) . clearly , the larger @xmath131 is , the sooner @xmath32 decays . ] [ [ how - is - the - singularity - approached ] ] how is the singularity approached ? + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    we find that in all numerical continuations @xmath305 , @xmath306 , @xmath82 and @xmath83 tend to finite values and @xmath307 as the singularity @xmath104 is approached . if we assume that @xmath305 and @xmath308 as @xmath309 then the @xmath310 and @xmath311 equations to leading order give @xmath312 making the ansatz that @xmath313 we find from the @xmath176 and @xmath314 equations to the leading order , @xmath315 that @xmath316 ( the expansion also holds in the special case @xmath101 . ) the next order , @xmath317 , gives @xmath318 the metric in the future patch contains a residual gauge freedom worth one periodic function of @xmath15 . near the singularity @xmath104 , we can fix this gauge freedom by setting @xmath319 to an arbitrary value . ( in our numerical evolutions starting at the cauchy horizon , the gauge has of course been fixed already by setting @xmath320 at the cauchy horizon . ) @xmath144 and @xmath321 are then physical free data which determine @xmath141 and all other coefficients of the expansion . ( [ constraint ] ) can be solved for @xmath322 only if the right - hand side has vanishing average , which means that @xmath323 can not be set completely freely . ) this expansion is therefore generic in the sense of depending on two free functions after the gauge has been fixed . the behavior just described is what we observe numerically for large values of the free data @xmath324 at the cauchy horizon . for small values of @xmath324 we find @xmath325 , that is @xmath326 . note that by our ansatz of exact dss with @xmath89 , @xmath82 must be an odd function of @xmath15 with zero average . we observe in fact that @xmath158 goes to a fundamental frequency square wave of amplitude @xmath327 , that is , in half of each @xmath15-period @xmath328 , and @xmath329 in the other half . this is shown in fig . [ choptuon_singular ] .    as the center @xmath104 is approached with @xmath330 , we observe empirically that the @xmath15-derivatives become dynamically negligible . this means that different values of @xmath15 effectively decouple , that at each point @xmath15 the evolution equations become an ode system in @xmath16 , while the constraint becomes algebraic , and that the spacetime becomes locally css . it also means that the evolution equations become `` velocity - dominated '' in the sense that all derivatives transversal to the singularity ( here in spherical symmetry , this is only the @xmath15-derivative ) become dynamically irrelevant compared to the one derivative running into the singularity ( here , the @xmath16-derivative ) . it is known that generic spacelike singularities in general relativity with massless scalar field matter are velocity - dominated @xcite . here we find this to be the case only in the limit of small data @xmath324 . as this class of continuations seems to be locally css near the singularity , it is interesting to study the exactly css solutions from the point of view of a dss ansatz . starting from a generic dss scalar field ( [ generalphi ] ) we introduce a new ( gauge - dependent ) variable @xmath331 which coincides with @xmath332 at the future lightcone @xmath333 , and obeys the equation @xmath334 in exactly css solutions of the system the metric functions and @xmath302 are independent of @xmath15 , and hence @xmath335 is constant with value @xmath336 . the css solution with @xmath89 was found in closed form by roberts @xcite and is described in appendix  [ appendix : roberts ] . the css solutions with @xmath337 were studied numerically by brady @xcite . we find that the small @xmath324 continuations locally approach a roberts solution ( with @xmath89 ) with a value of the parameter @xmath0 of the roberts solution that depends on @xmath15 , and not , as one might have expected , one of brady s solutions . this is discussed in appendix  [ appendix : roberts ] . figures [ 3dinnerplots ] , [ choptuon_outer ] , [ choptuon_future ] and [ choptuon_singular ] show in a very detailed way the structure of the choptuik spacetime . however it is difficult to get from them an idea of what it looks like globally . in this final section before the conclusions we present a number of additional figures that will fill this gap . we do it in two very different ways .      as shown in fig . [ figure : prech ] , our three @xmath338 patches match continuously , but we do not expect the resulting global coordinate system to be differentiable at the interfaces between the patches . the critical spacetime itself , however , is differentiable ( analytic at the past lightcone and @xmath171 at the ch ) , and it must be possible to construct global coordinate systems which are at least @xmath339 .        one simple possibility is synchronous slicing plus area gauge : @xmath340 we add the gauge condition @xmath341 on the past lightcone of the singularity . the coordinate transformation @xmath342 can be easily integrated and it is shown on the left panel of fig .  [ coordinates ] for a single @xmath9-period in @xmath15 . another simple possibility is double null coordinates : @xmath343 with gauge condition @xmath344 at the center , where @xmath345 is the time coordinate constructed in the previous paragraph . the coordinate change @xmath346 is shown on the right panel of fig . [ coordinates ] for a single @xmath9-period in @xmath15 . the fields @xmath32 , @xmath82 and @xmath83 in double null coordinates are shown in fig . [ fieldsuv ] . as these fields are spacetime scalars , they should be analytic in these coordinates at the past lightcone , and @xmath171 at the cauchy horizon . in the plot it looks as if they are only @xmath57 at the cauchy horizon , but that is due to a lack of resolution in the plot : the slopes on the two sides of the cauchy horizon are dominated by @xmath347 and @xmath348 , which are discontinuous at the resolution of the plot , but close enough to the cauchy horizon , the slope becomes @xmath349 , which is continuous . we shall now consider the spherical dss scalar field as an infinite - dimensional dynamical system where @xmath16 is the `` time '' coordinate . the dynamical variables in this system are @xmath350 , @xmath351 and @xmath149 ( or @xmath352 in the outer patch ) . the variable @xmath353 is not independent , but given by a constraint . however , many solutions of the dynamical system correspond to the same spacetime , namely all that are related by the coordinate transformations ( [ dssgaugefreedom ] ) . the pair of periodic functions @xmath350 and @xmath351 describes a closed , possibly self - intersecting curve in the @xmath354 plane . an entire evolution in @xmath16 gives a surface that is topologically a cylinder , which we may call a phase portrait of the solution . clearly the surface itself is invariant under the coordinate change @xmath355 . if we consider @xmath15 as our `` space '' coordinate in the usual `` 3 + 1 '' split , then by looking at the phase portrait we have eliminated the `` spatial '' gauge freedom . the `` slicing '' freedom @xmath356 however does change the shape of the phase portrait , so it is not completely gauge - invariant . the choptuik spacetime up to the ch is given in fig . [ dynbefore ] and its regular continuation is shown in fig . [ dynafter ] .         phase - portrait of the regular continuation of the choptuik spacetime from the ch at @xmath198 to the regular center at @xmath357 . again , from @xmath259 to @xmath358 we have distorted the axis to show the infinite number of oscillations that pile up there . from @xmath358 to @xmath359 we use a logarithmic axis and from this middle point to @xmath360 we use a different logarithmic axis . ]    imposing css means that the system is independent of @xmath15 and hence the phase - portrait in @xmath361 reduces to a line that can be easily projected on the @xmath354 plane . then the whole evolution of the system can be described by this curve in the @xmath354 plane , which is now completely gauge - invariant . this is essentially what brady has done @xcite .    in order to understand the singular continuations of the critical spacetime , we first look at phase portraits of the roberts spacetime . fig . [ roberts ] shows the phase - lines of the roberts solution for several values of @xmath0 in both branches ( see also appendix [ appendix : roberts ] ) . we see that the shape of the curves become constant for very small values of @xmath0 , but just translated along the @xmath362-axis . on the other hand fig . [ dynsingular ] gives the phase - portrait two of our singular evolutions , for small and large initial @xmath131 , respectively . in both cases as @xmath301 squeezes the phase cylinder into the diagonal , following a roberts solution in the former case , and not doing so in the latter . this figure demonstrates why it is not possible to obtain @xmath363 in the continuations of the choptuik solution . the roberts spacetime with @xmath364 does have @xmath363 on the singularity , @xmath365 ( minkowski ) has @xmath366 , and all other @xmath0 have @xmath278 . the @xmath364 roberts solution has the phase line @xmath367 line . to reach @xmath363 approach a solution must approach this line from the very beginning , because it is unstable . this requires @xmath83 of order @xmath327 at least and on the cauchy horizon we only have @xmath83 of order @xmath368 . phase - lines of the roberts spacetime . the colors encode different values of @xmath0 , for both branches : 1 ( red ) , 0.99 ( orange ) , 0.9 ( yellow ) , 0.5 ( light green ) , 0.1 ( dark green ) , 0.01 ( light blue ) , 0.001 ( dark blue ) and 0.0001 ( purple ) . the future ( past ) branch of the singularity is denoted with thick ( thin ) lines starting with @xmath369 ( @xmath370 ) . note that the vertical lines @xmath367 corresponding to @xmath364 are unstable . note also that the smaller @xmath0 is , the longer the curves stay close to the unstable flat line @xmath371 . ]     phase - portrait of two singular continuation of the choptuik spacetime ( upper : initial @xmath289 of order @xmath368 ; lower : initial @xmath131 of order @xmath232 ) , together with a few roberts solutions for @xmath372 . the @xmath373 solution almost coincide with the external envelope of the upper phase - portrait and sets a maximum for the values that are realized in this solution . the lower phase - portrait is shown only for @xmath374 . for @xmath375 it is essentially on the diagonal in the @xmath354 plane although the amplitude of the oscillations is not the same as in the roberts solution . note that for small @xmath16 this phase portrait has been truncated also in @xmath82 , because @xmath82 becomes large . for @xmath376 the phase portrait is essentially on the @xmath82-axis . ] in a companion paper @xcite we investigate the constraints that the kinematic assumptions of spherical symmetry and discrete self - similarity impose on the causal structure . the key elements of our analysis were the self - similarity horizons  radial null geodesics that are mapped onto themselves by the self - similarity . these come in two types : a `` fan '' connects a point singularity to a piece of null infinity , while a `` splash '' connects a piece of null singularity to a point at null infinity . all spherically symmetric and self - similar spacetimes can be enumerated in terms of a sequence of fans and splashes . we found that the singularity is central and consists of a middle segment which is either a point or null line , flanked by two segments which can be timelike , spacelike , or null lines , or which can be absent . there are never two or more disconnected singularities .    in this paper , we have focussed on scalar field matter and in particular choptuik s critical solution in gravitational collapse . this solution is found as an intermediate attractor in the evolution of regular asymptotically flat scalar field initial data near the black hole threshold @xcite . independently , it can be constructed as the solution of a nonlinear pde boundary value problem from the assumptions of spherical symmetry , discrete self - similarity , and analyticity at the past center and past lightcone of the singularity @xcite . it can then be uniquely continued up to the future lightcone ( cauchy horizon ) of the singularity . @xcite found that the scalar field oscillates roughly as @xmath377 as the cauchy horizon @xmath139 is approached . it was suggested on theoretical grounds that these oscillations are damped by a factor @xmath178 with @xmath169 , but the numerical value of @xmath165 was too small to distinguish it from zero numerically . we have repeated the numerical calculations of ref . @xcite from scratch and have increased the accuracy by more than four orders of magnitude . we have found the positive value for @xmath165 given in eq . ( [ epsilon_value ] ) . therefore the scalar field is continuous on the cauchy horizon , with null data given in eq . ( [ v_f_numerical ] ) . the structure of the fields near the cauchy horizon is quite complicated , and has been discussed in subsections [ section : singularch ] and [ section : expansion ] . all possible dss continuations beyond the cauchy horizon are determined by the ( fixed ) null data on the cauchy horizon , and one ( arbitrary ) periodic function @xmath161 that can be thought of as data emerging from the naked singularity . ( in the absence of the naked singularity , null data on the cauchy horizon would alone determine the continuation . ) there is a unique dss continuation that has a regular timelike center except for a naked point singularity . for all other values of the arbitrary function @xmath161 , the continuation has a timelike central singularity with infinite negative mass . in particular , neither fig . [ figure : spacelikesingular ] , nor the even more exotic spacetime diagrams found in self - similar perfect fluid solutions @xcite can arise as continuations of the choptuik solution ( although they can probably arise in other dss scalar field solutions ) . the global structure of the choptuik solution is of interest partly because of the connection between critical collapse and cosmic censorship . choptuik s work has established that ( assuming the validity of general relativity at arbitrarily high curvatures ) a naked singularity can be formed in the spherical einstein - scalar field system from generic regular and asymptotically flat initial data by fine - tuning any one parameter to the black hole threshold . solutions with a naked singularity therefore form a subset of codimension one ( the black hole threshold in the space of initial data ) in all solutions arising from regular data . the choptuik spacetime is an attractor on the critical surface and therefore , assuming that it is actually a global attractor on the surface , we can conclude that every naked singularity will have , at least in a neighborhood of the singularity , the structure of the choptuik spacetime . we can now summarize this structure as follows : the curvature at the cauchy horizon of the singularity is finite but not differentiable , with an infinite number of damped oscillations piling up as the cauchy horizon is approached . the continuation of the spacetime beyond the cauchy horizon is not unique , but we have shown that if it is dss then the naked singularity is either a single point or timelike with infinite negative mass . the spacetime near the timelike singularity is locally css and velocity - dominated . the question remains whether this structure is stable against perturbations which break self - similarity and/or spherical symmetry . because the background is spherically symmetric and periodic in @xmath15 , all such perturbations are of the form @xmath378 , where @xmath379 is periodic , times a suitable scalar , vector or tensor spherical harmonic . we have shown in the past @xcite that all but one of these modes decay , with @xmath380 . ( the one growing mode determines the critical exponent for the black hole mass . ) the functions @xmath379 were calculated explicitly only between the regular center and past lightcone , but as they are by construction analytic at the past lightcone , they can be analytically continued up to the cauchy horizon , with the same @xmath381 .    what needs to be done is to check that the functions @xmath379 remain bounded as @xmath16 approaches the cauchy horizon . we expect that this is true . nolan and waters @xcite have investigated a massless minimally coupled , non - spherically symmetric scalar field propagating on a class of spherically symmetric css spacetimes , and find that its gradient remains finite at the cauchy horizon . the structure of the full perturbation equations is similar .    as the functions @xmath379 can not be analytic at the cauchy horizon ( because the background is not ) , the spectrum of @xmath381 need not be the same in the continuation beyond the cauchy horizon , and is likely to depend on the choice of continuation . we have finally given , for the first time , a number of global images of the choptuik spacetime , trying to minimize the gauge content of the pictures . we would like to thank james vickers for helpful discussions . this research was supported by epsrc grant gr / n10172/01 . in this appendix we review a variety of coordinate systems defined by giving @xmath14 as a function of @xmath15 and @xmath16 . this class covers the coordinate systems used in refs . @xcite , but does not include the lagrangian fluid coordinates used in @xcite . we classify the gauges within this class in three stages . 1 ) we fix @xmath70 . 2 ) we impose an algebraic relation between @xmath12 , @xmath13 and @xmath4 . 3 ) we parameterize @xmath12 , @xmath13 and @xmath4 in terms of two functions , say @xmath32 and @xmath180 . ( it will turn out to be useful to always use the scalar @xmath32 defined above as one of the two parameters . )    in a dss spacetime , the unknowns are periodic in @xmath15 . we therefore solve the einstein equations by `` evolving '' in @xmath16 , with periodic boundary conditions in @xmath15 . when we use the einstein equations to eliminate the metric derivative @xmath382 , @xmath82 and @xmath83 obey a pair of transport equation of the form @xmath383 where the dots stand for a known function of @xmath82 , @xmath83 , @xmath32 and @xmath180 but not their derivatives . by a suitable parameterization of @xmath12 , @xmath13 and @xmath4 in terms of @xmath32 and @xmath180 , the einstein equations can always be brought into the form of two odes in @xmath16 , which for our purposes are evolution equations , and one ode in @xmath15 , which for our purposes is a constraint : @xmath384 this last equation can be made linear in @xmath124 ( that is , in @xmath31 ) by a suitable choice of @xmath180 , at least in all three coordinate systems we use . this inhomogeneous linear equation can then be solved uniquely for @xmath32 in terms of @xmath82 , @xmath83 and @xmath180 , using the fact that @xmath82 , @xmath83 and @xmath180 are periodic in @xmath15 and that we require @xmath32 to be periodic , too . in a css solution , where nothing depends on @xmath15 , the @xmath123 constraint becomes an algebraic equation linking @xmath82 , @xmath83 , @xmath32 and @xmath180 . the regular center is both timelike and an @xmath16 line , but @xmath15 is finite there . this requires @xmath35 at the center , and @xmath386 elsewhere . the simplest choice is @xmath385 . on a past patch , it is possible to choose the @xmath15 lines to be timelike , null or spacelike , or to change signature . we consider only the first two possibilities    [ [ bondi - past - patch ] ] bondi past patch + + + + + + + + + + + + + + + +    if we choose the @xmath15-lines to be null everywhere , this means @xmath189 . we then have @xmath49 between the past center and the past light cone , and @xmath45 beyond the past lightcone , so that @xmath12 and @xmath13 both change sign at the past lightcone . @xmath385 and @xmath189 has been used by brady @xcite to investigate css solutions with scalar field matter . he used the bondi metric coefficients that are traditionally called @xmath387 and @xmath388 as parameters . in a dss solution this gives a constraint for a combination of @xmath387 and @xmath388 , as well as evolution equations for @xmath387 and @xmath388 . the best parameterization we have found uses the mass function @xmath32 and the bondi metric coefficient @xmath389 : @xmath390 where the upper sign applies when @xmath15 is an outgoing null coordinate , and the lower sign applies when @xmath15 is an ingoing null coordinate , assuming that @xmath194 . this gives an ode evolution equation and a linear ode constraint for @xmath124 , and an ode evolution equation for @xmath180 . [ [ schwarzschild - past - patch ] ] schwarzschild past patch + + + + + + + + + + + + + + + + + + + + + + + +    if the @xmath15-lines are spacelike everywhere , this implies @xmath391 . a natural algebraic condition to impose is to make the @xmath15 lines orthogonal to the @xmath21 lines . this means @xmath392 . @xmath15 is then @xmath393 where @xmath394 is the schwarzschild time coordinate . a possible parameterization is in terms of the metric coefficients often called @xmath32 and @xmath395 in schwarzschild coordinates . ( @xmath32 is the scalar we have defined above ) . the best parameterization we have found is by @xmath32 and @xmath100 , and is given in ( [ pastpatch ] ) above . it gives ode evolution equations for @xmath32 and @xmath180 , and a linear ode constraint equation for @xmath124 . the transport equations for @xmath82 and @xmath83 are linear in @xmath82 and @xmath83 because @xmath382 reduces to a function of @xmath32 , @xmath15 and @xmath16 . the lightcones are at @xmath396 ( although of course only one lightcone can be covered at one time by a past patch ) .      on an outer patch that stretches from the past to the future lightcone the @xmath15 lines must be timelike at least somewhere , and it is possible to make them timelike everywhere . the simplest choice for @xmath14 compatible with all these possibilities is @xmath14=1 . [ [ schwarzschild - outer - patch ] ] schwarzschild outer patch + + + + + + + + + + + + + + + + + + + + + + + + +    assuming that the @xmath15 lines are timelike everywhere means @xmath398 . we have @xmath49 on the past lightcone and @xmath45 on the future lightcone , so that @xmath13 has to change sign somewhere between . this `` @xmath13-surface '' is where the @xmath15 and @xmath16 lines are orthogonal . in @xcite one of us used a coordinate system in the outer patch that was based on schwarzschild coordinates . as in the past patch based on schwarzschild coordinates , this meant @xmath399 . however , when we impose @xmath399 together with @xmath397 , we face an unexpected coordinate singularity at @xmath104 . in particular , if we again use @xmath32 and @xmath100 to parameterize the metric , we obtain an ode constraint and an ode evolution equation for @xmath32 as before on the past patch , but for @xmath180 we obtain an equation of the form @xmath400 in order to make the coordinates regular at @xmath104 , we need to impose the vanishing of the numerator of this equation there , and that is effectively what was done in @xcite . but this introduces an additional boundary condition just to keep the coordinate system regular on a surface where the solution itself is perfectly regular , and that is why we do not use these coordinates here . [ [ best - buy - outer - patch ] ] best buy outer patch + + + + + + + + + + + + + + + + + + + +    a better algebraic condition to combine with @xmath397 is @xmath401 . a workable parameterization is in terms of the mass function @xmath32 and @xmath13 itself . this gives ode evolution equations for @xmath32 and @xmath13 , and a ( nonlinear ) ode constraint for @xmath32 . the lightcones are at @xmath402 . if we replace @xmath13 by @xmath403 , we have the parameterization given above in ( [ outerpatch ] ) , and the ode constraint for @xmath32 becomes linear in @xmath124 , while the wave equation becomes linear in @xmath82 and @xmath83 . [ [ a - c - unworkable ] ] @xmath404 : unworkable + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    it is compatible with @xmath397 to assume that the @xmath15 lines `` bend round '' to become spacelike at large and small @xmath16 . means this that @xmath4 changes sign . a natural choice is that @xmath12 and @xmath4 change sign together , and we can impose this by setting @xmath404 . we have not found a way of parameterizing this choice in a way that brings the einstein equations into the standard form ( [ standardform ] ) . parameterizing with @xmath13 and @xmath31 gives two roots for @xmath12 . using @xmath12 and @xmath31 gives two roots for @xmath13 . [ [ x - and - tau - orthogonal - unworkable ] ] @xmath16 and @xmath15 orthogonal : unworkable + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    we can also force the `` bending round '' by making @xmath15 and @xmath16 lines orthogonal everywhere , or @xmath405 . this implies that , while they are independent , @xmath12 and @xmath4 change sign together . parameterizing by @xmath12 and @xmath4 we obtain an evolution equation and a constraint for @xmath12 and a constraint for @xmath4 . but the constraint for @xmath4 is homogeneous , @xmath406 , so that @xmath4 can not be periodic . therefore this coordinate system is not compatible with css or dss . [ [ forcing - the - lightcones - unworkable ] ] forcing the lightcones : unworkable + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    instead of finding a working outer patch and then forcing the two lightcones to fall on the lines @xmath43 and @xmath47 by imposing @xmath407 , we can make @xmath12 a given function of @xmath16 . the simplest such choice is @xmath408 , which vanishes at @xmath409 . ( the two lightcones will be distinguished by the sign of @xmath13 ) . a parameterization that brings the einstein equations into standard form is @xmath410 this gives an evolution equation and a constraint for @xmath127 and a linear constraint for @xmath411 . however , the right - hand sides of all three equations are of the form @xmath412 . furthermore , at @xmath413 the numerator of @xmath414 depends on the matter through @xmath415 , while the numerator of @xmath416 depends on the matter through @xmath417 . ( the numerator of @xmath418 is proportional to that of @xmath414 . ) that means that we would have to impose two independent regularity conditions at @xmath413 . together these would fix @xmath82 and @xmath83 completely . therefore this coordinate system is not sufficiently generic near the line @xmath405 . from ( [ expansion1],[expansion2 ] ) we see that @xmath419 the regular coefficients to @xmath186 are as follows : @xmath420 is the unique periodic solution of @xmath421 and @xmath422,\\ a_1&=&-{\xi\over 2}(u_0 ^ 2-v_0 ^ 2 ) , \\ b_1&=&-{\xi'\over a_0}-{\xi\over 2a_0}\left(-3+a_0 ^ 2 + 2v_0 ^ 2\right),\\ b_2&=&-{a_1\over 2 a_0}b_1 \nonumber \\ & & -{\xi\over 4a_0 } \left(2a_0a_1-b_1u_0 ^ 2+b_1v_0 ^ 2 + 4v_0v_1\right ) , \\ d_1&=&-{a_0 b_1\over \xi } .\end{aligned}\ ] ] the regular coefficients to higher orders continue in this style , with @xmath423 the solution of a linear inhomogeneous ode , while @xmath424 , @xmath425 and @xmath426 are given algebraically . the nonvanishing singular coefficients up to @xmath427 are @xmath428    we only give two examples of how these coefficients are derived . substituting the ansatz into the @xmath310 equation and isolating the terms of @xmath429 in the result , we obtain @xmath430 we solve this for all @xmath15 and @xmath182 by setting @xmath431 and by making @xmath432 the unique solution of the ode @xmath433    if we substitute the ansatz into the @xmath176 equation ( [ u_singular_eq ] ) and isolate the terms of @xmath434 , we find    @xmath435    because @xmath436 , the derivatives @xmath437 cancel out . taking also into account that @xmath438 ( this is an accident at this particular order ) , the equation can be rewritten as @xmath439\hat u_{1+\epsilon}^{(i ) }   \nonumber \\ + \left[-\epsilon d_2 - 2a_0a_2)\check u_\epsilon\right ] \hat u_\epsilon + \left[-kd_2\check u_\epsilon\right ] \hat u_\epsilon ' + \left[\check v_{1+\epsilon}-2a_0\check a_{1+\epsilon}u_0\right ]   \hat a_{1+\epsilon}=0.\end{aligned}\ ] ]    the terms in square brackets all depend on @xmath15 and each multiplies a different known function of @xmath182 . to solve this for all @xmath15 and @xmath182 , we assign one term @xmath440 to each function of @xmath182 : @xmath441 the corresponding coefficients @xmath442 are the unique solutions of the odes @xmath443 where the source terms @xmath444 can be read off directly from ( [ u1pe ] ) . the calculation of the other terms proceeds in a similar manner at all orders . note that the functions of type @xmath445 , where @xmath180 stands for @xmath83 , @xmath32 and @xmath127 , are given algebraically and @xmath446 obeys an ode . for @xmath82 it is the other way around . the limit @xmath447 of our expansion exists . the regular @xmath448 terms in the series vanish identically , and so do some of the @xmath449 terms but not all of them . @xmath349 and @xmath347 for example vanish because they are proportional to @xmath173 , but @xmath348 is proportional to @xmath450 , and so encodes the curvature to @xmath186 .    in the limit @xmath451 the curvature components proportional to @xmath287 are no longer continuous at the ch because they are now periodic in @xmath452 , but for the same reason they remain bounded . the components proportional to @xmath453 and @xmath454 are still continuous even then because @xmath83 is @xmath186 . @xmath32 is also still continuous , with @xmath282 on the cauchy horizon . in the notation of @xcite , the roberts solution @xcite is given by @xmath455 , \\ \label{roberts3 } \phi(u , v ) & = & \frac{1}{\sqrt{16\pi g } } \log \frac{(1-p)v - u}{(1+p)v - u},\end{aligned}\ ] ] with @xmath0 a constant parameter . @xmath365 is minkowski spacetime with zero scalar field , and without loss of generality @xmath456 . only the regions @xmath457 are physical and without loss of generality , we consider the right side ( @xmath458 ) of the spacetime .    for @xmath459 , the lines @xmath460 are central curvature singularities : the mass function @xmath461 diverges on them . the line @xmath462 , @xmath463 is timelike and has ( infinite ) negative mass , and for @xmath463 forms the past branch of the roberts singularity . the line @xmath464 is timelike with negative mass for @xmath465 , null with zero mass for @xmath364 , and spacelike with positive mass @xmath466 . for @xmath467 it forms the future segment of the singularity . some of our solutions become asymptotically css and have @xmath468 , therefore approaching the roberts spacetime . we know that the solutions develop a timelike , negative mass singularity , which must correspond to one of the two branches of the roberts singularity . in this appendix we study the issue of which of the two branches is actually approached in our numerical evolutions and for what values of @xmath0 . our scalars @xmath82 and @xmath83 are @xmath469 where both denominators are positive in the region of interest . it is clear that @xmath0 is covariantly defined by @xmath82 and @xmath83 on the light cones @xmath470 or @xmath471 , respectively . however , at the singularities @xmath460 the scalars take values @xmath472 for any @xmath473 in the roberts solution . if we use an expansion using a generic self - similar coordinate @xmath159 such that @xmath474 with the singularity at @xmath139 [ that is , @xmath475 : @xmath476 therefore @xmath0 is covariantly given by the ratio of the rates at which @xmath82 and @xmath83 approach their values at the singularity : @xmath477 this determines both @xmath0 and which branch of the singularity we approach locally : the upper sign applies for @xmath478 ( past branch , turned upside down ) and the lower sign for @xmath479 ( future branch of the singularity ) . we can relate the coordinates @xmath480 of ( [ roberts1]-[roberts3 ] ) to our coordinate @xmath16 through @xmath481 this expression applies only for negative @xmath16 and so only covers the future branch of the singularity . we obtain the expressions for the past branch exchanging the role of the functions @xmath82 and @xmath83 . that explains why exchanging @xmath82 and @xmath83 in eq . ( [ uvratio ] ) amounts to a change of branch . it also shows that we have to change branch four times per @xmath9-period as we move through the four quadrants determined by the lines @xmath482 and @xmath483 .        c. gundlach , living reviews in relativity * 2 * , 4 ( 1999 ) , published electronically at http://www.livingreviews.org . for an updated version see c. gundlach , critical phenomena in gravitational collapse , to be published in phys . , preprint gr - qc/0210101 .
Abstract:

Answer: at the threshold of black hole formation in the gravitational collapse of a scalar field a naked singularity is formed through a universal critical solution that is discretely self - similar . we study the global spacetime structure of this solution . it is spherically symmetric , discretely self - similar , regular at the center to the past of the singularity , and regular at the past lightcone of the singularity . at the future lightcone of the singularity , which is also a cauchy horizon , the curvature is finite and continuous but not differentiable . to the future of the cauchy horizon the solution is not unique , but depends on a free function ( the null data coming out of the naked singularity ) . there is a unique continuation with a regular center ( which is self - similar ) . all other self - similar continuations have a central timelike singularity with negative mass .