text:
# On Branched Galois Coverings (𝔹₁)^𝑛→ℙ^𝑛

## Abstract.

For any $`n>1`$, we construct examples branched Galois coverings $`M^n`$ where $`M`$ is one of $`(^1)^n`$, $`^n`$ and $`(𝔹_1)^n`$, and $`𝔹_1`$ is the 1-ball. In terms of orbifolds, this amounts to giving examples of orbifolds over $`^n`$ uniformized by $`M`$.

1. Introduction. In contrast with the considerable literature on the orbifolds over $`^2`$ uniformized by the 2-ball $`𝔹_2`$ (see  and references therein), not much is known about which orbifolds over $`^n`$ are uniformized by other symmetric spaces. In this article, we apply a simple orbifold-covering technique to construct some orbifolds over the projective space $`^n`$ uniformized either by $`(^1)^n`$, $`^n`$ or $`(𝔹_1)^n`$. Our main result is the following theorem.

###### Theorem 1.

Let $`(n,b)`$ be a pair of coprime integers with $`n2`$. There exists a Galois covering $`\left(D_{n,1}^{(b)}\right)^n^n`$ of degree $`n!b^{n^2n}`$ branched along an irreducible degree-$`2b(n1)`$ hypersurface $`D_n^{(b)}^n`$ where $`D_{n,1}^{(b)}D_n^{(b)}`$ is a curve of euler number $`e=b^{n1}(n+1+bnb)`$.

For $`b=1`$, the hypersurface $`D_n^{(1)}`$ is the discriminant hypersurface, and $`D_{n,1}^{(1)}^1`$ is a rational normal curve. In this case one obtains the well-known branched Galois covering $`(^1)^n^n`$. The subvarieties $`D_n^{(b)}`$ and $`D_{n,1}^{(b)}`$ are the liftings respectively of $`D_n^{(1)}`$ and $`D_{n,1}^{(1)}`$ by an abelian branched self-covering $`[Z_0,\mathrm{},Z_n]^n[Z_0^b,\mathrm{},Z_n^b]^n`$. For $`(n,b)\{(3,2),(2,3)\}`$ one has $`e\left(D_{n,1}^{(b)}\right)=0`$, and the universal covering of $`\left(D_{n,1}^{(b)}\right)^n`$ is $`^n`$. The curve $`D_2^{(3)}=D_{2,1}^{(3)}`$ is the nine-cuspidal sextic. For $`b>1`$ and $`(n,b)\{(3,2),(2,3)\}`$ one has $`e\left(D_{n,1}^{(b)}\right)<0`$, and the universal covering of $`\left(D_{n,1}^{(b)}\right)^n`$ is $`(𝔹_1)^n`$.

In case $`(n,b)=(2,3)`$, the claim of Theorem 1 was proved in . The case $`n=2`$ was established in . In this case, $`D_2^{(b)}`$ coincides with $`D_{n,1}^{(b)}`$, which is a curve of genus $`\frac{1}{2}(b^23b+2)`$ with $`3b`$ cusps of type $`x^2=y^b`$ and no other singularities. Irreducibility of $`D_n^{(b)}`$ is proved in Proposition 3. The remaining assertions of Theorem 1 are proved in Theorem 3. Our method leads naturally to the definition of braid groups of the orbifolds over $`^1`$, which we dicuss in Section 4 below. These groups were already introduced in  for some basic cases.

2. Orbifolds. Let $`M`$ be a connected complex manifold, $`G\text{Aut}(M)`$ a properly discontinuous subgroup and put $`N:=M/G`$. Then the projection $`\varphi :MN`$ is a branched Galois covering endowing $`N`$ with a map $`\beta _\varphi :N`$ defined by $`\beta _\varphi (p):=|G_q|`$ where $`q`$ is a point in $`\varphi ^1(p)`$ and $`G_q`$ is the isotropy subgroup of $`G`$ at $`q`$. In this setting, the pair $`(N,\beta _\varphi )`$ is said to be uniformized by $`\varphi :M(N,\beta _\varphi )`$. An orbifold is a pair $`(N,\beta )`$ of an irreducible normal analytic space $`N`$ with a function $`\beta :N`$ such that the pair $`(N,\beta )`$ is locally finitely uniformizable. A covering $`\varphi :(N^{},\beta ^{})(N,\beta )`$ of orbifolds is a branched Galois covering $`N^{}N`$ with $`\beta ^{}=(\beta \text{o}\varphi )/\beta _\varphi \text{o}\varphi `$. Note that the restriction $`(N^{},1)(N,\beta _\varphi )`$ is a uniformization of $`(N,\beta _\varphi )`$. Conversely, let $`(N,\beta )`$ and $`(N,\gamma )`$ be two orbifolds with $`\gamma |\beta `$, and let $`\varphi :(N^{},1)(N,\gamma )`$ be a uniformization of $`(N,\gamma )`$, e.g. $`\beta _\varphi =\gamma `$. Then $`\varphi :(N^{},\beta ^{})(N,\beta )`$ is a covering, where $`\beta ^{}:=\beta \text{o}\varphi /\gamma \text{o}\varphi `$. The orbifold $`(N^{},\beta ^{})`$ is called the lifting of $`(N,\beta )`$ to the uniformization $`N^{}`$ of $`(N,\gamma )`$.

Let $`(N,b)`$ be an orbifold, $`B_\beta :=\mathrm{supp}(\beta 1)`$ and let $`B_1,\mathrm{},B_n`$ be the irreducible components of $`B_\beta `$. Then $`\beta `$ is constant on $`B_i\backslash \mathrm{sing}(B_\beta )`$; so let $`b_i`$ be this number. The orbifold fundamental group $`\pi _1^{orb}(N,\beta )`$ of $`(N,\beta )`$ is the group defined by $`\pi _1^{orb}(N,\beta ):=\pi _1(N\backslash B_\beta )/\mu _1^{b_1},\mathrm{},\mu _n^{b_n}`$ where $`\mu _i^{b_i}`$ is a meridian of $`B_i`$ and $``$ denotes the normal closure. An orbifold $`(N,\beta )`$ is said to be smooth if $`N`$ is smooth. In case $`(N,\beta )`$ is a smooth orbifold the map $`\beta `$ is determined by the numbers $`b_i`$; in fact $`\beta (p)`$ is the order of the local orbifold fundamental group at $`p`$. Since the orbifolds to be considered in this article are exclusively smooth, we shall adopt the convention that such orbifolds are defined to be the pairs $`(N,B)`$ where $`B:=b_1B_1+\mathrm{}+b_nB_n`$ is a divisor with $`b_i1`$. We shall also allow $`b_i`$ to take infinite values, meaning that the corresponding hypersurface $`B_i`$ is removed from the base space $`N`$. If $`𝒪:=(N,B)`$ is an orbifold and $`C`$ a hypersurface in $`N`$, then we shall use the notation $`(𝒪,bC)`$ to denote the orbifold $`(N,B+bC)`$.

3. Discriminants. Let $`n1`$ be an integer and consider the action of the symmetric group $`\mathrm{\Sigma }_n`$ on $`(^1)^n`$. Let $`p_i=[u_i,v_i]^1`$ and let $`\sigma _j`$ ($`j[0,n]`$) be the homogeneous elementary symmetric polynomial

$$\sigma _j(p_1,\mathrm{},p_n):=\underset{A[1,n],|A|=j}{}\left(\underset{\alpha A}{}x_\alpha \underset{\beta [1,n]\backslash A}{}y_\beta \right)$$

It is well known that the map $`\varphi :(^1)^n^n`$ given by

$$\varphi :(p_1,\mathrm{},p_n):=[\sigma _0(p_1,\mathrm{},p_n):\mathrm{}:\sigma _n(p_1,\mathrm{},p_n)]$$

is $`\mathrm{\Sigma }_n`$\- invariant and gives an isomorphism $`(^1)^n/\mathrm{\Sigma }_n^n`$.

Let $`\pi _i:(^1)^n^1`$ be the $`ith`$ projection map, $`q`$ a point in $`^1`$, and put $`F_q^i:=\pi _i^1(q)`$. Let $`\tau _{ij}\mathrm{\Sigma }_n`$ be the transposition exchanging the $`i`$th and $`j`$th coordinates of $`(p_1,\mathrm{},p_n)(^1)^n`$. Since $`\tau _{1i}F_q^1=F_q^i`$, the hypersurface $`H_q:=\varphi _n(F_q^i)`$ does not depend on $`i`$.

###### Lemma 1.

For any $`q^1`$, the hypersurface $`H_q`$ is a hyperplane in $`^n`$. For any set $`\{q_0,\mathrm{},q_m\}^1`$ of distinct points, the hyperplanes $`H_{q_0},\mathrm{},H_{q_m}`$ are in general position.

Proof. Suppose without loss of generality that $`i=1`$. Then $`H_q`$ is parametrized as $`H_q=[X_0:X_1:\mathrm{}:X_n]^n`$, where $`X_j=\sigma _j(q,p_2,\mathrm{},p_n)`$ and $`p_i^1`$ ($`i[2,n]`$). If $`q=[u_1:v_1]=[x:y]`$ and $`p_i=[u_i:v_i]`$ ($`i[2,n]`$) then one has the identity

(1)
$$P(A,B):=\underset{j[0,n]}{}(1)^{nj}\sigma _j(q,p_2,\mathrm{},p_n)A^jB^{nj}=\underset{i[1,n]}{}(u_iAv_iB)$$

Substitute $`[A:B]=[y:x]`$ in (1). Since the right-hand side of (1) vanish at the point $`(q,p_2,\mathrm{},p_n)`$, so does the middle term, and thus $`H_q`$ satisfies the linear equation

(2)
$$\underset{j[0,n]}{}(1)^{nj}y^jx^{nj}X_j=0$$

Let $`\{q_i=[x_i:y_i]:i[0,n]\}`$ be a set of $`n+1`$ points. Since the determinant of the projective Vandermonde matrix $`𝒱an(q_0,\mathrm{},q_n)`$ given by

$$𝒱an_{i,j}(q_0,\mathrm{},q_n):=(1)^{nj}y_i^jx_i^{nj}i,j[0,n]$$

vanish if and only if $`q_i=q_j`$ for some $`i,j[0,n]`$, the hyperplanes $`H_{q_0},\mathrm{},H_{q_n}`$ are always in general position. $`\mathrm{}`$

The hypersurface

$$\mathrm{\Delta }_n:=\{(p_1,\mathrm{},p_n)(^1)^n:p_i=p_j\text{for some}\mathrm{\hspace{0.17em}\hspace{0.17em}1}ijn\}$$

of $`(^1)^n`$ consists of points fixed by an element of $`\mathrm{\Sigma }_n`$, so that the covering $`\varphi `$ is branched along the hypersurface $`D_n:=\varphi (\mathrm{\Delta }_n)`$, which is called the discriminant hypersurface since it is defined by the discriminant of the homogeneous polynomial $`P(A,B)`$.

In the orbifold terminology, one has an orbifold covering

(3)
$$((^1)^n,a\mathrm{\Delta }_n)(^n,2aD_n)$$

Let $`\{q_0,\mathrm{},q_m\}^1`$ be $`m+1`$ distinct points, $`b_0,\mathrm{},b_m`$ numbers in $`\{\mathrm{}\}`$ and consider the orbifold

$$(b_0,\mathrm{},b_m):=(^1,b_0q_0+\mathrm{}+b_mq_m)$$

Let $`n1`$ be an integer and consider the orbifold $`(b_0,\mathrm{},b_m)^n`$. Let $`𝒢_n`$ be the orbifold

$$𝒢_n(a,b_0,\mathrm{},b_m):=((b_0,\mathrm{},b_m)^n,a\mathrm{\Delta }_n)$$

and define the orbifold $`_n(a,b_0,\mathrm{},b_m)`$ as

$$_n(a,b_0,\mathrm{},b_m):=(^n,aD_n+b_0H_{q_0}+\mathrm{}+b_mH_{q_m})$$

By the covering in (3) and Lemma 1 one has the following fact

###### Lemma 2.

There is an orbifold covering of degree $`n!`$

$$\varphi :𝒢_n(a,b_0,\mathrm{},b_m)_n(2a,b_0,\mathrm{},b_m)$$

In particular, for $`a=1`$ one has the orbifold covering

$$\varphi :(b_0,\mathrm{},b_m)^n𝒢_n(1,b_0,\mathrm{},b_m)_n(2,b_0,\mathrm{},b_m)$$

The following facts are well known (see ):

###### Theorem 2.

