this paper is a detailed version of my preprint @xcite , which was written about five years ago . here
i do not discuss results that have appeared since then .
i plan to survey them soon in another paper .
the subject is now evolving into a _ real algebraic knot theory_.    this paper is dedicated to the memory of my teacher vladimir abramovich rokhlin .
it was v.  a.  rokhlin , who suggested to me , a long time ago , in 1977 , to develop a theory of real algebraic knots .
he suggested this as a topic for my second dissertation ( after phd , like habilitation ) . following this suggestion
, i moved then from knot theory and low - dimensional topology to the topology of real algebraic varieties .
however , in the topology of real algebraic varieties , problems on spatial surfaces and plane curves were more pressing than problems on spatial curves , and my second dissertation defended in 1983 was devoted to the constructions of real algebraic plane curves and spatial surfaces with prescribed topology .
the change in the topic occured mainly because i managed to obtain decent results in another direction , on plane curves .
there was also a less respectable reason : i failed to relate the traditional techniques of classical knot theory to real algebraic knots .
one of the obstacles was a phenomenon which became the initial point of this paper .
a large part of the traditional techniques in knot theory uses plane knot diagrams , i.e. , projections of knots to the plane .
the projection of an algebraic curve is algebraic , and one could try to apply results on plane real algebraic curves
. however , the projection contains extra real points , which do not correspond to real points of the knot .
these points are discussed below . in the seventies they ruined my weak attempts to study real algebraic knots .
now they allow us to detect crucial differences between topological and real algebraic knots .
i am grateful to alan durfee , tobias ekholm , and v.  m.  kharlamov for stimulating conversations .    the lengthy informal introduction , which follows , is intended to explain the matter prior to going into details .
i can not resist the temptation to write in the style of popular mathematics and apologize to the reader whom this style may irritate .      in classical knot theory , by a link one
means a smooth closed @xmath2-dimensional submanifold of the @xmath0-dimensional sphere @xmath3 , i.e. , the union of several disjoint circles smoothly embedded into @xmath3 .
a link may be equipped with various additional structures such as orientation or framing and considered up to various equivalence relations like smooth ( or ambient ) isotopy , pl - isotopy , cobordism or homotopy .
see , e.g. , @xcite or @xcite .    in algebraic geometry
classical links naturally appear as links of singular points of complex plane algebraic curves . given a singular point @xmath4 of a complex plane algebraic curve @xmath5 , the intersection of @xmath5 with the boundary of a sufficiently small ball centered at @xmath4 is called the _ link of the singularity_. it provides a base for a fruitful interaction between topology and algebraic geometry with a long history and lots of important results .    another obvious opportunity for interaction between algebraic geometry and knot theory
is based on the fact that a classical link may emerge as the set of real points of a real algebraic curve .
this opportunity was completely ignored , besides that a number of times it was proved that any classical link is approximated by ( and hence isotopic to ) the set of real points of a real algebraic curve .
there are two natural directions in which algebraic geometry and knot theory may interact in the study of real algebraic links : first , the study of relationships between invariants which are provided by link theory and algebraic geometry , second , developing a theory parallel to the classical link theory , but taking into account the algebraic nature of the objects .
from the viewpoint of this second direction , it is more natural to consider real algebraic links up to isotopy consisting of real algebraic links , which belong to the same continuous family of algebraic curves , rather than up to smooth isotopy in the class of classical links .
i call an isotopy of the former kind a _ rigid isotopy _ , following the terminology established by rokhlin @xcite in a similar study of real algebraic plane projective curves and the likes ( see , e.g. , the survey @xcite ) .
of course , there is a forgetting functor : any real algebraic link can be regarded as a classical link and a rigid isotopy as a smooth isotopy .
it is interesting to see how much is lost under that transition .    in this paper
i point out a real algebraic link invariant which is lost .
it is unexpectedly simple . in an obvious sense it is a nontrivial vassiliev invariant of degree @xmath2 on the class of real algebraic knots (
recall that a knot is a link consisting of one component ) . in classical knot theory
the lowest degree of a nontrivial vassiliev knot invariant is @xmath6 .
thus there is an essential difference between classical knot theory and the theory of real algebraic knots .
of course this difference has a simple topological explanation : a real algebraic link is more complicated topologically , besides its set of real points contains the set of complex points invariant under the complex conjugation and a rigid isotopy induces an equivariant smooth isotopy of this set .
the invariant of real algebraic links which is defined below is very similar to the self - linking number of a framed knot .
in @xcite i call it also the _ self - linking number_. its definition looks like a replacement of an elementary definition of the writhe of a knot diagram , but taking into consideration the imaginary part of the knot .      here
i propose to change this name ( i.e. , self - linking number ) to _ encomplexed writhe _ , and , in general , since many other characteristics can also be enhanced in a similar way , i suggest a new verb _
encomplex _ for similar enhancements by taking into consideration additional imaginary ingredients .
this would agree with the general usage of the prefix _
en- _ which is described in the oxford dictionary of current english as follows : `` en- _ prefix _
@xmath7 forming verbs @xmath7 * 1 * from nouns , meaning ` put into or on ' ( _ engulf ; entrust ; embed _ ) , * 2 * from nouns or adjectives , meaning ` bring into the condition of ' ( _ enslave _ ) @xmath7 '' .
complexification _ does not seem to be appropriate for what we do here with the writhe .
a complexification of the writhe should be a complex counterpart for the writhe , it should be a characteristic of complex objects , while our enhancement of the writhe is defined only for real objects possessing complexification .
the linking number is a well - known numerical characteristic of a pair of disjoint oriented circles embedded in three - dimensional euclidean space . roughly speaking
, it measures how many times one of the circles runs around the other .
it is one of the most classical topological invariants , introduced in the nineteenth century by gauss @xcite .    in the classical theory ,
a self - linking number of a knot is defined if the knot is equipped with an additional structure like a framing or just a vector field nowhere tangent to the knot .
the self - linking number is the linking number of the knot oriented somehow and its copy obtained by a small shift in the direction specified by the vector field .
it does not depend on the orientation , since reversing the orientation of the knot is compensated by reversing the induced orientation of its shifted copy .
of course , the self - linking number depends on the homotopy class of the vector field .
a knot has no natural preferable homotopy class of framings , which would allow us to speak about a self - linking number of the knot without a special care on the choice of the framing .
some framings appear naturally in geometric situations .
for example , if one fixes a generic projection of a knot to a plane , the vector field of directions of the projection appears . the corresponding self - linking number is called the _ writhe _ of the knot .
however , it depends on the choice of the projection and changes under isotopy .
the linking number is a vassiliev invariant of order @xmath2 of two - component oriented links .
this means that it changes by a constant ( in fact , by @xmath6 ) when the link experiences a homotopy with the generic appearance of an intersection point of the components . whether the linking number increases or decreases depends only on the local picture of orientations near the double point : when it passes from @xmath8{fo1.eps}}}$ ] through @xmath9{fo2.eps}}}$ ] to @xmath9{fo3.eps}}}$ ] , the linking number increases by @xmath6 .
generalities on vassiliev invariants see , e.g. , in @xcite .    in a sense
the linking number is the only vassiliev invariant of degree @xmath2 of two - component oriented links : any vassiliev invariant of degree @xmath2 of two - component oriented links is a linear function of the linking number .
similarly , the self - linking number is a vassiliev invariant of degree @xmath2 of framed knots ( it changes by @xmath6 when the knot experiences a homotopy with a generic appearance of a self - intersection point ) and it is the only vassiliev of degree @xmath2 of framed knots in the same sense . the necessity of a framing for the definition of self - linking number can now be formulated more rigorously : only constants are vassiliev invariants of degree @xmath2 of ( non - framed ) knots .
the diagrammatical definition of the writhe , which is imitated below , runs as follows : for each crossing point of the knot projection , one defines a _ local writhe _ equal to @xmath10 if near the point the knot diagram looks like @xmath9{fo3.eps}}}$ ] and @xmath1 if it looks like @xmath8{fo1.eps}}}$ ] .
then one sums the local writhes over all double points of the projection .
the sum is the writhe .
a continuous change of the projection may cause the vanishing of a crossing point .
this happens under the first reidemeister move shown in the left hand half of figure  [ f1 ] .
this move changes the writhe by @xmath11 .      if a link is algebraic , then its projection to a plane is algebraic , too .
a generic projection has only ordinary double points and the total number of its complex double points is constant .
the number of real double points can vary , but only by an even number .
a real double point can not turn alone into an imaginary one , as it seems to happen under the first reidemeister move . under an algebraic version of the first reidemeister move
, the double point stays in the real domain , but becomes solitary , like the only real point of the curve @xmath12 .
the algebraic version of the first reidemeister move is shown in the right hand half of figure  [ f1 ] .
it is not difficult to prove that the family of spatial curves that realizes this move can be transformed by a local diffeomorphism to the family of affine curves defined by the following system of equations @xmath13 where @xmath14 is the parameter of the deformation .
these are rational curves , admitting a rational parametrization @xmath15 the projection corresponds to the standard projection @xmath16 to the coordinate @xmath17-plane .
it maps these curves to the family of affine plane rational cubic curves defined by @xmath18 with @xmath19 .
a solitary double point of the projection is not the image of any real point of the link .
it is the image of two imaginary complex conjugate points of the complexification of the link .
the preimage of the point in the 3-space under the projection is a real line .
it is disjoint from the real part of the link , but intersects its complexification in a couple of complex conjugate imaginary points .    in the model of the first reidemeister move above
, @xmath20 is the double point of the projection for each @xmath21 .
if @xmath22 , it is a usual crossing point .
its preimage consists of two real points @xmath23 and @xmath24 .
if @xmath25 , it is a solitary double point .
its preimage consists of two imaginary conjugate points @xmath26 and @xmath27 , which lie on a real line @xmath28    below , in section [ s0.3 ] , with any solitary double point of the projection , a local writhe equal to @xmath29 is associated .
this is done in such a way that the local writhe of the crossing point vanishing in the first reidemeister move is equal to the local writhe of the new - born solitary double point . in the case of an algebraic knot ,
the sum of local writhes of all double points , both solitary and crossings , does not depend on the choice of projection and is invariant under rigid isotopy .
this sum is the encomplexed writhe .
a  construction similar to the construction of the encomplexed writhe number of an algebraic knot can be applied to an algebraic _ link_. however in this case there are two versions of the construction .    in the first of these
, we define an encomplexed writhe number generalizing the encomplexed writhe number defined above for knots .
we consider a link diagram and the sum of local writhes at solitary double points and crossing points where the branches belong the same connected component of the set of real points . at these crossing point , to define a local writhe , we need orientations of the branches . as
above , we choose an orientation on each of the components . if we make another choice , at a crossing point for which the branches belong the same component , either both orientations change or none .
hence the local writhe numbers at crossing points of this kind do not depend on the choice . in section [ s0.1 ]
below , we prove that the whole sum of local writhes over crossing points of this kind and solitary double points does not depend on the projection and is invariant under rigid isotopy .
we call this sum the _ encomplexed writhe number _ of the link @xmath30 and denote by @xmath31 .    in the second version of the construction , we consider a real algebraic link which is equipped with an orientation of the set of real points , use these orientations to define local writhe numbers at all crossing points and sum the local writhe numbers over all crossing points and all solitary double points .
the result is called the _
encomplexed writhe number of an oriented real algebraic link_. this encomplexed writhe number does not change when the orientation reverses .
an orientation considered up to reversing is called a _
semi - orientation_. thus the encomplexed writhe number depends only on the semi - orientation of the link .
the ( semi-)orientation may be an artificial extra structure , but it may also appear in a natural way , say , as a complex orientation , if the set of real points divides the set of real points , see @xcite .
in fact , the complex orientation is defined up to reversing , so it is indeed a semi - orientation .
another important class of semi - oriented algebraic links appears as transversal intersections of two real algebraic surfaces of degrees @xmath4 and @xmath32 with @xmath33 .
the encomplexed writhe number of ( semi-)oriented real algebraic link differs from the encomplexed writhe number of the same link without orientation by the sum of all pairwise linking numbers of the components multiplied by @xmath6 : let @xmath30 be a real algebraic link , let @xmath34 be the same link equipped with an orientation of its set of real points and @xmath35 the ( oriented ) connected components of this set , then @xmath36      in the case of a knot , the encomplexed writhe number defines a natural class of framings , since homotopy classes of framings are enumerated by their self - linking numbers and we can choose the framing having the self - linking number equal to the algebraic encomplexed writhe number .
i do not know any direct construction of this framing .
moreover , there seems to be a reason for the absence of such a construction . in the case of links
, the construction above gives a single number , while framings are enumerated by sequences of numbers with entries corresponding to components .
let @xmath30 be a nonsingular real algebraic curve in @xmath0-dimensional projective space .
then the set @xmath37 of its real points is a smooth closed @xmath2-dimensional submanifold of @xmath38 , i.e. , a smooth projective link .
the set @xmath39 of its complex points is a smooth complex @xmath2-dimensional submanifold of @xmath40 .
let @xmath41 be a point of @xmath38 .
consider the projection @xmath42 from @xmath41 .
assume that @xmath41 is such that the restriction to @xmath39 of @xmath43 is generic .
this means that it is an immersion without triple points and at each double point the images of the branches have distinct tangent lines .
it follows from well - known theorems that those @xmath41 s for which this is the case form an open dense subset of @xmath38 ( in fact , it is the complement of a @xmath6-dimensional subvariety ) .
the real part @xmath44 of the image consists of the image @xmath45 of the real part and , maybe , several solitary points , which are double points of @xmath46 .
there is a purely topological construction which assigns a local writhe equal to @xmath29 to a crossing belonging to the image of only one component of @xmath37 .
this construction is well - known in the case of classical knots .
here is its projective version .
i borrow it from drobotukhina s paper @xcite on the generalization of kauffman brackets to links in projective space .    , @xmath47 , @xmath48 .
]    let @xmath49 be a smooth connected one - dimensional submanifold of @xmath38 , and @xmath41 be a point of @xmath50 .
let @xmath51 be a generic double point of the projection @xmath52 and @xmath53 be the line which is the preimage of @xmath51 under the projection .
denote by @xmath54 and @xmath55 the points of @xmath56 .
the points @xmath54 and @xmath55 divide the line @xmath57 into two segments .
choose one of them and denote it by @xmath58 .
choose an orientation of @xmath49 .
let @xmath59 and @xmath60 be tangent vectors of @xmath49 at @xmath54 and @xmath55 respectively directed along the selected orientation of @xmath49 .
let @xmath47 be a vector tangent to @xmath57 at @xmath54 and directed inside @xmath58 .
let @xmath48 be a vector at @xmath54 such that it is tangent to the plane containing @xmath57 and @xmath60 and is directed to the same side of @xmath58 as @xmath60 ( in an affine part of the plane containing @xmath58 and @xmath60 ) .
see figure  [ f2 ] .
the triple @xmath59 , @xmath47 , @xmath48 is a base of the tangent space @xmath61 .
define the local writhe of @xmath51 to be the value taken by the orientation of @xmath38 on this frame .
the construction of the local writhe of @xmath51 contains several choices .
here is a proof that the result does not depend on them .
we have chosen an orientation of @xmath49 . had the opposite orientation been selected , then @xmath59 and @xmath48 would be replaced by the opposite vectors @xmath62 and @xmath63 .
this would not change the result , since @xmath62 , @xmath47 , @xmath63 defines the same orientation as @xmath59 , @xmath47 , @xmath48 .
we have chosen the segment @xmath58 .
if the other half of @xmath57 was selected , then @xmath47 and @xmath48 would be replaced by the opposite vectors .
but @xmath59 , @xmath64 , @xmath63 defines the same orientation as @xmath59 , @xmath47 , @xmath48 .
the construction depends on the order of points @xmath54 and @xmath55 .
the other choice ( with the same choice of the orientation of @xmath49 and segment @xmath58 ) gives a triple of vectors at @xmath55 .
it can be moved continuously without degeneration along @xmath58 into the triple @xmath48 , @xmath64 , @xmath59 , which defines the same orientation as @xmath59 , @xmath47 , @xmath48 .
let @xmath30 , @xmath41 , and @xmath43 be as in the beginning of section [ s0.1 ] and let @xmath65 be a solitary double point of @xmath43 .
here is a construction assigning @xmath29 to @xmath66 .
i will also call the result a _ local writhe _ of @xmath66 .
denote the preimage of @xmath66 under @xmath43 by @xmath57 .
this is a real line in @xmath38 connecting @xmath41 and @xmath66 .
it intersects @xmath39 in two imaginary complex conjugate points , say , @xmath54 and @xmath55 . since @xmath54 and @xmath55 are conjugate , they belong to different components of @xmath67
.    choose one of the common points of @xmath39 and @xmath68 , say , @xmath54 .
the natural orientation of the component of @xmath67 defined by the complex structure of @xmath68 induces an orientation on @xmath69 as on the boundary of its closure .
the image under @xmath43 of the local branch of @xmath39 passing through @xmath54 intersects the plane of the projection @xmath70 transversally at @xmath66 .
take the local orientation of the plane of projection such that the local intersection number of the plane and the image of the branch of @xmath39 is  @xmath10 .
thus the choice of one of two points of @xmath71 defines an orientation of @xmath69 and a local orientation of the plane of projection @xmath70 ( we can speak only of a local orientation of @xmath70 , since the whole @xmath70 is not orientable ) . the plane of projection intersects .
if you would like to think on it as on the set of lines of @xmath38 passing through @xmath41 , please , identify it in a natural way with any real projective plane contained in @xmath38 and disjoint from @xmath41 .
all such embeddings @xmath72 are isotopic . ]
transversally @xmath69 in @xmath66 .
the local orientation of the plane , the orientation of @xmath69 and the orientation of the ambient @xmath38 determine the intersection number .
this is the local writhe .
it does not depend on the choice of @xmath54 .
indeed , if one chooses @xmath55 instead , then both the orientation of @xmath69 and the local orientation of @xmath70 would be reversed .
the orientation of @xmath73 would be reversed , because @xmath69 inherits opposite orientations from the different halves of @xmath67 .
the local orientation of @xmath70 would be reversed , because the complex conjugation involution @xmath74 preserves the complex orientation of @xmath75 , preserves @xmath70 ( point - wise ) and maps one of the branches of @xmath46 at @xmath66 to the other reversing its complex orientation .
now for any real algebraic projective link @xmath30 , choose a point @xmath76 such that the projection of @xmath30 from @xmath41 is generic and sum the writhes of all crossing points of the projection belonging to the image of only one component of @xmath37 and the writhes of all solitary double points .
this sum is called the _
encomplexed writhe number of @xmath30_.    i have to show that it does not depend on the choice of projection . the proof given below proves more : the sum is invariant under _ rigid isotopy _ of @xmath30 . by rigid isotopy we mean an isotopy consisting of nonsingular real algebraic curves .
the effect of a movement of @xmath41 on the projection can be achieved by a rigid isotopy defined by a path in the group of projective transformations of @xmath38 .
therefore the following theorem implies both the independence of the encomplexed writhe number from the choice of projection and its invariance under rigid isotopy .
[ mainth ] for any two rigidly isotopic real algebraic projective links @xmath77 and @xmath78 whose projections from the same point @xmath76 are generic , the encomplexed writhe numbers of @xmath77 and @xmath78 defined via @xmath41 are equal .
this theorem is proved in section [ s0.6 ] .
[ cor1 ] the encomplexed writhe number of a real algebraic projective link does not depend on the choice of the projection involved in its definition .
a projection depends only on the center from which it is done .
the effect on the projection of a movement of the center can be achieved by a rigid isotopy defined by a path in the group of projective transformations of @xmath38 .
thus the encomplexed writhe number is a characteristic of a real algebraic link .
[ cor2 ] the encomplexed writhe number of a real algebraic projective link is invariant under rigid isotopy .          as in the purely topological situation of an isotopy of a classical link , a generic rigid isotopy of a real algebraic link may be decomposed into a composition of rigid isotopies , each of which involves a single local standard move of the projection .
there are @xmath79 local standard moves .
they are similar to the reidemeister moves .
the first of these @xmath79 moves is shown in the right hand half of figure  [ f1 ] .
the other moves are shown in figure [ o23 ] .
the first two of these coincide with the second and third reidemeister moves .
the fourth move is similar to the second reidemeister move : also two double points of projection come to each other and disappear .
however the double points are solitary .
the fifth move is similar to the third reidemeister move : a triple point also appears for a moment .
but at this triple point only one branch is real , the other two are imaginary conjugate to each other . in this move a solitary double point traverses a real branch .      to prove theorem [ mainth ] ,
first replace the rigid isotopy by a generic one and then decompose the latter into local moves described above , in section [ s0.5 ] .
only in the first , fourth and fifth moves solitary double points are involved .
the invariance under the second and the third move follows from the well - known fact of knot theory that the topological writhe is invariant under the second and third reidemeister moves .
thus the following three lemmas imply theorem [ mainth ] .
[ lem1 ] in the fifth move the writhe of the solitary point does not change .
[ lem2 ] in the fourth move the writhes of the vanishing solitary points are opposite .
[ lem3 ] in the first move the writhe of vanishing crossing point is equal to the writhe of the new - born solitary point .
proof of lemma [ lem3 ] is postponed to section [ s0.9 ] .
note that although lemma [ lem3 ] is the most difficult to prove , it is the least significant : here its only role is to justify the choice of sign made in the definition of local writhe in solitary double point of the projection .
it is clear that the writhes of vanishing double points involved in the first move are related , and if they were opposite to each other , then the definition of the encomplexed writhe number should be changed , but would not be destroyed irrecoverably .
this is obvious .
indeed , the real branch of the projection does not interact with the imaginary branches , it just passes through their intersection point .
at the moment of the fourth move take a small ball @xmath80 in the complex projective plane centered in the solitary self - tangency point of the projection of the curve .
its intersection with the projection of the complex point set of the curve consists of two smoothly embedded disks tangent to each other and to the disk @xmath81 . under the move each of the disks experiences a diffeotopy . before and after the move the intersection the curve with @xmath80 is the union of the two disks meeting each other transversally in two points , but
before the move the disks do not intersect @xmath70 , while after the move they intersect @xmath70 in their common points .    to calculate the writhe at both vanishing solitary double points ,
let us select the same imaginary branch of the projection of the curve passing through the points .
this means that we select one of the disks described above .
the sum of the local intersection numbers of this disk ( equipped with the complex orientation ) and @xmath81 ( equipped with some orientation ) is zero since under the fourth move the intersection disappears , while in the boundary of @xmath80 no intersection happens .
therefore the local orientations of the projective plane in the vanishing solitary double points defined by this branch define opposite orientations of @xmath81 .
( recall that the local orientations are distinguished by the condition that the local intersection numbers are positive . )    on the other hand , under the move the preimages of the vanishing solitary double points come to each other up to coincidence at the moment of the move and their orientations defined by the choice of the same imaginary branch are carried to the same orientation of the preimage of the point of solitary self - tangency .
indeed , the preimages are real lines and points of intersection of their complexifications with the selected imaginary branch of the curve also come to the same position .
therefore the halves of the complexifications containing the points come to coincidence , as well as the orientations defined by the halves on the real lines .
it follows that the intersection numbers of @xmath80 with the preimages of the vanishing solitary double points equipped with these orientations are equal .
since the local orientations of the projective plane in the vanishing solitary double points define distinct orientations of @xmath81 , the writhes are opposite to each other .
it is sufficient to consider the model family of curves described in section [ si3 ] .
recall that the curves of this family are defined by the following system of equations @xmath13 where @xmath14 is the parameter of the deformation .
these curves admit a rational parametrization @xmath82 the projection corresponds to the standard projection @xmath16 to the coordinate @xmath17-plane .
it maps these curves to the family of affine plane rational cubic curves defined by @xmath18 with @xmath19 .
we must prove that the local writhe at @xmath20 for @xmath22 coincides with the local writhe at @xmath20 for @xmath25 .
let us calculate the local writhe for @xmath22 .
denote @xmath83 by @xmath84 .
the preimage of @xmath20 consists of points @xmath85 and @xmath86 corresponding to the values @xmath87 and @xmath84 of @xmath88 , respectively , see figure [ f3 ] .
the tangent vectors to the curve at these points are @xmath89 and @xmath90 . the vector @xmath47 connecting @xmath54 and @xmath55 is @xmath91 . by definition ,
the writhe is the value taken by the orientation of @xmath92 on the frame @xmath59 , @xmath47 , @xmath48 .
this value is equal to the value of this orientation on the frame @xmath93 , @xmath94 , @xmath95 multiplied by the sign of @xmath96    let us calculate the local writhe for @xmath25 .
denote @xmath97 by @xmath84 .
the preimage of @xmath20 consists of points @xmath98 and @xmath99 corresponding to the values @xmath100 and @xmath101 of @xmath88 .
choose the branch which passes through @xmath102 .
it belongs to the upper half of the line @xmath28 , which induces the positive orientation of the real part directed along @xmath95 .
at @xmath102 the branch of the curve has tangent vector @xmath103 and the real basis consisting of @xmath59 and @xmath104 positively oriented with respect to the complex orientation of this branch .
the projection maps this basis to the positively oriented basis @xmath105 , @xmath106 of the projection of the branch .
the intersection number of this projection and @xmath107 in @xmath108 is the sign of @xmath109 hence the orientation of @xmath107 such that its local intersection number with the selected branch of the projection does not coincide with the orientation defined by the standard basis .
the intersection number of the line @xmath28 with the standard orientation and the @xmath17-plane with the standard orientation is the value of the orientation of the ambient space @xmath110 taken on the standard basis @xmath93 , @xmath94 , @xmath95 .
therefore the local writhe is opposite to this value .
[ rem ] there is a more conceptual proof of lemma [ lem3 ] . it is based on a local version of the rokhlin complex orientation formula , see @xcite and @xcite .
in fact , the original proof was done in that way . however , the complex orientation formula is more complicated than the calculation above .      to speak about vassiliev invariants , we need to fix a connected family of curves , in which links under consideration comprise the complement to a hypersurface . in the case of classical knots
one could include all knots in such a family by adjoining knots with self - intersections and other singularities .
a singular knot is a right equivalence class of a smooth map of the circle to the space ( recall that two maps from a circle are right equivalent if one of them is a composition of a self - diffeomorphism of the circle with the other one ) .    in the case of real algebraic knots ,
such a family including all real algebraic knots does not exist .
even the space of complex curves in the three - dimensional projective space consists of infinitely many components .
it is impossible to change the homology class realized by the set of complex points of an algebraic curve in @xmath40 by a continuous deformation .
recall that the homology class belongs to the group @xmath111 and is a positive multiple @xmath112 $ ] of the natural generator of @xmath113\in h_2(\c p^3)$ ] realized by a line .
the coefficient @xmath114 is called the _ order _ of the curve .
the genus is another numerical characteristic of a complex curve which takes the same value for all nonsingular curves in any irreducible family .
as is well known , the nonsingular complex curves of given order and genus in three - dimensional projective space are parametrized by a finite union of quasi - projective varieties . for each of these varieties , one can try to build a separate theory of vassiliev invariants on a class of nonsingular real algebraic curves whose complexifications are parametrized by points of this variety .
( a similar phenomenon takes place in topology : links with different numbers of components can not be included into a reasonable connected family , and therefore for each number of components there is a separate theory of vassiliev invariants . )    among the varieties of algebraic curves in three - dimensional projective space , there are two special families : for each natural number @xmath114 there is an irreducible variety of rational curves of order @xmath114 ( recall that a an algebraic curve is called rational if it admits an algebraic parametrization by a line ) , and for each pair of natural numbers @xmath4 and @xmath32 there is an irreducible variety of curves which can be presented as intersection of surfaces of degrees @xmath4 and @xmath32 .    in the class of real algebraic rational curves of order @xmath114 , singular curves comprise a _ discriminant hypersurface _ in which a generic point is a rational curve such that it has exactly one singular point and this point is an ordinary double point .
an ordinary double point may be of one of the following two types : either it is an intersection point of two real branches , or two imaginary conjugate branches .
any two real algebraic rational nonsingular curves of order @xmath114 can be connected by a path in the space of real rational curves of degree @xmath114 that intersects the discriminant hypersurface only transversally at a finite number of generic points .
such a path can be regarded as a deformation of a curve to the other one .
when it intersects the discriminant hypersurface at a point , which is a curve with singularity on real branches , the set of real points of the curve behaves as in classical knot theory : two pieces of the set of real points come to each other and pass through each other . as in classical knot theory , at the moment of intersection , the generic projection of the curve experiences an isotopy .
nothing happens besides that one crossing point becomes for a moment the image of a double point and then changes back into a crossing point , but with the opposite writhe .
when the path intersects the discriminant hypersurface at a point , which is a curve with singularity on imaginary branches , two complex conjugate imaginary branches pass through each other .
at the moment of passing , they intersect in a real isolated double point . at this moment
the set of real points of a generic projection experiences an isotopy .
no event happens besides that a solitary double point becomes for a moment the image of a solitary real double point of the curve and then changes back into an ordinary solitary double point of the projection ( which is not the image of a real point of the knot ) , but with the opposite writhe number .
it is clear that the encomplexed writhe number of an algebraic curve changes under a modification of each of these kinds by @xmath115 , with the sign depending only on the local structure of the modification near the double point .
this means that _ the encomplexed writhe number on the family of real rational curves under consideration is a vassiliev invariant of degree @xmath2_.    this is true also for any space of nonsingular real algebraic curves that can be included into a connected family of real algebraic curves by adjoining a hypersurface , penetration through which at a generic point looks as in the family of rational curves described above .
there are many families of this kind besides the families of rational knots .
however , in many families of algebraic curves a transversal penetration through the discriminant hypersurface in a generic point looks differently . in particular , for intersections of two surfaces it is a morse modification of the real part of the curve . at the moment , the old double points of the projection , both solitary and crossing , do not change .
an additional double point appears just for a moment .
however the division of crossing points to self - crossing points of a single component and crossing points of different components may change .
therefore the encomplexed writhe number changes in a complicated way .
if the degrees of the surfaces defining the curve are of the same parity , the real part of the curve has a natural semi - orientation .
the morse modification respects this semi - orientation .
therefore the encomplexed writhe number of the semi - oriented curve does not change .
[ vanishsl ] the encomplexed writhe number of any nonsingular semi - oriented real algebraic link which is a transversal intersection of two real algebraic surfaces whose degrees are of the same parity is zero .
any two nonsingular real curves of the type under consideration can be connected by a path as above . hence their self - linking numbers coincide . on the other hand , it is easy to construct , for any pair of natural numbers @xmath4 and @xmath32 of the same parity , a pair of nonsingular real algebraic surfaces of degrees @xmath4 and @xmath32 transversal to each other in three - dimensional projective space such that their intersection has zero self - linking number .
in contrast to this vanishing result , one can prove that the _ encomplexed writhe number of a real algebraic rational knots of degree @xmath114 can take any value in the interval between @xmath116 and @xmath117 including these limits and congruent to them modulo @xmath6_.
the same construction may be applied to real algebraic curves in @xmath38 having singular imaginary points , but no real singularities . in the construction
we can eliminate projections from the points such that some singular point is projected from them to a real point .
indeed , for any imaginary point there exists only one real line passing through it ( the line connecting the point with its complex conjugate ) , thus we have to exclude a finite number of real lines .
this gives a generalization of encomplexed writhe numbers with the same properties : it is invariant with respect to rigid isotopies ( i.e. , isotopies made of curves from this class ) , and is multiplied by @xmath1 under a mirror reflection .      the construction of this paper can be applied to algebraic links in the sphere @xmath3 .
although from the viewpoint of knot theory this is the most classical case , from the viewpoint of algebraic geometry the case of curves in the projective space is simpler .
the three - dimensional sphere @xmath3 is a real algebraic variety .
it is a quadric in four - dimensional real affine space .
the stereographic projection is a birational isomorphism of @xmath3 onto @xmath38 .
it defines a diffeomorphism between the complement of the center of the projection in @xmath3 and a real affine space .
given a real algebraic link in @xmath3 , one may choose a real point of @xmath3 from the complement of the link and project the link from this point to an affine space .
then include the affine space into the projective space and apply the construction above .
the image has no real singular points , therefore we can use the result of the previous section .
this construction blows up the center of projection , making a real projective plane out of it , and maps the complement to the center of the projection in the set of real points of the sphere isomorphically onto the complement of the projective plane . in the imaginary domain
, it contracts each generatrix of the cone which is the intersection of the sphere with its tangent plane at the center of projection .
the image of the cone is an imaginary quadric curve contained in the projective plane which appeared as the result of blowing up of the central point .
first , consider the most straightforward generalization .
let @xmath57 be a nonsingular real algebraic @xmath118-dimensional subvariety in the projective space of dimension @xmath119 .
its generic projection to @xmath120 has only ordinary double points . at each double point
either both branches of the image are real or they are imaginary complex conjugate . if the set of real points is orientable , then one can repeat everything with obvious changes and obtain a definition of a numeric invariant generalizing the encomplexed writhe number defined above .
let @xmath121 be a nonsingular three - dimensional real algebraic variety with oriented set of real points equipped with a real algebraic fibration over a real algebraic surface @xmath122 with fiber a projective line .
there is a construction which assigns to a real algebraic link ( i.e. , a nonsingular real algebraic curve in @xmath121 ) with a generic projection to @xmath122 an integer , which is invariant under rigid isotopy , is multiplied by @xmath1 under the orientation reversal in @xmath121 and is a vassiliev invariant of degree @xmath2 .
this construction is similar to the one presented above , but uses , instead of the projection to @xmath70 , an algebraic version of turaev s shadow descriptions of links @xcite .
here we discuss only one example .
however it can be easily generalized .
consider immersions of the sphere @xmath123 to @xmath124 .
up to regular homotopy ( i.e. , a homotopy consisting of immersions whose differentials also comprise a homotopy ) , an immersion @xmath125 is defined by its smale invariant @xcite , which is an element of @xmath126 . for a generic immersion
, it can be expressed as the sum of local self - intersection numbers over all double points of the immersion , see @xcite .
let us encomplex the smale invariant .
for this , first , we have to consider a real algebraic counterpart for the notion of generic immersion @xmath125 .
the identification is defined via the universal covering @xmath127 .
replace euclidean space @xmath124 by torus @xmath128 , which has the advantage of being compact .
the classification of immersions @xmath129 up to regular homotopy coincides with the smale classification of immersions @xmath125 .
the sphere @xmath123 is the real part of a quadric projective hypersurface .
the torus @xmath128 is the real part of a complex abelian variety .
consider real regular maps of the quadric to the abelian variety .
a generic map defines an immersion both for the complex and real parts .
the only singularities are transversal double points .
double points in the real part of the target variety are of two kinds . at a double point of the first kind
two sheets of the image of @xmath123 meet . at a double point of the second kind
the images two complex conjugate sheets of the complexification of @xmath123 meet .
the smale invariant is the sum of the local intersection numbers over the double points of the first kind .
one can extend the definition of the local intersection number to the double points of the second kind in such a way that the total sum of the local intersection numbers over double points of both kinds would be invariant under continuous deformations of regular maps .
this total sum is the _
encomplexed smale invariant_. notice that it is , in a sense , more invariant than the original smale invariant .
the smale invariant may change under homotopy , it is invariant only under regular homotopy .
the encomplexed smale invariant does not change under a homotopy in the class of regular maps , which corresponds to the class of all continuous maps .
drobotukhina , _ an analogue of the jones polynomial for links in @xmath130 and a generalization of the kauffman  murasugi theorem _
, algebra i analiz * 2 * ( 1990 ) , no .  3 , 171191 ; english transl .
leningrad math .
j. * 2 * ( 1991 ) , no .  3 , 613630 .                o.  ya .
viro , _ progress in the topology of real algebraic varieties over the last six years _ , uspekhi mat .
nauk * 41 * ( 1986 ) , no .  3 , 4567 ; english transl
russian math .
surveys * 41 * ( 1986 ) , no .  3 , 5582 .


