index coding , introduced in @xcite , considers the problem of efficiently broadcasting a number of messages available at a source , to receivers that already possess some prior knowledge of the messages . in the index coding framework ,
the source is allowed to encode the messages ( and thereby use the channel efficiently ) while satisfying the receiver demands . the general class of _ groupcast _
index coding problems consists of @xmath0 messages generated at a source , where each message is demanded by at least one receiver .
index coding problems where each receiver demands a unique message are called _ single unicast index coding _ problems and are the most widely studied class .
although index coding continues to be open in general , several researchers have made inroads into characterising the _ _ rate of index coding _ _ and presenting achievable schemes .
the landmark paper @xcite famously connected the scalar linear index coding problem to finding a quantity called _ minrank _ associated with the _ side - information graph _ related to the given single unicast index coding problem .
upper and lower bounds on the rate for single unicast index coding have been presented via graph theoretic ideas like clique cover , chromatic number @xcite , local chromatic number @xcite , fractional clique covering and hyperclique covering @xcite , and recently , the ` generalized interlinked cycle cover ' @xcite .
many of these papers naturally lead to constructions of ( scalar and vector ) linear index codes .
linear codes however are not always found to be optimal @xcite .
random coding approaches to index coding were studied in @xcite .
bounds on the rate of groupcast index coding were presented in @xcite .
interference alignment , well known as a powerful tool to study degrees of freedom in wireless interference networks , was employed to the linear index coding problem in @xcite , by modelling the unavailable side - information as interference .
the idea is to assign precoding matrices to the message vectors such that all receivers can decode even in the presence of interference ( thereby requiring some degree of linear independence between the precoding matrices , and hence reducing the rate ) .
however , at the same time the interference at the receivers must be as ` aligned ' as possible ( in order to reduce the amount of linear independence required , i.e. to increase the rate ) .
this technique was further explored in @xcite and several classes of index coding instances with certain feasible rates were identified based on the properties of the interference seen by the receivers .
this work builds primarily upon the results in @xcite . in @xcite , a necessary and sufficient condition for the feasibility of _ _ rate half _
_ ( some finite field ) symbols transmitted through the channel , one @xmath3 symbol of each message is conveyed . ] in a groupcast index coding problem was established based on the properties satisfied by two graphs obtained from the interference structure of the problem , called the _ conflict graph _ and the _ alignment graph_. also , a sufficient ( but not necessary ) condition on the structure of these graphs was given for rate @xmath2 feasibility .
prior work @xcite also gives a necessary ( but not sufficient ) condition for rate @xmath2 feasibility based on the interference structure of the problem .
the relevant details of the prior work are discussed elaborately in the forthcoming sections of this paper .
* firstly , we revise the definition of the conflict graph given in @xcite
. the conflict graph definition in @xcite does not capture the interference structure of the index coding problem sufficiently . in order to rectify this
, we define the _ conflict hypergraph _ , which is shown to capture the interference structure sufficiently .
( subsection [ conflicthypergraphsubsec ] ) * our main result shows that rate @xmath2 is achievable in a given index coding problem under certain conditions on the topology of the alignment graph and conflict hypergraph of the index coding problem .
the sufficient condition which we present for rate @xmath2 feasibility is looser than the sufficient condition shown in @xcite .
therefore , this class of index coding problems is bigger than the previously known class of rate @xmath2 feasible problems .
the achievability of rate @xmath2 in such problems is shown by presenting a construction of an index code by random generation of precoding vectors over a large field .
( subsection [ rateonethirdfeasibleproblems ] ) * towards obtaining our main result , we also obtain a necessary condition for rate @xmath2 feasibility , which is stricter than the prior condition from @xcite ( subsection [ necessaryconditionrateonethirds ] ) .
our feasibility conditions can thus be seen as being midway between those of @xcite and @xcite .    _
notations : _ throughout the paper , we use the following notations .
let @xmath4 $ ] denote @xmath5 .
for a set of vectors @xmath6 , @xmath7 denotes their span . for a vector space @xmath8
, @xmath9 denotes its dimension . an arbitrary finite field
is denoted by @xmath3 .
a vector from the @xmath10-dimensional vector space @xmath11 is said to be picked _ at random _ if it is selected according to the uniform distribution on @xmath11 .
formally , the index coding problem ( over some field @xmath12 ) consists of a broadcast channel which can carry symbols from @xmath3 , along with the following .    * a set of @xmath13 receivers * a source which has messages @xmath14\}$ ] , each of which is modelled as a vector over @xmath3 . * for each receiver @xmath15 , a set @xmath16 denoting the set of messages demanded by the receiver @xmath15 . * for each receiver @xmath15 , a set @xmath17 denoting the set of side - information messages available at the @xmath18 receiver .
this general class of index coding problems is known as _ groupcast _ index coding problems .
an index code of symmetric rate @xmath1 for a given index coding problem consists of an encoding function @xmath19 for some @xmath20 mapping the @xmath0 @xmath21-length message vectors ( @xmath22 ) to some @xmath23-length codeword which is broadcast through the channel , as well as decoding functions @xmath24 at the receivers @xmath25,$ ] mapping the received codeword and the side - information messages to the demanded messages @xmath26 , i.e. , @xmath27.\ ] ]    we could in general have different rates for different messages , but in this paper we restrict our attention to symmetric rates .
therefore any rate referred to in this paper is the symmetric rate .    for a given index coding problem , a rate @xmath1 is said to be _ achievable _ if there exists an index code of rate @xmath1 , and the index coding problem is said to be _ rate _ @xmath1
_ feasible_.    if a rate @xmath28 is achievable , the associated index code is a _ scalar index code of length
_ @xmath23 .
if the encoding and decoding functions are linear , then we have a _ linear index code_.    if we have a linear index code of rate @xmath1 , then we can represent the encoding function as follows .
@xmath29 where each @xmath30 is a @xmath31 matrix with elements from @xmath12 . in the case of scalar linear index coding
, we have @xmath32 . finding a scalar linear index code of length @xmath23 ( i.e. , with a feasible rate @xmath33 ) is equivalent to finding an assignment of these @xmath23-length vectors @xmath30s to the @xmath0 messages such that the receivers can all decode their demanded messages , i.e. , @xmath34.\ ] ] in the case of a scalar linear index code , the encoding function @xmath35 fixes the vectors assigned to the messages . however not every vector assignment is a valid index coding encoding function , as the decoding functions may not exist . in a number of proofs in this paper ,
we start with some vector assignment and show that it leads to an encoding function of a valid index code .
therefore , in such proofs we refer to the initial vector assignment as a valid encoding function @xmath35 by abusing the notation .
we restrict our attention to scalar linear index codes for the rest of this paper .
however we believe that our results can be extended to vector linear index codes as well .      for some receiver @xmath15 and for some message @xmath36 ,
let @xmath37 denote the set of messages ( except @xmath38 ) not available at the receiver @xmath15 .
the sets @xmath39 are called the _ interfering sets at receiver _ @xmath15 .
if receiver @xmath15 does not demand message @xmath38 , then we define @xmath40 .
if a message @xmath41 is not available at a receiver @xmath15 demanding at least one message @xmath42 , then @xmath41 is said to _ interfere at receiver _
@xmath15 , and @xmath41 and @xmath38 are said to be _ in conflict_.    for a set of vertices @xmath43 , let @xmath44 denote the vector space spanned by the vectors assigned to the messages in @xmath6 , under the specified encoding function @xmath35 .
if @xmath45 , we define @xmath44 as the zero vector .    for a given assignment of vectors to the messages ( or equivalently , for a given encoding function @xmath35 ) , we say that _ conflicts within a subset _ @xmath46 _ are resolved _ , if @xmath47,\end{aligned}\ ] ] where @xmath48 is the vector assigned to @xmath38 under the encoding function @xmath35 .
if ( [ resolvedconflicts ] ) holds for @xmath49 then _ all the conflicts _ in the given index coding problem _ are said to be resolved_.    we now state a simple lemma , rephrased from @xcite , which is easily proved .
[ lemmasuccessdecoding ] for any encoding function @xmath35 , successful decoding at the receivers is possible if and only if all the conflicts are resolved .
let @xmath48 be the vector assigned to @xmath38 under the encoding function @xmath35 .
_ if part : _ at all receivers @xmath15 , consider that we have @xmath50 $ ] .
what is received at a receiver @xmath15 is the codeword @xmath51 , from which it wants to obtain @xmath38 . as receiver @xmath15
can always subtract the contribution from the side - information messages from @xmath51 , it only remains to be shown that @xmath52 can be decoded from @xmath53 where the equality is because @xmath54 because of the assumption that @xmath55 , receiver @xmath15 can get @xmath38 from ( [ eqn1 ] ) .
_ only if part : _ consider now that there is some receiver @xmath15 and some message @xmath56 such that @xmath57 clearly , solving for @xmath38 from ( [ eqn1 ] ) does not lead to an unique solution . hence decoding fails .
this concludes the proof .    by lemma [ lemmasuccessdecoding ] ,
it should also be clear that if there is an assignment of @xmath23-length vectors @xmath30s to the messages @xmath41s such that the condition in lemma [ lemmasuccessdecoding ] is satisfied , then these vectors naturally define an index code of length @xmath23 for the given index coding problem .
in @xcite , the authors defined the notions of _ alignment graph _ and _ conflict graph _ whose properties were used to characterise index coding problems for which rate @xmath58 is feasible .
both of these graphs have the same vertex set , which is the set of messages @xmath59 .    in the alignment graph , the vertices @xmath41 and
@xmath52 are connected by an edge ( called an _ alignment edge _ , shown in our figures by a solid edge ) when the messages @xmath41 and @xmath52 are not available at a receiver demanding a message other than @xmath41 and @xmath52
. a connected component of the alignment graph is called an _ alignment set_.
it is easy to see that the alignment sets define a partition of the alignment graph . also , the messages in @xmath60 for all messages @xmath56 at all receivers @xmath15 are fully connected in the alignment graph .    in the conflict graph , @xmath41 and @xmath52
are connected by an edge ( called an _ conflict edge _ , shown by a dotted edge ) if @xmath41 is not available at a receiver demanding @xmath52 , or @xmath52 is not available at a receiver demanding @xmath41 .
the following example illustrates that the conflict graph definition does not capture the _ directionality _ of the conflicts .
[ conflictissueexm ] consider two single unicast index coding problems with four messages . in the first problem ,
the interfering sets are as follows @xmath61 . in the second problem ,
@xmath62 all other interfering sets corresponding to the receivers are empty .
these two problems have the same alignment and conflict graphs ( see fig.[fig : conflictgraphsame ] ) . in particular
, the conflict graphs are the same because the definition of the conflict graph does not model the directionality of the absent side - information messages .
note that they both have different solutions .
we leave it to the reader to check that the first problem is rate @xmath58 feasible , while the second problem is not rate @xmath58 feasible , but has a rate @xmath2 solution .    to overcome this issue
, we define the _ conflict hypergraph _ as follows .
the conflict hypergraph is an undirected hypergraph with vertex set @xmath59 ( the set of messages ) , and its hyperedge set defined as follows .
* for any receiver @xmath15 demanding any message @xmath38 , @xmath38 and @xmath63 are connected by a hyperedge , which is denoted by @xmath64 .
for example , the two problems presented in example [ conflictissueexm ] are now represented using different conflict hypergraphs in fig .
[ fig : conflicthyp1 ] and fig .
[ fig : conflicthyp2 ] ( the alignment graphs remains the same ) .
note that even though this definition for the conflict hypergraph does not explicitly contain direction , the directionality of interference seen by any receiver is modelled correctly whenever the number of interfering messages is more than one .
this is sufficiently general as we see in the following lemma .
suppose two index coding problems , denoted by @xmath65 and @xmath66 are modelled by the same conflict hypergraph .
then any index coding solution for @xmath65 is an index coding solution for @xmath67 .
let @xmath68 be the encoding function of the given index code for @xmath65 .
let @xmath48 be the vector assigned to @xmath38 .
by lemma [ lemmasuccessdecoding ] , we must have that for @xmath50,~\forall j\in[1:t]$ ] in @xmath65 .
now assume that the same index code is used for @xmath67 .
for any receiver @xmath15 and message @xmath69 with @xmath70 in @xmath67 , the conflict hyperedge @xmath71 present in the conflict hypergraph of @xmath67 is present also in that of @xmath65 as both @xmath65 and @xmath67 have the same conflict hypergraph . in other words ,
message @xmath72 is an interfering set of a receiver @xmath15 in @xmath65 also and hence we must have @xmath73 , which means that @xmath69 is recoverable at @xmath15 in @xmath67 also .
the only case left to check is when @xmath74 in @xmath67 for some receiver @xmath15 and message @xmath69 .
let us assume that @xmath75 because @xmath65 and @xmath67 share the same conflict hypergraphs , this conflict ( in some direction ) is present in @xmath67 also . by assignment @xmath35
, we must have that @xmath30 and @xmath48 are linearly independent , which ensures that this conflict is resolved in @xmath67 also , irrespective of its directionality .      towards showing a necessary and sufficient condition for rate @xmath58 feasibility
, the following definition for _ internal conflicts _ was given in @xcite .
a conflict between two messages within an alignment set is called an _ internal conflict_.    the following theorem
was proved in @xcite on rate @xmath58 feasible index coding problems .
[ ratehalfthm ] an index coding problem is rate @xmath58 feasible if and only if there are no internal conflicts .
the following lemma plays a crucial role in our proof of theorem [ ratehalfthm ] , which is included for the sake of completeness as it is not available in its complete form in prior literature .
we will also see in subsection [ rateonethirdfeasibleproblems ] that it is central to proving our main theorem which deals with the feasibility of rate @xmath2 .
[ intersectingkdimspaces ] let @xmath76 be @xmath77 sets of vectors , such that @xmath78 $ ] . then the space spanned by @xmath79 has dimension @xmath80 if and only if each @xmath81 spans a vector space of dimension @xmath80 .    _
if part : _ suppose each @xmath81 spans a vector space of dimension @xmath80 , and @xmath78.$ ] clearly , all the vector spaces @xmath82s are exactly equal , and thus @xmath83    _ only if part : _ suppose @xmath84 and @xmath85.\ ] ] by ( [ eqn3 ] ) , @xmath86 . if some @xmath87 , then @xmath88 , and we have a contradiction . thus , we must have that @xmath89 .
we now use lemma [ intersectingkdimspaces ] and lemma [ lemmasuccessdecoding ] to prove theorem [ ratehalfthm ] .
corresponding to any vertex @xmath38 in the alignment graph , let @xmath90 denote the alignment set it belongs to ( this is unique as the alignment sets partition the alignment graph ) .
we first note that in any ( scalar linear ) index coding scheme for the given problem , all the vertices must be assigned some non - zero vectors ( zero vector can not be assigned to any message as this means that the message can not be decoded by any receiver ) .
_ if part : _ suppose that there are no internal conflicts .
we assume a large field @xmath12 . for each alignment set ,
we independently generate a random @xmath91 vector over @xmath3 and assign it to the vertices of the alignment set . because of random generation
, we can assume that any assigned vector is non - zero and any two assigned vectors are linearly independent with high probability ( whp ) .
let @xmath35 denote the associated encoding function and @xmath48 denote the vector assigned to vertex @xmath38 .
since there are no internal conflicts , we only have to check conflicts between alignment sets . for any vertex @xmath38 ,
the set @xmath63 for any receiver @xmath15 which demands @xmath38 belongs to a @xmath92 different and @xmath93 unique alignment set than @xmath90 ( @xmath92 is because there are no internal conflicts , @xmath93 is because all the messages in @xmath63 must be in the same alignment set ) .
since any two alignment sets get independent vectors ( whp ) , we have that @xmath55 , and the same argument is true for all receivers @xmath15 and all messages @xmath38 .
hence this assignment of vectors ensures successful decoding by lemma [ lemmasuccessdecoding ] .    _ only if part : _ suppose that there is some internal conflict ( represented in the conflict graph as an edge between node @xmath94 and node @xmath56 ) in an alignment set @xmath90 .
because @xmath95 are part of the same alignment set ( connected component ) @xmath90 , there lies a path from @xmath94 to @xmath56 , given by an ordered set @xmath96 , such that every adjacent pair of elements ( @xmath97 , etc . )
belong to the interfering set of some receiver .    in some assignment corresponding to a rate @xmath58 solution ,
let @xmath98 be the non - zero vectors assigned to the vertices @xmath96 .
we define the sets , @xmath99 $ ] and @xmath100 .     and @xmath56 are in conflict
, they must be assigned linearly independent vectors .
this requires that at least one of the sets @xmath101 has to be two dimensional.,width=211 ]    suppose some @xmath102 for some @xmath103 .
then a receiver @xmath15 ( at which the messages corresponding to @xmath104 are unavailable ) ` sees ' an interfering space of dimension @xmath105 .
thus we need to assign a vector linearly independent from @xmath104 ( which itself has dimension @xmath105 ) to the corresponding demanded message of receiver @xmath15 .
this can be possible only if the assigned vectors are of length at least @xmath106 , i.e. the rate can be at most @xmath2 .
therefore , for a rate @xmath58 index coding assignment , all @xmath104 should spanning a space of dimension @xmath107 .
thus we have @xmath108 and @xmath109.$ ] by lemma [ intersectingkdimspaces ] , we should thus have @xmath110    however @xmath94 and @xmath56 are in conflict , which means that they should be assigned linearly independent vectors , i.e. , @xmath111 which means that @xmath112 . thus there is a contradiction and thus any internal conflicts forces the rate to be less than @xmath58 .
this concludes the proof .
from section [ sec3 ] , the following is clear .    *
if there are no conflicts ( not even conflicts between two alignment sets ) in the alignment graph , rate @xmath107 is achievable ( this is the case when any receiver demanding a message has all the other messages as side - information ) . * for rate @xmath107 infeasible index coding problems ,
rate @xmath58 is feasible if and only if there are no internal conflicts .    towards obtaining our main result , which characterises a class of index coding problems which are rate @xmath2 feasible
, we first give a prior known necessary condition for feasibility of rate @xmath2 .
suppose there are @xmath113 messages @xmath114 such that message @xmath115 is demanded by some receiver ( say receiver @xmath116 ) and @xmath117 following @xcite , we call such a set of messages @xmath114 as an _ acyclic subset of messages of size _ @xmath113
. the following theorem can be obtained from the results in @xcite .
[ thmmais ] an index coding problem which is rate @xmath2 feasible can not have an acyclic subset of messages of size @xmath113
.     assigned to the four messages must be linearly independent.,width=192 ]    the following example however shows that theorem [ thmmais ] is not a sufficient condition for rate @xmath2 feasibility .    [ ex : onethirdinfeasible ]    consider a single unicast index coding problem with six messages .
the interfering sets are as follows : @xmath118 . the alignment graph and the conflict hypergraph corresponding to the problem
are given in fig .
[ fig : infeasible ] .
it can be seen that the conflict hypergraph of this problem does not have an acyclic subset of messages of size @xmath113 .
if we assume that the problem is rate @xmath2 feasible , then it has to be necessarily true that both @xmath119 and @xmath120 are two dimensional . we leave it to the reader to check that @xmath121 must also be two dimensional .
however , the conflict between the messages @xmath122 and @xmath123 are resolved only if @xmath121 is three dimensional .
thus , the problem is rate @xmath2 infeasible .      in @xcite , the following theorem about the achievability of rate
@xmath2 was proved ( this is a special case of corollary 9 of @xcite ) .
[ thmnocyclesforks ] consider a rate @xmath58 infeasible index coding problem with no acyclic subset of size @xmath124 if none of its alignment sets have both forks ( a _ fork _ is a vertex connected by three or more edges ) and cycles , then the index coding problem is rate @xmath2 feasible .
as mentioned in @xcite , the condition that there are no alignment sets with both forks and cycles of length @xmath106 means that there is no interfering set @xmath63 at any receiver of size @xmath125 , as such a set would mean that there is both a cycle and fork within an alignment set ( since the messages in @xmath63 are fully connected in the alignment set ) . therefore theorem [ thmnocyclesforks ]
characterises a rather limited class of index coding problems which are @xmath2 feasible .
example [ ex : onethirdfeasible ] shows an index coding problem which does not satisfy the conditions of theorem [ thmnocyclesforks ] but is rate @xmath2 feasible .
[ ex : onethirdfeasible ] consider a single unicast index coding problem with six messages .
the interfering sets of the problem are as follows : @xmath126 . the alignment graph and the conflict hypergraph corresponding to the problem
are given in fig .
[ fig : feasible ] . it can be seen that the alignment graph of this problem has both forks and a cycle .
consider @xmath106 linearly independent @xmath127 vectors @xmath128 .
we note that the following assignment of vectors resolves all conflicts : ( i ) vector @xmath129 to messages @xmath130 and @xmath123 , ( ii ) vector @xmath131 to messages @xmath132 and @xmath133 , ( iii ) vector @xmath134 to messages @xmath135 and @xmath136 .
thus , the problem is rate @xmath2 feasible .      in this subsection
, we develop a new framework for studying the rate @xmath2 feasibility of groupcast index coding problems . towards this
, we define the notions of a triangular interfering set and a type-2 alignment set .    a subset @xmath137 of size three of the set of messages @xmath138
is said to be a _ triangular interfering set _
if all the messages in @xmath137 interfere simultaneously at some receiver , and at least two of the messages in @xmath137 are in conflict .
two distinct triangular interfering sets @xmath139 and @xmath140 are said to be _ adjacent _ if they ` meet ' at a conflicting edge , i.e. , @xmath141 such that @xmath41 and @xmath52 are in conflict .
[ type2align ] two triangular interfering sets @xmath139 and @xmath140 are said to be _ connected _ if there exists a _
( i.e. , a sequence ) of adjacent triangular interfering sets starting from @xmath139 and ending at @xmath140 .
type-2 alignment set _ is a maximal set of triangular interfering sets which are connected to each other .
examples to illustrate the above definitions are shown in fig .
[ fig : triangleinterfset ] .
note that the maximality in definition [ type2align ] means that we can not add another triangular interfering set to a type-2 alignment set and still maintain connectivity ( as in definition [ type2align ] ) .    by definition
, every type-2 alignment set must be a subset of a ( regular ) alignment set , and there could be many type-2 alignment sets within any alignment set . given an index coding problem , we can identify type-2 alignment sets as follows . within any alignment set
, we identify a triangular interfering set of messages ( if there is no such set then there is no type-2 alignment set inside that alignment set )
. then we repeatedly add triangular interfering sets which are adjacent to the existing connected triangular interfering sets .
when we can no longer add such adjacent triangular interfering sets , then we have our type-2 alignment set .
we now prove a necessary condition for rate @xmath2 feasibility based on the vectors assigned to type-2 alignment sets .
this theorem is another application of the key lemma , lemma [ intersectingkdimspaces ] .
[ type2sets2dim ] in any rate @xmath2 solution to a given index coding problem , all the messages in any type-2 alignment set must be assigned vectors from a vector space of dimension two .
let @xmath35 be the encoding function of a rate @xmath2 solution . consider a type-2 alignment set @xmath142 with triangular interfering sets @xmath143.$
]    suppose for any triangular interfering set @xmath144 of @xmath142 , we have @xmath145 .
since all the vertices in @xmath144 interfere at some receiver ( say , a receiver which requests some other message @xmath52 ) , the message @xmath52 must be assigned a vector which is linearly independent from those assigned to the messages in @xmath144 , and thus we need at least @xmath113 linearly independent vectors , and hence the rate has to be @xmath146 .
thus no triangular interfering set @xmath144 of @xmath142 has @xmath145 .
however , any triangular interfering set @xmath144 of @xmath142 _ must _ have @xmath147 , as @xmath144 has a conflict .
suppose @xmath148 .
consider three messages @xmath149 and @xmath150 in @xmath151 that have been assigned three linearly independent vectors , belonging to some three triangular interfering sets ( not necessarily different ) , @xmath152 and @xmath153 respectively .
we have already shown that we can not have @xmath154 .
so at least two of the three triangular interfering sets are different .
suppose all three sets @xmath152 and @xmath153 are different . because the three triangular interfering sets are within the same type-2 alignment set , it must be the case that there exists a path consisting of adjacent triangular interfering sets starting from @xmath155 through @xmath156 and upto @xmath157 let @xmath77 be the number of triangular interfering sets on this path ( counted as we go along the path ; repetitions are counted separately )
. for @xmath158 $ ] , let @xmath159 denote the set of @xmath106 vectors assigned to the @xmath160 triangular interfering set in this path .
[ fig : type2assignment ] illustrates this scenario for the type-2 alignment set example shown in fig . [
fig : type2ex ] .    .
to determine the dimension of the span of these messages , consider the following path of triangular interfering sets : @xmath161 and @xmath162 are the sets of @xmath106 vectors assigned to @xmath160 triangular interfering set.,width=249 ]    by the previous arguments , we have that @xmath163 also , @xmath164,$ ] as the @xmath160 and the @xmath165 triangular interfering sets are adjacent by construction of the path .
therefore , by lemma [ intersectingkdimspaces ] , it must be the case that @xmath166 .
however , the vector assigned to the message @xmath167 belongs to @xmath168 for @xmath169 , and according to our assumption the vectors assigned to these three messages are linearly independent vectors .
thus there is a contradiction , which means that we can not have three messages in three different triangular interfering sets which have been assigned linearly independent vectors .
a similar claim can be proved if the three messages come from two different triangular interfering sets .
thus , no three messages in a type-2 alignment set can be assigned linearly independent vectors .
in other words , any type-2 alignment set @xmath142 in a rate @xmath2 solution must have @xmath170 .
theorem [ type2sets2dim ] is stricter than theorem [ thmmais ] , as theorem [ thmmais ] applies only to an acyclic subset of messages of size @xmath171 which basically is equivalent to a triangular interfering set .
theorem [ type2sets2dim ] on the other hand considers a ` connected component ' of such triangular interfering sets , and is therefore more strict
. we leave it to the reader to verify that the problem in example [ ex : onethirdinfeasible ] , while ` passing ' the condition of theorem [ thmmais ] , ` fails ' the condition of theorem [ type2sets2dim ] .
theorem [ type2sets2dim ] prescribes that type-2 alignment sets must be ` two - dimensional ' in a rate @xmath2 code . in this subsection
, we give a necessary and sufficient condition for achieving this two - dimensionality .
for this purpose , we require the notion of a restricted index coding problem .
let @xmath172 denote an index coding problem with message set @xmath59 .
for some @xmath173 , a @xmath142__-restricted index coding problem _
_ is defined as the index coding problem @xmath174 consisting of    * the messages @xmath142 . *
the subset @xmath175 ( of size @xmath176 ) of the receivers of @xmath172 which demand messages in @xmath142 .
* for each @xmath177 the demand sets @xmath178 and the side - information sets @xmath179 are restricted within @xmath142 , i.e. , @xmath180    the alignment graph and the alignment sets of the restricted index coding problem @xmath174 are called the @xmath142-_restricted alignment graph _ and @xmath142-_restricted alignment sets _ respectively . a @xmath142-_restricted internal conflict _
is a conflict between any two messages within a restricted alignment set of @xmath142 .
the proof of the following theorem is a direct application of theorem [ ratehalfthm ] , and hence is skipped .
[ restrictedratehalfthm ] the restricted index coding problem @xmath174 is rate @xmath58 feasible if and only if there are no @xmath142-restricted internal conflicts .
[ 3cross1restrictedcorr ] for a given index coding problem @xmath172 , there exists an assignment of @xmath181 vectors from a two dimensional vector space ( over a large enough field @xmath3 ) to the messages in a subset @xmath142 such that the conflicts within @xmath142 are resolved , if and only if the restricted index coding problem @xmath174 is rate @xmath58 feasible .
we first recall that conflicts within @xmath142 are said to be resolved if ( [ resolvedconflicts ] ) is satisfied .    _ if part : _ the proof for the if - part follows the achievability scheme shown in the proof of theorem [ ratehalfthm ] , with the difference that we assign to the messages in @xmath142 randomly generated @xmath181 vectors ( not @xmath91 vectors as in theorem [ ratehalfthm ] ) from a two dimensional space over a suitably large field @xmath3 . because @xmath174 is rate @xmath58 feasible , such an assignment resolves the conflicts within @xmath142 .    _ only if part : _ suppose there is an assignment of @xmath181 vectors from a two dimensional vector space to messages in @xmath142 such that all conflicts within @xmath142 are resolved .
then we can always obtain a @xmath182 matrix @xmath6 such that premultiplying all the vectors assigned to @xmath142 with @xmath6 gives us a rate @xmath58 ( length 2 ) index coding solution for @xmath174 .
the following theorem gives a necessary and sufficient condition for assigning vectors from a two dimensional space to the type-2 alignment sets ( i.e. , for satisfying the necessary condition of theorem [ type2sets2dim ] ) .
[ norestrconflictstype2 ] let @xmath142 be a type-2 alignment set of the given index coding problem @xmath172 .
if @xmath172 is rate @xmath2 feasible , then @xmath174 must be rate @xmath58 feasible which holds if and only if there are no @xmath142-restricted internal conflicts .
the proof follows by combining the claims of theorem [ type2sets2dim ] , corollary [ 3cross1restrictedcorr ] and theorem [ restrictedratehalfthm ] .
we now prove the main result of this paper , which connects all the previously proved results and widens the class of index coding problems for which rate @xmath2 is achievable .
because of the framework we have developed in the previous subsections , the proof of this theorem is simpler than theorem [ thmnocyclesforks ] , while also subsuming that result .
[ thm : main ] a rate @xmath58 infeasible index coding problem @xmath172 is rate @xmath2 feasible if every alignment set of @xmath172 satisfies _ either _ of the following conditions .    1 .   it does not have both forks and cycles .
it is a type-2 alignment set with no restricted internal conflicts .
we first make some observations before giving the achievable index coding scheme .    _ observation 1 : _ consider an alignment set @xmath183 of @xmath184 which does not contain both forks and cycles and is also not a triangular interfering set .
we claim that there is no triangular interfering set within @xmath183 .
this is because if there is a triangular interfering set ( say @xmath139 ) then there should be at least one more message which is not included within @xmath139 .
however this would mean that there is both a fork and a cycle within @xmath183 , which is not allowed .
_ observation 2 : _ now suppose there are three messages in @xmath183 such that all of them interfere with a particular receiver .
then there must necessarily be no conflicts in - between the three messages .
this is because if there were conflicts in between the three messages , then these three messages will form a triangular interfering set .
thus no three messages from @xmath183 having at least one conflict in - between themselves can interfere at the same receiver .
_ observation 3 : _ by a similar argument as in observation @xmath107 , if @xmath183 has three messages interfering at a receiver , then it can not have any other message than these three messages ( because if it did , then we would have both cycles and forks within @xmath183 .
_ observation 4 : _ consider an alignment set @xmath185 of @xmath184 which does not contain both forks and cycles , but is also a triangular interfering set .
then , by definition @xmath185 must be a type-2 alignment set .
for the sake of this proof , we consider an alignment set such as @xmath186 as being under the class of type-2 alignment sets .
we thus have three kinds of alignment sets in @xmath172 .    1 .
alignment sets which have no three messages interfering at any receiver .
alignment sets which consists only of three messages , all three interfering at some receiver , without any conflicts in - between .
( these three messages may interfere at other receivers also , but at least one common receiver where they all interfere exists ) .
alignment sets which are also type-2 alignment sets without restricted internal conflicts .     are shown . in the type-2 alignment
set , the receivers at which the triangular interfering sets interfere are not explicitly shown in the figure . since there should not be any restricted internal conflict within a type-2 alignment set , no triangular interfering set interferes at a receiver inside the same type-2 alignment set.,width=249 ]    we now give the achievability index coding scheme by assigning vectors independently to each alignment set of @xmath172 , and the technique of assignment , which follows , depends on the type of the alignment set .
_ alignment set which has no three messages interfering at any receiver : _ for each message in such an alignment set , we assign an independently generated random @xmath181 vector ( over a large field @xmath3 ) .
_ alignment set which consists only of three messages interfering at any receiver without any conflicts in - between : _ we randomly generate a @xmath127 vector and assign it to all the messages in such an alignment set .
_ alignment set which is a type-2 alignment set without restricted internal conflicts : _ let @xmath142 be the type-2 alignment set under concern . for each @xmath142-restricted alignment set , we assign an independently generated random @xmath181 vector from a two dimensional space .
note that this resolves all the conflicts between the messages within @xmath142 , by theorem [ restrictedratehalfthm ] and corollary [ 3cross1restrictedcorr ] .
all the messages fall under one of these alignment sets , and hence all of the messages have been assigned vectors at this point .
let @xmath35 denote encoding function corresponding to this assignment and @xmath48 denote the vector assigned to message @xmath38 .
we now show that this assignment resolves all the conflicts in @xmath172 .
consider a receiver @xmath15 which requests a message @xmath38 .
we proceed on a case by case basis , depending on the size of @xmath63 . for each case
we check whether the condition , @xmath55 ( whp ) , is met .
we call this condition as the _ no conflict condition _ for the sake of this proof .
_ case 1 : _
@xmath187 there are two cases here , either @xmath38 and @xmath63 are in the same alignment set or they are in different alignment sets . if @xmath38 and @xmath63 are in different alignment sets , then @xmath38 and the messages @xmath63 have been assigned independently and randomly generated @xmath181 vectors .
thus the no conflict condition is met .
now , @xmath38 and @xmath63 are in the same alignment set .
this can only be an alignment set where no three messages interfere at any receiver , or a type-2 alignment set . in the former case
, the no conflict condition is met as any three messages in such an alignment set are assigned independent vectors ( whp ) . now if it is a type-2 alignment set , the conflict(s ) between @xmath38 and @xmath63 is(are ) within that type-2 alignment set .
because there are no restricted internal conflicts in any type-2 alignment set , it must be the case that @xmath38 and @xmath63 are in different restricted alignment sets . by our scheme , such conflicts should therefore be resolved .
hence , the no conflict condition is met in this case too .    _
case 2 : _ @xmath188
then we have two cases .
the first case is that no two messages in @xmath63 are in conflict .
this means that @xmath63 must be an alignment set which consists only of three messages without any conflicts in - between , and @xmath38 necessarily belongs to another alignment set . by our scheme , @xmath189 is an one - dimensional space generated by a random vector which is independently generated from the vector assigned to @xmath38 .
hence , the no conflict condition is met .
finally we consider the case when @xmath190 and at least one conflict exists within @xmath63
. then @xmath63 must be within some type-2 alignment set ( say @xmath142 ) .
again we have two sub - cases here , i.e. , @xmath38 and @xmath63 are within the same ( type-2 ) alignment set , or @xmath38 is in a different alignment set than @xmath63 which is within a type-2 alignment set . in the former case ,
the conflicts between @xmath38 and @xmath63 are within the type-2 alignment scheme , which is resolved by our scheme ( by the same arguments as in the last subcase of case 1 ) . in the latter case
, we must have that @xmath189 lies within a two dimensional space ( as it is within a type-2 alignment set ) which is generated independently from the vector assigned to @xmath38 .
hence again the no conflict condition is met .    by the previous arguments
, the no conflict condition is met for any receiver @xmath15 and any demand @xmath38 at @xmath15 .
thus , all the conflicts in @xmath172 are resolved .
this proves the theorem .
in this work , we presented a class of index coding problems for which rate @xmath2 is feasible .
this class of problems is larger than what was previously known .
we believe that the framework developed in this work in order to obtain our results can be leveraged to settle the rate @xmath2 feasibility completely .
in particular , we conjecture that the the necessary condition for rate @xmath2 feasibility of theorem [ norestrconflictstype2 ] is also sufficient .
developing conditions for feasibility of rates of the form @xmath191 would also be an interesting area of future study .
the connection to topological interference management problems follows from @xcite . from @xcite , it is known that the length of scalar linear index codes for the single unicast index coding problem is known to be equal to _ minrank _ of the side - information graph .
thus , our results also imply a class of graphs whose minrank is equal to @xmath106 .
this is a promising result as the general minrank problem is known to be np hard @xcite .
160 y. birk and t. kol ,  coding on demand by an informed source ( iscod ) for efficient broadcast of different supplemental data to caching clients  , ieee transactions on information theory , vol .
52 , no . 6 , june , 2006 , pp .
2825 - 2830 .
a. blasiak , r. kleinberg , and e. lubetzky ,  broadcasting with side information : bounding and approximating the broadcast rate  , ieee transactions on information theory , vol .
2013 , pp .
5811 - 5823 .


the study of ultra high energy cosmic rays ( uhecr ) and neutrinos ( uhe@xmath0s ) is currently a high priority in astroparticle physics with many experimental efforts being dedicated to these two related areas of research .
uhecrs are routinely detected through the extensive atmospheric showers ( eas ) they produce when interacting in the atmosphere . despite the recent advances in the measurement of the flux of uhecrs  @xcite , their primary composition remains unknown  @xcite , and this is one of the main obstacles to extract precise conclusions on their origin .
there are strong reasons to believe that uhe@xmath0s should be produced in the interactions of uhecrs with the material surrounding the sources , and/or in their propagation through the observed cosmic microwave background radiation @xcite .
however their detection has not yet been achieved .
efforts are being made to improve the experimental situation in both fields .
the challenge is to instrument sufficiently large target volumes to compensate for the low cross section and low fluxes for uhe@xmath0 detection , and to find measurements that provide large aperture at the highest energies and help to constrain the primary composition of uhecrs .
the radio technique is being explored in both fields .    as early as in 1962 , g.  askaryan proposed to detect uhecrs and uhe@xmath0s by observing the coherent radio pulse from the excess of electrons in a shower developing in a dense , dielectric and nonabsorptive to radiowaves medium  @xcite .
soon after , pulses were observed in coincidence with air shower arrays  @xcite .
the emission is coherent in wavelengths which are large compared to the characteristic size of the electric charge and current distributions associated to the induced showers .
the technique has been receiving a lot of attention in the field of astroparticle physics in the last decade because of the relatively low cost of the antennas needed for the detection systems .
also , coherence implies that the pulse energy scales with the square of the primary energy which favours long range detection , a requirement to achieve large areas and volumes both for uhecr and uhe@xmath0 detection .
indeed , quadratic scaling has been confirmed in accelerator experiments   @xcite leading to very strong pulses associated to uhe energy showers .
a large number of initiatives have been made or are currently in development or planning stages .
in dense media neutrino - induced showers are less than a meter in width and full coherence is expected to be mantained up to the ghz range ( see for instance @xcite ) .
a variety of past and present experiments search for pulses at those frequencies produced by neutrinos in ice  @xcite , the moon regolith  @xcite or salt domes  @xcite .
the largest and most promising experiments are in planning stages looking at ice  @xcite .
electromagnetic pulses emitted in the mhz - ghz frequency range by eas induced by uhecr are being measured in the hope of using this complementary information to constrain its composition and/or to develop new cost effective detection systems  @xcite .
in addition the anita balloon flown antennas , initially devised to search for neutrinos , has recorded coherent pulses up to the ghz range  @xcite consistent with eas induced by uhecrs and there are plans for future developments  @xcite .
the success of the radio technique requires an accurate and computationally efficient calculation of the radio emission properties of uhe showers .
it is necessary to perform this in an efficient way since the number of particles in a shower at eev energies is @xmath1 .
simulation techniques for evaluating pulses in dense media have been used for more than 20 years  @xcite .
the shower is simulated to obtain , for every particle track , the information needed to calculate its contribution to the electric field .
the formula used for this purpose stems from an approximate solution of maxwell s equations in the fraunhofer limit  @xcite adapted for simulation purposes in  @xcite .
the electric field due to a short charged particle track , assumed to be travelling at a constant speed , can be approximated by two terms which correspond to the start and end of the track .
the resulting electric field depends on the particle speed and the angle between the line of sight from the track to the observer  @xcite . to calculate the emission in a shower , tracks
are chosen so that the particle velocity can be approximated to be constant , and all track contributions are added taking into account interference effects .
although the approximate formula is sufficient for many practical applications , its range of validity is limited in frequency and position of the observer with respect to the track .
it is possible to extend its range of validity by subdividing each track in sub - intervals .
the resulting algorithm ( often referred to as  zhs algorithm " ) is easy to implement and fast enough for the simulation of particle showers .
this is convenient since it has allowed the simulation of pulses in the fresnel region for neutrino detection  @xcite and in measurements of eas  @xcite .
although the range of applicability of the algorithm is enhanced when used in this way , it is not obvious that the sum of the subcontributions correclty accounts for all the radiation in regions close to the emission source .
the object of this article is to study the conditions for the zhs algorithm to be valid and to establish its accuracy . for this purpose
we obtain in section ii an exact solution to the problem of a charged particle instantaneously accelerating to a constant speed and stopping abruptly after a discrete time interval @xcite .
this allows us to establish the precise conditions necessary to turn this solution into the basic formula of the zhs algorithm . in section iii
we compare the exact solutions for an infinite and a finite track to the result of the zhs algorithm , and stress the compatible interpretation of the radiation regime in terms of cherenkov radiation for both cases . in section
iv we compare the exact solution of the single track problem in nearby regions with the results of applying the zhs algorithm with track subdivisions to test its validity and accuracy .
section v is devoted to discussing the zhs algorithm in relation to other approximations made to calculate the pulses emitted from showers .
section vi presents the summary and conclusions of our work .
let us assume an electron ejected from an atom at time @xmath2 that travels at constant speed @xmath3 through a medium along a finite track until it is absorbed by another atom at time @xmath4 . neglecting the movement of the atoms
, we can model the electric current associated to the electron as , @xmath5 where @xmath6 is the charge of a positron , @xmath7 is its position and @xmath8 an arbitrary reference position
. the step @xmath9-functions account for the fact that the electron only moves in the time interval ( @xmath10 , @xmath4 ) .    in a dielectric medium with permitivity @xmath11 and magnetic susceptibility @xmath12 , maxwell s equations for the vector potential in the frequency domain @xmath13
can be written as  @xcite : @xmath14 = -\mu\mathbf{j}(\mathbf{x},\omega ) , \end{split } \label{eq : avec}\ ] ] where @xmath15 is the fourier transform of the scalar potential and we use the following convention for the fourier - transform of a function @xmath16 : @xmath17
. in principle @xmath11 and @xmath12 can depend on frequency and our results below would be equally valid , but we drop the explicit dependence of @xmath11 and @xmath12 with @xmath18 for simplicity .
we use the lorenz gauge condition which implies that eq .
( [ eq : avec ] ) for the vector potential becomes : @xmath19 with @xmath20 and @xmath21 the refractive index of the medium .
the lorenz gauge condition in the frequency domain  @xcite : @xmath22 implies that the scalar potential for non - zero frequencies is entirely determined by the divergence of the vector and only the vector potential is needed to calculate the field .      the solution of the helmholtz equation in eq .
( [ eq : helmholtz ] ) is standard physics and can be obtained , using green s method , as an integral over source positions  @xcite : @xmath23 where from now on , @xmath24 will denote the position of the observer and @xmath25 the position of the source .
for simplicity we rewrite the current in eq .
( [ eq : current ] ) as : @xmath26 with @xmath27 ; @xmath28 and @xmath29 , where we assume without loss of generality that the electron travels along the @xmath30-axis which is parallel to @xmath31 .
the equations that follow below are clearly valid for any continuous and differentiable functions @xmath32 and @xmath33 .
they are also valid for the case considered in eq .
( [ eq : current ] ) which can be obtained as a limiting case of suitable continuous and differentiable functions .    transforming this definition of the current to the frequency domain and substituting into eq .
( [ eq : avec_w ] ) for the vector potential we obtain : @xmath34    in the lorenz gauge , the only non - zero component of the vector potential for a charge moving along @xmath30 is the @xmath30-component . the scalar potential @xmath35 is obtained through eq .
( [ eq : gauge ] ) . for this purpose
we need to obtain @xmath36 : @xmath37 \frac{(z - z')}{|\mathbf{x } - \mathbf{x'}| } \nonumber\end{aligned}\ ] ]    the electric field @xmath38 can be obtained from the scalar and vector potentials @xcite as : @xmath39    due to the cylindrical symmetry of a charged particle track we calculate only the radial ( @xmath40 ) and @xmath30 ( @xmath41 ) components of the field . since @xmath42 only has @xmath30-component , then the radial component of the field is given by :    @xmath43    where @xmath44 is the radial coordinate , and @xmath45 denotes @xmath46 .
the @xmath47component of the field is given by : @xmath48 with @xmath49 denoting @xmath50 .
performing the derivatives in eqs .
( [ field_rho ] ) and ( [ field_z ] ) we obtain : @xmath51 \end{split } \label{eq : rhofield}\ ] ] and ,    @xmath52 ~+~   \\ & i\omega~\frac{\mu}{4\pi}~qv~\int_{t_1}^{t_2 }   \mathrm{d}t'~e^{i\omega t'}~\frac{e^{ikr}}{r } \end{split } \label{eq : zfield}\ ] ]    where @xmath53 and @xmath54 are both functions of the source time @xmath55 which is defined as , @xmath56 and @xmath57    eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) provide an exact solution for the electric field of a finite track .
in general they do not have an analytical form and a numerical integration needs to be performed to obtain the field . for the calculations in this paper , we have divided the integration interval and applied simpson s rule , increasing the number of divisions until the integral converged . under certain conditions ,
however , we can give analytical approximations for relevant physical situations .    in the following
we show that the basic expression used in the zhs algorithm is a particular case of eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) under certain approximations .
a simple expression for the approximate calculation of the electric field from a single charged particle track moving at constant speed was found in @xcite . in this section
we derive the expression used for the zhs algorithm from the exact solution and compare the electric field as obtained in both the exact calculation and with the zhs expression .
this allows us to establish under which circumstances the formula gives a good account of the electric field
.    the zhs algorithm can be obtained from eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) if the following set of conditions are fulfilled :    1 .
the observer is in the  far field " zone i.e. @xmath58 2 .
the fraunhofer approximation holds .
this can be stated as a condition for the phase factor to be approximated as : @xmath59   \label{eq : fraunhofer}\ ] ] where @xmath60 is the distance from the observation point to a reference point along the track where the particle is located at a reference time @xmath61 , and @xmath62 is the angle between the particle track and the direction from the reference point to the observer .
this approximation holds provided the parameter @xmath63 with @xmath64 defined as : @xmath65 ^ 2}{r}\sin^2\theta ,   \label{eq : cond3}\ ] ] this condition should be fulfilled at any time @xmath66 from @xmath10 to @xmath4 .
a more commonly used and nearly equivalent form of this condition  @xcite is : @xmath67 where @xmath68 is the length of the track .
this condition is necessary to ensure that the second and higher order terms for the phases @xmath69 in eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) have no relevance , even when the sum of the leading and first order terms in the taylor expansion of the phases is zero , as it occurs for observation at the cherenkov angle @xcite defined as @xmath70 with @xmath71 .
3 .   finally
the distance to the observer appearing in the denominators of several terms of eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) must be approximated as : @xmath72 over the length @xmath73 of the track , where @xmath60 is the distance to a reference point along the track ( in the zhs algorithm the mid - point of the track is selected ) .
the error when making this approximation is of order @xmath74 .    in particular , the condition in eq .
( [ eq : cond1 ] ) implies : @xmath75    with these approximations the radial component of the field in eq .
( [ eq : rhofield ] ) becomes : @xmath76 where we have used : @xmath77 and @xmath78    eq .  ( [ eq : rhozhs ] ) can be easily integrated yielding : @xmath79 .
\end{split}\ ] ] if we make @xmath80 this becomes the expression for the radial field as used in the zhs algorithm @xcite except for a factor 2 due to the fourier transform convention used in @xcite .    similarly applying the approximations in eqs .
( [ eq : fraunhofer ] ) , ( [ eq : cond4 ] ) and ( [ eq : b ] ) to eq .
( [ eq : zfield ] ) for the @xmath30-component of the field , and using that @xmath81 , it is straightforward to show that : @xmath82 which can be cast as : @xmath83    performing the integral , and taking @xmath80 , the zhs formula is recovered :    @xmath84~. \end{split}\ ] ]      to calculate the electric field of the pulse emitted from a current distribution , such as that produced in a high energy shower , the zhs algorithm uses the zhs expressions in eqs .
( [ eq : rhofield_zhs ] ) and ( [ eq : zfield_zhs ] ) to calculate the emission from all the charged particle tracks .
the final result is obtained adding up all the contributions .
the value of @xmath60 used for the phase factor in each particle track is the distance from the first point of the track to the observation point .
this definition is consistent with the convention to account for the phase change between emission arising from the start and end points of the track .
the actual value of @xmath60 used for the denominator is the distance between the midpoint of the track and the observer which is for all practical purposes the same value when the approximation is valid .
the algorithm used in alternative simulation programs is similar but different in these technical details  @xcite .
naturally it is possible to divide any charged particle track into arbitrarily small subtracks in order to make the computation more accurate .
this will extend the range of validity of the approximation .
it is interesting to note that it makes absolutely no difference to subdivide the track of a uniformly moving charge when observed in the fraunhofer limit because the term associated to the end of one subtrack cancels the term associated to the beginning of the next subtrack .
this is because in this limit the differences in @xmath60 between adjacent subtracks are arbitrarily small .
however in practical situations this cancelation is not exact if @xmath60 is allowed to change from a track to the next one according to geometry . in the original zhs program  @xcite the tracks were not subdivided but soon it was realized that a more accurate result was obtained by subdividing all tracks at every point there was a discrete interaction in the simulation program  @xcite .
this was the main modification that was ever made to the original zhs monte carlo @xcite and has been effective since then . with this subdivision
the distribution of tracks for a shower in ice has a peak at about @xmath85 mm .
the exact electric field for a charged particle track derived in eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) must account for every single feature of the electric field , since no approximations have been made .
in particular , it must reproduce cherenkov radiation which is the only radiation emitted by a charged particle moving at constant speed , when @xmath86 in the limit of an infinite track .
an analytical solution for the electric field produced by such particle can be obtained and it is given in @xcite(chapter 4 ) .
the radial @xmath87 and @xmath47components of the fields given in @xcite converted to si units , and using the fourier transform convention adopted in this work , can be written as : @xmath88 @xmath89 where @xmath87 is the radial coordinate , @xmath90 and @xmath91 are the modified bessel functions of the second kind , and @xmath92 is a function that can take two different values depending on the magnitude of the particle speed , @xmath93 or @xmath86 ( subluminal or superluminal regime ) . in our convention : @xmath94 @xmath95 if the particle travels below the speed of light in the medium , the argument @xmath96 of the bessel functions is real and the particle does not radiate as shown in @xcite . on the contrary , if the speed of the particle is larger than the speed of light , then @xmath96 is imaginary and the particle radiates .
the latter case corresponds to pure cherenkov radiation @xcite .    with the help of the asymptotic forms for the bessel functions ,
we can obtain the limits of eqs .
( [ eq : rhofield_inf ] ) and ( [ eq : zfield_inf ] ) when @xmath97 i.e. for small distances to the track compared to the radiation wavelength .
conversely we can also obtain the limits when @xmath98 i.e. for large distances compared to the wavelength . if @xmath97 the @xmath91 bessel function dominates over the @xmath90 and only the radial component of the field @xmath40 matters . in this case ,
@xmath99 and a @xmath100 dependence with distances is obtained , as well as no dependence with frequency .    if @xmath98 and keeping in mind that when @xmath86 the argument @xmath101 is imaginary the fields can be written as : @xmath102 @xmath103 in this case the field is proportional to @xmath104 .
this is in agreement with @xcite where the same behavior is deduced using simple arguments of energy conservation through a cylindrical surface surrounding the track .    in fig .
[ fig : inf_vs_finzhs ] the fourier components of the modulus of the electric field for an infinite track as obtained from eqs .
( [ eq : rhofield_inf ] ) and ( [ eq : zfield_inf ] ) are shown , for different frequencies , as a function of @xmath87 , the radial distance to the track . the particle speed is @xmath105 travelling in homogeneous ice with refractive index @xmath106 . under these circumstances
the quantity @xmath107 can be approximated as : latexmath:[\[\label{eq : u }     at large distances to the track when @xmath98 the fields shown in fig .
[ fig : inf_vs_finzhs ] scale with distance as @xmath109 and with frequency as @xmath110 , in agreement with the aysmptotic field components in eqs .
( [ eq : limit_rho_urho_large ] ) and ( [ eq : limit_z_urho_large ] ) .
this behavior takes place when @xmath111 , 1 and 10 m for frequencies @xmath112 ghz , 100 mhz and 10 mhz respectively , in agreement with eq .
( [ eq : u ] ) , as can be seen in fig .
[ fig : inf_vs_finzhs ] .    as the distance to the track decreases and the condition @xmath97 starts to be valid
, the field behaves as @xmath100 and does not depend on frequency as expected from eq .
( [ eq : limit_urho_small ] ) . as can be clearly seen in fig .
[ fig : inf_vs_finzhs ] , the transition from the @xmath100 behaviour to @xmath104 occurs at a distance that depends on frequency because @xmath107 involves the frequency ( eq .  [ eq : u ] ) .
for instance at distances @xmath113 m the condition @xmath97 applies for both @xmath114 and 10 mhz .
the fourier component of the field scales with @xmath100 and has the same value for the two frequencies as seen in fig .
[ fig : inf_vs_finzhs ] , while this is not the case for a frequency of @xmath112 ghz .    in fig .
[ fig : inf_vs_fin ] the modulus of the field for a charged particle in an infinite track - as obtained from eqs .
( [ eq : rhofield_inf ] ) and ( [ eq : zfield_inf ] ) - is shown as a function of frequency for an observer at a fixed radial distance . at large enough frequencies
so that the condition @xmath115 applies , the field scales as @xmath110 as expected from eqs .
( [ eq : limit_rho_urho_large ] ) and ( [ eq : limit_z_urho_large ] ) , while it is constant with frequency for small enough frequencies so that @xmath116 as predicted from the asymptotic eq .
( [ eq : limit_urho_small ] ) .
more quantitatively , since the observer in fig .
[ fig : inf_vs_fin ] is located at @xmath117 m , the field should behave as @xmath110 for @xmath118 mhz ( applying eq .
( [ eq : u ] ) ) .
this is approximately the case as can be seen in fig .
[ fig : inf_vs_fin ] .
the result of the exact calculation for a track of length @xmath119 km is also shown in figs .
[ fig : inf_vs_finzhs ] and [ fig : inf_vs_fin ] .
the agreement between both calculations is excellent because the finite track is long compared to observation distance .
this confirms that eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) must also account for cherenkov radiation . as can be appreciated in fig .
[ fig : inf_vs_fin ] , the exact results for finite and infinite tracks differ for wavelengths larger than the length of the track - frequencies typically below @xmath120 with @xmath121 which gives @xmath122 mhz .
we have numerically evaluated the exact expressions for the @xmath30 and @xmath87 components of the electric field in eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) at different frequencies and observer distances , for a single tracks of different lenghts , and for the tracks constituting a shower in ice as obtained in full simulations performed with the zhs monte carlo code @xcite . in this section
we present the results of this comparison .
the applicability of the zhs expressions in eqs .
( [ eq : rhofield_zhs ] ) and ( [ eq : zfield_zhs ] ) relies on the conditions @xmath123 in section  [ sec : zhs ] .
( [ eq : cond3 ] ) and ( [ eq : cond4 ] ) are easily fulfilled for any wavenumber @xmath124 simply by dividing the particle track in a sufficiently large number of subtracks .
once this is guaranteed , the zhs expression is applied to every sub - track and the electric field is obtained adding the corresponding contributions .
the validity of this procedure will be numerically confirmed below when comparing the exact calculation with the electric field obtained using the zhs algorithm in the manner just described .
however the condition @xmath125 , eq .
( [ eq : cond1 ] ) , does not depend on the size of the track and can not be enforced by applying the procedure outlined above .
as a consequence @xmath125 is an intrinsic limit to the range of observing frequencies and distances in which the zhs algorithm gives accurate results as will be shown in the following .    in fig .
[ fig : single_track_1 mm ] we compare the fourier components of the electric field modulus for a single particle track as obtained with the exact calculation , eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) , to that obtained with the zhs algorithm .
the length of the track is chosen to be small @xmath126 m ( close to the peak value of the distribution of track lengths in the standard zhs code ) . to test the validity of the zhs algorithm , we have calculated the spectra for observers at different distances ( @xmath60 ) measured with respect to the center of the track and placed at the cherenkov angle
.    the condition @xmath127 in ice with refractive index @xmath106 can be cast as : @xmath128 for observers at distances @xmath129 and @xmath130 m from the particle track , the condition in eq .
( [ eq : kr_ice ] ) is fulfilled as long as @xmath131 mhz and 1 ghz respectively .
the zhs algorithm is expected to reproduce the results of the exact calculation in this range .
this can be seen in fig .
[ fig : single_track_1 mm ] .
the relative difference between the zhs and exact calculations is less than @xmath132 at the frontier of the validity range in the explored frequency and distance space .
the accuracy can be however orders of magnitude better .
if @xmath133 is enforced for instance the corresponding relative difference is below @xmath134 .
this condition is satisfied for @xmath135 m and @xmath136 mhz what corresponds to a range of frequencies and distances typically encountered in experiments that search for neutrino induced radio transients .
what is striking is that these conclusions also apply in the `` near field '' provided the track is subdivided in sufficiently small sub - tracks
. this can be seen in fig .
[ fig : single_track_1 m ] comparing the exact and zhs results for a @xmath119  m track .
the same range of validity is obtained because conditions  ( [ eq : fraunhofer ] ) , ( [ eq : cond3 ] ) and ( [ eq : cond4 ] ) are guaranteed by reducing the length of the sub - tracks .    for fixed frequencies the condition in eq .
( [ eq : kr_ice ] ) only applies at sufficiently large distances @xmath60 to the track .
this is illustrated in fig .
[ fig : single_track_distance ] for a track of length @xmath119  m. the fourier components at frequencies of @xmath137 , 100 mhz and 1  ghz are in agreement with the exact calculation respectively at @xmath138 , 1 and 0.1  m as expected .
since the typical distance between antennas in experiments such as the askaryan radio array ( ara ) @xcite is @xmath139 m , we expect the results of the zhs algorithm obtained through the procedure outlined above , to be accurate enough in most practical situations .
it has been questioned whether the zhs algorithm reproduces cherenkov radiation from a single charged particle track  @xcite . in figs .
[ fig : inf_vs_finzhs ] and [ fig : inf_vs_fin ] the algorithm is shown to be in very good agreement with the exact solution for a 1.2 km track as long as @xmath127 as explained above .
the emission from such track is in turn practically equivalent to the cherenkov emission from an infinite track ( also displayed ) and it follows that the zhs calculation must account for cherenkov radiation .      as can be seen in fig .
[ fig : single_track_1 mm ] ( solid lines ) the exact solution for the modulus of the electric field scales linearly with frequency provided that @xmath140 . as a result the electric field can be approximated with the zhs expressions , eqs .
( [ eq : rhofield_zhs ] ) and ( [ eq : zfield_zhs ] ) . for an observer close to the cherenkov angle as in fig .
[ fig : single_track_1 mm ] , the factor @xmath141 and the term in brackets in eqs .
( [ eq : rhofield_zhs ] ) and ( [ eq : zfield_zhs ] ) can be approximated by : @xmath142 \approx \\ &   \frac{1 + i(\omega-\mathbf{k}\cdot\mathbf{v})t_2 - 1 - i(\omega-\mathbf{k}\cdot\mathbf{v})t_1}{i(\omega-\mathbf{k}\cdot\mathbf{v } ) }   = t_2 - t_1   \end{split}\ ] ] making the dependence of the field with @xmath18 apparent .
clearly the field can not grow indefinetely with frequency .
this apparent  ultraviolet divergence " is only an artifact of considering the unrealistic medium in which the permittivity @xmath143 is constant with frequency . in a physical medium absorption at high frequencies will tame the growth with frequency of the electric field .    when @xmath144 the field behaves with frequency as @xmath145 as can be also seen in fig .
[ fig : single_track_1 mm ] . in the model of a charged particle at rest for @xmath146 , moving with a speed @xmath147 between @xmath2 and @xmath148 , and becoming again at rest for @xmath149
, the coulomb field dominates at small distances to the track and/or low frequencies .
the field can be modeled as : @xmath150 for @xmath151 and @xmath152 for @xmath153 , where @xmath154 and @xmath155 are respectively the distances from the charge to the observer at times @xmath10 and @xmath4 , and @xmath156 and @xmath157 denote the instants of time at which the coulomb field arrives at the observer .
the fourier transform of a heaviside function at non - zero frequency is proportional to @xmath145 , and the two coulomb fields interfere coherently at low frequencies , explaining the frequency dependence of the electric field .
naturally the zhs algorithm does not reproduce this behavior which is not associated to radiation .
the growth of the field at low frequencies is an artifact of not accounting for screening of the field by the atoms in the medium .      in fig .
[ fig : single_track_distance ] the dependence of the fourier components of the field modulus with distance is shown for several frequencies . at sufficiently large distances to the track the electric field behaves as @xmath158 for all frequencies .
this is the radiation zone , the field is expected to behave as @xmath158 as explained in conventional radiation theory @xcite .
if the observer is placed at very small distances compared to the length of the track , the situation resembles that of an infinite track . in this case the discussion in section  [ sec : cherenkov ] applies .
the field behaves as @xmath159 at distances much smaller than the length of the track provided the frequency is high enough to satisfy @xmath115 - see eq .
( [ eq : u ] ) . this can be seen in the curve of the fourier component at @xmath112 ghz for @xmath160  m. at small distances and sufficiently low frequencies , when @xmath116 , the field becomes proportional to @xmath158 and independent of frequency - see eq .
( [ eq : zfield_inf ] ) .
this feature can be appreciated in fig .
[ fig : single_track_distance ] at distances below 0.1  m for the calculations at 10 and 100 mhz .
the zhs algorithm reproduces the calculation provided the @xmath127 condition is satisfied as could be expected , reproducing both the @xmath158 and @xmath159 behaviors .
it is possible to test the zhs algorithm in a more realistic situation . in this section
we compare the fourier components of the electric field predicted by the zhs algorithm with those obtained using the exact calculation in a full simulation of electromagnetic showers .    for this purpose
we have applied the exact solutions of the field of a track given in eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) in the zhs monte carlo code @xcite for the simulation of electron and photon - induced showers in ice .
the zhs monte carlo calculates the start and end points of small sub - tracks of all charged particles ( electrons and positrons ) in an electromagnetic shower down to a kinetic energy threshold of @xmath161 kev . with these we can calculate the exact electric field produced by each single sub - track and add the fields up accounting for interference between different tracks . since eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) are only valid for a charged particle travelling along the @xmath30 axis ( parallel to the shower axis ) , we perform the necessary rotations of eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) to obtain the field for a particle track moving along an arbitrary direction .
simultaneously with the exact calculation , we also obtain the field as predicted by the zhs algorithm for exactly the same shower ( i.e. the same set of tracks and sub - tracks ) . as explained above the subdivisions are such that the conditions in eqs .
( [ eq : fraunhofer ] ) , ( [ eq : cond3 ] ) and ( [ eq : cond4 ] ) are fulfilled for all the sub - tracks in the shower .
the result is qualitatively the same as in the case of single tracks .
as long as the condition @xmath127 is fulfilled , the zhs algorithm gives an accurate prediction for the fourier components of the electric field with a difference of less than a few percent relative to those obtained with the exact calculation .
as can be seen in fig .
[ fig : zhs_10tev ] this occurs for distances to the shower axis as small as @xmath162 m and frequencies above @xmath163 mhz , well in the distance and frequency ranges relevant for experiments looking for particle shower induced radio pulses in dense media @xcite .
we stress here that the accuracy reported above refers to the approximation of using the zhs formula applied to the standard subdivision of tracks in the zhs code , instead of the exact expression for the radiation emitted by the same particle sub - tracks . by comparing the results obtained in the fraunhofer limit with the standard subdivision of tracks to those obtained with a much finer subdivison
, it was determined that the accuracy of the zhs code is @xmath164 at frequencies @xmath165  ghz , improving significantly at lower frequencies .
we do not further address this uncertainty in this paper , nor the uncertainty due to the shower simulation itself .
it is also worth remarking that in terms of computing time the exact calculation is roughly a factor @xmath165 slower than the calcuation performed with the zhs algorithm .
since the zhs algorithm can only be applied in a limited range of frequencies , an accurate representation of the electric field in the time - domain can not be obtained with an inverse fourier transform of eqs .
( [ eq : rhofield_zhs ] ) and ( [ eq : zfield_zhs ] ) . the low frequency components which do not satisfy the condition
@xmath127 are not accurately described by the zhs algorithm as shown before . also at very high frequencies
the number of steps in which the tracks have to be divided in order to fulfill eqs .
( [ eq : fraunhofer ] ) , ( [ eq : cond3 ] ) and ( [ eq : cond4 ] ) can become prohibitively large from the computational point of view what can compromise the calculation of the pulses below the 100 picosecond scale . in practice however these techniques , although not exact , provide fast and accurate calculations in the region of interest to uhe neutrino detection . for experiments with typical time resolutions of the order of 1 ns , and which are only sensitivity to frequencies from @xmath166 mhz to few ghz ,
the zhs algorithm has been shown to give a very accurate representation of the fourier components of the electric field ( fig .
[ fig : zhs_10tev ] ) .
several calculations of the field emitted in showers developing in dense media can be found in the literature . in @xcite
the finite diference time domain method is used for calculating the field of a pancake - like shower with a gaussian longitudinal development and gaussian radial profile in the time - domain which is then transformed to the frequency - domain . in  @xcite using the saddle - point approximation , an analytic equation for the calculation of the electric field of a charge distribution exhibiting a longitudinal profile with a well - pronounced maximum
is derived .
the result is factorized into an integral accounting for the longitudinal variation of the charge and a form factor that accounts for the lateral spread of the shower , a procedure revisited in @xcite for realistic showers . assuming a gaussian longitudinal and lateral development for the charge distribution both results
were directly compared and turned out to be in good overall agreement as shown in @xcite .
minor differences could be attributed to the form factors used .    with the exact calculation of the electric field performed in this work , the field due to a gaussian profile
can also be obtained and compared to the calculations mentioned before .
the electric current for a shower with a gaussian profile is given by eq .
( [ eq : current ] ) with the following replacements : @xmath167 and @xmath168 the shower develops in the longitudinal direction parallel to the @xmath169 coordinate ( shower axis ) , and radially along the @xmath170 and @xmath171 coordinates .
@xmath172 is a normalization constant .
@xmath173 characterizes the width of the shower along shower axis and @xmath174 the corresponding lateral width . after substituting this current in eq .
( [ eq : vectorpotential ] ) , and following the same steps as in section  [ sec : exact ] , the expression for the field is the same as in eqs .
( [ eq : rhofield ] ) and  ( [ eq : zfield ] ) with the following change : @xmath175 eqs .
( [ eq : rhofield ] ) and ( [ eq : zfield ] ) after the changes in eq .
( [ eq : field_gaussian ] ) can be solved numerically .
the zhs algorithm can also be applied to this situation as long as the condition @xmath125 for all distances @xmath176 to the shower .
the procedure consists on slicing the volume occupied by the bulk of the shower in small cubes and approximating each as a track with constant charge given by the gaussian distributions in eq .
( [ eq : current_gaussian ] ) .
comparison of the result of the exact calculation in this work ( or the zhs algorithm ) with the saddle - point calculation in @xcite requires knowing the form factor @xmath177 for a gaussian profile .
@xmath177 is defined in @xcite as : @xmath178 with @xmath179 with @xmath180 the radial position of the observer .
also , @xmath181 with @xmath182 .
@xmath60 is the distance from the maximum of the shower to the observer .
the function @xmath183 represents a normalized charge density of the travelling pancake . assuming a gaussian for @xmath183 of the form : @xmath184 and substituting @xmath183 into eq .
( [ eq : form_factor ] ) , the form factor reads , @xmath185 setting @xmath186 , @xmath187 m , @xmath188 m and @xmath189 m with the refractive index of ice @xmath106 , the electric fields for a gaussian charge profile were calculated using the saddle - point approach as described in @xcite , and compared to the exact formula given in this work .
the results are shown in fig .
[ fig : comparison ] .
since @xmath60 is large the condition @xmath127 is satisfied for frequencies above @xmath190 mhz ( see eq .
[ eq : kr_ice ] ) and the agreement between the zhs and the exact calculation presented in this work is very good ( not shown in fig .
[ fig : comparison ] for clarity ) .
the saddle - point approach is also in very good agreement with both the exact and zhs calculations .
the results of the field for a gaussian charge profile are also in good agreement with those obtained in  @xcite using the finite difference time domain method .
we have obtained the results of an exact calculation of the electric field produced by a charged particle that accelerates instantaneously moves at constant speed and intantaneously deccelerates .
the results are used to obtain the approximate expression used in the zhs algorithm to calculate the radio emission from showers in dense media using shower simulations .
this allows the precise determination of the conditions necessary for this approximation to be valid , namely , the observer must be in the far field zone @xmath127 and the frauhofer approximation must apply , @xmath191 for each track .
the exact and simulated results are compared for tracks in a variety of circumstances to illustrate both the behavior of the fields and the validity of the approximation made in different frequency and distance ranges .
we have shown that the range of validity of the expressions can be greatly enlarged by subdividing long charged particle tracks in smaller sub - tracks and adding up the conributions , as it is done in the zhs algorithm .
nevertheless the far field condition @xmath127 has been shown to be an intrinsic limit of the zhs algorithm .
the zhs algorithm has been tested for long tracks in regions where the fraunhofer approximation does not hold comparing it to the exact solution .
the results clearly indicate that the zhs algorithm reproduces the exact behavior provided the far field condition is satisfied and the length of the sub - tracks is small enough for the fraunhofer approximation to be valid . by comparing the emission with the exact emission from an infinite track it is shown that the emission calculated with the zhs algorithm does indeed contain what is conventionally described as cherenkov radiation .
the precision of the zhs algorithm is shown to be below the @xmath192 level provided that @xmath193 , corresponding to @xmath135  m and @xmath194  mhz . by enforcing @xmath133
the precision improves to better than @xmath195 , what corresponds to @xmath135  m and @xmath196  mhz , which are conditions met in most experimental arrangements trying to detect radio pulses form neutrinos interacting in dense media .
the zhs algorithm is tested for completitude when applied to a shower simulation .
this is done comparing the result of the zhs algorithm to that obtained when the zhs expression is replaced by the exact calculation for every charged particle sub - track for exactly the same shower .
the results indeed confirm that the accuracies reported for individual sub - tracks are approximately mantained in the final zhs result .
finally , in order to compare to alternative calculations the results are compared to the saddle point approximation .
this approximation has been used to test solutions using the method of finite differences to solve maxwell s equations directly in the time domain ( fdtd ) using a simplified shower front based on gaussian distributions in the shower plane and in time .
the comparison of the saddle point approximation to both the exact solution and the zhs algorithm give compatible results confirming that both approaches reproduce the radiation emitted from showers .
the results presented here in summary confirm that the zhs algorithm can be used to describe most practical applications to detect pulses emitted from high energy showers produced in dense media by neutrinos .
the approach only begins to show significant discrepancies when the observer is at distances comparable to the lateral dimensions of the shower ( @xmath197  1  m in ice ) .
since the typical distance between antennas in experiments such as the askaryan radio array ( ara ) @xcite is @xmath139 m , we expect the results to be accurate enough in most practical situations .
the results obtained are also of interest for the application of the zhs algorithm to calculate pulses from eas .
indeed the application of the method of track subdivision has been applied in zhaires @xcite , a recent code developed for this purpose , by making track subdivisions that are forced to satisfy the fraunhofer condition .
j.a - m , w.r.c . , d.g .- f . and e.z .
thank xunta de galicia ( incite09 206 336 pr ) and consellera de educacin ( grupos de referencia competitivos  consolider xunta de galicia 2006/51 ) ; ministerio de educacin , cultura y deporte ( fpa 2010 - 18410 , fpa2012 - 39489 and consolider cpan - ingenio 2010 ) ; aspera ( pri - pimasp-2011 - 1154 ) and feder funds , spain .
we thank cesga ( centro de supercomputacin de galicia ) for computing resources .
part of this research was carried out at the jet propulsion laboratory , california institute of technology , under a contract with the national aeronautics and space administration .
[ hires collaboration ] phys .
* 100 * , 101101 ( 2008 ) j. abraham _
[ pierre auger collaboration ] , phys .
* 101 * , 061101 ( 2008 ) ; j. abraham _ et al . _
[ pierre auger collaboration ] , physics letters b * 685 * , 239 ( 2010 ) .
j. abraham _
[ pierre auger collaboration ] , phys .
lett . * 104 * , 091101 ( 2010 ) r.u .
[ hires collaboration ] phys .
lett . * 104 * , 161101 ( 2010 ) k. kotera , a. olinto , jcap * 10 * , 013 ( 2010 ) .
askaryan , _
soviet physics _
jetp * 14,2 * 441443 ( 1962 ) ; * 48 * 988990 ( 1965 ) .
jelley , nuovo cimento x * 46 * , 649 ( 1966 ) .
h. r. allan , in : j. g. wilson , s. a. wouthuysen ( eds . ) , progress in elementary particle and cosmic ray physics , north holland , ( 1971 ) , p. 169 .
d.  saltzberg _ et al .
_ , phys .  rev
.  lett .  * 86 * , 2802 ( 2001 ) .
p.  miocinovic _ et al .
d * 74 * , 043002 ( 2006 ) .
d * 72 * , 023002 ( 2005 ) .
* 99 * , 171101 ( 2007 ) .
e. zas , f. halzen , t. stanev , phys .
d * 45 * , 362 ( 1992 ) .
i. kravchenko _
_ , astropart .
* 19 * , 15 ( 2003 ) ; i. kravchenko _ et al . _
[ rice collaboration ] , phys .
d * 85 * , ( 2012 ) 062004 .
gorham _ et al . _
[ anita collaboration ] , astropart .
phys . * 32 * , ( 2009 ) 10 ; p.w .
gorham _ et al .
* 103 * , 051103 ( 2009 ) ; p.w .
gorham _ et al . _
[ anita collaboration ] , phys .
d * 82 * , ( 2010 ) 022004 ; phys .
d * 85 * , ( 2012 ) 049901(e ) .
t. hankins _ et al .
_ mon .  not .
royal astron .
soc .  * 283 * , 1027 ( 1996 ) .
gorham _ et al .
lett .  * 93 * , 041101 ( 2004 ) .
o. scholten _
astropart .
* 26 * , 219 ( 2006 ) .
a.r . beresnyak _
_ astronomy reports * 49 * 2 , 127 ( 2005 ) .
james _ et al .
d * 81 * , 042003 ( 2010 ) .
jaeger , r.l .
mutel , k.g .
gayley , astropart . phys . * 34 * , 293 ( 2010 ) .
p. allison _
[ ara collaboration ] astropart
. phys . * 35 * , ( 2012 ) 457 s.w .
barwick et al
. j. phys .
60 , 276 - 286 ( 2007 ) .
s. fliescher for the pierre auger collaboration , nucl .
a * 662 * , s124 ( 2012 ) d. ardouin _ et al .
_ , astropart .
phys . * 31 * , 192 ( 2009 ) .
s.  hoover _ et al . _
[ anita collaboration ] , phys .
* 105 * , 151101 ( 2010 ) .
gorham , f.e .
baginski , p. allison , k.m .
liewer , c. miki , b. hill , g.s .
varner , astropart .
35 * 242 , ( 2011 ) f.  halzen , e.  zas , t.  stanev , phys .
b * 257 * , 432 ( 1991 ) .
j. alvarez - muiz , e. zas , phys .
b * 411 * , 218 ( 1997 ) .
j. alvarez - muiz , e. zas , phys .
b * 434 * , 396 ( 1998 ) .
j. alvarez - muiz , r.a .
vzquez and e. zas , phys .
d * 61 * , 023001 ( 1999 ) .
j.  alvarez - muiz , r.a .
vzquez , e.  zas , phys .
d * 62 * , 063001 ( 2000 ) .
s.  razzaque _ et al .
d * 65 * , 103002 ( 2002 ) ; j.  alvarez - muiz , e.  marqus , r.a .
vzquez , e.  zas , phys .
d * 67 * , 101303 ( 2003 ) .
s.  razzaque _ et al .
d * 69 * , 047101 ( 2004 ) .
s.  hussain , d.w .
mckay , phys .
d * 70 * , 103003 ( 2004 ) .
j.  alvarez - muiz , e.  marqus , r.a .
vzquez , e.  zas , phys .
d * 74 * , 023007 ( 2006 ) .
j. alvarez - muiz , c.w .
james , r.j .
protheroe and e. zas , astropart .
32 * , 100 ( 2009 ) .
j. alvarez - muiz , a. romero - wolf , e. zas , phys .
d * 81 * , 123009 ( 2010 ) .
j. alvarez - muiz , w. rodrigues , m. tueros and e. zas , astroparticle physics * 35 * , 287 ( 2012 ) j. alvarez - muiz , a. romero - wolf , e. zas , phys .
d * 84 * , 103003 ( 2011 ) j. alvarez - muiz , w. rodrigues , m. tueros and e. zas , astroparticle physics * 35 * , 325 ( 2012 ) i.  tamm , journal of physics * i * , no . 5 - 6 , 439 ( 1939 ) j.d .
jackson , _ classical electrodynamics 3rd ed .
_ wiley , new york , ( 1998 ) .
james , h. falcke , t. huege , and m. ludwig phys . rev .
e * 84 * , 056602 ( 2011 ) r.v .
buniy , j.p .
ralston , phys .
d * 65 * , 016003 ( 2002 ) .
hu , c .- c .
chen , p. chen , astroparticle physics * 35 * , 421 ( 2012 ) .
j.  alvarez - muiz , g. parente , e. zas , in procs . of the 24@xmath198 icrc , rome ( italy ) , vol . 1 , 1023 ( 1995 ) .
afanasiev , _ vavilov - cherenkov and synchrotron radiation : foundations and applications _ , kluwer academic publishers , ( 2004 )