\[Bundgaard-Nielsen,Fox\] The orbifold $`(b_0,\mathrm{},b_m)`$ admits a finite uniformization if $`n>1`$, $`b_i<\mathrm{}`$ ($`1im`$) and if $`n=2`$, then $`b:=b_0=b_1`$. Let $`R(b_0,\mathrm{},b_m)`$ be a finite uniformization.
(i) $`R^1`$ if $`n=2,b_0=b_1<\mathrm{}`$ or $`n=3`$, $`b_0^1+b_1^1+b_2^1>1`$. In this case , $`^1`$ is also the universal uniformization. The groups $`\pi _1^{orb}\left((b,b)\right)`$ and $`\pi _1^{orb}\left((b_0,b_1,b_2)\right)`$ are finite of orders $`b`$ and $`2\left[b_0^1+b_1^1+b_2^11\right]^1`$ respectively.
(ii) $`R`$ is of genus 1 if $`n=3`$, $`b_0^1+b_1^1+b_2^1=1`$ or $`n=4`$, $`b_0=b_1=b_2=b_3=2`$. Hence, $``$ is the universal uniformization of these orbifolds. Moreover, $`(\mathrm{},\mathrm{})`$ and $`(2,2,\mathrm{})`$ are uniformized by $``$. The corresponding orbifold fundamental groups are infinite solvable.
(iii) $`R`$ is of genus $`>1`$ otherwise, and the universal uniformization is $`(𝔹_1)^n`$, where $`𝔹_1`$ is the unit disc in $``$. The corresponding orbifold fundamental groups are big (i.e. they contain non-abelian free subgroups).

In virtue of the covering $`\varphi :(b_0,\mathrm{},b_m)^n_n(2,b_0,\mathrm{},b_m)`$ one has the following corollary.

###### Corollary 1.

Let $`n>1`$, $`b_i<\mathrm{}`$ ($`1im`$) and if $`n=2`$, then $`b_0=b_1`$. Then the orbifold $`_n(2,b_0,\mathrm{},b_m)`$ admits a finite uniformization by $`R^n`$, where $`R`$ is the uniformization of $`(b_0,\mathrm{},b_m)`$ given in Theorem 2. The orbifolds $`(2,\mathrm{},\mathrm{})`$ and $`(2,2,2,\mathrm{})`$ are uniformized by $`^n`$. Moreover, $`\pi _1^{orb}\left((b,b)\right)`$ is a finite group of order $`n!b^n`$ and $`\pi _1^{orb}\left((b_0,b_1,b_2)\right)`$ is a finite group of order $`n!2^n\left[b_0+b_1+b_21\right]^n`$ if $`b_0^1+b_1^1+b_2^1>1`$.

4. Braid groups. Following and generalizing , let us call the groups

$$𝐏_n(a,b_0,\mathrm{},b_m):=\pi _1^{orb}\left(𝒢_n(a,b_0,\mathrm{},b_m)\right)$$

the pure braid groups of $`(b_0,\mathrm{},b_m)`$ on $`n`$ strands, and the groups

$$𝐁_n(a,b_0,\mathrm{},b_m):=\pi _1^{orb}\left(_n(a,b_0,\mathrm{},b_m)\right)$$

the braid groups of $`(b_0,\mathrm{},b_m)`$ on $`n`$ strands. Obviously, the group $`𝐁_n(a,b_0,\mathrm{},b_m)`$ is a quotient of $`𝐁_n(a^{},b_0^{},\mathrm{},b_m^{})`$ provided $`a|a^{}`$ and $`b_i|b_i^{}`$ for $`0in`$. The group $`𝐁_n(a,b_0,\mathrm{},b_m)`$ is a subgroup of $`𝐁_{n+k}(a,b_k,\mathrm{},b_m)`$ in case the equality $`a=b_0=\mathrm{}=b_{k1}`$ holds. The group $`𝐁_n(2a,b_0,\mathrm{},b_m)`$ is a normal subgroup of index $`n!`$ in the group $`𝐏_n(a,b_0,\mathrm{},b_m)`$. The group $`𝐁_n(a,b_0,\mathrm{},b_m)`$ admits the presentation (see  for the case $`n=2`$ and 111The projective relation was kindly communicated by Paolo Bellingeri., ,  for the general case )

generators
$`\sigma _1,\mathrm{},\sigma _{n1}`$, $`\tau _0,\mathrm{},\tau _m`$
braid relations
$`[\sigma _i,\sigma _j]=1`$, $`|ij|>1`$
$`\sigma _i\sigma _{i+1}\sigma _i=\sigma _{i+1}\sigma _i\sigma _{i+1}`$, $`1in1`$
mixed relations
$`(\sigma _1\tau _i)^2=(\tau _i\sigma _1)^2`$, $`1im`$,
$`[\tau _i,\sigma _j]=1`$, $`j1`$, $`1im`$
$`[\sigma _1\tau _i\sigma _1^1,\tau _j]=1`$, $`1i<jm`$,
projective relation $`^1`$
$`\sigma _1\sigma _2\mathrm{}\sigma _{n1}\tau _0\mathrm{}\tau _m\sigma _{n1}\mathrm{}\sigma _2\sigma _1=1`$
orbifold relations
$`\tau _0^{b_0}=\mathrm{}=\tau _m^{b_m}=\sigma _1^a=1`$

In particular, the group $`𝐁_n(\mathrm{},\mathrm{})`$ is the usual braid group of $``$ introduced by Artin . The group $`𝐁_n(\mathrm{})𝐁_n(\mathrm{},1)`$ is the braid group of the sphere, see . On the other hand, one has

$$𝐁_1(b_0,\mathrm{},b_m)\tau _0,\mathrm{},\tau _m|\tau _0^{b_0}=\mathrm{}=\tau _m^{b_m}=\tau _0\mathrm{}\tau _m=1$$

In case $`n=2`$, the discriminant hypersurface $`D_2^{(1)}`$ is a smooth quadric, and the lines $`H_{q_i}`$ are tangent to $`D_2^{(1)}`$ (see ). In particular, the groups $`𝐁_2(a,b)`$ are abelian. The group $`𝐁_2(a,b,c)`$ admits the presentation

(4)
$$𝐁_2(a,b,c)\tau ,\sigma |(\tau \sigma )^2=(\sigma \tau )^2,\tau ^b=(\tau \sigma ^2)^c=\sigma ^a=1$$

###### Proposition 1.

For $`b,c<\mathrm{}`$, the group $`𝐁_2(a,b,c)`$ is a finite central extension of the triangle group $`T_{2,a,d}:\tau ,\sigma |(\tau \sigma )^2=\tau ^d=\sigma ^a=1`$, where $`d:=\mathrm{gcd}(b,c)`$. Hence, $`𝐁_2(a,b,c)`$ is finite $`1/a+1/b>1/2`$, infinite almost solvable if $`1/d+1/a=1/2`$, and big otherwise (i.e. it contains non-abelian free subgroups). The group $`𝐁_2(a,b,b)`$ is of order $`2b\left[a^1+b^12^1\right]^1`$ if $`1/a+1/b>1/2`$.

Proof. Note that $`\delta :=(\tau \sigma )^2`$ is central in $`𝐁_2(a,b,c)`$, so that $`(\tau \sigma ^2)^c=1`$ $``$ $`(\sigma \tau \sigma )^c=1`$ $``$ $`(\tau ^1\delta )^c=\tau ^c\delta ^c=1`$. The element $`\delta `$ is of finite order. Adding the relation $`\delta =1`$ to the presentation (4) yields the triangle group $`T_{2,a,d}`$, which is finite if $`1/a+1/d>1/2`$, infinite solvable if $`1/a+1/d=1/2`$, and big otherwise. In case $`c=b`$, one has $`d=b`$ and the triangle group is of order $`2\left[a^1+b^12^1\right]^1`$ if $`1/a+1/b>1/2`$, which shows that $`𝐁(a,b,b)`$ is of order $`2b\left[a^1+b^12^1\right]^1`$. ∎

Let $`R^n`$ be a uniformization of the orbifold $`\left((b_0,\mathrm{},b_m)\right)^n`$. If $`km`$ then any orbifold $`𝒢_n(2a,c_0b_0,\mathrm{}c_mb_m,c_{m+1},\mathrm{},c_k)`$ can be lifted to $`R^n`$. In case $`R^1`$ or $`R`$ one obtains some arrangements associated to reflection groups as follows. Suppose that $`q_0=[0:1]`$ and $`q_1=[1:0]`$. Lifting $`𝒢_n(2a,cb,\mathrm{})`$ to the uniformization of $`𝒢_n(2,b,b)`$ yields the orbifold $`(^n,a\mathrm{\Delta }_n^{(b)}+cF)`$ where $`F:=\{(X_1,\mathrm{},X_n)^n:X_1\mathrm{}X_n=0\}`$ and $`\mathrm{\Delta }_n^{(b)}`$ is the lifting of the superdiagonal

$$\mathrm{\Delta }_n^{(b)}:=\{(X_1,\mathrm{},X_n)^n:\psi _b(p_i)=\psi _b(p_j)\text{for some}\mathrm{\hspace{0.17em}\hspace{0.17em}1}ijn\}$$

with $`\psi _b(X)=X^b`$ if $`b<\mathrm{}`$ and $`\psi _{\mathrm{}}(X)=\mathrm{exp}(2\pi \mathrm{i}X)`$. Setting $`b=2`$ in this construction identifies the group $`𝐁_n(\mathrm{},\mathrm{},\mathrm{})`$ with the Artin group corresponding to the diagrams $`B_n`$ (see ).

The groups $`𝐁_2(a,b,c,d)`$ admits the simplified presentation (see )

$$𝐁_2(a,b,c,d)\tau ,\rho ,\sigma |\begin{array}{c}(\tau \sigma )^2=(\sigma \tau )^2,(\rho \sigma )^2=(\sigma \rho )^2,[\rho ,\tau ]=1,\hfill \\ \tau ^b=(\sigma \tau \sigma \rho )^d=\rho ^c=\sigma ^a=1\hfill \end{array}$$

We summarized the known information about the orbifolds $``$ and the corresponding braid groups in Table 1 below. Suppose that if $`(n,m)=(2,1)`$ then $`b_0=b_1`$. We believe that the group $`𝐁_n(a,b_0,\mathrm{},b_m)`$ is finite if

$$\frac{2(n1)}{a}+\underset{i[0,m]}{}\frac{1}{b_i}>n+m2,$$

infinite solvable if the equality holds, and big otherwise.

| Orbifold | Uniformization | Braid group | Reference |
| --- | --- | --- | --- |
| $`_n\left(2\right)`$ | $`\left(^1\right)^n`$ | $`n!`$ | Cor. 1 |
| $`_n(2,b,b)`$ | $`\left(^1\right)^n`$ | $`n!b^n`$ | Cor. 1 |
| $`\begin{array}{c}_n(2,b,c,d)\hfill \\ \left(1/b+1/c+1/d>1\right)\hfill \end{array}`$ | $`\left(^1\right)^n`$ | $`n!2^n\left[\frac{1}{b}+\frac{1}{c}+\frac{1}{d}1\right]^n`$ | Cor. 1 |
| $`\begin{array}{c}_n(2,b,c,d)\hfill \\ \left(1/b+1/c+1/d=1\right)\hfill \end{array}`$ | $`^n`$ | Crystallographic | Cor. 1 |
| $`_n(2,2,2,2,2)`$ | $`^n`$ | Crystallographic | Cor. 1 |
| $`\begin{array}{c}_n(2,b_0,\mathrm{},b_m)\hfill \\ \left(\mathrm{otherwise}\right)\hfill \end{array}`$ | $`\left(𝔹_1\right)^n`$ | Linear | Cor. 1 |
| $`_n(\mathrm{},\mathrm{},\mathrm{})`$ | - | $`B_n`$-Artinian |  |
| $`_2(\mathrm{},\mathrm{},\mathrm{},\mathrm{})`$ | - | $`\stackrel{~}{C}_2`$-Artinian |  |
| $`\begin{array}{c}_2(a,b,b)\hfill \\ \left(1/a+1/b>1/2\right)\hfill \end{array}`$ | - | $`2b\left[\frac{1}{a}+\frac{1}{b}\frac{1}{2}\right]^1`$ | Prop. 1 |
| $`\begin{array}{c}_2(a,b,b)\hfill \\ \left(1/a+1/b=1/2\right)\hfill \end{array}`$ | - | $`\mathrm{}`$ almost solvable | Prop. 1,  |
| $`_3(a,\mathrm{})\left(a=3,4,5\right)`$ | - | 24, 96, 600 |  |
| $`_n(3,\mathrm{})\left(n=4,5\right)`$ | - | 648, 155520 |  |
| $`_3(\mathrm{},2)`$ | - | 192 | Maple |
| $`_4\left(a\right)\left(a=4,5\right)`$ | - | 192, 60 | Maple |
| $`_5\left(4\right)`$ | - | 120 | Maple |
| $`_2(a,2,2,2)`$ | K3 $`\left(a=4\right)`$ | $`4a^3`$ |  |
| $`_2(3,3,2,2)`$ | K3 | 576 |  |
| $`\begin{array}{cc}_2(3,3,4,4)\hfill & _2(4,4,4,4)\hfill \\ _2(3,6,6,2)\hfill & _2(3,3,3,6)\hfill \end{array}`$ | $`𝔹_2`$ | Picard Modular | , |
| $`_2(3,3,4,2)`$ | $`𝔹_1\times 𝔹_1\mathrm{?}`$- | Unknown |  |