gamma - ray bursts are the most powerful explosions in the universe .
( for reviews see piran 1999 and mszros 2002 . ) a crucial advance in understanding gamma - ray bursts began with the discovery of `` afterglows '' , starting with _ bepposax _ observations in the soft x - ray band of grb 970228 ( costa et al .
1997 ; wijers , rees , & mszros 1997 ) .
( for a review of grb afterglows see van paradijs , kouveliotou , & wijers 2000 . )
if grbs were isotropic , then the measured redshifts would imply total explosion energies of @xmath5 ergs ( frail et al .
. theoretical work on relativistic jet expansion , however , shows that one expects a steepening in the decay light curve if one is looking down the axis of a jet as the flow decelerates from a bulk lorentz factor @xmath6 to @xmath7 , where @xmath1 is the jet beaming angle ( e.g. , rhoads 1997 , sari , piran , & halpern 1999 , panaitescu & mszros 1999 , panaitescu & kumar 2001ab , 2002 ) .
one does in fact see such steepenings in the light curves ( e.g. , stanek et al .
1999 ; harrison et al .
prior to the time when @xmath8 the expansion is effectively `` spherical '' from the observer s viewpoint because the relativistic beaming is narrower than the jet itself . in other words
, if the grb emission were coming from one spot on a large , relativistically expanding sphere , aimed directly at the observer , the observer would not see any emission from the other parts of the sphere .
after the time when @xmath8 the observer can `` see '' the entire jet , and a faster rate of decline in the luminosity is predicted .
a separate issue that we will address in this work is the sideways or lateral expansion of the jet as the increasing solid angle of the jet enables a larger fraction of the circumstellar medium ( csm ) surrounding the progenitor star to be intercepted and provide decelerating gas .
several groups have claimed that this leads to a faster ( exponential ) decrease in @xmath2 , which acts as an additional agent to diminish the amplitude of the relativistic beaming .
the concept of a `` break '' in the afterglow light curves occurring when @xmath8 has been used to infer the presence of strong beaming in grbs ( frail et al .
2001 ; panaitescu & kumar 2001ab , 2002 ; panaitescu , mszros , & rees 1998 ; piran , kumar , panaitescu , & piro 2001 ) .
frail et al .
( 2001 ) utilize the theoretical framework of sari , piran , & halpern ( 1999 ) which takes the jet evolution to be spherical adiabatic expansion to show that , after correcting the `` isotropic '' energies to account for the specific beaming factor for each burst , the total burst energy reduces to a narrow range centered on @xmath9 ergs .
frail et al .
( 2001 ) give a tabulation of 17 afterglows for which redshifts had been measured , up to 2001 january .
the redshifts range from 0.433 to 4.5 . of these
, 10 also have break times known to within @xmath10% , ranging from 1 d to 25 d. ( in addition , three grbs have listed lower limits on the break time , and two have listed upper limits . )
the combination of redshifts , fluences , and break times lead to estimates of the jet angles @xmath1 ranging from @xmath11 to @xmath12 , with an aggregation near @xmath13 .
( the number of grb s with redshift determinations is currently @xmath14 . for an update of the work described in frail et al .
2001 , see berger , kulkarni , & frail 2003 . )
prior investigations of grb jet expansion have been largely analytical , or involved computational models with some imposed symmetry , typically spherical or axial ( e.g. , rhoads 1997 , 1999 ; sari 1997 , sari , piran , & halpern 1999 ; panaitescu & mszros 1999 ; granot , piran , & sari 1999a ; kumar & panaitescu 2000ab ; panaitescu & kumar 2001ab , 2002 ; granot et al .
these works have tended for the most part to confirm estimates and scalings based on the analytical formalism of blandford & mckee ( 1976=bm76 ) , although for many of the studies the agreement is ( to some extent ) circular , given that they rely at least in part on the bm76 formalism .
kobayashi , piran , & sari ( 1999 ) present results based on a spherically symmetric relativistic lagrangian code .
they identify three regimes during the evolution of the grb jet : ( i ) an acceleration phase during which @xmath15 and @xmath2 increases rapidly to a large value @xmath16 , ( ii ) a coasting phase during which @xmath2 is relatively constant as the mass accumulated from the csm is small compared to that in the jet , and ( iii ) a deceleration phase during which the accumulated mass forces a rapid decrease in @xmath2 .
previous studies using analytical methods divided the evolution into regimes defined by some ordering of distance , velocity , or energy scales ( see table 3 of piran 1999 and figure 7 of mszros 2002 ) .
the main findings that seem to be common to all studies are that the break in the decaying light curve @xmath17 occurs roughly when the deceleration has decreased @xmath2 to roughly the reciprocal of the jet beaming angle , and that the subsequent decrease in @xmath2 is roughly exponential with distance . also , panaitescu & mszros ( 1999 ) present axisymmetric calculations to study the combined effects of the transition from @xmath6 to @xmath7 and the lateral jet expansion .    in other subdisciplines of astronomy
the use of relativistic hydrodynamics codes has been standard for some time .
for example , workers studying extragalactic jets have used such codes to continuously inject a collimated supersonic beam into a surrounding medium , usually under the assumption of pressure equilibrium .
( see , e.g. , norman et al .
1982 for a thorough discussion ) . in the context of grbs
, work has been done using 2d and 3d relativistic hydro codes to consider the evolution of the grb jet as it propagates through the envelope of the progenitor star , up to the point where it breaks out of the stellar surface and produces the prompt grb emission ( zhang , woosley , & macfadyen 2003 , zhang , woosley , & heger 2003 ) . in this work
we consider the evolution covering the afterglow time ( i.e , after the period considered by zhang et al .
we utilize a three dimensional relativistic hydrodynamical code to study the propagation of an initially ultrarelativistic blob into a dense csm .
we study the spatial spreading of the blob both along the direction of propagation and orthogonal to it , as well as the evolution of @xmath2 in space and time .
we also calculate afterglow light curves , taking a simple prescription in which the local emissivity scales with the local pressure .
the model used is that of del zanna & bucciantini ( 2002=db02 ) .
these authors present a simple and efficient numerical scheme for special relativistic hydrodynamics ( srhd ) that does not rely on computationally expensive spectral decomposition ( and the accompanying matrix inversions ) .
the semi - discrete form of the srhd equations is solved , so that time integration can be carried out using a standard runge - kutta method .
unlike many previous multidimensional calculations run on supercomputers , characteristic decomposition and riemann solvers are not required .
the only local pieces of information required are the highest characteristic speeds .
the simplicity and efficiency of the model allow for runs to be carried out on a modern pc with @xmath18 ghz cpu speed and @xmath19 gb ram .
we refer the reader to db02 and to londrillo & del zanna ( 2000 ) for the details .
( also , del zanna & bucciantini 2003 add mhd to the srhd formalism of db02 . )
the basic idea is to calculate a vector of `` fluxes '' * f * that is then used within each time step to advance a vector of `` conserved '' variables * u * by computing the spatial divergences of the * f * s .
the time integration is carried out with a third order runge - kutta scheme .
the vector of `` primitive '' or physical variables @xmath20^t$ ] must be recovered from the conserved variables @xmath21^t$ ] in each runge - kutta step , and at each grid point .
primitive variables are reconstructed at the left ( l ) and right ( r ) cell faces using the `` convex essentially non - oscillatory '' interpolation described in db02 .
this method gives third order accuracy by adaptively adjusting the `` stencil '' , or spatial interval , used for computing differences in the vicinity of a given cell .
db02 present two versions of the fluxes , * f*@xmath22 , where `` hll '' denotes harten , lax , & van leer ( 1983 ) , that are precise but at the same time potentially prone to numerical instabilities , and * f*@xmath23 , where `` llf '' denotes local lax - friedrichs ( cf . lax & liu 1998 ) , that represent a smoother , numerically dissipative flux . for the main application we consider in this work we will utilize the llf fluxes .
db02 present a suite of test results for 1 , 2 , and 3 dimensions .
we have reproduced these tests , and now show the results of a 3d spherical expansion that results in higher @xmath2 values than the tests discussed in db02 .
we utilize a cubical grid of @xmath24 points , with each side having length unity .
minutes per runge - kutta integration for a 3d problem with @xmath25 nodes , utilizing the third order reconstruction for the l and r values of the primitives , and using a 1 ghz linux pc ; our run time of @xmath26 s is consistent with the 1.7 ghz linux pc used in this work . ]
initial conditions are that everywhere @xmath27 , @xmath28 , and @xmath29 within a radius 0.1 of the center of the cube , and @xmath30 outside .
this is the 3d analog of the 1d `` piston '' problem .
( an identical test with less extreme initial conditions is shown and discussed by hughes , miller , & duncan 2002 , see their fig .
3 . ) the enormous over - pressure launches a relativistic blast wave at @xmath31 .
figure 1 shows the evolution of @xmath32 , @xmath33 , @xmath34 , and @xmath2 for a slice taken along the @xmath35axis .
the density inside the sphere becomes small for late times . to avoid numerical instabilities in this region
we utilize the smoother llf fluxes for small densities ( @xmath36 ) , and the hll fluxes elsewhere .
we show the initial configuration plus 10 time slices from a run with 360 total time steps , taking a courant number of @xmath37 . after the run begins one can see the development of a strong spherical outflow from the center of the cube .
the lorentz @xmath2 factor has increased to @xmath38 in the outer parts of the expanding shell by the end of the run .
the run was halted before the expansion reached the edge of the grid . for the last three times steps , a kink develops in the @xmath33 profile at the center of the sphere because of the prolonged spherical outflow from localized point .
figure 2 shows the corresponding evolution for pure llf fluxes .
all variables now show a smoother evolution and @xmath39 at large radii for late times @xmath40 somewhat less than for the hll fluxes .
conservation of rest - mass energy @xmath41 and energy @xmath42 ( where the specific energy @xmath43 ) integrated over the computational grid is as follows : by the end of the run shown in fig . 1 , @xmath44 and @xmath45
have increased by about 1 part in @xmath46 ; by the end of the run shown in fig .
2 , @xmath44 has increased by @xmath473 parts in @xmath48 , and @xmath45 is constant ( to within machine accuracy ) .
we now present results for the case of interest , namely a relativistic blob of material that expands roughly axisymmetrically ( orthogonal to the direction of motion ) as it decelerates from ultrarelativistic to relativistic speeds .
we set up our initial problem as follows : our blob is made to propagate along the positive @xmath35axis through a 3d rectangular grid of @xmath49 cells .
the box measures 5 units along the @xmath50 axis , and 1 unit each along the @xmath51 and @xmath52 axes .
therefore the grid spacing along each axis is 0.01 units .
we explore a variety of configurations for our initial state .
we performed runs in which the blob initially is either a small cone , sphere , or plate , with symmetry about the propagation vector .
the dimensions of the initial blob are optimized , according to its shape , so that @xmath26 grid points out of the @xmath53 in the computational domain initially comprise the blob .
an initial spread is imparted to each fluid element by setting the @xmath51 and @xmath52 components of velocity such that @xmath54 .
the motivation for taking @xmath55 is the study of frail et al .
( 2001 ) that finds a peak in the frequency histogram distribution of inferred spreading angles @xmath1 from observed grbs of @xmath13 or 0.07 radian .
we divide this value by two to obtain the half - angle spread @xmath56 radian .
the blob is given a lorentz factor @xmath57 so that the relativistic beaming angle @xmath3 lies within the physical spreading angle @xmath1 of the jet .
we set @xmath58 ( inside the blob ) , @xmath59 ( outside the blob ) , and @xmath60 .
we also perform one trial with @xmath61 , i.e. , @xmath62 at all points .
this is the `` null hypothesis '' run . the effective density contrast perceived in the lab ( i.e. , csm ) frame is @xmath63
. soon after the computation starts , the initial profile relaxes to one in which @xmath64 , therefore the effective density contrast for the early run is @xmath65 .
this is less than expected astrophysically , but required in our computations in order to see significant deceleration of the blob by the time it reaches the end of the grid @xmath66 .
in other words , we must telescope the evolution into the finite dimensions of our grid . from conservation of momentum
, the condition for significant deceleration is that the total swept - up mass - energy roughly equals that in the initial blob .
experimentation using hll fluxes shows that strong internal shocks almost immediately create large - amplitude sawtooth @xmath32 variations within the blob , leading to noisy results .
therefore we utilize the smoother llf fluxes in this work .
we follow the evolution of the blob in terms of both its motion in @xmath50 and its spreading in @xmath51 and @xmath52 . taking a courant number of @xmath67 necessitates @xmath68 time steps for the blob to reach the end of the grid , traveling at @xmath69 .
this is a simple consequence of the box length and grid spacing @xmath70 , from which it follows that the time to traverse the box is @xmath71 units traveling at @xmath72 .
a density contrast of @xmath73 is sufficient to see the desired deceleration from @xmath74 to @xmath75 during a run .    in order to avoid undue complexity in these experiments , that currently are purely hydrodynamic and do not yet contain proper prescriptions for emission from bremsstrahlung and synchrotron processes , we calculate a simple measure of the emissivity by taking the local emission to scale as @xmath34 , which would be expected roughly for optically thin synchrotron emission characteristic of frequencies significantly above the self - absorption frequency .
we avoid the issues of synchrotron self - absorption ( granot , piran , & sari 1999b , 2000 ; panaitescu & kumar 2000 ; sari & esin 2001 ) , and of whether the evolution is adiabatic or radiative ( panaitescu & mszros 1998b ; sari , piran , & narayan 1998 ) @xmath40 our calculations are adiabatic .
we calculate effective `` light curves '' for observers at various viewing angles between @xmath76 and @xmath77 from the center of the jet .
the amplification of the photon energy flux is dictated by three factors : ( 1 ) the rate of photon emission is increased by a factor of @xmath2 as one goes from the blob frame to the lab frame .
( 2 ) the individual photons are doppler - shifted so that their lab frame energy is increased by a factor of @xmath78 , where @xmath79 is the angle between the line of sight and blob s direction of motion .
( 3 ) the photons are focused in the direction of motion , so that an angle of emission ( relative to the direction of motion ) in the blob s frame , @xmath80 , is related to the lab frame angle @xmath79 by @xmath81 .
the differential solid angle @xmath82 $ ] .
thus , the net amplification of the photon energy flux for each fluid element is @xmath83 , similar expressions have been derived previously ( e.g. , granot , piran , & sari 1999a , eqns . [ 3 ] , [ 4 ] ; wang , & loeb 2001 , eqns . [ 13 ] , [ 31 ] ) .
the formalisms of granot et al ( 1999a ) and wang & loeb ( 2001 ) also contain a provision for time delays based on the relative positions and velocities of different fluid elements .
the `` time '' of interest in summing the relative emission contributions is the retarded time @xmath84 .
we compute light curves by keeping running totals of the contributions to the flux from 100 slices through the computational domain ( straddling the position of local disturbance in the fluid caused by the advancing blob ) that are orthogonal to the blob velocity vector ( i.e. , cuts in the @xmath85 plane ) and advance in time toward the observer at @xmath86 . over the course of a run
we build up a light curve based on constant retarded times with respect to the observer .
this approach is essential insofar as the photons emitted by a given fluid element outpace the blob itself only slightly .
figure 3 shows the evolution of physical quantities for the run which begins with a small cone along a cut through the center of the computational box , along the direction of propagation of the blob .
the initial conditions quickly disappear and give way to the deceleration phase noted by kobayashi et al .
( 1999 ) . during the entire evolution
there is a local maximum in @xmath34 in advance of the propagating blob that accompanies the leading pile - up in density .
the region of highest @xmath2 lies within the density minimum that lies just behind this shocked region .
although the initial pressure is small , after the bow shock is fully developed we have @xmath87 throughout the subsequent evolution .
the evolution proceeds in a roughly self - similar manner and the forward / reverse shock system is stable , in accord with analytic estimates ( wang , loeb , & waxman 2002 ) .
figures 4 and 5 show a time series of contour plots in @xmath2 and @xmath32 that follow the blob propagation and radial spreading .
the full 500 grid points along the @xmath50 direction are shown , and each frame represents 240 time steps ( i.e. , @xmath88 ) .
the leading contours show the density enhancement associated with the shock .
the @xmath32 contours indicate the values 0.2 , 0.3 , 0.4 , and 0.5 , while the ( trailing ) @xmath2 contours , indicate values of 1.25 , 2.5 , 5 , 10 , 15 , and 20 .
the contours are formed by taking a cut at @xmath89 through the @xmath90 plane .
the initial velocity vector of the blob points toward the top of the plot , with @xmath57 .
one sees a strong bow shock associated with the maximum in pressure at the point where the blob encounters the csm .
the lack of a bow shock in the @xmath2 contours suggests that little or no material in advance of the blob is accelerated to a significant bulk lorentz factor .
as the evolution progresses the deceleration causes the higher @xmath2 contours to disappear , and those that remain become increasingly distorted .
the bilateral symmetry of the contours evident in figures 4 and 5 attests to the power of the third order differencing scheme given in db02 , insofar as the 3d model has no enforced symmetry .
figure 6 shows the evolution of the total rest - mass energy @xmath44 ( dotted ) and total energy @xmath45 ( dashed ) within the grid , over the course of the conical run .
the values of @xmath44 and @xmath45 have been normalized to their initial values .
the curves that extend significantly below an abscissa value of unity indicate the values summed over the grid . at the six faces of the computational box we continually reset the values of all variables to their initial values , and keep track of the differences between those values and the initial ones .
as the run progresses and more high velocity gas reaches the edge of the computational domain and is extracted , there is a net extraction of positive @xmath44 and @xmath45 from the grid .
in addition , there is a low density wake that trails the blob , and therefore at the contact points between the wake and the sides of the grid there is a net extraction of negative @xmath44 and @xmath45 .
from a practical standpoint , the density inside the wake is so small that the net extraction of negative @xmath44 and @xmath45 does not significantly affect the bookkeeping .
the curves in figure 6 that lie close to an abscissa value of unity indicate the summed values for @xmath44 and @xmath45 corrected for both net extracted positive and negative rest - mass energies .
the second panel of figure 6 gives an expanded view of the first panel around the abscissa value of unity .
after doing the bookkeeping on the @xmath44 and @xmath45 in the extracted gas , conservation of total rest - mass energy and total energy are good to within @xmath91 parts in @xmath92 .
figure 7 shows values of total rest - mass energy within selected high @xmath2 cuts through the computational domain for the initially conical run , and the evolution of the @xmath35component of momentum @xmath93 for five different runs @xmath40 the initial blob as a plate ( `` p1 '' and `` p2 '' ) , sphere ( `` s '' ) , cone ( `` c '' ) , and cone with zero spreading ( `` @xmath94 '' ) .
the weighting function used in evaluating @xmath93 is @xmath95 .
the run p2 uses fewer grid points ( @xmath96 ) than the other runs , thus material shunted aside by the bow shock reaches the edges of the computational domain earlier . in spite of this , the evolution of p1 and p2 is quite similar , i.e. , there is no dynamically significant back reaction on the material inside the computational domain induced by the departing fluid .
as we detail in the discussion section , the drag force on the blob increases as the effective area presented by the blob , which for a roughly constant lateral spreading goes as @xmath97 .
therefore the integral of the relativistic impulse equation @xmath98 provides a solution of the form @xmath99 , where our fitting for the conical run ( @xmath100 , @xmath101 , and @xmath102 ) is shown by the curve labeled `` @xmath103 '' .
the curves in figure 7 showing the variation of the rest - mass energy contained within different @xmath2 cuts follow the values @xmath104 1.001 , 1.01 , 1.1 , 1.5 , 2 , 4 , 6 , and 8 . the higher @xmath2 cuts @xmath105 reveal that a negligible fraction of csm matter gets accelerated to significant values .
this is because the curved bow shock shunts material laterally in front of the advancing blob , rather than accelerating it up to a significant fraction of the blob s bulk lorentz factor @xmath106 .
the rest - mass energy curves are relatively constant in time up to @xmath107 . at roughly this point
the total accumulated csm gas becomes comparable to that in the initial blob .
if the blob did not spread laterally , this time would occur when the blob had plowed through a density @xmath108 , where @xmath109 and @xmath110 .
this means that the blob would have shed of order its initial momentum by the time it had traveled @xmath65 times its initial length @xmath111 , or @xmath112 units . in practice
, the initially imposed tangential velocity @xmath113 leads to the lateral expansion that effectively increases the cross section for interaction of the blob .
the negligible decrease in @xmath93 for the @xmath94 run demonstrates the importance of this effect .
figure 8 shows the lateral expansion of the blob for the conical run .
the blob edge is computed in each time step by first finding the local maximum in either @xmath32 , @xmath2 , or @xmath114 along the @xmath35axis , and then stepping laterally to the position at which the background ( csm ) value of the relevant quantity has been increased by 10% due to expansion of the blob .
the transient lateral expansion @xmath115 lasting until @xmath116 is unphysical insofar as it occurs during the period of adjustment to the initial density and velocity profiles .
the value of the later spreading rate @xmath117 is basically dictated by the initial @xmath113 value given to the gas .
figure 9 shows the evolution of @xmath93 and @xmath1 as the blob propagates .
the solid curve in the top panel indicates the same weighted value of @xmath93 as shown in fig .
after the transient physical conditions associated with the initial state have vanished , one sees a period of deceleration associated with the increasing drag force .
the dashed curves indicate the reciprocal of the full width spreading angle determined by summing and averaging the local values of @xmath118 . due to the complicated pattern of spreading and the mismatch between contours of constant @xmath32 and constant @xmath2 , the specific value for @xmath119 depends on where one cuts off the outer edge of the blob .
the three dashed curves use limiting lower values of 4 , 5 and 6 , respectively , for the @xmath2 value , in determining the volume averages that enter into @xmath119 .
the spherical / jetlike transition occurs when @xmath120 drops below the @xmath119 value computed using the lower limit @xmath121 = 6 .
the bottom panel shows the number of cells used in computing the @xmath119 values . at late times
the disappearance of higher velocity matter makes problematic a calculation of @xmath119 based on cuts in @xmath2 .
figure 10 shows the evolution of our canonical `` luminosity '' measure , namely @xmath122 , for various viewing angles @xmath79 .
as noted previously , the light curves are calculated by summing the emission from all points with the computational domain that lie , in a given time step , on 100 @xmath85 planes which straddle the position of the blob and advance toward the observer at @xmath86 .
one sees a slight break in the `` light curve '' for the @xmath123 observer at @xmath124 due to the transition from spherical to jetlike expansion .
although our model is too simple to allow any detailed quantitative comparison with observations at this point , it is worth noting that many of the observed afterglow transitions are also smooth ( e.g. , fig . 2 of stanek et al .
1999 ; fig . 1 of harrison et al .
1999 ) , rather than the abrupt breaks evident in some of the semi - analytical models .
the transition is not so visible for off - axis viewing angles .
the general pattern in which increasingly off - axis observers see significantly lower emission until well past maximum light curve can be seen in the light curves shown in fig .
4b of panaitescu & mszros ( 1999 ) .
afterglows coming from those grbs for which the observer does not lie within the initial cone of grb emission are referred to `` orphan''afterglows .
future surveys of orphan afterglows might provide constraints on the beaming that are different from those of frail et al .
( 2001 , e.g. , totani & panaitescu 2002 ; levinson , ofek , waxman , & gal - yam 2002 ) .
utilizing 3d relativistic hydrodynamical calculations , we have examined the evolution of an expanding relativistic blob of gas intended to be representative of a jet associated with ejecta from an extremely energetic event such as a hypernova , that produces a gamma - ray burst ( aloy et al .
2000 ; tan , matzner , & mckee 2001 ; macfadyen , woosley , & heger 2001 , zhang , woosley , & heger 2003 , zhang , woosley , & macfadyen 2003 ) . since these are the first such calculations applied to the blob during the time in which the afterglow radiation is produced , we have purposely kept them simple in an effort to concentrate on the most fundamental aspects of the physics .
we restrict our attention to the transition from spherical to jetlike expansion that occurs during the time that the lorentz factor becomes less than the reciprocal of the jet spreading angle .
we have not yet attached specific numbers to our results . from the srhd equations
, one sees that the relevant quantities are the ratios of pressure to density , and of distance to time . if we specify either one of these two sets of numbers , the other one is also determined .
the column giving the observed afterglow break time @xmath125 in table 1 of frail et al .
( 2001 ) indicates @xmath126 d as being representative .
for an observer directly on the velocity vector of the blob , the time @xmath127 between the grb and afterglow    @xmath128 \simeq { 1\over 2c }      \int \delta x \
\gamma(x)^{-2},\ ] ]    where the dominant contribution to the integral comes from later times .
thus the light travel time of 1 day is multiplied by @xmath129 , assuming the spherical - to - jetlike transition giving the break in the afterglow light curve happens at @xmath130 . for the conical run ,
the break in the light curve occurs at at @xmath131 .
if we designate this point as corresponding to a time 2 d , then @xmath132 translates to @xmath133 cm @xmath134 cm .
since we have 100 grid points per unit length , this means each grid point spacing corresponds to @xmath135 cm .
zhang , woosley , & macfadyen ( 2003 ) utilize a 2d relativistic hydrodynamics code to follow the evolution of the blob through the envelope of the progenitor star and past the point of `` break - out '' through the stellar surface .
they find that a blob launched at @xmath136 km with an initial spreading angle of @xmath137 and a bulk lorentz factor of 50 experiences strong shock heating and also lateral confinement , so that when it emerges from the progenitor it has a large internal energy @xmath138 , a bulk lorentz factor @xmath139 , and is more confined than initially ( @xmath140 ) . the subsequent expansion due to the large @xmath141 value leads to an effective @xmath142 for a distant observer .
in the standard model for grbs , the gamma radiation is produced by strong internal shocks in the expanding fireball at the point where it becomes optically thin to its own radiation , at @xmath143 s , where @xmath144 is the total energy carried by the jet , corrected for beaming , in units of @xmath145 ergs , and @xmath146 is the csm density in units of @xmath147 @xmath148 . during the time of gamma ray emission , the bulk lorentz factor of the ejecta @xmath149 .
zhang , woosley , & heger ( 2003 ) compute the evolution after break - out covering the time of the prompt grb emission .
they estimate it should occur at about @xmath150 to @xmath151 cm and have a duration @xmath152 s. in this work we focus on the time after this event , i.e. , subsequent to the evolution considered by zhang et al .
( 2003ab ) , when the `` afterglow '' is produced .
this emission is thought to arise primarily in the bow shock where the strong heating leads to @xmath153 ( wang & loeb 2001 ) .
in essence , this energy production arises from `` external shocks '' , as opposed to the `` internal shocks '' that generate the gamma radiation ( rees & mszros 1994 ; panaitescu , & mszros 1998a ; spada , panaitescu , & mszros 2000 , kobayashi & sari 2001 ) .
after the blob is ejected from the envelope of the star that was the hypernova precursor and continues to propagate , the medium through which it travels should be dominated by the density profile left from the precursor s stellar wind ( chevalier & li 1999 , 2000 , li & chevalier 2001 ) .
if the density varies as @xmath154 away from the star , the mass loss of the wind @xmath155 , and the wind velocity @xmath156 km s@xmath157 , the circumstellar density at the point where the afterglow radiation is emitted @xmath158 cm will be @xmath159 @xmath148 . in the calculation of zhang et al ( 2003 ) , the density inside the jet is @xmath160 @xmath148 at @xmath161 cm , and @xmath162 . if the spreading angle of the jet remains roughly constant from @xmath161 cm to @xmath158 cm ( lithwick & sari 2001 )
, then the density inside the jet at @xmath158 cm should be lower by @xmath163 , or about @xmath164 @xmath148 . at this point @xmath34 inside the jet
should be negligible compared to @xmath32 .
the density contrast between the blob and csm is greater than what we have assumed in the results shown previously .
a separate run taking an initial density contrast of @xmath165 shows the same basic effects as the previous runs , however , namely a lateral expansion rate of @xmath166 and the non - accretion of csm gas .
our computations lie within the `` deceleration phase '' discussed by kobayashi et al .
( 1999 , see also kobayashi & sari 2000 , 2001 ) .
we find a change in the form of the luminosity decrease corresponding to the transition between spherical and jetlike expansion .
the determination of the average spreading angle @xmath119 is nontrivial because it depends on how one does the averaging , and how much of the diffuse , sideways - expanding jet material is included in the computation . in figure 9
we presented cuts for gas possessing @xmath167 4 , 5 , and 6 as representative of material in the flow that partakes most strongly in producing the observed radiation .
we do not see a dramatic increase in @xmath119 during the deceleration phase ; the @xmath119 value basically reflects the ballistic motion of material following its initial @xmath113 value . also ,
because the deceleration in our problem is forced not by the accumulation of gas from the csm but rather the drag force of the csm on the blob , we do not see an exponential decrease in @xmath2 with distance during the deceleration phase , but rather a decrease well - approximated by a cubic dependence that is expected for a roughly constant lateral spreading rate .
many of the previous studies , utilizing the formalism of bm76 in which matter accreted from the csm is added onto the grb shell and facilitates the deceleration , have discussed the sideways expansion of the jet that supposedly occurs near the time that most of the initial kinetic energy of the blob has been shed .
( a qualitative depiction of this effect is shown in fig . 1 of piran 2002 . )
our results do not appear to support this viewpoint .
we find that the csm gas hitting the bow shock does not add significantly to the mass of the jet , but rather is swept back into its wake as the blob passes a given location . in this sense , the relevant factor in determining the deceleration is simply the projected csm surface density relative to that in the blob . as a result , there is no runaway phase of rapid lateral expansion @xmath40 a physical process that has been `` built in '' to some previous theories .
we do not see an abrupt jump in the lateral spreading coincident with the deceleration phase ; the lateral spreading rate is mandated basically by the initial @xmath113 value imparted to the gas .
the pressure - driven lateral expansion is negligible until such time as @xmath168 , which for @xmath169 is so late as to be uninteresting . in the aforementioned works there is a close connection between the transition @xmath170 and the change from spherical to jetlike emission , whereas in our study
there is not .
the break in the light curve produced by emission from the bow shock is due solely to the lessening of relativistic beaming brought about by a slow , uniformly increasing drag from the csm , and the fact that once @xmath3 becomes less than @xmath1 , the observer perceives the finite width of the jet .
rhoads ( 1999 ) estimates that the lateral spreading speed of the jet should be @xmath171 .
sari , piran , & halpern ( 1999 ) argue for a faster spreading rate @xmath172 .
rhoads ( 1999 ) divides the dynamical evolution of the blob into two regimes , characterized by a power law decay ( in time ) of the bulk lorentz factor , followed by an exponential decay . during the latter stage , the swept - up mass increases exponentially in time .
panaitescu & mszros ( 1999 ) calculate light curves for observers at varying angles from the jet axis , and calculate separately the effects of including and excluding the lateral jet expansion .
they find that the maxima in the light curves occur substantially later in runs which do not take into account the jet broadening ( see their fig .
4 ) . in the analytical models of both rhoads ( 1999 ) and
panaitescu & mszros ( 1999 ) , the physics of mass accumulation from the csm is an integral component of the formalism ; all mass within the solid angle of the expanding shell is assumed to accrete .
workers have applied the results of rhoads ( 1999 ) and sari , piran , & halpern ( 1999 ) to the afterglow evolution , however the results of zhang et al .
( 2003 ) cast doubt on the validity of this exercise , because at the time corresponding to the afterglow emission one anticipates that @xmath173 and the lateral spreading rate would not be governed by the internal sound speed but rather the ballistic motions of the ejecta comprising the blob as they leave the vicinity of the progenitor star .
the blob has a large internal thermal energy @xmath174 as it emerges from the progenitor star , and even by the time the expanding ejecta have become optically thin to their own emission from internal shocks ( producing the grb ) , one still expects @xmath175 . by the ( much later ) time of the afterglow emission ,
however , the blob would have cooled to the point that @xmath176 , which provides the impetus for our initial condition @xmath177 . following the arguments of rhoads ( 1999 ) and sari , piran , & halpern ( 1999 ) ,
if the lateral spreading rate were mandated by the internal sound speed , then in our calculations it should be @xmath178 , whereas we find it to be @xmath139 times larger .
the spreading rate basically reflects our initial @xmath113 value .
thus the physics of the lateral expansion is different than in rhoads ( 1999 ) and sari , piran , & halpern ( 1999 ) .
how might this result be influenced by systematic effects present in our model ?
one obvious potential shortcoming is the absence of cooling . in this work
we have assumed an adiabatic gas , whereas in reality one might envision the presence of cooling within the shock .
this might then reduce the ability of the shock to deflect the gas as efficiently as it now does , which in our model prevents the acceleration of csm material to lorentz factors approaching that of the blob .
we adopt a relatively small density contrast between material inside the blob and the exterior csm in order to see significant evolution of the blob during the course of the simulation within our eulerian grid .
we have run additional models using @xmath179 , and although these could not be run for sufficient time to see the deceleration phase begin , we find in these runs also a lateral spreading rate @xmath180 , for the same initial @xmath113 values , and negligible accretion of gas from the csm .
we may understand the deceleration of the blob with a simple analytical model .
if we assume that the only significant contribution to @xmath181 is through forward motion , that the drag force varies as @xmath182 , and that @xmath183 , then @xmath184 implies    @xmath185    where @xmath186 , @xmath187 is the cross sectional area for interaction @xmath188 , @xmath189 is the ( constant ) lateral spreading speed ( @xmath180 in our calculations ) , and @xmath190 is a dimensionless number ( typically of order 0.2 in laboratory applications ) . taking @xmath191 and @xmath192 gives    @xmath193    where the blob mass @xmath194 .
the blob volume is given by the number of cells comprising the blob initially ( @xmath26 ) times the volume of an element @xmath195 .
therefore @xmath196 .
integrating gives    @xmath197    where @xmath198 is the initial value of @xmath199 and @xmath200 . from the fitting to @xmath201
presented in figure 7 , we infer that @xmath202 . the specific numerical value for @xmath190 is probably influenced by our numerical resolution .
a comparison of the curves labelled `` c '' and `` @xmath103 '' in fig .
7 shows that the functional decrease in @xmath201 with @xmath50 is reasonably described by a cubic , as expected if the cross sectional area increases quadratically with @xmath50 , or equivalently @xmath203 .
note that we assumed @xmath204 in this exercise , which is a good approximation for the evolution of interest .
the calculations we present are the first 3d relativistic hydrodynamical calculations of grb jet evolution pertinent to the afterglow phase that do not enforce any special symmetry ( e.g. , spherical or axial ) .
we find that ( i ) the csm gas does not accrete onto the advancing blob , but rather is shunted aside by the bow shock , ( ii ) the decay light curve steepens roughly when one first `` sees '' the edge of the jet @xmath4 , with this effect being strongest for `` face - on '' observers ( confirming previous studies ) , and ( iii ) the rate of decrease of the @xmath35component of momentum @xmath93 is well - characterized by a simple model in which the cross sectional area of the blob increases quadratically with laboratory time ( or distance ) .
the primary impetus for the built - in assumption of accretion of matter in previous studies was the influential work of bm76 in which spherical relativistic expansion was considered .
accretion of gas onto the relativistically expanding shell is obviously justified for spherical expansion , but subsequent grb workers applied the results to the case of the grb jet , in which a thin wedge of material propagates through a low density medium .
in such a situation the natural tendency of material in front of the jet is to be pushed aside and to form a `` channel flow '' around the jet , rather than to accrete .
a separate issue is that workers used the results of rhoads ( 1999 ) and sari et al .
( 1999 ) that give lateral expansion rates of order @xmath86 to explain the phase of afterglow evolution , while it now turns out from the work of zhang et al .
( 2003ab ) that the much smaller internal energy in the afterglow - producing blob would be expected to give a much lower spreading rate , were that the only relevant spreading mechanism .
the evolution of the blob prior to the afterglow epoch , however , in particular its emergence from the progenitor envelope , did encompass a period of much greater @xmath141 during which a large @xmath189 was imparted to the gas . in summary , there is nothing in our results to suggest that bm76 , rhoads ( 1999 ) , and sari et al . ( 1999 ) are not internally consistent , rather it appears that the subsequent application of their results to afterglow evolution may have been inappropriate .    recent semi - analytical work has centered on `` structured '' relativistic jet modeling , wherein one replaces the `` top - hat '' assumption of uniform jet properties , e.g. , lorentz factor and @xmath141 constant within the initial cone angle , with more physically motivated forms in which the physical parameters vary with angle away from the jet symmetry axis @xmath79 ( e.g. , kumar & granot 2003 , granot & kumar 2003 ) .
structured jets represent a refinement in the semi - analytical and analytical work in the sense that unphysical artifacts of the top - hat models are avoided . in this work we allow naturally occurring gradients within the flow
determine the evolution .
in comparing the runs with different starting conditions , for instance , we see that only the gross initial characteristics ( bulk lorentz factor and spreading rate ) are of relevance .
one could impose a structured initial state , in terms of having an angle - dependent relation between , for instance lorentz factor and @xmath79 , and/or @xmath141 and @xmath79 , but one suspects that the initial intricacies would be washed out , just as the initial blob shape is .
one might in fact imagine inverting the problem and using our results to determine physically motivated functional forms for @xmath205 and @xmath206 as input to a structured jet formalism , but in our models the detailed evolution of the physical parameters and the interaction of the blob with the bow shock reveal gross changes along the jet axis .
therefore the motivation for trying to characterize the jet properties as simple functions of @xmath79 , averaged along the axis of the jet , seems questionable . in addition
, one would still be missing important physical effects , such as the non - accretion of the forward csm material .
two obvious refinements , currently being carried out , are to ( i ) treat the problem on a lagrangian grid in which the mesh points follow the blob and are adaptively inserted in regions with strong gradients , so as to be able to explore regimes in which the density contrast between the blob and csm is much larger , and ( ii ) include provisions for realistic bremsstrahlung and synchrotron physics , in order to produce light curves that can be compared directly with observations so as to test different aspects of the theory and thereby constrain the allowed parameter space .
our sincerest thanks go to luca del zanna who provided key insights into the workings of db02 .
we also thank john baker , david band , tom cline , chris fragile , chris fryer , markos georganopoulos , peter goldreich , demos kazanas , pawan kumar , zhi - yun li , andrew macfadyen , peter mszros , ewald mueller , jay norris , alin panaitescu , sterl phinney , steve reynolds , james rhoads , steve ruden , sabrina savage , craig wheeler , and weiqun zhang for helpful comments .
the evolution of a spherical relativistic expansion due to a large over - pressure inside a sphere of radius 0.1 , adopting hll fluxes .
the evolution encompasses 360 time steps , taking a courant number of 0.05 .
shown are the initial conditions plus ten equally spaced time steps taken from a slice along the @xmath35axis depicting the evolution of ( i ) pressure @xmath34 ( _ top panel _ ) , ( ii ) density @xmath32 ( _ second panel _ ) , ( iii ) @xmath207 ( _ third panel _ ) , and ( iv ) lorentz factor @xmath2 ( _ bottom panel _ ) .
the small numbers beside each curve indicate the evolution .
one sees an expansion of the central over - pressurized sphere into the surrounding medium .
the strong relativistic outflow peaks at @xmath208 near the end of the evolution .
the general trend in which @xmath209 within the expansion can be derived from fundamental principles and is well known from classical solutions . at late times
the strong decrease in @xmath32 at the center of the sphere , which is a singularity in this test , begins to cause numerical instability .
we utilize the llf fluxes for @xmath36 and the hll fluxes for @xmath210 .
the evolution of a spherical relativistic expansion due to a large over - pressure inside a sphere of radius 0.1 , adopting pure llf fluxes .
quantities shown are the same as in fig .
the evolution is nearly identical to that in fig .
1 , except the variations in the physical variables are smoother , and @xmath211 at larger radii by the end of the run @xmath40 somewhat smaller than in fig . 1 .
the evolution of a relativistic blob launched as a small cone traveling along the @xmath212-axis , initially confined to between @xmath213 and @xmath214 and a maximum radius 0.025 .
the small numbers beside each curve indicate the time step .
the panels are the same as shown in fig .
1 . snapshots represent conditions along a slice through the center of the rectangular grid , taken every 120 time steps ( @xmath215 ) . for this trial and all that follow
we utilize a courant number of 0.25 and pure lff fluxes .
the initial bulb lorentz factor @xmath2 was set to 25 , and the initial full width flaring angle to 0.07 radian as measured by @xmath216 .
the spreading in @xmath32 in the @xmath217 plane .
the eight panels depict the initial state , and seven subsequent snapshots taken every 240 time steps ( i.e. , @xmath88 ) .
the evolution represents a slice through the midplane of the grid , and the full length and width of the box are shown .
contours indicate values of 0.2 , 0.3 , 0.4 , and 0.5 ( in the same dimensionless units as shown in fig .
3 ) .    figure 5 . the spreading in @xmath2 in the @xmath217 plane .
the conventions are the same as in the previous figure .
contours indicate values of 1.25 , 2.5 , 5 , 10 , 15 , and 20 . at later times
the higher @xmath2 contours disappear due to deceleration .
the absence of a bow shock in this depiction shows that material in front of the advancing blob is not accelerated to @xmath2 values much above unity .
conservation of total rest - mass energy @xmath44 and total energy @xmath45 . at @xmath218
material starts to reach the edges of the grid and flow off , therefore one sees the start of a linear decrease in @xmath44 and @xmath45 .
the curves which closely follow @xmath219 show the values for @xmath44 and @xmath45 after correcting for material extracted from the grid .
the bottom panel shows an expanded view of the region near @xmath219 in the first panel .
the evolution of rest mass energy contained within various mass cuts , for the conical run ( solid lines ) , and the evolution of the weighted @xmath35component of the momentum @xmath93 , divided by 5 so as to be on a common scale with the rest mass energy curves ( dashed lines ) .
the numbers associated with the solid lines indicate the @xmath2 value used in each rest mass energy cut .
we show the @xmath93 evolution for the conical run ( c ) , the spherical run ( s ) and the two plate runs ( p1 : @xmath220 , p2 : @xmath221 ) . for these four runs the initial spread in velocities
is such that @xmath222 , whereas for the run labelled @xmath94 there is no tangential component initially to the velocities .
consequently the reduced lateral spreading rate leads to a slower @xmath93 decrease .
the function @xmath103 represents a fitting to the conical run given by @xmath223 .
although the rest mass energy curves for @xmath224 and 1.01 show a sharp increase ( up until @xmath225 when material starts to leave the edges of the grid ) , the near constancy of the @xmath226 curve and the slight decline for higher @xmath227 curves indicate that there is negligible acceleration of high@xmath228 material . also , the functional form of the decrease in @xmath93 is similar for all runs in which @xmath222 initially .
the location of the edge of the jet , determined by first finding the location along the direction of propagation of the maximum in either @xmath2 ( _ top curve _ ) , @xmath114 ( _ middle curve _ ) , or @xmath32 ( _ bottom curve _ ) , then stepping laterally to the point at which the background value has increased by 10% of its original value .
( the precise value of this constant does not affect the results . )
the initial maximum radius of the blob is 0.05 .
after @xmath229 the two lower curves coincide .
the evolution of the weighted mean of the @xmath35 component of momentum and a measure of the spreading angle @xmath119 , taken to be the ratio of tangential to axial speeds , in the observer s frame . shown
are ( i ) @xmath93 ( solid line ) and the reciprocal of @xmath119 , measured in radians ( dashed lines ) , where the lower limiting @xmath2 value used in the averaging for @xmath119 is taken to be either 4 , 5 , or 6 ( _ top panel _ ) .
( ii ) the @xmath119 values whose reciprocal values are indicated in the first panel ( _ middle panel _ ) , and ( iii ) the number of cells entering into the averages for the three @xmath121 values in the first two panels ( _ bottom panel _ ) . at late times
the deceleration makes the @xmath119 calculations problematic because the number of high @xmath2 cells drops .
this progression in the loss of high @xmath2 cells can be seen at @xmath230 in the bottom panel .
figure 10 .
light curves constructed by summing the quantity @xmath231 over the grid .
we take this global sum as being a measure of the luminosity seen by an observer looking down the jet .
the light curves are built up during the course of the run by summing the emission from 100 slices moving toward the observer at @xmath86 that straddle the location of the disturbance induced by the blob .
thus the light curve consists of a maximum of 100 points , where @xmath232 is measured with respect to the arrival time of the first comoving slice .
the five values of the angle @xmath79 between the observer and the local velocity vector of a fluid element are ( going to smaller emission values ) @xmath76 , @xmath233 , @xmath234 , @xmath235 , and @xmath77 .
for the `` face - on '' viewer ( inclination @xmath236 ) , there is a mild change in the decay slope , i.e. , break , between the regimes of `` spherical '' and `` jetlike '' expansion .