Table 1.

5. Another covering of $`(a,b_0,\mathrm{},b_m)`$. Let $`b`$ be an integer and consider the orbifold $`𝒦_n(b):=(^n,bH_{q_0}+\mathrm{}+bH_{q_n})`$. By Lemma 1, the hyperplanes $`H_{q_0},\mathrm{},H_{q_n}`$ are in general position. It is well known that the universal uniformization of this orbifold is $`^n`$. Applying a projective transformation one may assume that the hyperplanes $`H_{q_i}`$ are given by the equations $`Y_i=0`$ where $`[Y_0:\mathrm{}:Y_n]^n`$. In this case the uniformization $`\psi _b:^n𝒦_n(b)`$ is nothing but the map

$$[Y_0:\mathrm{}:Y_n]=\psi _b([Z_1:\mathrm{}:Z_n])=[Z_1^b:\mathrm{}:Z_n^b]$$

It is clear that the orbifold $`_n(a,bb_0,\mathrm{},bb_n)`$ lifts to the uniformization of $`𝒦_n(b)`$. Put $`D_n^{(b)}:=\psi _b^1(D_n)`$, denote $`M_{q_i}:=\psi ^1(H_{q_i})`$ and define the orbifold

$$_n^{(b)}(a,b_0,\mathrm{},b_m):=(^n,aD_n^{(b)}+b_0M_{q_0}+\mathrm{}b_nM_{q_n})$$

to be this lifting. In case $`n=2`$ these liftings were studied in . For $`n>2`$ the following proposition is valid:

###### Proposition 2.

For $`n>2`$ and $`b2`$ the orbifolds $`_n^{(b)}(2,b_0,\mathrm{},b_m)`$ are uniformized by $`(𝔹_1)^n`$ except the orbifold $`_3^{(2)}(2)`$, which is uniformized by $`^3`$.

Proof. There is an orbifold covering

$$_n^{(b)}(2,b_0,\mathrm{},b_m)_n(a,bb_0,\mathrm{},bb_n)$$

The claim follows, since by Corollary 1 the latter orbifold is uniformized by $`^3`$ if $`b=2`$, $`n=3`$, $`b_0=\mathrm{}=b_n=1`$ and by $`(𝔹_1)^n`$ otherwise. $`\mathrm{}`$

For $`k[1,n]`$, define the $`k`$-dimensional subvarietiy $`\mathrm{\Delta }_{n,k}`$ of $`\mathrm{\Delta }_n`$ by

$$\mathrm{\Delta }_{n,k}:=\{(p_1,p_2,\mathrm{},p_n)(^1)^n:p_k=p_{k+1}=\mathrm{}=p_n\}(^1)^k$$

Thus, $`\mathrm{\Delta }_{n,n1}`$ is an irreducible component of $`\mathrm{\Delta }_n`$ and $`\mathrm{\Delta }_{n,1}`$ is the diagonal in $`(^1)^n`$. The subgroup of $`\mathrm{\Sigma }_n`$ acting on $`\mathrm{\Delta }_{n,k}`$ is a symmetric group $`\mathrm{\Sigma }_{k1}`$, so that $`D_{n,k}:=^1\times ^{k1}`$. These varieties admits the parametrizations

(5)
$$D_{n,k}:[X_0:\mathrm{}:X_n]^nX_j=\sigma _j(p_1,\mathrm{},p_n),p_k=\mathrm{}=p_n$$

In particular, the curve $`D_{n,1}`$ is a rational normal curve parametrized as

$$[(\begin{array}{c}n\\ 0\end{array})v^n,(\begin{array}{c}n\\ 1\end{array})uv^{n1},\mathrm{},(\begin{array}{c}n\\ n\end{array})u^n]([u:v]^1)$$

Applying the projective transformation $`𝒱an(q_0,\mathrm{},q_n)`$ to the parametrizations (5) gives the parametrization $`D_{n,k}:[Y_0:\mathrm{}:Y_n]^n`$, where

(6)
$$\underset{j[0,n]}{}(1)^{nj}y_i^jx_i^{nj}\sigma _j(p_1,\mathrm{},p_n),p_k=\mathrm{}=p_n$$

Let $`p_i=[u_i:v_i]`$ and let $`[u:v]=[u_k:v_k]=\mathrm{}=[u_n:v_n]`$. In virtue of the identity (1) one has the parametrizations $`D_{n,k}:[Y_0:\mathrm{}:Y_n]^n`$ where

(7)
$$Y_j=(uy_jvx_j)^{nk+1}\underset{i[1,k1]}{}(u_iy_jv_ix_j)$$

In particular, the curve $`D_{n,1}`$ is parametrized as

(8)
$$D_{n,1}:[(uy_0vx_0)^n:\mathrm{}:(uy_nvx_n)^n]$$

The varieties $`D_{n,k}^{(b)}`$ are parametrized as

(9)
$$D_{n,k}^{(b)}:[Z_0:\mathrm{}:Z_n]Z_j^b=(uy_jvx_j)^{nk+1}\underset{i[1,k1]}{}(u_iy_jv_ix_j)$$

Note that the parametrizations (7) and (9) are not generically one-to-one unless $`k2`$, since (7) is a map $`(^1)^kD_{n,k}`$.

###### Proposition 3.

(i) The curve $`D_{n,1}^{(b)}`$ is irreducible if and only if $`\mathrm{gcd}(n,b)=1`$. Hence, the subvarieties $`D_{n,k}^{(b)}`$ are irreducible if $`\mathrm{gcd}(n,b)=1`$.

Definition. Let $`t`$ and $`\psi _t`$ be the map

$$\psi _t:[Z_0:\mathrm{}:Z_n]^n[Z_0^t:\mathrm{}Z_n^t]^n$$

Let $`V^n`$ be a subvariety and $`r,s`$ such that $`s>1`$. Then $`V^{(r/s)}`$ is the subvariety of $`^n`$ defined as

$$V^{(r/s)}:=(\psi _r^1\text{o}\psi _s)(V)$$

In particular, $`V^{(r/r)}`$ is the orbit of $`V`$ under the $`(/(r))^n`$-action on $`^n`$.

Proof of the Proposition. The parametrization (8) shows that $`D_{n,1}L^{1/n}`$, where $`L`$ is a line $`^n`$ in general position with respect to $`\psi _n`$, in other words $`L`$ intersects the hyperplane arrangement $`Z_0\mathrm{}Z_n=0`$ transversally at smooth points. Hence there is a surjection of fundamental groups

(10)
$$\pi _1(L\backslash \{\stackrel{~}{q}_0,\mathrm{},\stackrel{~}{q}_n\})\pi _1(^n\backslash \{Z_0,\mathrm{},Z_n\})$$

where $`\stackrel{~}{q}_i:=Z_iL`$. Let $`(b)`$, $`𝒦(b)`$ be the orbifolds

$$(b):=(L,b\stackrel{~}{q}_0+\mathrm{}+b\stackrel{~}{q}_n),𝒦(b):=(^n,bZ_0+\mathrm{}+bZ_n)$$

Then (10) induce a surjection of orbifold fundamental groups

$$\pi _1^{orb}\left((b)\right)\pi _1^{orb}\left(𝒦(b)\right)$$

(one may say: $`(b)`$ is a sub-orbifold of $`𝒦(b)`$). This shows that the curve $`L^{(b)}`$ is irreducible and is a uniformization of $`(b)`$. Since $`\mathrm{gcd}(n,b)=1`$, one has $`D_{n,1}^{(b)}=L^{(b/n)}`$, showing that $`D_{n,1}^{(b)}`$ is irreducible. Note that $`D_{n,1}^{(b)}`$ is the maximal abelian orbifold covering of $`(b)`$. Irreducibility of $`D_{n,k}^{(b)}`$ follows since $`D_{n,1}^{(b)}`$ is a subvariety of $`D_{n,k}^{(b)}`$. ∎

Let $`𝒪(b)`$ be the orbifold

$$𝒪(b):=(D_{n,1},b\overline{q}_0+\mathrm{}+b\overline{q}_n),$$

where $`\overline{q}_i:=Y_iD_{n,1}`$. The orbifold $`𝒪(b)`$ is identified via the covering $`\varphi `$ with the orbifold

$$𝒫(b):=(\mathrm{\Delta }_{n,1},bq_0^{}+\mathrm{}+bq_n^{}),$$

where this time $`q_i^{}:=\varphi ^1(\overline{q_i})`$. In turn, $`𝒪(b)`$ is identified with the orbifold $`(b,\mathrm{},b)`$ via the coordinate projection. By the proof of Proposition 3, these orbifolds are identified with the orbifold $`(b)`$ in case $`(n,b)=1`$.

###### Theorem 3.

Let $`\mathrm{gcd}(n,b)=1`$. Then there is a finite uniformization $`\xi _n:(D_{n,1}^{(b)})^n_n^{(b)}(2)`$ which is of degree $`n!b^{n^2n}`$.

Proof. One has the diagram

$$\begin{array}{ccc}& \xi _n& \\ _n^{(b)}(2)& & (D_{n,1}^{(b)})^n\\ & & \\ \psi _b& & \zeta _b\\ & \varphi _n& \\ _n(2,b,\mathrm{},b)& & 𝒪(b)^n\end{array}$$

where $`\zeta _b:(D_{n,1}^{(b)})^n𝒪(b)^n`$ is the maximal abelian orbifold covering and $`\xi _n`$ is to be shown to be a branched Galois covering of degree $`n!b^{n^2n}`$. It suffices to show that the group $`H:=(\varphi _n\text{o}\zeta _b)_{}\pi _1\left((D_{n,1}^{(b)})^n\right)`$ is a normal subgroup of $`K:=(\psi _b)_{}\pi _1^{orb}\left(_n^{(b)}(2)\right)`$. Let $`\sigma `$ be a meridian of $`D_n`$. Then since $`\pi _1^{orb}\left(_n(2,b,\mathrm{},b)\right)/\sigma \pi _1^{orb}\left(𝒦_n(b)\right)\left(/(b)\right)^n`$ is the Galois group of $`\psi _b`$, the group $`K`$ is the normal subgroup of $`\pi _1^{orb}\left(_n(2,b,\mathrm{},b)\right)`$ generated by $`\sigma `$, i.e. $`K\sigma `$. The group $`\pi _1^{orb}\left(_n(2,b,\mathrm{},b)\right))/K`$ being abelian, one has $`[\tau _i,\tau _j]K`$ for $`i,j[0,n]`$. On the other hand one has

$$\pi _1^{orb}\left(_n(2,b,\mathrm{},b)\right)/\tau _0,\mathrm{},\tau _n\pi _1^{orb}\left(_n(2)\right)\mathrm{\Sigma }_n$$

Since $`\mathrm{\Sigma }_n`$ is the Galois group of $`\varphi _n`$, one has $`\varphi _{}𝒪(b)^n\tau _0,\mathrm{},\tau _n.`$ Since $`\zeta _n`$ is the maximal abelian orbifold covering, one has $`H[\tau _i,\tau _j]`$. This shows that $`H`$ is a normal subgroup of $`K`$. Since $`\mathrm{deg}(\zeta _b)=b^{n^2}`$, $`\mathrm{deg}(\mathrm{\Phi }_n)=n!`$ and $`\mathrm{deg}(\psi _b)=b^n`$, one has

$$\mathrm{deg}(\xi _n)=\frac{\mathrm{deg}(\zeta _b)\mathrm{deg}(\varphi _n)}{\mathrm{deg}(\psi _b)}=n!b^{n^2n}$$

The euler number of $`D_{n,1}^{(b)}`$ is easily computed by Riemann-Hurwitz formula. $`\mathrm{}`$

6. Remarks. Consider the restriction of $`D_{n,k}`$ to the $`nk+1`$ dimensional linear subspace $`M_{nk+1}:=\{[Y_0:\mathrm{}:Y_n]^n|Y_{nk+2}=\mathrm{}=Y_n=0\}`$ of $`^n`$. Setting $`[u:v]=[x_n:y_n]`$ and $`[u_i:v_i]=[x_{ni}:y_{ni}]`$ for $`i[1,k2]`$ in (7) we see that $`D_{n,k}`$ has a 1-dimensional linear component $`L`$ in $`M_{nk+1}^{nk+1}`$, parametrized as $`[Y_0:\mathrm{}:Y_{nk+1}]M_{nk+1}`$ where

$$Y_l=(u_ny_lv_nx_l)(x_ny_ly_nx_l)^{nk+1}\underset{i[2,k1]}{}(x_{ni}y_ly_{ni}x_l)$$

for $`l[0,nk+1]`$ and $`[u_n:v_n]^1`$. It is readily seen that there are $`k1`$ such lines. In case $`[u_i:v_i]=0`$ for $`i[1,k1]`$, one has the curve $`C`$ in $`D_{n,k}M_{nk+1}`$ parametrized as $`[Y_0:\mathrm{}:Y_{nk+1}]M_{nk+1}`$ where

$$Y_l=(uy_lvx_l)^{nk+1}\underset{i[1,k1]}{}(x_{ni}y_ly_{ni}x_l)$$

for $`l[0,nk+1]`$ and $`[u:v]^1`$, which shows that $`C`$ is the curve $`E^{(1/nk+1)}`$ for some line $`E`$ in $`^{nk+1}`$. The lines $`L`$ are tangent to $`C`$ with multiplicity $`nk+1`$. In case $`k=n1`$, one has $`M_{nk+1}^2`$, and one obtains an arrangement of a quadric $`C`$ with $`n2`$ tangent lines. The lines $`Y_0=0`$, $`Y_1=0`$ and $`Y_2=0`$ are also tangent to this quadric.

From these considerations it is easy to obtain a description of the intersection of $`D_{n,k}^{(b)}`$ with $`^{nk+1}Z_{nk+2}=\mathrm{}=Z_n=0`$. For $`D_{3,2}^{(2)}`$, this is the arrangement of a quadric with four tangent lines.

Let $`H^n`$ be a hyperplane. The intersection $`H^{(1/2)}M_2`$ is a quadric, tangent to the lines $`Y_0=0`$, $`Y_1=0`$ and $`Y_2=0`$, which is very similar to the intersections $`D_nM_2`$. In contrast with this, there is the following fact: In a recent article , it was proved that the dual of $`D_n`$ is one dimensional (we believe that $`D_{n,k}`$ and $`D_{n,nk}`$ are duals), whereas it is easy to show that $`H^{(r/s)}`$ and $`H^{(r/rs)}`$ are duals, so that the dual of $`H^{(1/2)}`$ is the cubic hypersurface $`H^{(1)}`$. Note also that $`D_n`$ is of degree $`2(n1)`$, whereas $`H^{(1/2)}`$ is of degree $`2^{n1}`$. It is of interest to know more about the varieties $`D_{n,k}^{(r/s)}`$ and their duals.

Appendix: The curves $`L^{(r/s)}`$. In $`^2`$, many interesting curves appears as $`L^{r/s}`$. For example, $`L^{1/2}`$ is the curve $`D_{2,1}`$, a quadric tangent to the coordinate lines, $`L^{3/2}D_{2,1}^3`$ is a nine cuspidal sextic, $`L^{2/3}`$ is a Zariski sextic with 4 nodes and 6 cusps, $`L^{1/2}D_{2,1}^1`$ is a three cuspidal quartic, $`L^{(1)}`$ is a quadric passing through the intersection points of the coordinate lines.

###### Proposition 4.

If $`r,s0`$ are coprime integers , then $`L^{r/s}`$ is an irreducible curve of degree $`sr`$ and genus $`(r1)(r2)/2`$, with $`3r`$ points of type $`x^r=y^s`$ and $`r^2(s1)(s2)/2`$ nodes.

Proof. We begin by proving that the curves $`L^{1/s}`$ are nodal. For this, it suffices to show that the orbit of $`L`$ under the action of the group $`/(s)/(s)`$ has only double points on $`^2\backslash \{xyz=0\}`$. If $`\omega :=e^{2\pi i/s}`$, then the orbit of $`L`$ consists of the lines $`L_{ij}:=a\omega ^ix+b\omega ^jy+cz=0`$ for $`1i,js`$. Suppose that no pairs of lines among the lines $`L_{i,j}`$, $`L_{k,l}`$, $`L_{p,q}`$ meet on $`xyz=0`$. Then they meet at a point $`\{xyz=0\}`$ only if the determinant of the matrix

| $`a\omega ^i`$ | $`b\omega ^j`$ | $`c`$ |
| --- | --- | --- |
| $`a\omega ^k`$ | $`b\omega ^l`$ | $`c`$ |
| $`a\omega ^p`$ | $`b\omega ^q`$ | $`c`$ |

vanish. Since $`abc0`$, this is equivalent to the vanishing of

det    $`\omega ^\alpha 1`$ $`\omega ^\beta 1`$  $`\omega ^\gamma 1`$ $`\omega ^\theta 1`$

where $`\alpha :=ki`$, $`\beta :=lj`$, $`\gamma :=pi`$ and $`\theta :=qj`$. The integers $`\alpha `$, $`\beta `$, $`\gamma `$, $`\theta `$ are not multiples of $`s`$ by hypothesis. Then vanishing of the determinant implies

$$\frac{(\omega ^\alpha 1)(\omega ^\theta 1)}{(\omega ^\beta 1)(\omega ^\gamma 1)}=1\frac{(\omega ^{\alpha /2}\omega ^{\alpha /2})(\omega ^{\theta /2}\omega ^{\theta /2})}{(\omega ^{\beta /2}\omega ^{\beta /2})(\omega ^{\gamma /2}\omega ^{\gamma /2})}=\omega ^{(\beta +\gamma \alpha \theta )/2}$$

Since the left-hand side of the latter expression is real, so must be the right-hand side. Therefore

$$\text{Im}(e^{\pi i(\beta +\gamma \alpha \theta )/s})=0s|\beta +\gamma \alpha \theta .$$

But this means that there is a pair of lines meeting at $`z=0`$, contradiction. This shows that the curves $`L^{1/s}`$ are nodal.

Since $`L^{1/s}`$ is a rational curve of degree $`s`$, it must have $`(s1)(s2)/2`$ nodes. Since $`L^{r/s}=\varphi _r^1(L^{1/s})`$, the number of nodes of $`L^{r/s}`$ is $`r^2(s1)(s2)/2`$. Obviously, three flex points of $`L^{1/s}`$ are lifted as $`3r`$ cusps of type $`x^r=y^s`$. The genus of $`L^{r/s}`$ can be calculated by the genus formula, or by noting that the curves $`L^{r/s}`$ are coverings of $`L^{1/s}`$ branched at these three flex points, with the branching index $`r`$. $`\mathrm{}`$

Galatasaray University, Department of Mathematics, Ortaköy/İstanbul, Turkey
e-mail adress: muludag@gsu.edu.tr


format: text

text:
# On the generation of UHECRs in GRBs: a reappraisal

## 1 Introduction

We suggested some time ago (Vietri 1995) that Gamma Ray Bursts provided the ideal accelerating sites for UHECRs, supporting this view with a computation of the highest energy attainable by protons accelerated around the fireball’s external shock (which was found to exceed $`10^{20}eV`$), and of the average energy deposited by bursts per unit time and volume in gamma ray photons (which came close to that inferred for UHECRs, if these have a cosmological origin); assuming comparable efficiencies in the production of gamma ray photons and non–thermal particles, we argued that this provided an excellent explanation, without fitting any parameter, of the UHECRs’ flux levels at Earth. Since then, our knowledge of bursts has increased dramatically, and the fireball theory (Rees and Meszaros 1992) has received a spectacular confirmation with the accurate description of the properties of afterglows (see Piran 1999, 2000 for reviews). Independently of this, the recently reported data from HiRes (Abu–Zayyad et al., 2002) appear once again to support the view that the origin of the Ultra High Energy Cosmic Rays (UHECRs), at least as observed at Earth, is to be found among astrophysical sources. These two facts together make it worthwhile to re-examine the objections that were raised against the above–mentioned scenario, in the wake of its first appearance. Also, an alternative model for the production of UHECRs in GRBs has been proposed (Waxman 1995), and thus it seems pertinent to discuss some differences between these two models, which will be subject to different observational tests.

The plan of the paper is as follows: in Section 2, we address the criticisms levied against the model in the past years which concern the acceleration mechanism. In Section 3, we deal with global energetics. We discuss in Section 4 two caveats on the previous discussion on energetics. Doublets and triplet in the data reported by AGASA, are discussed in Section 5; we compare our work to that of other authors in Section 6, then we summarize our conclusions.

## 2 Acceleration

The details of the model for the acceleration of non–thermal particles remain essentially those outlined in Vietri (1995), and in Dermer (2002a): particles are accelerated around the fireball’s external shock, which is highly relativistic for most of the relevant time. The major reason for this is that the basic fireball model employed (Meszaros, Laguna, and Rees, 1993, where a detailed model for the external shock evolution is presented, together with the physical characteristics of the plasma behind the shock) has proved itself able to explain the characteristics of the afterglow (Piran 1999, 2000), and it thus appears in need of no major revision. The model was used mostly in determining the largest energy which could bounce off the ejecta shell, without crossing it unscathed (and also, in showing that no radiative losses where incurred into, by the highest energy particles). We derived for the highest energy

$$E_{max}=10^{20}eVE_{51}^{1/3}\mathrm{\Gamma }_2^{1/3}n_1^{5/6}ϵ_B^{1/2}\frac{1}{1\mathrm{cos}\theta }$$
(1)
which differs from the original (Eq. 34 in Vietri 1995) because a more modern initial Lorentz factor $`\mathrm{\Gamma }=100`$ has been used, the true energy release $`E`$ and beaming angle $`\theta `$ have been introduced. The above equation defines the largest energy which can be turned back in a typical wind which produces GRBs, as seen in the upstream frame (the laboratory); the very important question of why we should see particles leaving the shock region through the upstream frame rather than the more customary downstream frame will be addressed shortly; let us just remark that particles leaving through the downstream region will be seen in the laboratory frame as having an energy roughly $`\mathrm{\Gamma }^2`$ smaller than the above.

The average length for deflecting back a non–thermal particle has been taken as $`40`$ times its Larmor radius, but the computation holds only for $`\mathrm{\Gamma }1`$. It was shown in Vietri (1995) that losses due to all radiative processes are negligible. Furthermore, no particle loss may occur through the sides of the ejecta (which would be conceivable, because the emission is beamed along a cone of semiopening angle $`\theta `$); in fact, in the comoving frame, the ejecta have a (radial) width given by $`r/4\mathrm{\Gamma }`$, where $`r`$ is the shock’s distance from the explosion site, while the ejecta width in the transverse direction is $`r\theta `$, and, since $`1/\mathrm{\Gamma }<\theta `$ (Panaitescu and Kumar 2002) it will be easier for the particles to leave through the bottom end than through the sides. Recently, Panaitescu and Kumar (2002) have obtained values for the circumburst’s density $`n`$, for the departures of the magnetic field from its equipartition value $`ϵ_B`$, for the total energy $`E`$ release and beam opening angle $`\theta `$, for a variety of bursts’s afterglows for which sufficient time and frequency coverage were available; using the data in their Tables 2 and 3, we compile our Table 1, illustrating the largest energies achievable within this scenario, which are conveniently large with respect to observational data.

A critique against this argument was levied by Gallant and Achterberg (1999), who pointed out that the limiting factor for the highest energy achievable is due to the ability of the insterstellar medium to deflect high–energy particles toward the shock, and not the vice versa, which is what was computed above. Their analysis showed that this highest energy was

$$E_{max}10^{15}eV\frac{B}{1\mu G}.$$
(2)
This argument however assumes a definite value for the magnetic field around the burst, which is very model–dependent: there is no compelling reason to assume for $`B`$ the value of the average ISM. In particular, two of the leading models for GRBs, the binary neutron star merger (Narayan, Paczynski and Piran 1992), and the SupraNova (Vietri and Stella 1998) are clearly at odds with this value. In fact, these two models make very similar predictions, because in both cases the GRB goes off inside a pulsar wind bubble (PWB, Königl and Granot 2002; Inoue, Guetta and Pacini 2001, Guetta and Granot 2002, but the same, bare proposal was made also in Vietri 1995, and in Gallant and Achterberg 1999), filled with an electron/positron gas, and its associated magnetic field. In either case, a pulsar–type dipolar magnetic field ($`B=10^{12}10^{13}G`$) has a light cylinder at $`r_{lc}10^7cm`$, and decreases like $`Br^1`$ (at least in its toroidal component), from there on, as dictated by standard pulsar electrodynamics (Contopoulos, Kazanas and Fendt, 1999). This results in fields of order $`B=0.110G`$ at the radii appropriate for afterglows ($`r=10^{16}10^{17}cm`$). These estimates are in agreement with determinations of the magnetic field of the (much slower rotating!) Vela bubble (Helfand, Gotthelf and Halpern 2001, but see Arons 2002 for a different estimate of the magnetic field).

In order to adapt Eq. 2 to our needs, we repeat Gallant and Achterberg’s computations for $`Br^1`$. We consider a particle leaving the shock at radius $`R_0`$ with Lorentz factor $`\gamma `$, and moving initially perpendicular to the shock, in the observer’s frame; the shock has Lorentz factor $`\mathrm{\Gamma }_s=\mathrm{\Gamma }_s(R)1`$, so that most particles will be emitted perpendicular to the shock (Bednarz and Ostrowski 1998). We also have allowed for the shock deceleration. Because of a component of the magnetic field $`B_{}(R)=B_1R_1/R`$ perpendicular to the shock normal, and thus to the particle initial direction of motion, the particle speed perpendicular to the shock $`v_z`$ changes according to

$$\frac{v_z}{c}1\frac{e}{m\gamma c^2}_{R_0}^RB_{}𝑑R=1\frac{eB_1R_1}{m\gamma c^2}\mathrm{ln}(R/R_0).$$
(3)
The shock moves instead with speed $`v_{shock}/c11/2\mathrm{\Gamma }_s^2`$. It is necessary that the shock overcomes the particle promptly; only in this case, in fact, the returning particle will see a shock with the same Lorentz factor as when it left the post–shock region, and conventional Fermi–Type acceleration will take place; if this were not the case, we would have a strongly evolutionary shock, and the ordinary predictions of steady–state theory would not follow (and, most likely, the shock would not be an efficient particle accelerator). We must thus have $`v_{shock}dR>v_zdR`$, or, using the equations above,

$$\mathrm{ln}(R/R_0)>\frac{m\gamma c^2}{2\mathrm{\Gamma }_s^2eB_1R_1}$$
(4)
where $`R`$ is the radius at which the particle recrosses the shock. Demanding again that the shock promptly overcomes the particle means that $`\mathrm{ln}(R/R_0)\delta R/R_01`$, and thus we find that only particles with

$$E_p=m\gamma c^2<2\mathrm{\Gamma }_s^2eB_1R_1E_{mb}$$
(5)
are accelerated by the usual Fermi mechanism. Particles with $`E>E_{mb}`$ will find a weaker, slower shock if they come back, and will not be significantly accelerated any more. Because of Eq. 10, we also see from the above that the condition that particles never come back (assuming that the PWB merges into the ISM field for $`R10^{18}cm`$) is only marginally stronger than the above: we find

$$E_pE_{mb}=2\mathrm{\Gamma }_s^2eB_1R_1.$$
(6)
Gallant and Achterberg (1999) used $`R_1=R_d=10^{16}cm`$, the approximate radius at which the afterglow begins, the typical ISM field $`B=1\mu G`$, and the initial shock Lorentz factor to obtain Eq. 2. But for a PWB, following the discussion above we may take for the field outside the light cylinder radius (Eq. 31 of Goldreich and Julian 1969):

$$B(R_1)R_1=10^{17}Gcm\frac{B}{5\times 10^{12}G}\left(\frac{1ms}{P}\right)^2\mathrm{sin}\theta \frac{\mathrm{\Psi }(\theta )}{I}.$$
(7)
Here $`P`$ is the pulsar’s rotation period, which we have scaled to $`1msec`$, which is about right for both neutron star mergers and SupraNovae; we took for the dipolar field of the neutron star a standard value, $`\theta `$ is the beams’ semiopening angle and $`\mathrm{\Psi }(\theta )/I`$ is a normalized fudge factor which Gooldreich and Julian (1969) do not estimate. However, inspection of Fig. 3 of Contopoulos, Kazanas and Fendt (1999) shows that $`\mathrm{\Psi }(\theta )/I(\mathrm{sin}\theta )^2`$; using $`\theta 0.1`$ from Frail et al. (2001) and Panaitescu and Kumar (2002), we find

$$B(R_1)R_110^{14}Gcm,$$
(8)
independent of $`R_1`$. Strictly speaking, we have $`100\mathrm{\Gamma }_s1`$, but for a Newtonian shock particle acceleration proceeds on a slower, diffusive timescale (Meli and Quenby 2002; see also the discussion later in this section). From Fig. 2 of Kirk et al. (2000), we see that the particle spectral index departs from it asymptotic value for very large shock Lorentz factors when $`\mathrm{\Gamma }_s2`$. So we shall use $`100\mathrm{\Gamma }_s1`$, but keeping this caveat in mind:

$$10^{17}eVE_{mb}6\times 10^{20}eV,$$
(9)
depending on the instantaneous value of $`\mathrm{\Gamma }`$. Since we know that (Piran 2000)

$$\mathrm{\Gamma }_s^2=\mathrm{\Gamma }^2\left(\frac{R_d}{R}\right)^3$$
(10)
where $`\mathrm{\Gamma }_s100`$ (Panaitescu and Kumar 2002) is the initial shock Lorentz factor, we see that $`E_{mb}`$ spans continuously the range above, as $`R`$, the instantaneous shock position, changes.

We thus have the following picture. So long as $`E_{mb}>E_{max}`$, particles cannot escape in the upstream direction: they are promptly overtaken by the shock because they are effectively kicked backwards by the environment magnetic field, and continue their acceleration cycle until they reach an energy $`E_{max}`$. When this occurs, the particles cannot be effectively deflected by the ejecta’s magnetic field: they will cross unhindered the hyperrelativistic shell, and be lost through the back of the shell. When this occurs, their energies in the laboratory (upstream) frame are a full factor $`1/\mathrm{\Gamma }^2`$ lower than given by Eq. 1, which holds only for particles leaving the shell in the forward direction. However, when $`E_{mb}<E_{max}`$, particles which have been accelerated up to $`E_{mb}`$ (in the laboratory=upstream frame) cannot be deflected backwards toward the shock any longer, and will be able to escape in the upstream direction without further acceleration. The critical moment occurs when $`E_{mb}=E_{max}`$ (where now in Eq. 1 the shell Lorentz factor must be taken to vary as in Eq. 10): equating the two, we find that $`\mathrm{\Gamma }30`$, and thus that $`E_{max}`$ at this moment is a factor $`1.5`$ lower than Eq. 1 (and thus, within the accuracy of these estimates, it is identical to Eq. 1). We thus have that the energy range for particles escaping in the upstream direction is

$$10^{17}eVEE_{max}.$$
(11)
It should be pointed out that the above picture differs from Gallant and Achterberg’s (1999).

We now have to derive the spectrum of these particles. Here there are two distinct possibilities. On the one hand, one may take Vietri’s (1995) and Gallant and Achterberg’s (1999) view that the particles which end up as UHECRs are the relativistic particles pre–existing the burst which must necessarily exist in a PWB. Under this assumption, Gallant and Achterberg showed that the spectrum will be of the form

$$dn(E)E^2dE.$$
(12)
On the other hand, one may assume that the total number of energetic particles is dwarfed by the (formerly thermal) particles injected at the shock, in which case the spectrum of particles can be derived as follows. Assume that a fraction $`q`$ of the total fireball energy $`E_f`$ is converted on the timescale $`t_{acc}`$ in particles of energy given by Eq. 6. There will then be a total amount of energy

$$\frac{qE_f}{t_{acc}}dt=\frac{qE_f}{2}\frac{dR}{\mathrm{\Gamma }_sR}$$
(13)
in particles of energy $``$ given by Eq. 6. Eliminating $`\mathrm{\Gamma }_s`$ and $`R`$, and using (Gallant and Achterberg 1999)

$$t_{acc}t_u=\frac{E_p}{\mathrm{\Gamma }_seBc}$$
(14)
where $`t_u`$ is the upstream residence time, we find a particle spectrum

$$dnE^{5/2}dE,$$
(15)
Putting together all of the above, we shall take a spectrum of the form

$$dnE^\gamma ;\mathrm{\hspace{0.33em}\hspace{0.33em}2}\gamma 2.5.$$
(16)
We shall then let $`\gamma `$ be a free a parameter in the range above, and use it to fit observations, in the next section.

It is worth remarking at this point that the above results depends upon our having assumed that the GRB goes off inside a PWB (Königl and Granot 2002), a natural assumption when GRBs are due to netron star mergers or to SupraNovae. While currently neutron star binary mergers are not favored as an explanation for long–duration GRBs, they appear attractive as a model for short GRBs; while SupraNovae are one of the leading contenders for the long bursts (see Dermer 2002b for a review). But the above discussion applies, at least qualitatively as a lower limit, to a wide class of models, including Usov’s (1992) model of a fast pulsar with exceptionaly strong magnetic field, models based upon ADIOS flows (Blandford and Begelamn 1998,Blandford 2002), Spruit’s models based on MHD winds (Drenkhahn and Spruit 2002, and references therein), or several distinct models all involving magnetized neutron stars (Duncan and Thompson 1992, Dai and Lu 1998, Blackman and Yi 1998).

One may wonder whether losses in the somewhat more photon- and magnetic field- rich environment provided by a PWB make the estimate of Eq. 1 unrealistically optimistic. This was studied by Inoue, Guetta and Pacini (2002), who proved this to be the case only under exceptional circumstances. It was in fact their aim to explain the high–energy ($`>1GeV`$) emission occasionally seen by EGRET in GRBs, and this could be achieved only with rather extreme assumptions on the environment parameters. This is not at all in contradiction with the above discussion: $`GeV`$ emission in GRBs is a rare phenomenon, so it makes sense that models to explain it require extreme values for the parameters in play.

A different critique was levied (Ostrowski and Bednarz 2002) concerning the spectral slope for particles accelerated at highly relativistic shocks. These authors pointed out that the claims of universality for the spectral index $`s4.2`$ (Kirk and Schneider 1987, Heavens and Drury 1988, Bednarz and Ostrowski 1998, Kirk et al., 2000) possibly overstate their case, because they neglect the fact that, as the regular magnetic field in the downstream region becomes stronger than the field associated with the turbulent component, the spectral slope must increase, the spectrum softens, and the claim to universality evaporates. This criticism, however, can be countered in two (orthogonal) ways. On the one hand, it is possible that it applies to a somewhat idealized situation, whereby post–shock material is incapable of generating a small–scale magnetic field, and thus the field compressed at the shock dominates the turbulent component. This is the view of Medvedev and Loeb (1999) or Thompson and Madau (2000), who proposed ways in which magnetic field of non–negligible amplitude could be generated around the shock, and would be intrinsically turbulent. On the other hand, one may take the view of Königl and Granot (2002), according to whom the very large values of $`ϵ_B`$ derived for some bursts (in particular, GRB 970508 has $`ϵ_B=0.1`$ as determined by several independent studies: Wijers and Galama 1999, Granot, Piran and Sari 1999, Chevalier and Li 2000, Panaitescu and Kumar 2002) cannot possibly be accounted for in this way, but are instead easily accounted for by the compression of a magnetic field surrounding the burst explosion site, provided this occurs in a PWB, which again occurs naturally if bursts are due to binary neutron star mergers, or SupraNovae. The dominant component of this magnetic field is of course toroidal (Goldreich and Julian 1969), and its sweep up and compression would give a post–shock magnetic field mostly oriented along the shock surface, which would increase, not decrease, the shock effectiveness in the acceleration of particles.

Another criticism was that the acceleration time scale was shorter than hypothesized in Vietri (1995), because the average energy gain $`<G>`$ does not scale as $`\mathrm{\Gamma }^2`$, where $`\mathrm{\Gamma }`$ ($`100`$, Panaitescu and Kumar 2002) is the shock Lorentz factor. It was pointed out that this energy gain applies to the first shock crossing, but the following crossings would provide a more mundane $`<G>2`$ (Gallant and Achterberg 1999). This criticism is surely correct, but does not spoil the gist of the original argument, that the acceleration time scale is much shorter than for Newtonian shocks. Indeed, the fact that the energy gain per shock crossing is no longer infinitesimal, as it is for Newtonian shocks (where $`G1=v_{shock}/c1`$), implies that the required number of shock crossings for acceleration to the EeV range becomes finite: from $`m\mathrm{\Gamma }^2c^2`$ to the EeV range only $`20`$ shock crossings are required, in agreement with the cosiderable reduction of the acceleration time scale seen in numerical simulations (Meli and Quenby 2002). This very fact has led some of the severest critics of the scenario to propose themselves that UHECRs are accelerated at highly relativistic shocks (Kirk, Achterberg, Guthmann and Gallant 2000).

There is a deep reason for this: we showed elsewhere (Vietri 2002) that, for particle acceleration around shocks of any speed, the following relation holds:

$$P<G^{3s}>=1.$$
(17)
Here $`P`$ is the probability of returning to the shock, $`G`$ the energy amplification for one cycle across the shock, $`s`$ the distribution function spectral index, and the average is to be taken after the raising to the $`3s`$ power; the above generalizes to relativistic shocks Bell’s (1978) relation for the spectral index. The above relationship holds irrespective of the assumed scattering law. In the case of Newtonian shocks, both $`P`$ and $`G`$ differ from unity by quantities $`𝒪(v/c)1`$, and the two combine to give $`s=4`$. In the hyperrelativistic case, both $`P`$ and $`G`$ differ significantly from unity, so that, since $`P1/2`$ (the post–shock fluid runs away from the shock with speed $`c/3`$, so that it is easy for particles to be advected away), also $`G1`$ must grow away from its infinitesimal Newtonian values, and tends to $`G2`$. This, of course, shortens the acceleration time scale, and is a robust prediction of the theory of acceleration around hyperrelativistic shocks, independent of the scattering law.

## 3 Energetics

Another criticsm which has been levied against this scenario is that, contrary to our original claims, there is a mismatch between the observations of the energy injection rate in $`\gamma `$–ray photons due to GRBs, and in UHECRs. This claim has been made independently by Stecker (2000) and Berezinsky et al. (2002), who pointed out that the realization that the redshift distribution of GRBs should follow the star formation history in the Universe must lower the local energy injection rate. In particular, Stecker found the evolving GRB population to produce a $`\gamma `$ ray photon flux which is an order of magnitude lower than that of just the fraction of UHECRs beyond $`10^{20}eV`$, let alone the rest of the spectrum.. Waxman (2002) made some quite reasonable criticisms of the paper by Berezinsky et al. (2002). We consider here Stecker’s work. We shall first show that Stecker’s claims are not supported by currently available data, then derive the spectrum of UHECRs if these are emitted by GRBs, and then comment on the (rather large) uncertainties still surrounding the overall energetics of GRBs.

Stecker’s (2000) estimates are based upon old data: he used the AGASA flux at $`E>10^{20}eV`$ (Takeda et al., 1998, Hayashida et al., 1999) which he took to be an unbroken power–law with index $`2.78`$ beyond $`10^{20}eV`$, while Abu–Zayyad et al. (2002) explicitly state that their data beyond $`10^{19.8}eV`$ depart from such a law, and state that the number of events they observe is a full factor of $`4/19.10.21`$ below the extrapolation of such law. Thus, while he took the local energy release rate in UHECRs to be

$$\dot{ϵ}_{(AGASA)}=2\times 10^{44}ergMpc^3yr^1,$$
(18)
HiRes data suggest a value a factor of $`0.21`$ lower:

$$\dot{ϵ}_{(HiRes)}=4\times 10^{43}ergMpc^3yr^1.$$
(19)
which compares already much better with the value he took for the energy release rate in $`\gamma `$–ray photons, $`2.2\times 10^{43}ergMpc^3yr^1`$. But even this value is superseded, since it is based on Schmidt’s (1999) estimate of the GRB rate. In his new work, Schmidt (2001) gives the burst’s rate (for isotropic bursts) as

$$\dot{n}_{GRB}=0.5Gpc^3yr^1.$$
(20)
As correctly pointed out by Stecker (private communication), the above value includes short bursts, about which, strictly speaking, we know very little. In particular, we ignore what their overall energy budget is, and, most important to this paper, whether they have relativistic afterglows at all. So to play it safe, we reduce the above rate to exclude short bursts: According to Fishman and Meegan (1995) short bursts account for $`30\%`$ of all bursts, and we consequently take the local rate of long GRBs as

$$\dot{n}_{lGRB}=0.35Gpc^3yr^1$$
(21)
Schmidt does not give the average (isotropic) burst energy release, but we can get this from Frail et al. (2001), who, among many other things, give the isotropic burst energy release for all bursts with known redshifts. The reason why we prefer this over using Schmidt’s work is that their total energy releases (taken from Bloom, Frail and Sari 2001) cover the spectral range $`0.22000keV`$, while Schmidt restricts himself to the range $`50300keV`$: there seems no obvious reason why we should adopt Schmidt’s values, when more suitable ones are available. From Frail et al.’s Table I, we find

$$E_{av}=3.3\times 10^{53}erg.$$
(22)
We thus find the energy release rate in $`\gamma `$–ray photons by GRBs as

$$\dot{ϵ}_{GRB}=1.1\times 10^{44}ergMpc^3yr^1.$$
(23)
Strictly similar results are obtained in an analogous way by using data from Panaitescu and Kumar (2000) and from Piran et al. (2001). Please notice that, while not all data from these three sets of authors are different (many bursts are in common) they fitted different quantities with different methods, so the results are, for the most part, independent. This is larger than $`\dot{ϵ}_{(HiRes)}`$ by a hefty factor of $`\dot{ϵ}_{GRB}/\dot{ϵ}_{(HiRes)}=3`$, and allows a full explanation of UHECRs by GRBs, even in the case in which the efficiency for particle acceleration is smaller than for photon production.

We thus consider Stecker’s (2000) objections superseded by new data. However, we are still left with the question of how much energy is produced by GRBs in the form of UHECRs in the whole energy range, Eq. 11, and not just for $`E>10^{20}eV`$, and whether the observed spectrum of UHECRs at Earth can be fitted by a cosmological distribution of GRBs. In order to answer both of these questions, we repeat here a computation of the spectra of UHECRs as observed at Earth,for the evolving population discussed above. This does not  mimick Waxman’s (1995b) computations for several reasons. First, he uses a constant comoving density for GRBs, while, following Schmidt (2001), we use the second star formation history (SFR2) proposed by Porciani and Madau (2001); second, he uses a continuous energy loss approximation, while we take into account the discreteness of the photopion production events (but not of pair production); and third, what we present here are simulations of observations (De Marco, Blasi, Olinto 2003) in order to stress that the fits to observations contain some amount of noise due to the finite number of events now available. For a full discussion of the method, simulations, and approximations, we refer the reader to De Marco, Blasi, Olinto (2003). Following Eq. 16,11, the particle spectrum is taken as

$$n(E)dE=E^\gamma dE;E_{min}EE_{max}$$
(24)
where $`\gamma `$ is a free parameter that is to be varied to allow fitting of Earth–based observations; $`E_{min}`$ is from Eq. 11, while $`E_{max}`$ is irrelevant (provided it exceeds $`10^{20}eV`$, which it does, Eqs. 1,11) because most energy is in the lowest energy particles (for $`\gamma >2`$). The spectrum we obtain is shown in Figs. 1 and 2 is obtained for the two values of $`\gamma =2.2`$ and $`2.5`$. A smaller value of $`\gamma `$ yields only marginally different results. The steeper spectrum, like $`\gamma =2.5`$ fits data below $`10^{19}eV`$ much better, but, because of the obvious possibility of a Galactic contamination, we do not regard this accurate fit as either desirable or likely. We shall thus take the value

$$\gamma =2.2$$
(25)
as our canonical result, For both values of $`\gamma `$ we under-reproduce the AGASA counts beyond $`10^{20}eV`$, but not the HiRes data.

The local energy release rate in UHECRs, for the whole range of Eq. 11 is

$$\dot{ϵ}_{UH}=1.5\times 10^{45}\left(\frac{E_{min}}{10^{17}eV}\right)^{0.2}ergMpc^3yr^1$$
(26)
which exceeds Eq. 23 by a factor of $`13`$. This overestimate is reduced when we consider the suggestion (De Marco, Blasi and Olinto 2003) that the energies of AGASA events are overestimated by a full $`15\%`$, which brings their data in much better agreement with HiRes’. In this case, the total required energy release becomes

$$\dot{ϵ}_{UH}=1\times 10^{45}\left(\frac{E_{min}}{10^{17}eV}\right)^{0.2}ergMpc^3yr^1$$
(27)
which now exceeds Eq. 23 by a factor of $`9`$ only.

The estimate above is much lower than the claims by Scully and Stecker (2002), who hypothesized an excess by a factor of $`10^{23}`$, over the much smaller range $`E>10^{19}eV`$, while the above equation applies to $`E>10^{17}eV`$. The origin for the disagreement between our paper and Scully and Stecker’s is twofold. First, we fit a different, flatter injection spectrum ($`2.2`$ rather than $`2.75`$), but we still obtain excellent fits. Second, we use a $`\dot{ϵ}_{GRB}`$ (Eq. 23) which differs from theirs. However, they do not compute it anew, but they refer to Stecker (2000). We have already shown that this estimate is based upon old data, and will thus not discuss their results any further.

The slight overabundance of energy required to explain UHECRs as observed at Earth with respect to the energy injection rate in photons by GRBs, which we estimated above as included in the range $`1013`$, is not an embarrassment for the present model. Numerical simulations of radiative efficiencies at internal shocks, the currently favoured model for the prompt phase of GRB emission (Spada, Panaitescu and Meszaros 2000) show that these are $`0.005`$, so that even the overproduction of UHECRs by a factor of $`9`$ amounts to a total efficiency for UHECRs of a few percent, surely in line with conventional thinking (Draine and McKee 1993) about emission efficiencies for non–thermal particles around non–relativistic shocks. For relativistic shocks, given the shortening of the acceleration time–scale (Meli and Quenby 2002), acceleration should be, if anything, more efficient.

For this reason, and for the reasons to be discussed in the next paragraph, we are not worried by this slight discrepancy.

## 4 Two caveats

There are two important caveats which need to be made when discussing bursts’ energetics.

Several authors have stressed that GRBs could have a redshift distribution strongly skewed toward large redshifts, but while very likely, this idea has not been confirmed yet. In order to stress this point, we have taken the set of all GRB redshifts known to date, to see whether the low redshift distribution agrees better with a Porciani and Madau–like evolution law (Porciani and Madau 2001), or with a non–evolving one. We restrict our attention to small redshifts ($`z<1`$) for two reasons: first, the Porciani and Madau evolution law is very distinctive at such low redshifts, and is not disputed by different authors; second, because it seems reasonable to assume that this sample is unbiased by the selection effects which obviously mar the higher end of the distribution.

In Fig. 3 we show the observed cumulative distribution in redshift of known GRBs; data are taken from the compilation publicly available at http://www.mpe.mpg.de/j̃cg/grbgen.html, which also gives the relevant references. We compare it with the theoretical expectation that the GRBs explosion rate follows the star formation rate. We consider the three different star formation histories SF1, SF2 and SF3 given by Porciani and Madau (2001). We use in the figure $`\mathrm{\Omega }_M=0.3`$ and $`\mathrm{\Omega }_\mathrm{\Lambda }=0.7`$ and a Hubble constant $`H_0=65h_{65}`$km s<sup>-1</sup> Mpc<sup>-1</sup>. For completeness we compare the observed distribution also with a distribution $`(1+z)^3`$. A KS test shows that the observed distribution without GRB 980425 differs from the theoretical ones SF1, SF2, SF3 and $`(1+z)^3`$ with a significance of 24%, 25%, 28% and 32% respectively, while inclusion of GRB 980425 changes these values to 10%, 11%, 13% and 16%.

Fig. 3 clearly shows that we have right now no compelling reason to believe that GRBs are characterized by very strong evolution with redshift, a uniform distribution providing an equally good, if not better, fit to current data. But furthermore, even if larger samples like the one to be collected by SWIFT should in the future provide a confutation of this (which we very much believe), still the point would remain that a perfectly normal subsample of a dozen or so objects may depart greatly from the average distribution, and appear like a uniform one. In this sense, then, the argument made by Stecker (2000) and Berezinsky et al. (2002) appears as an oversimplification of reality: in discussing sources of UHECRs, one cannot forget that the total number of sources one is dealing with is statistically very small, and likely to stay that way, and thus a proper discussion of the model on an energy basis should include the effects of fluctuations.

Second, there is an important uncertainty in the global energetics of bursts, as we now discuss. Comparing the late–time evolution of bursts (which fixes the total energy in the kinetic form) with the energy released in $`\gamma `$–ray photons, it is possible to deduce the radiative efficiency for several individual bursts (Panaitescu and Kumar 2000). The radiative efficiency thusly derived is very high, always eceeding $`66\%`$ (except for one case) and occasionally even approaching an incredible $`90\%`$ (see Table III of Panaitescu and Kumar). We remark that these derived efficiencies always exceed the largest efficiency achievable in converting rest mass into photons around maximally rotating black holes ($`43\%`$). Although this is not a strict theoretical upper limit (models approaching $`100\%`$ have been proposed, Lazzati et al., 2000), still values derived from observations are so large to throw doubts on the validity of the present version of internal shock model. In fact, the maximum efficiency that can be reached with this model is of the order of 20% (Guetta, Spada and Waxman 2001a) but only under ad hoc assumptions: for random distributions in the wind properties, the efficiency cannot exceed $`.005`$ (Spada, Panaitescu, Meszaros 2000). Perfectly efficient models without internal shocks can be conceived (Lazzati, Ghisellini, Celotti and Rees 2000), but they cannot account for bursts’ spectra (Ghisellini, Lazzati, Celotti, Rees 2000).

For this reason, Kumar and Piran (2000) have proposed an ingenious modification of the fireball spectra, whereby hot spots in the emitting surface occasionally dominate the emission from the whole GRB. These hot spots are not typical of the whole surface; when our line of sight chances through one hot spot, we claim to see a burst, but when otherwise, we wouldn’t recognize the event as such because too dim. In this way, they proposed that there be a significant number of objects (dissimulated GRBs) escaping our detection, and redressing the statistics. The amusing consequence of this is that, while of course this implies we would have to decrease the value of Eq. 23 (thus apparently going in Stecker’s direction) there would be a substantially higher number of afterglows acting as acceleration sites. From Table II of Panaitescu and Kumar, we see that the apparently isotropic kinetic energy of bursts’ afterglows is

$$E_{ag}=2.5\times 10^{53}erg.$$
(28)
This of course is strictly relevant to our acceleration mechanism, because as stated above, we are proposing that UHECRs are accelerated at external (= afterglow) shocks, so that this is the relevant source of energy to tap. Comparison with Eq. 22 shows indeed that the radiative efficiency, defined as $`E_{av}/(E_{av}+E_{ag})0.6`$, very large indeed. Combining this with Schmidt’s rate estimate, we see that there is a kinetic energy injection rate of

$$\dot{ϵ}_{ag}=1.3\times 10^{44}ergMpc^3yr^1,$$
(29)
comparable with Eq. 23. If the argument by Kumar and Piran (2000) is right, Eq. 29 is to be increased. In fact, the kinetic energy per burst in the afterglow is derived from very late observations, when the hot spots are not working any longer (Kumar and Piran 2000), but the rate is to be increased to account for all bursts for which our line of sight does not meet a hot spot. This factor is estimated by Kumar and Piran as roughly a factor of $`100`$, bringing Eq. 29 in agreement with Eq. 27, with a vengeance. This incompleteness factor can be estimated also as follows: conventional models (Spada, Panaitescu and Meszaros 2000) find a radiative efficiency of $`0.005`$; observations give $`0.6`$, so the hot spots bias our statistics by a factor of $`100`$.

## 5 Doublets and triplet

One argument often used against GRBs as sources of UHECRs (and, in general, against all explosive sources) is the presence of doublets and triplet in the data from AGASA (Takeda et al., 1999, Nagano and Watson 2000, Hayashida et al., 1997), where, occasionally, the lower energy event precedes the higher energy one. The argument goes that, since the time delay with respect to photons is a monotonically decreasing function of particle energy, no impulsive source can have produced these clustered events, while obvioulsy no such problem exists with continuous sources. In the light of the present disagreement between the observational results of HiRes and AGASA, one may wonder whether the energy estimates even of the events in the clusters are correct. We shall not take this view, but accept the existence of the doublets and triplet as established, despite the caveat above, and show that the purported inconsistency of the explosive models is not water tight.The argument in fact is an oversimplification, and neglects the fact that the time delay is a statistically distributed quantity, with both a mean and a variance, which we now compute to allow a further discussion.

We take as a model for particle motion a Langevin–like equation

$$\frac{d\stackrel{}{p}}{dt}=m\gamma \frac{d\stackrel{}{v}}{dt}=e\frac{\stackrel{}{v}}{c}\stackrel{}{B}$$
(30)
where, that is, we assume $`\stackrel{}{B}`$ to be a Gaussian field with zero mean and known self–correlation function. Initially, we behave as if $`\stackrel{}{B}`$ were known exactly; we also use the assumption that the particle Larmor radius $`R_L`$ is large with respect to the typical distances over which the field is correlated, say $`r_c`$: $`r_c/R_L1`$, and return later to discuss this approximation. We take $`z`$ as the axis of the initial motion, and decompose in approximate form the equations of motion as

$$m\gamma \frac{dv_z}{dt}e\frac{v_x}{c}B_yeB_x\frac{v_y}{c};m\gamma \frac{d\stackrel{}{v}_{}}{dt}eB_x\widehat{y}eB_y\widehat{x}.$$
(31)
Here we have assumed that the particle’s deflections from the unperturbed trajectory are small (which means that $`v_zc,v_xv_yc`$) and thus the equations are valid to smallest significant order only; obviously, the fields are to be computed along the unperturbed trajectory, $`z`$. For a source located at distance $`D`$ from the observer, we find

$$y=\frac{e}{m\gamma c}_0^D𝑑Z_0^ZB_x(z)𝑑z;x=\frac{e}{m\gamma c}_0^D𝑑Z_0^ZB_y(z)𝑑z$$
(32)
and inserting this into the equation for the $`z`$ component we find

$`z={\displaystyle \frac{e^2}{m^2\gamma ^2c^4}}{\displaystyle _0^D}𝑑Z{\displaystyle _0^Z}𝑑z{\displaystyle _0^z}𝑑y\left(B_x(z)B_x(z+y)+B_y(z)B_y(z+y)\right)=`$
$`{\displaystyle \frac{e^2}{m^2\gamma ^2c^4}}{\displaystyle _0^D}(Dz)𝑑z{\displaystyle _0^z}𝑑y\left(B_x(z)B_x(z+y)+B_y(z)B_y(z+y)\right)`$ (33)
where the last identity has been obtained by means of an integration by parts. We see that $`y,x`$, are linear in the fields, and thus are distributed with zero mean, while the time delay $`z`$ (which is, by definition, the departure from an unperturbed motion with $`\gamma 1`$, and thus equals the time delay with respect to photons), is quadratic in the fields and in the particle charge. Because of this, the time delay has a non–zero mean, which we now compute by means of the above equation: we obtain

$$<z>=\frac{e^2}{m^2\gamma ^2c^4}_0^D(Dz)𝑑z_0^z2B_{xx}(y)𝑑y.$$
(34)
Here $`B_{xx}(s)=<B_x(\stackrel{}{x})B_x(\stackrel{}{x}+s\widehat{z})>`$ is the rms (since we assumed it to have zero mean) of one component of the field pependicular to the line of separation, averaged over all field configurations; we also used isotropy to set $`B_{xx}=B_{yy}`$, and homogeneity to drop the dependence of $`B_{xx}`$ on $`\stackrel{}{x}`$. We show in the Appendix that, under very reasonable physical assumptions,

$$B_{xx}(s)=\frac{1}{3}B_0^2g(\frac{s}{r_c},\alpha )$$
(35)
where

$$g(x,\alpha )\frac{3_0^{\mathrm{}}y^\alpha f(y)𝑑y(\frac{\mathrm{sin}yx}{yx}+\frac{\mathrm{cos}yx}{(yx)^2}\frac{\mathrm{sin}yx}{(yx)^3})}{2_0^{\mathrm{}}y^\alpha f(y)𝑑y}$$
(36)
For $`\alpha =5/3`$ we would have the well–known Kolmogorov–Okubo law.

When we insert the above into Eq. 34, we see that we have to compute the integral

$$_0^z𝑑yg(y/r_c,\alpha ).$$
(37)
In the above, $`z`$ is the ($``$ cosmological) distance travelled by the particle; thus we expect $`zr_c`$, since, as remarked in the Appendix, $`z`$ is of the order of the radius of the Univers today, while $`r_c`$ is of the order of the radius of the Universe at the time of the generation of the magnetic field. We can then approximate

$$_0^z𝑑yg(y/r_c,\alpha )_0^{\mathrm{}}𝑑yg(y/r_c,\alpha )=\frac{\pi r_c}{4}\frac{3_0^{\mathrm{}}y^{\alpha 1}f(y)𝑑y}{2_0^{\mathrm{}}y^\alpha f(y)𝑑y}=\frac{3\pi }{8}(\alpha 1)r_cr_cq(\alpha ).$$
(38)
where we used Eq. 51. $`q(\alpha )`$ is illustrated in Fig. 3. Plugging the above into Eq. 34 we find

$$<z>\frac{q(\alpha )e^2B_0^2D^2r_c}{3m^2\gamma ^2c^4}.$$
(39)
We are interested not just in the mean values, but also in the rms deviations from the mean. For the directions perpendicular to the $`z`$ axis, we easily find with a computation analogous to the above

$$\sqrt{<(y)^2>+<(x)^2>}=\sqrt{\frac{8q(\alpha )}{3}}\frac{eB_0r_cD}{m\gamma c^2}\left(\frac{D}{r_c}\right)^{1/2}.$$
(40)
Two particles arriving on the same experimental apparatus will, in average, be separated by $`\sqrt{2}`$ of the above, and this can be expressed in units of $`r_c`$:

$$\delta n=\sqrt{\frac{q(\alpha )}{3}}\frac{4eB_0D}{m\gamma c^2}\left(\frac{D}{r_c}\right)^{1/2}.$$
(41)
This adimensional quantity $`\delta n`$ is the number of average correlation lengths which separate two average particles reaching the same experiment. It essentially measures whether two particles of the same energy, have travelled through the same patches of magnetic field: this is so, whenever $`\delta n1`$, or otherwise when $`\delta n1`$.

In the case $`\delta n1`$, we may reasonably assume that the particles of the same energy emitted simulataneously in an explosive event, will arrive at the experiment also simultaneously, with the time–delay derived above. In the case $`\delta n1`$ or greater, since these particles will have travelled through patches of uncorrelated magnetic field, we cannot assume any longer that the particles will arrive simultaneously; a more reasonable assumption is that their time delays be distributed according to the statistics implicit in Eq. 5. We will present the full statistics elsewhere, and will limit ourselves here to derive the rms deviation of the time delay around its mean. We proceed as follows: using Eq. 5 we find

$`(z)^2={\displaystyle \frac{e^4}{m^4\gamma ^4c^8}}{\displaystyle _0^D}𝑑z{\displaystyle _0^D}𝑑z^{}(Dz)(Dz^{}){\displaystyle _0^z}𝑑s{\displaystyle _0^z^{}}𝑑s^{}`$
$`(B_y(z)B_y(z+s)+B_x(z)B_x(z+s))(B_y(z^{})B_y(z^{}+s^{})+B_x(z^{})B_x(z^{}+s^{})).`$ (42)
We now average the above over all field configurations; to do so, we use Wick’s theorem since we assumed $`\stackrel{}{B}`$ to be a Gaussian field:

$`<(B_y(z)B_y(z+s)+B_x(z)B_x(z+s))(B_y(z^{})B_y(z^{}+s^{})+B_x(z^{})B_x(z^{}+s^{}))>=`$
$`2(2B_{xx}(s)B_{xx}(s^{})+B_{xx}(|z^{}z|)B_{xx}(|ss^{}|)+B_{xx}(|zz^{}s^{}|)B_{xx}(|z+sz^{}|).`$ (43)
We recognize in the first term to the right of the equality the term which leads to $`(<z>)^2`$; since it is our aim to derive the rms in the time delay, we drop it from now on. The second term on the rhs of the above can be seen to vanish (or more precisely, to give exponentially small terms when $`B_{tt}\mathrm{exp}(D/r_c)`$), and will thus be neglected. Putting together all of the above, we find

$$\mathrm{\Delta }_{rms}<(z)^2>(<z>)^2=\frac{q^2(\alpha )e^4r_c^2B_0^4D^4}{18m^4\gamma ^4c^8}$$
(44)
and we now define a convenient quantity as

$$\mathrm{\Delta }\frac{\sqrt{\mathrm{\Delta }_{rms}}}{<z>}=\frac{1}{\sqrt{2}},$$
(45)
which is the result we were searching for. It is worthwhile remarking that this is independent of all quantities, so long as our approximations (Eq. 38) are satisfied.

We now may begin our discussion by remarking that propagation in the IGM clearly satisifies the hypothesis under which these results were derived, i.e., $`r_c/R_L1`$ and homogeneity, isotropy, and gaussianity (HIG, from now on) of the field (see Grasso and Rubinstein 2001 for a review), so that the computation is immediately applicable. We find for Eq. 41:

$$\delta n=1\frac{B_0}{2\times 10^{11}G}\frac{D}{100Mpc}\frac{3\times 10^{10}}{\gamma }\left(\frac{D}{100r_c}\right)^{1/2}$$
(46)
showing that, for very reasonable values of the parameters involved, particles which reach Earth experiments cross regions with uncorrelated magnetic field. In this case, these particles’ arrival times are spread out around their mean value (Eq. 5) with an rms of order $`\mathrm{\Delta }`$; if doublets and triplets include particles with energies different by a factor $`p`$, their mean arrival times differ by a factor $`11/p^2`$, and this corresponds to a splitting of $`(11/p^2)/\mathrm{\Delta }`$ times the rms arrival time, for the lower energy ones. For typical observed values of $`p2`$, we find that the separation of the means in arrival times is only $`3\sqrt{2}/41`$, implying a separation of $`1\sigma `$ only. This shows that there must be a superposition between the higher energy ones (assumed here for sake of argument to be travelling without dispersion in arrival times) and the lower energy ones (assumed instead to be spread out). Also, it is worthwhile to remark that the flux suppression at the lower energy due to statistics ($`e^2`$, as will be shown in aforthcoming paper) is nearly exactly compensated for by larger number of low energy particles: for a spectrum going like $`E^3`$ at Earth, the product of the two, $`e^2p^32`$, again in agreement with observations.

However, besides crossing the IGM, any particle reaching the Earth must also cross the Local Supercluster, and its associated magnetic field. Since the Local SuperCluster is obviously not a virialized structure, we may reasonably assume that HIG holds again. The current value of $`B`$ is rather uncertain. If we take $`D=10Mpc`$, $`B3\times 10^8G`$ and $`r_c=1Mpc`$ as discussed by Blasi and Olinto (1999), we find $`r_c/R_L=1/3`$, so that the above computations are still (marginally) applicable; also $`\delta n10`$. A more recent determination (Vallee 2002) has $`B1\mu G`$, and $`r_c=100kpc`$, so that again $`r_c/R_L=1`$, and the computation is marginally applicable; for these parameter values, $`\delta n100`$.

The bottom line of this argument is that, either in the IGM (for reasonable parameter values) or in the Local Supercluster (for currently favored, but still very uncertain parameter values), the protons which reach the same experimental apparatus on Earth will have travelled through disconnected patches of magnetic field, and their time of arrivals will be spread around the mean by their large (Eq. 45) variance, and will mostly wash out the strict time ordering which is to be expected in the case of particles of (formally) infinite energy.

A different, and important question, is whether the statistics implicit in Eq. 5 is consistent with the several instances of doublets, and triplet, observed, and the number of observed inversion in the ordering of arrivals. A full answer to this question requires working out the full statistics implicit in Eq. 5, and the comparison with the real data; the details of this work lie outside the scope of this work, and will be presented elsewhere. But it is nonetheless worth remarking that, on the basis of the arguments presented above, we expect roughly equal occurrences of multiplets with lower (higher) energy particles preceding higher (lower) energy particles. From Table 2 of Takeda et al. (1999) we see that this is precisely the case: in multiplets C1,C4,C5 the higher energy particle precedes the lower energy one, and the reverse occurs in multiplets C2 and C3.

## 6 Comparison with other work

Simultaneously with, and independently of Vietri (1995), two other papers (Milgrom and Usov 1995, Waxman 1995) proposed that UHECRs were accelerated in GRBs, but, while Milgrom and Usov did not propose a specific acceleration mechanism, Waxman (1995) proposed that Type II Fermi acceleration at the subrelativistic internal shocks could account for the acceleration of UHECRs.

It is however still controversial whether the internal shocks can account for the GRB emission in the prompt phase. Spectra of GRBs display both a universal break feature at $`200keV`$ (Preece et al., 2000) and a wildly variable low energy spectral index (Crider et al., 1997, Preece et al., 1998, Ghirlanda, Celotti and Ghisellini 2002), which have proved difficult to reproduce theoretically: to wit, look at Fig. 1 of Meszaros and Rees (2000). In particular, it seems that pure synchrotron emission is incapable of explaining at least some (but not all! Tavani 1996, 1997) of the bursts’ spectral properties (Ghisellini, Celotti, Lazzati 2000). Possible alternatives have been discussed (Pilla and Loeb 1998, Ghisellini and Celotti, 1999, Panaitescu and Meszaros 2000) but none of them is widely credited with being able to explain the above–mentioned properties of the observed bursts. These difficulties should be contrasted with the success of the external shock model, which is widely credited with the explanation of the properties of the afterglow (Piran 2000).

Since internal shocks are Newtonian (Waxman 2002), we do not expect the acceleration time scales to be short enough: the large average energy gain discussed above ($`G2`$, Gallant and Achterberg 1999) only applies to relativistic shocks, and the ability of internal shocks to provide UHECRs in the short times allotted remains to be proved. This is especially so, since Waxman (1995, 2002) is appealing to the much less efficient Type II acceleration, rather than the standard Type I. But, even if TypeII acceleration at a Newtonian shock should prove rather fast, still numerical simulations (see e.g. Meli and Quenby 2002) show explicitly that acceleration at relativistic shocks is much faster than Type I acceleration at Newtonian shocks, and this, in its own turn, is much faster than Type II acceleration (again at a Newtonian shock). Thus we may always expect Type I acceleration at external shocks to dominate Waxman’s assumed mechanism.

There is no obvious reason to suspect that particle acceleration at internal shocks should be particularly effective: as emphasized by Waxman (2002), internal shocks are Newtonian, and thus a strong counterpressure due to particles will result in a number of instabilities (Ryu, Kang and Jones 1993, and references therein) and in the smoothing out of the shock. The major consequence of this would be the suppression of the lower energy electrons which are supposed to radiate away the winds’ internal energy, and thus the muting of the burst. Acceleration of UHECRs at external shocks, instead, suffers no such problem <sup>1</sup><sup>1</sup>1We are indebted to P. Blasi for pointing this out to us.: since it takes only a deflection by an angle $`𝒪(1/\mathrm{\Gamma })`$ to allow the shock to overtake a particle in the upstream region, as $`\mathrm{\Gamma }\mathrm{}`$ the linear momentum communicated to the upstream section vanishes as $`1/\mathrm{\Gamma }^2`$, and no counterpressure/smoothing of the shock/instabilities ensue.

The simple energy budget favors acceleration at external shocks. The apparently large radiative efficiency of bursts, contrary to all models for internal shocks (Spada, Panaitescu, Meszaros 2000), has prompted Kumar and Piran (2000) to propose an ingenious model to salvage internal shocks. According to them, hot spots in the shock surface strongly bias our detection abilities in their favor, and make it look like the overall radiative efficiency approaches unity; the many bursts without hot spots aiming at us, which are necessary to redress the statistic, would not be detected. In this model, the true radiative efficiency in the internal shock phase would be about $`0.01`$ of what is observed, and there should be about $`1/0.01=100`$ more GRBs than we detect. But what this means is that there would be $`100`$ times as many external shocks from which to draw UHECRs, because external shocks, for which beaming effects are smaller, (and thus also the hot spot mechanism would not be operational, Kumar and Piran 2000) suffer through no such criticism (Vietri 1998a). In other words: either one believes the largish ($`0.18`$ for GRB 970508) radiative efficiencies, in which case it is hard to believe in the internal shock model; or, one tries to salvage the internal shock model, in which case the energy balance is dominated by a large margin by the kinetic energy of the afterglow phase, which can be tapped by acceleration of UHECRs at external shocks.

Apart from these objections, it is not easy to conceive at this point an observational test capable of distinguishing between the two models. The only hope appears to be the production of high energy neutrinos which must accompany the in situ acceleration of particles; occasionally, in fact, UHECRs will produce pions (and then neutrinos) through collisions with photons, in the moderately photon rich environment provided by the post–shock shells. If UHECRs are accelerated at internal shocks, the neutrinos thusly produced will arrive at Earth simultaneously with the photons of the burst proper and will have an energy $`10^{15}eV`$ (Waxman and Bachall 1997, Guetta Spada Waxman 2001b). If UHECRs are accelerated at external shocks they will arrive at Earth simultaneously with the photons of the afterglow and will have a higher energy, $`10^{17}`$ (Vietri 1998a, 1998b).

## 7 Conclusions

In short, an analysis of the objections raised against acceleration of UHECRs at the external shocks of GRBs shows that none of them is, at the time of writing, dangerous for the model.

The strongest objection levied so far is that due to Gallant and Achterberg (1999) that the highest energies of particles accelerated at external shocks is smaller than surmised in Vietri (1995); while in fact we pointed out that this objection is model–dependent, so is the reply, and ultimately the success of the scenario depends upon whether models in which a GRB explodes in a PWB (as discussed by Königl and Granot 2002) or some magnetic–field–rich environment (Blandford and Begelman 1998, Usov 1992, Drenkhahn and Spruit 2002, Duncan and Thompson 1992, Dai and Lu 1998, Blackman and Yi 1998) are correct, or not.

The strongest success of the model is once again the close agreement between the energy release in $`\gamma `$–ray photons and that necessary to account for UHECRs at Earth: to wit, compare Eqs. 27 and 29 (including the whole discussion, in Section 4). This agreement requires slightly different production efficiencies for photons and high energy particles, but is otherwise unique to this model.

It is also worth remarking that much information is contained in the observed properties of UHECRs: as noted above, in fact, multiplets contain information on the properties of the magnetic field along the line of sight which is peculiar to strongly time–varying sources, and which, coupled with other observed features such as the spectrum and the overall number of sources on the plane of the sky, can put the model through a severe test. A detailed analysis of this kind will be presented elsewhere.

It is a pleasure to acknowledge extremely helpful conversations with P. Blasi.

## Appendix

We wish to determine here the correlator for the Gaussian field:

$$B_{ij}(\stackrel{}{s})<B_i(\stackrel{}{x})B_j(\stackrel{}{x}+\stackrel{}{s})>;.$$
(47)
We start out by considering the Fourier transform of $`B_{ij}`$:

$$b_{ij}(\stackrel{}{k})B_{ij}(\stackrel{}{s})e^{ı\stackrel{}{k}\stackrel{}{s}}d^3s.$$
(48)
Because of homogeneity, isotropy, and Maxwell’s equation $`\stackrel{}{B}=0`$ one finds (Landau and Lifshitz 1987)

$$b_{ij}(\stackrel{}{k})=\frac{2\pi ^2E(k)}{k^2}(\delta _{ij}\frac{k_ik_j}{k^2})$$
(49)
where $`E(k)`$ is the mode energy distribution. If $`E(k)`$ were to follow a pure Kolmogorov–Okubo law, we would have $`E(k)=Ak^{5/3}`$, but we have to modify this in two ways. First, the magnetic field outside of galaxy clusters cannot be due to the usual small scale turbulence that we are accustomed to in the laboratory, and thus the exponent in the above law need not equal $`5/3`$. Second, if phase transitions at an early epoch (Grasso and Rubinstein 2001) are responsible for the generation of the magnetic field outside galaxy clusters, the field cannot be correlated over lengthscales exceeding the horizon scale at field generation. So we need to introduce a cutoff $`f(r_ck)`$ for all modes with $`k<1/r_c`$, for some $`r_c`$ which we shall assume known from now on. We thus take

$$E(k)=Ak^\alpha f(r_ck).$$
(50)
Strictly speaking, we do not know $`f(x)`$ either; for computational purposes, in the following we shall take

$$f(x)e^{1/x}.$$
(51)
Taking $`\stackrel{}{s}=s\widehat{z}`$, we now compute $`B_{xx}(s)=B_{yy}(s)`$. From the inverse of Eq. 48 we find

$$B_{xx}(s)=b_{xx}e^{ı\stackrel{}{k}\stackrel{}{s}}\frac{d^3k}{(2\pi )^3}$$
(52)
which, using Eqs. 50 and 49 yields

$`B_{xx}(s)={\displaystyle \frac{\pi A}{4}}{\displaystyle k^\alpha f(r_ck)(\mathrm{cos}^2\theta +\mathrm{sin}^2\theta \mathrm{sin}^2\varphi )e^{ıks\mathrm{cos}\theta }𝑑k𝑑\varphi d\mathrm{cos}\theta }=`$ (53)
$`{\displaystyle \frac{\pi ^2A}{2}}{\displaystyle _0^{\mathrm{}}}k^\alpha f(r_ck)𝑑k{\displaystyle _1^1}{\displaystyle \frac{1+\mathrm{cos}^2\theta }{2}}e^{ıks\mathrm{cos}\theta }d\mathrm{cos}\theta .`$
We find

$$B_{xx}(s)=\pi ^2A_0^{\mathrm{}}k^\alpha f(r_ck)𝑑k\left(\frac{\mathrm{sin}ks}{ks}+\frac{\mathrm{cos}ks}{(ks)^2}\frac{\mathrm{sin}ks}{(ks)^3}\right).$$
(54)
When $`s0`$, because of isotropy obviously $`B_{xx}B_0^2/3`$, where $`B_0`$ is the rms field. Letting $`s0`$ in the above equation, we find

$$B_0^2=2\pi ^2A_0^{\mathrm{}}k^\alpha f(r_ck)𝑑k$$
(55)
which allows us to rewrite

$$B_{xx}(s)=\frac{1}{3}B_0^2g(\frac{s}{r_c},\alpha );g(x,\alpha )\frac{3_0^{\mathrm{}}y^\alpha f(y)𝑑y(\frac{\mathrm{sin}yx}{yx}+\frac{\mathrm{cos}yx}{(yx)^2}\frac{\mathrm{sin}yx}{(yx)^3})}{2_0^{\mathrm{}}y^\alpha f(y)𝑑y}$$
(56)
and similarly for $`B_{yy}`$.