Squaring the Circle and Cubing the Sphere: Circular and Spherical Copulas

Abstract: Do there exist circular and spherical copulas in $R^d$? That is, do there exist circularly symmetric distributions on the unit disk in $R^2$ and spherically symmetric distributions on the unit ball in $R^d$, $d\ge3$, whose one-dimensional marginal distributions are uniform? The answer is yes for $d=2$ and 3, where the circular and spherical copulas are unique and can be determined explicitly, but no for $d\ge4$. A one-parameter family of elliptical bivariate copulas is obtained from the unique circular copula in $R^2$ by oblique coordinate transformations. Copulas obtained by a non-linear transformation of a uniform distribution on the unit ball in $R^d$ are also described, and determined explicitly for $d=2$.

\section{Introduction}
\label{sec:intro}

\noindent  Do there exist spherically symmetric distributions on the closed unit ball
$B_d$ in $\mathbb{R}^d$ that have uniform one-dimensional marginal distributions on $[-1,1] $?
A distribution on $B_d$ with this property may be said to ``square the circle"
when $d=2$ and to ``cube the sphere" when $d\ge3$.

The cumulative distribution function (cdf) of a multivariate distribution
on the unit cube $[0,1]^d$ whose marginal distributions are uniform[0,1] is
commonly called a {\it copula}; see
\citet{MR2197664}
for an accessible introduction to this topic. However, although it is
customary to confine attention to distributions on the {\it unit} cube,
our interest is in \emph{spherically symmetric (= orthogonally invariant)
distributions on $B_d$ with uniform marginal distributions.}
Therefore we take ``copula" to mean a multivariate cdf on the
{\it centered} cube $C_d:=[-1,1]^d$ with uniform[-1,1] marginals.

For $d=2$ (resp., $d\ge3$), such a copula, if it exists, will be called a \emph{circular copula (resp., spherical copula)} if it is the cdf of a circularly symmetric (resp., spherically symmetric) distribution on the unit disk $B_2$ (resp., unit ball $B_d$).

It will be noted in Sections~\ref{sec:uniqueexist} and ~\ref{sec:circularcopula}
that circular and spherical copulas are unique if they exist, but exist only for
dimensions $d=2$ and $d=3$. The proof of non-existence for $d\ge4$ is remarkably simple.
Explicit expressions for these copulas are given in Sections \ref{sec:circularcopula}
and \ref{sec:sphericalcopula} respectively.

In Section \ref{sec:ellipticalcopula}, a new one-parameter family of bivariate copulas
called {\sl elliptical copulas} is
obtained from the unique circular copula in $\mathbb{R}^2$ by oblique coordinate transformations.
Finally, in Section \ref{sec:nonlinearcopula}, copulas obtained by a non-linear
transformation of a uniform distribution on the unit ball in $\mathbb{R}^d$ are described,
and determined explicitly for $d=2$.

\section{Uniqueness and Existence of Circular and Spherical Copulas}
\label{sec:uniqueexist}

\begin{prop}\label{uniqueness}
Circular and spherical copulas are unique if they exist.
\footnote{This result is well-known (e.g.,
\citet{MR0270403},
pp.31-33,
who uses ``random direction" to indicate the uniform distribution of
$U\in\partial B_3$), and reappears frequently (e.g.
\citet{MR1857936},
Theorem 3.1).
The essence of the result goes back at least to
\citet{MR1503439}.}
\end{prop}
\begin{proof}
If a circular or spherical copula exists on $C_d$, it is the cdf of a random
vector $Z\equiv(Z_1,\dots,Z_d)$  with a spherically symmetric distribution
on $B_d$ and with each $Z_i\sim\mathrm{uniform}[-1,1]$.
The latter implies that $Z$ has no atom at the origin, i.e., $P[Z=0]=0$,
so we may consider the  ``polar coordinates" representation $Z = R\cdot U$,
where $R=\|Z\|\le1$ and $U=Z/\|Z\|$. It is well known (e.g.,
\citet{MR629795},
Lemmas 1 and 2) that the random unit vector $U\equiv(U_1,\dots,U_d)$ is independent
of $R$ and is uniformly distributed on the unit sphere $\partial B_d$,
which implies that each $U_i^2\sim\mathrm{Beta}(1/2,\,(d-1)/2)$.
Since $Z_i=RU_i$, we have that
\begin{equation}\label{logZRU}
\log(Z_i^2)=\log(R^2)+\log(U_i^2).
\end{equation}
Because $R$ and $U_i$ are independent, it follows that the characteristic function
of $\log(R^2)$ is the quotient of the characteristic functions of $\log(Z_i^2)$ and $\log(U_1^2)$.
Thus the distribution of $\log(R^2)$, and therefore that of $R$,  is uniquely
determined by the distributions of $Z_i^2$ and $U_i^2$, which are already specified above.
Thus the the joint distribution of $(R,U)$ is uniquely determined,
hence so is the distribution of $Z$, hence so its cdf = copula.
\end{proof}

The existence of spherical copulas is easy to determine in three or more dimensions:

\begin{prop}\label{threeormoredim}
Spherical copulas do not exist for $d\ge4$. For $d=3$, the unique spherical copula
is generated by the uniform distribution on the unit sphere
$\partial B_3:=\{(x_1,x_2,x_3)\mid x_1^2+x_2^2+x_3^2=1\}$.
\end{prop}
\begin{proof} Let  $Z$ be as in the proof of Proposition \ref{uniqueness}. Then
\begin{equation}\label{EasyMomentInequality}
\frac{1}{3} = E(Z_i^2) = E (R^2 ) E(U_i^2) \le  \frac{1}{d}.
\end{equation}
since $Z_i\sim\mathrm{uniform}[-1,1]$, $0\le R\le1$,
and $U_i^2\sim\mathrm{Beta}(1/2,\,(d-1)/2)$. Thus $d\le3$,
so a spherical copula cannot exist when $d\ge4$.

Furthermore, if a spherical copula is to exist for $d=3$, it follows from
\eqref{EasyMomentInequality} that its generating random vector $Z\in B_3$
must satisfy $E(R^2)=1$, hence $R=1$ with probability one.
This can occur only if $Z$ is uniformly distributed on the unit sphere
$\partial B_3$. But it is well known\footnote{This follows from the fact that the area of a spherical
zone is proportion to its altitude -- cf.
\citet{MR0270403},
 Proposition (i), p. 30. }
that this distribution does indeed have uniform marginal distributions
on $[-1,1]$, hence generates  the unique spherical copula for $d=3$.
\end{proof}

\section{The Bivariate Case: the Unique Circular Copula}
\label{sec:circularcopula}

The following three questions constitute an engaging classroom exercise.
\medskip

\noindent {\bf Question 1.} Let $(X,Y)$ be a random vector uniformly distributed on the unit
disk (= ball) $B_2$ in $\mathbb{R}^2$. Find the marginal probability distributions of $X$ and $Y$.
\smallskip

\noindent {\bf Answer 1.} One can easily show that $X$ has
the ``semi-circular" probability density function (pdf) given
by
\begin{equation}
\label{pdf1}
f(x)=\frac{2}{\pi}\sqrt{1-x^2},\quad -1\le x\le1.
\end{equation}
(See Figure~\ref{fig3.1b}.)
By symmetry, $Y$ has the same pdf as $X$.
\medskip

\noindent  {\bf Question 2.} Let $(X,Y)$ be a random vector uniformly distributed on the unit
circle $\partial B_2$ in $\mathbb{R}^2$. Find the marginal probability distributions of $X$ and $Y$.
\smallskip

\noindent {\bf Answer 2.} We can represent $(X,Y)$ as $(\cos\Theta,\,\sin\Theta)$ where
 $\Theta\sim{\rm uniform}[0,2\pi)$. It follows readily that $X$ has pdf
\begin{equation}\label{pdf2}
f(x)=\frac{1}{\pi\sqrt{1-x^2}},\quad -1<x<1.
\end{equation}
(See Figure~\ref{fig3.1b}.)  By symmetry, $Y$ has the same pdf as $X$.

\begin{figure}[h]
\centering
\includegraphics[scale=1]{TwoElementaryDensities.pdf}
\caption{The densities (3.1) (lower, blue) and (3.2) (upper, purple).}
\label{fig3.1b}
\end{figure}

\smallskip

In both cases, the joint distribution of $(X,Y)$ is circularly symmetric, that
is, invariant under all orthogonal transformations of $\mathbb{R}^2$. A comparison of the shapes of
the pdfs in Figure~\ref{fig3.1b} suggest a third question:
\medskip

\noindent  {\bf Question 3.} Does a circularly symmetric bivariate
distribution with uniform$[-1,1]$  marginals exist on $B_2$? If so, it
determines a circular copula on $C_2$, which is unique by Proposition \ref{uniqueness}.
This also follows from uniqueness results for the Abel transform; see e.g. \citet{MR924577}.
\medskip
\smallskip

\noindent {\bf Answer 3.}
Optimistically, let's seek an absolutely continuous solution.
That is, we seek a bivariate pdf on $B_2$ of the form
\begin{equation}
f(x,y)=g(x^2+y^2)
\nonumber
\end{equation}
such that the marginal pdf
\begin{equation}
f(x)\equiv\int_{-\sqrt{1-x^2}}^{\sqrt{1-x^2}}f(x,y)dy,\qquad -1<x<1,
\nonumber
\end{equation}
is constant in $x$. Here $g$ is a nonnegative function on $(0,1)$ that must satisfy
\begin{equation}\label{gcondition}
 2\pi\int_0^1 rg(r^2)dr=1
\end{equation}
in order that $\int\!\!\int_{B_2} f(x,y)dxdy=1$ (transform to polar coordinates:
$(x,y)\to(r,\theta)$).

To determine a suitable $g$, first set $h(t)=g(1-t)$, then let $u=\frac{y}{\sqrt{1-x^2}}$ to
obtain
\begin{eqnarray*}
f(x)&=&\int_{-\sqrt{1-x^2}}^{\sqrt{1-x^2}}h(1-x^2-y^2)dy\label{marginalx2}\\
 &=&\sqrt{1-x^2}\int_{-1}^{1}h((1-u^2)(1-x^2))du\label{marginalx3}\\
  &=&2\sqrt{1-x^2}\int_{0}^{1}h((1-u^2)(1-x^2))du\label{marginalx4}. \label{hRepresentOfMarginalDensity}
\end{eqnarray*}
If we take $h(t)=c\, t^{-1/2}$ then clearly $f(x)$ does not depend on $x$,
and choosing $c=1/2\pi$ satisfies (\ref{gcondition}). Thus the
bivariate pdf
\begin{equation}\label{fcirc}
f(x,y)=\frac{1}{2\pi\sqrt{1-x^2-y^2}},\qquad x^2+y^2<1,
\end{equation}
determines a circularly symmetric bivariate distribution on $B_2$ and yields the desired circular copula.
\medskip

\begin{figure}[h]
\centering
\includegraphics[scale=1]{Density-2Dsymm.pdf}
\caption{Circularly symmetric bivariate density on $B_2$}
\label{fig3.1}
\end{figure}

\noindent{\bf Question 4.} Having determined the unique circularly
symmetric distribution \eqref{fcirc} on $B_2$ with uniform marginals,
what is the corresponding cdf $F(x,y)$, that is, what is the corresponding circular copula?
\smallskip

\noindent{\bf Answer 4.}
The circular symmetry of $(X,Y)$ implies that its distribution
is invariant under sign changes, i.e., $(X,Y)\buildrel d\over=(\pm X,\,\pm Y)$.
By  the following lemma, the cdf $F(x,y)\equiv P[X\le x,\,Y\le y]$ on $C_2\equiv[-1,1]^2$
can be expressed in terms of  $F_0(x,y)$, its truncation to the first quadrant:
\begin{equation}\label{F0}
F_0(x,y)\equiv P[0\le X\le x,\,0\le Y\le y]
\end{equation}
for $0\le x,y\le1$, and also in terms of the complementary cdf
$\bar F(x,y)\equiv P[X>x,\,Y>y]$ for $0\le x,y\le1$.
Because $(X,Y)\buildrel d\over=(\pm X,\,\pm Y)$ and  has uniform$[-1,1]$ marginals,
\begin{eqnarray}
F_0(x,y)&=&P[0\le X\le1,\,0\le Y\le1]-P[X>x,\,0\le Y\le 1]\nonumber\\
 &&-P[0\le X\le1,\,Y>y]+P[X>x,\,Y>y]\nonumber\\
 &=&\frac{1}{4}-\left(\frac{1-x}{4}\right)-\left(\frac{1-y}{4}\right)+\bar F(x,y)\nonumber\\
 &=&\frac{x+y-1}{4}+\bar F(x,y),\quad 0\le x,y\le 1.\label{F0Fbar}
\end{eqnarray}

\begin{lem}\label{absvalues}
Let $(X,Y)$ be a bivariate random vector on $C_2$ with uniform$[-1,1]$
marginal distributions and sign-change invariance, i.e., $(X,Y)\buildrel d\over=(\pm X,\,\pm Y)$. Then for $(x,y)\in C_2$,
\begin{eqnarray}
F(x,y)&=&\frac{x+y+1}{4}+\sigma(xy)\,F_0(|x|,|y|)\label{XYabsXY}\\
&=&\frac{x+y+1}{4}+\sigma(xy)\left[\frac{|x|+|y|-1}{4}+\bar F(|x|,|y|)\right],\label{XYabsXYcomp}
\end{eqnarray}
where $\sigma(w)=\mathrm{sign}(w)$ if $w\ne0$ and $\sigma(0)=0$.
\end{lem}
\begin{proof} To obtain \eqref{XYabsXY}, consider four cases:
\medskip

\noindent{\it Case 1:} $0\le x,y\le1$.
Because $(X,Y)$ is sign-change invariant and has uniform$[-1,1]$  marginals,
\begin{eqnarray*}
F(x,y)&=&P[0< X\le x,\,0< Y\le y]+P[0< X\le x,\,Y\le0]\\
 &&+P[X\le0,\,0< Y\le y]+P[X\le0,\,Y\le0]\\
&=& F_0(x,y)+\frac{x}{4}+\frac{y}{4}+\frac{1}{4}\\
 &=&\frac{x+y+1}{4}+\sigma(xy)\,F_0(|x|,|y|).
\end{eqnarray*}
\smallskip

\noindent{\it Case 2:} $-1\le x\le0\le y\le1$. Similarly,
\begin{eqnarray*}
F(x,y)&=&\ \ P[X\le0,\,0< Y\le y]-P[x< X\le 0,\,0\le Y\le y]\\
 &&+P[X\le0,\,Y\le0]-P[x< X\le 0,\,Y\le 0]\\
&=&\frac{y}{4}- F_0(-x,y)+\frac{1}{4}-\frac{(-x)}{4}\\
 &=&\frac{x+y+1}{4}+\sigma(xy)\,F_0(|x|,|y|).
\end{eqnarray*}
\smallskip

\noindent{\it Case 3:} $-1\le y\le0\le x\le1$. Similarly,
\begin{eqnarray*}
F(x,y)&=&P[0< X\le x,\,Y\le0]-P[0\le X\le x,\,y< Y\le 0]\\
 &&+P[X\le0,\,Y\le0]-P[X\le 0,\,y< Y\le 0]\\
&=&\frac{x}{4}-F_0(x,-y)+\frac{1}{4}-\frac{(-y)}{4}\\
&=&\frac{x+y+1}{4}+\sigma(xy)\,F_0(|x|,|y|).
\end{eqnarray*}
\smallskip

\noindent{\it Case 4:} $-1\le x,y\le0$. Similarly,
\begin{eqnarray*}
F(x,y&)=&P[X\le 0,\,Y\le 0]-P[x<X\le0,\,Y\le 0]\\
 &&-P[X\le 0,\,y<Y\le0]+P[x<X\le0,\,y<Y\le0]\\
&=&\frac{1}{4}-\frac{(-x)}{4}-\frac{(-y)}{4}+ F_0(-x,-y)\\
&=&\frac{x+y+1}{4}+\sigma(xy)\,F_0(|x|,|y|).
\end{eqnarray*}
\medskip

Finally, \eqref{XYabsXYcomp} follows from \eqref{XYabsXY} by \eqref{F0Fbar}.
\end{proof}
\medskip

Thus, to determine the circular copula $F(x,y)$ for the pdf \eqref{fcirc},
it suffices to determine the complementary cdf $\bar F(x,y)$ for $0\le x,y\le1$
and apply \eqref{XYabsXYcomp}. Because $\bar F(x,y)=0$ when  $x^2+y^2\ge1$,
we need only consider the case where $x^2+y^2<1$.
\medskip

\noindent{\it First approach:}
When $0\le x,y\le 1$ and $x^2+y^2<1$, $\bar F(x,y)$ can be expressed as follows.
By using Figure~\ref{fig3.2} we find that

\begin{figure}[h]
\centering
\includegraphics[scale=.9]{Integr-fig1-new.pdf}
\caption{Region of integration, 2-dimensional case}
\label{fig3.2}
\end{figure}

\begin{eqnarray}
\bar F(x,y)
&=&\frac{1}{2\pi}\int_{x}^{\sqrt{1-y^2}} \left \{ \int_{y}^{\sqrt{1-s^2}}\frac{1}{\sqrt{1-s^2-t^2}}dt\right \}ds \nonumber\\
  &=&\frac{1}{2\pi}\int_{x}^{\sqrt{1-y^2}} \left \{ \int_{y}^{\sqrt{1-s^2}}\frac{1}{\sqrt{(1-\frac{t^2}{1-s^2})}}\frac{dt}{\sqrt{1-s^2}}
                         \right \} ds \nonumber\\
  &=&\frac{1}{2\pi}\int_x^{\sqrt{1-y^2}} \left \{ \int_{\frac{y}{\sqrt{1-s^2}}}^1 \frac{dv}{\sqrt{1-v^2}}\right \} ds \nonumber\\
   &=&\frac{1}{2\pi}\int_x^{\sqrt{1-y^2}}\left[\frac{\pi}{2}-\arcsin\left(\frac{y}{\sqrt{1-s^2}}\right)\right]ds.\label{Fbarxysin}
\end{eqnarray}
However, we were unable to evaluate this integral directly.
\medskip

\noindent{\it Second approach:}
Fortunately, we have found a solution in the molecular biology and optics literatures,
where the problem of finding the area of the intersection of two spherical caps on the
unit sphere $\partial B_3$ has been addressed.
The following general result is due to \citet{TovchiVakser:01} and also appears in
\citet{DBLP:conf/si3d/2007}.

\begin{lem}\label{surfacearea}
Let $S_1$ and $S_2$ be spherical caps on $\partial B_3$. Let $r_1$ and $r_2$
denote their angular radii and let $d$ denote the angular distance between their
centers ($0<d\le\pi$). Assume that $0< r_1, r_2\le \pi/2$ and $d\le  r_1+ r_2$,
so that the intersection $S_1\cap S_2\ne\emptyset$ and consists of a single ``diangle";
(see Figures~\ref{fig3.4} and ~\ref{fig3.5}.)
Then Area$(S_1\cap S_2)$ is given by
\begin{eqnarray}
A(r_1,r_2;d)
& = & 2 \pi - 2\pi \cos(r_1) - 2\pi \cos (r_2)
           - 2 \arccos \left (\frac{\cos(d)-\cos(r_1)\cos(r_2)}{\sin(r_1) \sin (r_2)} \right )\nonumber \\
&& + \ 2 \cos(r_1) \arccos \left ( \frac{\cos(d)\cos(r_1)-\cos(r_2)}{\sin(d)\sin(r_1)} \right )\nonumber\\
&&  + \ 2 \cos(r_2) \arccos \left ( \frac{\cos(d)\cos(r_2)-\cos(r_1)}{\sin(d)\sin(r_2)} \right ).\label{areacapintersect}
\end{eqnarray}
\end{lem}

This result can be applied to obtain our desired circular copula as follows.
\medskip

If $(X,Y,Z)$ is uniformly distributed on $\partial B_3$, then the event $\{X>x,\,Y> y\}$
corresponds to the intersection of the two spherical caps
$\{X>x\}$ and $\{Y>y\}$, so $P[X>x,\,Y> y]$ is given by the area $A(x,y)$ of this
intersection divided by the total area of $\partial B_3$, i.e., by $4\pi$.
(See Figures~\ref{fig3.4} and ~\ref{fig3.5}.)
Also, the joint distribution of $(X,Y)$ is circularly symmetric on the unit disk $B_2$
and has uniform marginals, so must be the unique such bivariate distribution,
namely the distribution with pdf \eqref{fcirc}.

\begin{figure}
\centering
\includegraphics[scale=1]{Intersection2sphericalCaps-new.pdf}
\caption{Intersection of two spherical caps}
\label{fig3.4}
\end{figure}

\begin{figure}
\centering
\includegraphics[scale=1]{CircularView2Caps-new4.pdf}
\caption{Intersection of two spherical caps,
circular representation (modified from \citet{TovchiVakser:01})}
\label{fig3.5}
\end{figure}

Thus, for $0\le x,y\le 1$ and $x^2+y^2<1$, our desired complementary cdf is given by
\begin{eqnarray}
\bar F(x,y)&=&\frac{1}{4\pi}A(x,y)\label{FbarAxy}\\
&=& \frac{1}{4\pi}A(\arccos(x),\,\arccos(y);\,\pi/2)\label{AA}\\
& = & \frac{1}{2} - \frac{x}{2} - \frac{y}{2}
          - \frac{1}{2\pi} \arccos \left ( - \frac{xy}{ \sqrt{(1-x^2)(1-y^2)}} \right )\nonumber\\
&&  + \ \frac{x}{2\pi} \arccos \left ( \frac{- y}{\sqrt{1-x^2}} \right )
       + \frac{y}{2\pi} \arccos \left ( \frac{-x}{\sqrt{1-y^2}} \right )\nonumber\\
 &\equiv& \frac{1-x-y}{4}+\alpha(x,y),\label{psidef}
\end{eqnarray}
where for $0\le x,y\le1$  and $x^2+y^2<1$,
\begin{eqnarray}
\alpha(x,y)
&=& \frac{1}{2\pi}\left[x\arcsin \left ( \frac{y}{\sqrt{1-x^2}} \right )
           +y\arcsin\left ( \frac{x}{\sqrt{1-y^2}} \right )\right.\nonumber\\
&& \left.\qquad-\arcsin\left ( \frac{x y}{ \sqrt{(1-x^2)(1-y^2)}} \right )\right].\label{psi}
\end{eqnarray}
\begin{thm}\label{circularcopula}
The unique circular copula on $C_2$ is given by
\begin{equation}\label{circcopformula}
F(x,y)=\frac{x+y+1}{4}+\alpha(x,y),
\end{equation}
where $\alpha(x,y)$ is defined by \eqref{psi} for $x^2+y^2<1$ and by
\begin{equation}\label{psiextension}
\alpha(x,y)=\sigma(xy)\cdot\left(\frac{|x|+|y|-1}{4}\right)
\end{equation}
 for $x^2+y^2\ge1$. Note that \eqref{psi} and \eqref{psiextension}
 agree when $x^2+y^2=1$ and both are sign-change equivariant on
 $C_2$: for all $(x,y)\in C_2$ and all $\epsilon,\delta=\pm1$,
\begin{equation}\label{psiantisymmetry}
\alpha(\epsilon x,\delta y)=\epsilon \delta\cdot\alpha(x,y).
\end{equation}
\end{thm}
\begin{proof}
From \eqref{XYabsXYcomp} and \eqref{psidef}, when $x^2+y^2<1$ we have
\begin{eqnarray}
F(x,y)&=&\frac{x+y+1}{4}+\sigma(xy)\,\alpha(|x|,|y|)\\
 &=&\frac{x+y+1}{4}+\alpha(x,y)
\end{eqnarray}
by \eqref{psiantisymmetry}. When $x^2+y^2\ge1$, $\bar F(|x|,|y|)=0$ so
\eqref{circcopformula} again holds by \eqref{XYabsXYcomp} and \eqref{psiextension}.
\end{proof}

See Figure~\ref{fig3.3} for
a plot of the resulting copula (on $[-1,1]^2$).

\begin{figure}[h]
\centering
\includegraphics[scale=1]{CDF-circularCopula-2D-try2.pdf}
\caption{Copula (joint distribution function), Theorem 3.3}
\label{fig3.3}
\end{figure}

\medskip

\section{The Trivariate Case: the Unique Spherical Copula}
\label{sec:sphericalcopula}

\noindent{\bf Question 5.}
Having determined the unique spherically symmetric distribution on $B_3$
with uniform marginals, namely, the uniform distribution on the unit sphere $\partial B_3$,
what is the corresponding cdf $F(x,y,z)$ on $C_3$, i.e,  the unique spherical copula?
\medskip

\noindent{\bf Answer 5.}
As in Section \ref{sec:circularcopula},
let $(X,Y,Z)$ be uniformly distributed on $\partial B_3$, so that $F(x,y,z)= P[X\le x,\,Y\le y,\,Z\le z]$.
Again we first determine the complementary cdf
$\bar F(x,y,z)\equiv P[X>x,\,Y>y,\,Z>z]$ for $0\le x,y,z\le1$ and $x^2+y^2+z^2<1$,
the intersection of the first octant of $C_3$ with the interior of $B_3$.
Here the event $\{X>x,\,Y>y,\,Z>z\}$ corresponds to the intersection of the three
spherical caps $\{X>x\}$, $\{Y>y\}$, and $\{Z>z\}$ on $\partial B_3$,
so $\bar F(x,y,z)$ is the area $A(x,y,z)$ of this intersection divided by the total area $4\pi$ of $\partial B_3$.
\smallskip

Recall that two approaches were proposed in Section \ref{sec:circularcopula}
to obtain the area $A(x,y)$ of the intersection of {\it two} circular caps $\{X>x\}$ and $\{Y>y\}$.
The first approach led to the integral \eqref{Fbarxysin} that we were unable to evaluate explicitly,
so we adopted a second approach based on the geometric Lemma \ref{surfacearea} of \citet{TovchiVakser:01}.
Andrey Tovchigrechko has kindly suggested a method for extending Lemma \ref{surfacearea}
to the case of three spherical caps in general position, which if carried out would yield an explicit expression for $A(x,y,z)$.
However, we have found that because the axes of our three caps are mutually orthogonal,
the two approaches just mentioned for the bivariate case can be combined
to obtain $\bar F(x,y,z)\equiv\frac{1}{4\pi}A(x,y,z)$ directly for the trivariate case, as now described.
\smallskip

We begin by extending \eqref{Fbarxysin} to obtain an integral expression for
$\bar F(x,y,z)$ when $0\le x,y,z\le1$ and $x^2+y^2+z^2<1$. We require the fact that
\begin{equation}\label{arcsinab}
0 \le a,b \le 1\ \  \mathrm{and}\ \ a^2+b^2=1  \ \ \mbox{implies} \ \  \arcsin(a) + \arcsin(b) = \pi/2.
\end{equation}

\begin{lem} If $0\le x,y,z\le1$ and $x^2+y^2+z^2<1$, then
\begin{equation}\label{integralFbarxyz}
\bar F(x,y,z)
   =\frac{1}{4\pi}\int_x^{\sqrt{1-y^2-z^2}}\left[\frac{\pi}{2}
     -\arcsin\left(\frac{y}{\sqrt{1-s^2}}\right)-\arcsin\left(\frac{z}{\sqrt{1-s^2}}\right)\right]ds.
\end{equation}
Because $(X,Y,Z)$ is exchangeable, \eqref{integralFbarxyz} remains valid
under any permutation of $x,y,z$ on the right-hand side.
\end{lem}
\begin{proof}
Since $X^2+Y^2+Z^2=1$ and $(X,Y,Z)\buildrel d\over=(X,Y,-Z)$,
it follows from \eqref{fcirc} by using Figure~\ref{fig4.1} that

\begin{figure}[h]
\centering
\includegraphics[scale=1.2]{Integr-fig2-new.pdf}
\caption{Region of integration, 3-dimensional case, Lemma 4.1}
\label{fig4.1}
\end{figure}

\begin{eqnarray*}
\lefteqn{P[X>x,\,Y>y,\,Z>z]}\\
&=&\frac{1}{2}P[X>x,\,Y>y,\,X^2+Y^2<1-z^2]\\
&=&\frac{1}{4\pi}\int_{x}^{\sqrt{1-y^2-z^2}} \left \{ \int_{y}^{\sqrt{1-s^2-z^2}}\frac{1}{\sqrt{1-s^2-t^2}}dt\right \} ds \nonumber\\
  &=&\frac{1}{4\pi}\int_{x}^{\sqrt{1-y^2-z^2}}
            \left \{ \int_{y}^{\sqrt{1-s^2-z^2}}\frac{1}{\sqrt{(1-\frac{t^2}{1-s^2})}}\frac{dt}{\sqrt{1-s^2}}\right \} ds \nonumber\\
  &=&\frac{1}{4\pi}\int_x^{\sqrt{1-y^2-z^2}}
           \left \{ \int_{\frac{y}{\sqrt{1-s^2}}}^{\sqrt{\frac{1-s^2-z^2}{1-s^2}}} \frac{dv}{\sqrt{1-v^2}}\right \} ds \nonumber\\
   &=&\frac{1}{4\pi}\int_x^{\sqrt{1-y^2-z^2}}
           \left[\arcsin\left(\sqrt{\frac{1-s^2-z^2}{1-s^2}}\right)-\arcsin\left(\frac{y}{\sqrt{1-s^2}}\right)\right]ds.
   \end{eqnarray*}
Now apply \eqref{arcsinab} to obtain \eqref{integralFbarxyz}.
\end{proof}

As noted above, the integral in \eqref{integralFbarxyz} appears difficult to evaluate explicitly,
but the following indirect argument succeeds.
Recall from \eqref{Fbarxysin} and \eqref{psidef} that when $0\le x,y\le 1$ and $x^2+y^2<1$,
\begin{eqnarray*}
\bar{F} (x,y) & = & \frac{1}{2\pi} \int_x^{\sqrt{1-y^2}} \left[ \frac{\pi}{2} - \arcsin\left ( \frac{y}{\sqrt{1-s^2}} \right ) \right] ds \\
& = & \frac{1-x-y}{4}  + \alpha (x,y),
\end{eqnarray*}
where $\alpha(x,y)$ is given by \eqref{psi}. Because $z\le\sqrt{1-x^2-y^2}\le\sqrt{1-y^2}$ when $0\le x,y,z\le1$ and $x^2+y^2+z^2<1$, it follows that
\begin{eqnarray}
\lefteqn{ \frac{1}{2\pi} \int_z^{\sqrt{1-x^2 - y^2}}  \left[ \frac{\pi}{2} - \arcsin\left ( \frac{y}{\sqrt{1-s^2}} \right ) \right]ds}\nonumber \\
& = & \frac{\sqrt{1-x^2 - y^2}-z}{4} + \alpha (z,y)
  - \alpha (\sqrt{1-x^2 - y^2},\,y).\label{badifference}
\end{eqnarray}
Therefore from \eqref{integralFbarxyz} and \eqref{badifference}, if $0\le x,y,z\le1$ and $x^2+y^2+z^2<1$ then
\begin{eqnarray*}
\lefteqn{4\pi \bar F(x,y,z)} \\
& = &  \int_z^{\sqrt{1-x^2 - y^2}} \left[ \frac{\pi}{2} - \arcsin \left ( \frac{x}{\sqrt{1-s^2}} \right ) \right] ds \\
&&+  \int_z^{\sqrt{1-x^2 - y^2}} \left[ \frac{\pi}{2} - \arcsin \left ( \frac{y}{\sqrt{1-s^2}} \right )  \right] ds
 - \frac{\pi}{2} \left[ \sqrt{1-x^2-y^2} -z \right] \\
& = & \frac{\pi}{2} ( \sqrt{1-x^2-y^2} -z) +2\pi \left[  \alpha (z,x)-\alpha( \sqrt{1-x^2-y^2},\,x)  \right] \\
&&+ \frac{\pi}{2} ( \sqrt{1-x^2-y^2} -z) +2\pi \left[  \alpha (z,y)-\alpha( \sqrt{1-x^2-y^2},\,y)  \right]  \\
&& - \frac{\pi}{2} ( \sqrt{1-x^2-y^2} -z)   \\
& = & \frac{\pi}{2} (\sqrt{1-x^2-y^2} -z) \\
 &&+ 2\pi \left[ \alpha (z,x) + \alpha (z,y) - \alpha (\sqrt{1-x^2-y^2},\, x)   - \alpha (\sqrt{1-x^2-y^2},\, y) \right]\\
& = & \frac{\pi}{2} ( \sqrt{1-x^2-y^2} - z)  \\
&&+\,x \arcsin \left ( \frac{z}{\sqrt{1-x^2}} \right )
           + z \arcsin \left ( \frac{x}{\sqrt{1-z^2}} \right )
           - \arcsin \left ( \frac{xz}{ \sqrt{(1-x^2)(1-z^2)}} \right ) \\
&& +\,y \arcsin \left ( \frac{z}{\sqrt{1-y^2}} \right )
           + z \arcsin \left ( \frac{y}{\sqrt{1-z^2}} \right )
           - \arcsin \left ( \frac{yz}{ \sqrt{(1-y^2)(1-z^2)}} \right ) \\
&& -\,x \arcsin \left ( \frac{\sqrt{1-x^2-y^2}}{\sqrt{1-x^2}} \right )
          - \sqrt{1-x^2-y^2} \arcsin \left ( \frac{x}{\sqrt{x^2+y^2}} \right ) \\
 && + \arcsin \left ( \frac{x\sqrt{1-x^2-y^2}}{ \sqrt{(1-x^2)(x^2+y^2)}} \right )
  -y \arcsin \left ( \frac{\sqrt{1-x^2-y^2}}{\sqrt{1-y^2}} \right ) \\
  &&  - \sqrt{1-x^2-y^2} \arcsin \left ( \frac{y}{\sqrt{x^2+y^2}} \right )
    +\arcsin \left ( \frac{y\sqrt{1-x^2-y^2}}{ \sqrt{(1-y^2)(x^2+y^2)}} \right ) .
\end{eqnarray*}
By \eqref{arcsinab}, however,
\begin{eqnarray*}
\sqrt{1-x^2-y^2} \arcsin \left ( \frac{x}{\sqrt{x^2+y^2}} \right ) &+& \sqrt{1-x^2-y^2} \arcsin \left ( \frac{y}{\sqrt{x^2+y^2}} \right ) \\
&=&\sqrt{1-x^2-y^2} \left(\frac{\pi}{2}\right),
\end{eqnarray*}
so if we define $h(x,y)$ by
\begin{eqnarray}
h(x,y)
&=&\arcsin\left ( \frac{x y}{ \sqrt{(1-x^2)(1-y^2)}} \right )
          + \arcsin \left ( \frac{x\sqrt{1-x^2-y^2}}{\sqrt{1-x^2}\sqrt{x^2+y^2}} \right )
           \nonumber \\
 && + \arcsin \left ( \frac{y\sqrt{1-x^2-y^2}}{\sqrt{1-y^2}\sqrt{x^2+y^2}} \right )  \label{hxydef}
\end{eqnarray}
for $0\le x,y\le1$ and $x^2+y^2<1$, then
\begin{eqnarray*}
\lefteqn{4 \pi \bar F(x,y,z)}\\
& = & -\frac{\pi}{2}  z + h(x,y)  - \arcsin\left ( \frac{xy}{\sqrt{(1-x^2)(1-y^2)}} \right )  \\
&& +\ x \arcsin \left ( \frac{z}{\sqrt{1-x^2}} \right )
           + z \arcsin \left ( \frac{x}{\sqrt{1-z^2}} \right )
           - \arcsin \left ( \frac{xz}{ \sqrt{(1-x^2)(1-z^2)}} \right ) \\
&&+\ y \arcsin \left ( \frac{z}{\sqrt{1-y^2}} \right )
           + z \arcsin \left ( \frac{y}{\sqrt{1-z^2}} \right )
           - \arcsin \left ( \frac{yz}{ \sqrt{(1-y^2)(1-z^2)}} \right ) \\
&&-\ x \arcsin \left ( \frac{\sqrt{1-x^2-y^2}}{\sqrt{1-x^2}} \right )
      -y \arcsin \left ( \frac{\sqrt{1-x^2-y^2}}{\sqrt{1-y^2}} \right ) \\
& = &  - \frac{\pi}{2}  z + h(x,y) + \alpha (x,y) \\
&&  - \ x \arcsin \left ( \frac{y}{\sqrt{1-x^2}} \right ) - y \arcsin \left ( \frac{x}{\sqrt{1-y^2}} \right )\\
&& + \ x \arcsin \left ( \frac{z}{\sqrt{1-x^2}} \right )
           + z \arcsin \left ( \frac{x}{\sqrt{1-z^2}} \right )
           - \arcsin \left ( \frac{xz}{ \sqrt{(1-x^2)(1-z^2)}} \right ) \\
&& + \ y \arcsin \left ( \frac{z}{\sqrt{1-y^2}} \right )
           + z \arcsin \left ( \frac{y}{\sqrt{1-z^2}} \right )
           - \arcsin \left ( \frac{yz}{ \sqrt{(1-y^2)(1-z^2)}} \right ) \\
&& - \ x \arcsin \left ( \frac{\sqrt{1-x^2-y^2}}{\sqrt{1-x^2}} \right )
       - y \arcsin \left ( \frac{\sqrt{1-x^2-y^2}}{\sqrt{1-y^2}} \right),
\end{eqnarray*}
where $\alpha(x,y)$ is given by \eqref{psi}. Now \eqref{arcsinab} gives
\begin{eqnarray*}
&&x \arcsin \left ( \frac{y}{\sqrt{1-x^2}} \right ) + x \arcsin \left ( \frac{\sqrt{1-x^2-y^2}}{\sqrt{1-x^2}} \right ) = x \left(\frac{\pi}{2}\right) \\
&& y \arcsin \left ( \frac{x}{\sqrt{1-y^2}} \right ) + y \arcsin \left ( \frac{\sqrt{1-x^2-y^2}}{\sqrt{1-y^2}} \right )  =  y \left(\frac{\pi}{2}\right),
\end{eqnarray*}
so the above simplifies to
\begin{equation}\label{FBarFirstFormula}
4 \pi \bar F (x,y,z)= - \frac{\pi}{2}  (x+y+z) + h(x,y) + 2\pi\Delta(x,y,z),
\end{equation}
where
\begin{equation}\label{Delta}
\Delta(x,y,z)=\alpha (x,y) + \alpha (x,z) + \alpha (y,z),
\end{equation}
a symmetric function of $(x,y,z)$. By \eqref{arcsinab}, however,
\begin{eqnarray*}
h(x,y)
&=&\frac{\pi}{2}-\arcsin\left ( \frac{\sqrt{1-x^2-y^2}}{ \sqrt{(1-x^2)(1-y^2)}} \right )
          +\arcsin \left ( \frac{x\sqrt{1-x^2-y^2}}{\sqrt{1-x^2}\sqrt{x^2+y^2}} \right )\\
 &&  +\arcsin \left ( \frac{y\sqrt{1-x^2-y^2}}{\sqrt{1-y^2}\sqrt{x^2+y^2}} \right ) \\
 &\equiv&\frac{\pi}{2}-\gamma+\alpha+\beta,
\end{eqnarray*}
and
\begin{eqnarray*}
\lefteqn{\sin(\alpha+\beta)}\\
 &=&\sin\alpha\cos\beta+\cos\alpha\sin\beta\\
 &=& \frac{x\sqrt{1-x^2-y^2}}{\sqrt{1-x^2}\sqrt{x^2+y^2}}
         \frac{x}{\sqrt{1-y^2}\sqrt{x^2+y^2}}+\frac{y}{\sqrt{1-x^2}\sqrt{x^2+y^2}}
         \frac{y\sqrt{1-x^2-y^2}}{\sqrt{1-y^2}\sqrt{x^2+y^2}}\\
  &=&\frac{\sqrt{1-x^2-y^2}}{ \sqrt{(1-x^2)(1-y^2)}} \\
  &=&\sin\gamma,
\end{eqnarray*}
so $\alpha+\beta=\gamma$, hence $h(x,y)\equiv\frac{\pi}{2}$ identically in $(x,y)$. Therefore we conclude that
\begin{equation}\label{Fbarxyznice}
\bar F(x,y,z)=\frac{1- x-y-z}{8}  + \frac{\Delta(x,y,z)}{2}
\end{equation}
for $0\le x,y,z\le1$ and $x^2+y^2+z^2<1$.
\smallskip

We now apply \eqref{Fbarxyznice} to obtain the cdf $F(x,y,z)$ for all $(x,y,z)\in C_3$.
For this, extend the definition of $\Delta$ in \eqref{Delta} to all $(x,y,z)\in C_3$
by means of \eqref{psi} and \eqref{psiextension}.
\begin{thm}\label{uniquesphericalcopula}
The unique spherical copula $F(x,y,z)$ on $C_3$ is given as follows:
\smallskip

\noindent for $x^2+y^2+z^2<1$,
\begin{equation*}
F(x,y,z)=
\begin{cases}
\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2},&\mathrm{if}\ x^2+y^2+z^2<1;\\
\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}&\\
+\ \sigma(xyz)\left[\frac{1-|x|-|y|-|z|}{8}+\frac{\Delta(|x|,|y|,|z|)}{2}\right],&\mathrm{if}\ x^2+y^2+z^2\ge1.
\end{cases}
\end{equation*}
\end{thm}

\begin{proof}
We use repeatedly the facts that
$(X,Y,Z)\buildrel d\over=(\pm X,\pm Y, \pm Z)$ and $(X,Y,Z)$ has uniform$[-1,1]$ marginals.
First, by \eqref{circcopformula} and \eqref{psiantisymmetry}, the complementary
two-dimensional marginal cdf is given on all of $C_2$ by
\begin{equation}\label{compcdfB2}
\bar F(x,y)=F(-x,-y)=\frac{1-x-y}{4}+\alpha(x,y),\qquad (x,y)\in C_2.
\end{equation}
\noindent{\it Case 1a:} $0\le x,y,z$, $x^2+y^2+z^2<1$.
By inclusion-exclusion, \eqref{compcdfB2}, and \eqref{Fbarxyznice},
\begin{eqnarray*}
\lefteqn{F(x,y,z)}\\
 &=&1-P[X>x]-P[Y>y]-P[Z>z]+P[X>x,\,Y>y]+P[X>x,\,Z>z]\\
  &&+\ P[Y>y,\,Z>z]-P[X>x,\,Y>y,\,Z>z]\\
 &=&1-\frac{1-x}{2}-\frac{1-y}{2}-\frac{1-z}{2}+\frac{1-x-y}{4}+\alpha(x,y)+\frac{1-x-z}{4}+\alpha(x,z)\\
  &&+\ \frac{1-y-z}{4}+\alpha(y,z)-\frac{1- x-y-z}{8}  - \frac{\Delta(x,y,z)}{2}\\
   &=&\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}.
\end{eqnarray*}
\noindent{\it Case 1b:} $0\le x,y,z$, $x^2+y^2+z^2\ge1$. Here $P[X>x,\,Y>y,\,Z>z]=0$, so
\begin{eqnarray*}
\lefteqn{F(x,y,z)}\\
 &=&1-P[X>x]-P[Y>y]-P[Z>z]+P[X>x,\,Y>y]+P[X>x,\,Z>z]\\
  &&+\ P[Y>y,\,Z>z]\\
 &=&1-\frac{1-x}{2}-\frac{1-y}{2}-\frac{1-z}{2}+\frac{1-x-y}{4}+\alpha(x,y)+\frac{1-x-z}{4}+\alpha(x,z)\\
  &&+\ \frac{1-y-z}{4}+\alpha(y,z)\\
   &=&\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}+\ \sigma(xyz)\left[\frac{1-|x|-|y|-|z|}{8}+\frac{\Delta(|x|,|y|,|z|)}{2}\right].
\end{eqnarray*}
\noindent{\it Case 2a:} $x\le0\le y,z$, $x^2+y^2+z^2<1$.
By inclusion-exclusion, \eqref{compcdfB2}, and \eqref{Fbarxyznice},
\begin{eqnarray*}
\lefteqn{F(x,y,z)}\\
&=&P[X\le x]-P[X\le x,\,Y>y]-P[X\le x,\,Z>z]+P[X\le x,\,Y>y,\,Z>z]\\
&=&\frac{x+1}{2}-P[-X\le x,\,Y>y]-P[-X\le x,\,Z>z]+P[-X\le x,\,Y>y,\,Z>z]\\
&=&\frac{x+1}{2}-P[X\ge-x,\,Y>y]-P[X\ge-x,\,Z>z]+P[X\ge-x,\,Y>y,\,Z>z]\\
  &=&\frac{x+1}{2} -\frac{1+x-y}{4}-\alpha(-x,y) -\frac{1+x-z}{4}-\alpha(-x,z)\\
   &&+\ \frac{1+x-y-z}{8}  + \frac{\Delta(-x,y,z)}{2}\\
   &=&\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}.
\end{eqnarray*}
\noindent{\it Case 2b:} $x\le0\le y,z$, $x^2+y^2+z^2\ge1$.  Here $P[X\ge-x,\,Y>y,\,Z>z]=0$, so
\begin{eqnarray*}
\lefteqn{F(x,y,z)}\\
 &=&P[X\le x]-P[X\le x,\,Y>y]-P[X\le x,\,Z>z]\\
&=&\frac{x+1}{2}-P[X\ge-x,\,Y>y]-P[X\ge-x,\,Z>z]\\
  &=&\frac{x+1}{2} -\frac{1+x-y}{4}-\alpha(-x,y) -\frac{1+x-z}{4}-\alpha(-x,z)\\
   &=&\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}+\ \sigma(xyz)\left[\frac{1-|x|-|y|-|z|}{8}+\frac{\Delta(|x|,|y|,|z|)}{2}\right].
\end{eqnarray*}
\noindent{\it Case 3a:} $y\le0\le x,z$, $x^2+y^2+z^2<1$. Similar to Case 2a.
\smallskip

\noindent{\it Case 3b:} $y\le0\le x,z$, $x^2+y^2+z^2\ge1$. Similar to Case 2b.
\smallskip

\noindent{\it Case 4a:} $z\le0\le x,y$, $x^2+y^2+z^2<1$. Similar to Case 2a.
\smallskip

\noindent{\it Case 4b:} $z\le0\le x,y$, $x^2+y^2+z^2\ge1$. Similar to Case 2b.
\medskip

\noindent{\it Case 5a:} $x,y\le0\le z$, $x^2+y^2+z^2<1$.
By inclusion-exclusion, \eqref{compcdfB2}, and \eqref{Fbarxyznice},
\begin{eqnarray*}
F(x,y,z)
&=&P[X\le x,\,Y\le y]-P[X\le x,\,Y\le y,\,Z>z]\\
&=&P[X\ge -x,\,Y\ge -y]-P[X\ge -x,\,Y\ge -y,\,Z>z]\\
&=&\frac{1+x+y}{4}+\alpha(-x,-y)-\frac{1+x+y-z}{8} -\frac{\Delta(-x,-y,z)}{2}\\
   &=&\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}.
\end{eqnarray*}
\noindent{\it Case 5b:} $x,y\le0\le z$, $x^2+y^2+z^2\ge1$. Here $P[X\ge-x,\,Y\ge-y,\,Z>z]=0$, so
\begin{eqnarray*}
\lefteqn{F(x,y,z)}\\
&=&P[X\le x,\,Y\le y]\\
&=&P[X\ge -x,\,Y\ge -y]\\
&=&\frac{1+x+y}{4}+\alpha(-x,-y)\\
&=&\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}+\ \sigma(xyz)\left[\frac{1-|x|-|y|-|z|}{8}+\frac{\Delta(|x|,|y|,|z|)}{2}\right].
\end{eqnarray*}
\noindent{\it Case 6a:} $x,z\le0\le y$, $x^2+y^2+z^2<1$. Similar to Case 5a.
\smallskip

\noindent{\it Case 6b:} $x,z\le0\le y$, $x^2+y^2+z^2\ge1$. Similar to Case 5b.
\smallskip

\noindent{\it Case 7a:} $y,z\le0\le x$, $x^2+y^2+z^2<1$. Similar to Case 5a.
\smallskip

\noindent{\it Case 7b:} $y,z\le0\le x$, $x^2+y^2+z^2\ge1$. Similar to Case 5b.
\medskip

\noindent{\it Case 8a:} $x,y,z\le0\le z$, $x^2+y^2+z^2<1$.
\begin{eqnarray*}
F(x,y,z)
&=&P[X\le x,\,Y\le y,\,Z\le z]\\
&=&P[X\ge -x,\,Y\ge -y,\,Z\ge-z]\\
&=&\frac{1+x+y+z}{8} +\frac{\Delta(-x,-y,-z)}{2}\\
   &=&\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}.
\end{eqnarray*}
\noindent{\it Case 8b:} $x,y,z\le0$, $x^2+y^2+z^2\ge1$.
\begin{eqnarray*}
\lefteqn{F(x,y,z)=0}\\
&=&\frac{1+x+y+z}{8}+\frac{\Delta(x,y,z)}{2}+\ \sigma(xyz)\left[\frac{1-|x|-|y|-|z|}{8}+\frac{\Delta(|x|,|y|,|z|)}{2}\right].
\end{eqnarray*}
\end{proof}

\section{A One-parameter Family of Elliptical Copulas}
\label{sec:ellipticalcopula}

Let $(X,Y)\sim f(x,y)$ in \eqref{fcirc}, the unique circularly symmetric
distribution on the unit disk $B_2$ with uniform$[-1,1]$ marginals.
For any angle $\gamma\in(-\pi/2, \pi/2)$, consider the transformed variables
\begin{equation}\label{gammarotation}
U=X,\quad V_\gamma=X\sin\gamma + Y\cos\gamma.
\end{equation}
By the circular symmetry of $(X,Y)$, $V_\gamma\buildrel d\over= Y\sim\mathrm{uniform}[-1,1]$,
so the random vector $(U,V_\gamma)$ again generates a copula on the centered square $C_2$.
Denote the pdf and cdf of $(U,V_\gamma)$ by $f_\gamma(u,v)$ and $F_\gamma(u,v)$ respectively.
Then $\{F_\gamma\mid\gamma\in(-\pi/2, \pi/2)\}$ is a one-parameter family of
\emph{elliptical copulas}, so-called because the support of $(U,V_\gamma)$ is the ellipse
\begin{equation}\label{Egamma}
E_\gamma:=\{(u,v)\mid u^2 + v^2 - 2uv \sin\gamma\le \cos^2 \gamma \}.
\end{equation}
(Note that $E_0=B_2$.) From \eqref{gammarotation}, the
correlation coefficient of $U$ and $V_\gamma$ is given simply by
\begin{equation}\label{gammacorr}
\rho(U,V_\gamma)=\sin\gamma,
\end{equation}
so $\gamma$ indicates the degree of linear dependence between $U$ and $V_\gamma$.

\begin{prop}\label{fgammapdf} The pdf  of $(U,V_\gamma)$ is given by
\begin{equation}\label{fgam}
f_\gamma(u,v)= \frac{1}{2 \pi \sqrt{ \cos^2 \gamma - u^2 - v^2 + 2uv \sin \gamma } }1_{E_\gamma}(u,v).
\end{equation}
\end{prop}
\begin{proof} The pdf can be obtained by a standard Jacobian computation. From \eqref{gammarotation},
\begin{equation}\label{xuvyuv}
x(u,v) = u,\quad y(u,v) = \frac{v - u\sin \gamma}{\cos(\gamma)},
\end{equation}
so
\begin{eqnarray*}
\begin{array}{l l}
 \frac{\partial x(u,v)}{\partial u}  = 1, \qquad  & \frac{\partial x(u,v)}{\partial v} = 0, \\
 \frac{\partial y(u,v)}{\partial u} = -\tan \gamma,  & \frac{\partial y(u,v)}{\partial v} = \frac{1}{\cos \gamma} .
 \end{array}
\end{eqnarray*}
Thus the Jacobian of the transformation is $J=1/\cos \gamma $,  so from \eqref{fcirc} we obtain
\begin{eqnarray*}
f_\gamma (u,v) & = & f (x(u,v), y(u,v)) \cdot J \\
& = & \frac{1}{2 \pi \sqrt{1 - u^2 - (v-u\sin\gamma)^2/\cos^2 \gamma} } 1_{B_2} (x(u,v), y(u,v)) \cdot \frac{1}{\cos \gamma}\\
& = & \frac{\cos \gamma}{2 \pi \sqrt{(1-u^2)\cos^2 \gamma  - (v - u \sin \gamma)^2}} \frac{1}{\cos \gamma} \cdot
          1_{B_2} (x(u,v), y(u,v)) \\
& = & \frac{1}{2 \pi \sqrt{ \cos^2 \gamma - u^2 - v^2 +2uv \sin \gamma } }
          1_{E_\gamma}(u,v).
\end{eqnarray*}
\end{proof}

Figure~\ref{fig5.1}
shows the density $f_{\gamma}(u,v)$ with $\gamma = \pi/4$.

\begin{figure}[h]
\centering
\includegraphics[scale=1]{Density-fgamma1.pdf}
\caption{The density $f_{\gamma}(u,v)$ (with $\gamma = \pi/4$)}
\label{fig5.1}
\end{figure}

To describe the family of elliptical copulas $F_\gamma$,
we extend the definitions \eqref{psi} and \eqref{psiextension} as follows.
First, for $(u,v)\in E_\gamma$ define
\begin{eqnarray}
\alpha_\gamma(u,v)&=&
\frac{1}{2\pi}\left[u\arcsin \left ( \frac{v-u\sin\gamma}{\cos\gamma\sqrt{1-u^2}} \right )
         +v\arcsin\left ( \frac{u-v\sin\gamma}{\cos\gamma\sqrt{1-v^2}} \right )\right.\nonumber\\
&& \left.\qquad-\arcsin\left ( \frac{uv-\sin\gamma}{ \sqrt{(1-u^2)(1-v^2)}} \right )\right].\label{psigamma}
\end{eqnarray}
Note that $\alpha_\gamma$ reduces to $\alpha$ in \eqref{psi} when $\gamma=0$, i.e., when $V_\gamma = Y$. From \eqref{areacapintersect},
\begin{equation}\label{Apsigamma}
A(\arccos u,\,\arccos v;\,\frac{\pi}{2}-\gamma)=(1-u-v)\pi+4\pi\alpha_\gamma(u,v).
\end{equation}

Next, extend the definition of $\alpha_\gamma(u,v)$ to $C_2\setminus E_\gamma$ as follows (see Figure~\ref{fig5.2}):
\begin{equation}\label{psigammaextension}
\alpha_\gamma(u,v)=
\begin{cases}
\frac{u+v-1}{4}& \mathrm{if}\ (u,v)\in R_5(\gamma):=(C_2\setminus E_\gamma) \cap\{(u,v)\mid u+v>1+\sin\gamma\},\\
\frac{u-v+1}{4}& \mathrm{if}\ (u,v)\in R_6(\gamma):=(C_2\setminus E_\gamma) \cap\{(u,v)\mid v-u>1-\sin\gamma\},\\
\frac{-u+v+1}{4}& \mathrm{if}\ (u,v)\in R_7(\gamma):=(C_2\setminus E_\gamma) \cap\{(u,v)\mid v-u<\sin\gamma-1\},\\
\frac{-u-v-1}{4}& \mathrm{if}\ (u,v)\in R_8(\gamma):=(C_2\setminus E_\gamma) \cap\{(u,v)\mid u+v<-\sin\gamma-1\}.
\end{cases}
\end{equation}

Note that \eqref{psigamma} and \eqref{psigammaextension} agree on
$\partial E_\gamma$, i.e., when $u^2 + v^2 - 2uv \sin\gamma=\cos^2 \gamma$.
Also note that \eqref{psigammaextension} reduces to $\alpha$ in \eqref{psiextension}
when $\gamma=0$. The following lemma will be useful for the proof of Theorem \ref{Fgamcdf}.

\begin{lem}\label{Fgammasymm}
Let $(U,V)$ be a bivariate random vector in $C_2$ with uniform$[-1,1]$
marginals that satisfies $(U,V)\buildrel d\over=(-U,-V)$. Then the cdf $F(u,v)$ satisfies
\begin{equation}\label{FuvFminusuv}
F(u,v)=\frac{u+v}{2}+F(-u,-v),\quad (u,v)\in C_2.
\end{equation}
\end{lem}
\begin{proof} By the symmetry condition,
\begin{eqnarray*}
F(u,v)
&=&P[-U\le u,\,-V\le v]
        =P[U\ge-u,\,V\ge-v]\\
 &=&1-P[U<-u]-P[V<-v]+P[U<-u,\,V<-v]\\
  &=&1-\left(\frac{-u+1}{2}\right)+\left(\frac{-v+1}{2}\right)+P[U\le-u,\,V\le-v]\\
   &=&\frac{u+v}{2}+F(-u,\,-v).
\end{eqnarray*}
\end{proof}

\begin{figure}[h]
\centering
\includegraphics[scale=1]{Regionplots-sec5.pdf}
\caption{Eight regions $R_1 (\gamma) - R_8 (\gamma)$ for an elliptical copula (with $\gamma = \pi/8$)}
\label{fig5.2}
\end{figure}

\begin{thm}\label{Fgamcdf} The  cdf $\equiv$ copula of $(U,V_\gamma)$ is given by
\begin{equation}\label{Fgam}
F_\gamma(u,v)=\frac{u+v+1}{4}+\alpha_\gamma(u,v),\qquad\qquad(u,v)\in C_2.
\end{equation}
\end{thm}

\begin{proof}
 To find $F_\gamma(u,v)$ we again use the formula \eqref{areacapintersect} for the
 area of the intersection of two spherical caps on $\partial B_3$.
 Here, unlike \eqref{AA}, the axes of the two caps are not necessarily perpendicular.
 The single formula \eqref{Fgam} is obtained by considering the partition
 $C_2=\cup_{i=1}^8R_i(\gamma)$, where $R_5(\gamma)-R_8(\gamma)$ are defined in
 \eqref{psigammaextension} and (see Figure~\ref{fig5.2})
 \begin{eqnarray*}
 R_1(\gamma)&=&E_\gamma\cap\{(u,v)\mid 0\le u,v\le1\},\\
  R_2(\gamma)&=&E_\gamma\cap\{(u,v)\mid-1\le u\le0\le v\le1\},\\
   R_3(\gamma)&=&E_\gamma\cap\{(u,v)\mid-1\le v\le0\le u\le1\},\\
    R_4(\gamma)&=&E_\gamma\cap\{(u,v)\mid-1\le u,v\le0\}.
 \end{eqnarray*}

\noindent {\it Case 1:} $(u,v)\in R_1(\gamma)$.   By using Figure~\ref{fig5.3},

\begin{figure}[h]
\centering
\includegraphics[scale=1]{Regionplot2-sec5-new.pdf}
\caption{The region $[X>u, X\sin(\gamma) + Y\cos(\gamma) >v]$ for Case 1 (with $\gamma = \pi/8$)}
\label{fig5.3}
\end{figure}

\begin{eqnarray*}
F_\gamma(u,v)&=&P[U\le u,\,V_\gamma\le v]\\
&=&1-P[U>u]-P[V_\gamma>v]+P[U>u,\,V_\gamma>v]\\
&=&1-\left(\frac{1-u}{2}\right)-\left(\frac{1-v}{2}\right)+P[X>u,\,X\sin\gamma+Y\cos\gamma>v]\\
 &=&\frac{u+v}{2}+\frac{1}{4\pi}A(\arccos u,\arccos v;\,\pi/2-\gamma)\\
 &=& \frac{u+v+1}{4}+\alpha_\gamma(u,v)\qquad\qquad\qquad\qquad\qquad\mathrm{[by\ \eqref{Apsigamma}]}.
\end{eqnarray*}

\noindent {\it Case 2:}   $(u,v))\in R_2(\gamma)$.  Because $(X,Y)\buildrel d\over=(-X,Y)$ and using Figure~\ref{fig5.3-2}
\begin{eqnarray*}
F_\gamma(u,v)&=&P[U\le u]-P[U\le u,\,V_\gamma>v]\\
&=&\frac{u+1}{2}-P[X\le u,\,X\sin\gamma+Y\cos\gamma>v]\\
&=&\frac{u+1}{2}-P[-X\le u,\,-X\sin\gamma+Y\cos\gamma>v]\\
&=&\frac{u+1}{2}-P[X\ge -u,\,X\sin(-\gamma)+Y\cos(-\gamma)>v]\\
&=&\frac{u+1}{2}-\frac{1}{4\pi}A(\arccos(-u),\arccos v;\,\pi/2+\gamma)\\
&=& \frac{u+v+1}{4}-\alpha_{-\gamma}(-u,v)\ \ \qquad\qquad\qquad\qquad\mathrm{[by\ \eqref{Apsigamma}]}\\
&=& \frac{u+v+1}{4}+\alpha_\gamma(u,v)\qquad\qquad\qquad\qquad\qquad\mathrm{[by\ \eqref{psigamma}]}.
\end{eqnarray*}

\begin{figure}[h]
\centering
\includegraphics[scale=1]{Regionplot3-sec5-try2.pdf}
\caption{The region $[X<u, X\sin(\gamma) + Y\cos(\gamma) >v]$ for Case 2 (with $\gamma = \pi/8$)}
\label{fig5.3-2}
\end{figure}

\noindent {\it Case 3:}   $(u,v)\in R_3(\gamma)$. Then $(-u,-v)\in R_2(\gamma)$, so by Lemma \ref{Fgammasymm} and Case 2,
\begin{eqnarray*}
F_\gamma(u,v)&=&\frac{u+v}{2}+F_\gamma(-u,-v)\\
 &=&\frac{u+v}{2}+\frac{-u-v+1}{4}+\alpha_\gamma(-u,-v)\\
 &=& \frac{u+v+1}{4}+\alpha_\gamma(u,v)\qquad\qquad\qquad\qquad\mathrm{[by\ \eqref{psigamma}]}.
\end{eqnarray*}

\noindent {\it Case 4:}  $(u,v)\in R_4(\gamma)$. Then $(-u,-v)\in R_1(\gamma)$, so by Lemma \ref{Fgammasymm} and Case 1,
the argument for Case 3 applies verbatim.
\bigskip

\noindent {\it Case 5:}  $(u,v)\in R_5(\gamma)$.
\begin{eqnarray*}
F_\gamma(u,v)&=&1-P[U>u]-P[V_\gamma>v]\\
&=&1-\left(\frac{1-u}{2}\right)-\left(\frac{1-v}{2}\right)\\
&=&\frac{u+v+1}{4}+\alpha_\gamma(u,v)\qquad\qquad \qquad[\mathrm{by\ \eqref{psigammaextension}}].
\end{eqnarray*}

\noindent {\it Case 6:}  $(u,v)\in R_6(\gamma)$.
\begin{eqnarray*}
F_\gamma(u,v)&=&P[U\le u]\\
 &=&\frac{u+1}{2}\\
&=&\frac{u+v+1}{4}+\alpha_\gamma(u,v)\qquad\qquad \qquad[\mathrm{by\ \eqref{psigammaextension}}].
\end{eqnarray*}

\noindent {\it Case 7:}  $(u,v)\in R_7(\gamma)$. Then $(-u,-v)\in R_6(\gamma)$,
so by Lemma \ref{Fgammasymm} and Case 6, the argument for Case 3 applies verbatim.
\bigskip

\noindent {\it Case 8:} $(u,v)\in R_8(\gamma)$. Then $(-u,-v)\in R_5(\gamma)$,
so by Lemma \ref{Fgammasymm} and Case 5, the argument for Case 3 applies verbatim.
\end{proof}

\begin{figure}[h]
\centering
\includegraphics[scale=1]{CDF-fgamma1-new.pdf}
\caption{The copula $F_{\gamma}(u,v)$ (with $\gamma = \pi/4$)}
\label{fig5.4}
\end{figure}

\bigskip

\section{Copulas Derived from the Uniform Distribution on the Unit Ball}
\label{sec:nonlinearcopula}

Up to now we have addressed the question of whether copulas can be
generated by means of linear functions of  a circularly symmetric or
spherically symmetric random vector. Now we ask whether non-linear
functions of such random vectors can generate copulas.
We shall restrict attention to random vectors uniformly distributed
over the unit ball $B_d$ and produce relatively simple non-linear functions that generate copulas on $C_d$.

We begin with the bivariate case. Suppose that $(X,Y) $ is distributed
uniformly on the unit disk $B_2 = \{ (x,y) \in \mathbb R^2\mid \ x^2 + y^2 \le 1 \}$.   Because
\begin{eqnarray*}
X\mid Y&\sim&\mathrm{uniform}[-\sqrt{1-Y^2},\,\sqrt{1-Y^2}]\\
\mathrm{and}\qquad Y\mid X&\sim&\mathrm{uniform}[-\sqrt{1-X^2},\,\sqrt{1-X^2}],
\end{eqnarray*}
it follows that the random variables
\begin{equation}\label{UVdef}
U := \frac{X}{\sqrt{1- Y^2} },  \qquad V := \frac{Y}{\sqrt{1-X^2}}
\end{equation}
satisfy
\begin{eqnarray*}
U\mid Y&\sim&\mathrm{uniform}[-1,1]\\
 V\mid X&\sim&\mathrm{uniform}[-1,1].
\end{eqnarray*}
Thus, $U$ and $Y$ are independent, $V$ and $X$ are independent, and unconditionally,
\begin{eqnarray*}
U&\sim&\mathrm{uniform}[-1,1]\\
\qquad V&\sim&\mathrm{uniform}[-1,1],
\end{eqnarray*}
so the joint distribution of $(U,V)$ generates a copula $F_(u,v)$
on the centered cube $C_2\equiv[-1,1]^2$.
Note that $U$ and $V$ are not linear functions of $(X,Y)$.
\medskip

\noindent{\bf Question 6:}
Are $U$ and $V$ independent, and if not, what is the nature of their dependence?
\smallskip

\noindent {\bf Answer 6.} Clearly $U$ and $V$ are uncorrelated, since $E(U)=E(V)=0$ and
\begin{equation}
E(UV)=E\left(\frac{XY}{\sqrt{(1-X^2)(1-Y^2)}}\right)=0,
\nonumber
\end{equation}
all by the circular symmetry of $(X,Y)$. However, the joint pdf and cdf of $(U,V)$
derived below show that they are not independent.
\medskip

\begin{prop}\label{Prop1}
 The joint density of $(U,V)$  is given by
\begin{eqnarray}
f(u,v) = \frac{1}{\pi} \frac{\sqrt{(1-u^2)(1-v^2)}}{(1- u^2 v^2)^2} 1_{C_2} (u,v) .
\label{CopulaDensityOnMinusOneToOneSquare}
\end{eqnarray}
\end{prop}

\begin{figure}[h]
\centering
\includegraphics[scale=1]{DensityPlotSec6.pdf}
\caption{Joint density $f(u,v)$ of $(U, V)$ in Proposition 6.1}
\label{fig6.1}
\end{figure}

\begin{proof} This pdf is again obtained via the Jacobian method. It follows from \eqref{UVdef} that
\begin{eqnarray*}
&& u^2 (1-y^2) = x^2, \qquad \mbox{and} \qquad  v^2 (1-x^2) = y^2 .
\end{eqnarray*}
Substitution of the second expression for $y^2$ into the left side of the first relation and vice-versa
yields
\begin{eqnarray*}
x^2 = \frac{u^2 (1-v^2)}{1 - u^2 v^2}, \qquad y^2 = \frac{v^2(1-u^2)}{1 - u^2 v^2} ,
\end{eqnarray*}
so, since $x$ and $u$ ($y$ and $v$) have the same signs by \eqref{UVdef}, we obtain
\begin{equation}\label{xyuv}
x\equiv x(u,v)= \frac{u \sqrt{1-v^2}}{\sqrt{1 - u^2 v^2}}, \qquad  y\equiv y(u,v) = \frac{v \sqrt{1-u^2}}{\sqrt{1 - u^2 v^2}} .
\end{equation}
Thus
\begin{eqnarray*}
\frac{\partial x}{\partial u}
& = & \sqrt{1-v^2}\left[\frac{1}{\sqrt{1-u^2v^2}} + u  (1- u^2 v^2)^{-3/2} (u v^2)\right] \\
& = & \frac{\sqrt{1-v^2}}{\sqrt{1-u^2 v^2}} \left \{ 1 + \frac{u^2 v^2}{1 - u^2 v^2} \right \} \\
& = & \frac{\sqrt{1-v^2}}{(1-u^2 v^2)^{3/2}},
\end{eqnarray*}
\begin{eqnarray*}
\frac{\partial x}{\partial v}
& = & u\left[\frac{  (1-v^2)^{-1/2} (-v)}{\sqrt{1-u^2v^2}} +  \sqrt{1-v^2}(1- u^2 v^2)^{-3/2} (u^2 v) \right]\\
& = & \frac{uv}{\sqrt{1-v^2} \sqrt{1-u^2 v^2}} \left \{ -1 + \frac{u^2 (1-v^2)}{1 - u^2 v^2} \right \}  \\
& = & - \frac{uv (1-u^2)}{\sqrt{1-v^2} (1 - u^2 v^2)^{3/2} } .
\end{eqnarray*}
By symmetry it follows that the Jacobian is given by
\begin{eqnarray*}
J = \left ( \begin{array}{ c c}   \frac{\sqrt{1-v^2}}{(1-u^2 v^2)^{3/2}} & - \frac{uv (1-u^2)}{\sqrt{1-v^2} (1 - u^2 v^2)^{3/2} } \\
- \frac{uv (1-v^2)}{\sqrt{1-u^2} (1 - u^2 v^2)^{3/2} } & \frac{\sqrt{1-u^2}}{(1-u^2 v^2)^{3/2}}
\end{array} \right ) ,
\end{eqnarray*}
and hence the determinant of $J$ is given by
\begin{eqnarray*}
| J | = \frac{\sqrt{(1-u^2)(1-v^2)}}{(1 - u^2 v^2)^2} .
\end{eqnarray*}
Because the pdf of $(X,Y)$ is $f(x,y)=\frac{1}{\pi}1_{B_2}(x,y)$, the result  \eqref{CopulaDensityOnMinusOneToOneSquare} follows.
\end{proof}

For $0\le u,v\le1$, $(u,v)\ne(1,1)$, let $E_1(u)$ and $E_2(v)$ be the ellipses
\begin{eqnarray}\label{twoellipses}
E_1(u)&=&\left\{(x,y)\Bigm| \frac{x^2}{u^2}+y^2\le1 \right\},\\
E_2(v)&=&\left\{(x,y)\Bigm| x^2+\frac{y^2}{v^2}\le1 \right\}.
\end{eqnarray}
The next  lemma leads to the cdf $F(u,v)$ corresponding to the pdf \eqref{CopulaDensityOnMinusOneToOneSquare}.
\begin{lem}\label{AreaE1E2}
\begin{equation*}
 \mathrm{Area}(E_1(u)\cap E_2(v))
 =2u\arcsin\left( \frac{v \sqrt{1-u^2}}{\sqrt{1 - u^2 v^2}}\right)+2v\arcsin\left(\frac{u \sqrt{1-v^2}}{\sqrt{1 - u^2 v^2}}\right).
\end{equation*}
\end{lem}
\begin{proof}Define the points $o,a,b,c,d,d,f,g$ as follows:
see Figure~\ref{fig6.2}

\begin{figure}[h]
\centering
\includegraphics[scale=1]{Integr-fig3.pdf}
\caption{Integration regions for Lemma 6.2}
\label{fig6.2}
\end{figure}

\begin{eqnarray*}
o&=&(0,0),\\
a&=&(x(u,v),\,y(u,v))=\left(\frac{u \sqrt{1-v^2}}{\sqrt{1 - u^2 v^2}},\,\frac{v \sqrt{1-u^2}}{\sqrt{1 - u^2 v^2}}\right),\\
b&=&(u,0),\\
c&=&(0,v),\\
d&=&(\sqrt{1-y^2(u,v)},\,y(u,v)),\\
e&=&(x(u,v),\,\sqrt{1-x^2(u,v)}),\\
f&=&(1,0),\\
g&=&(0,1).
\end{eqnarray*}
Then
\begin{eqnarray*}
\frac{1}{4}\mathrm{Area}(E_1(u)\cap E_2(v))&=&\mathrm{Area}(oab)+\mathrm{Area}(oac)\\
 &=&u\,\mathrm{Area}(odf)+v\,\mathrm{Area}(oeg)\\
 &=&\frac{u}{2}\arcsin(y(u,v))+ \frac{v}{2}\arcsin(x(u,v)),
\end{eqnarray*}
from which the result follows.
\end{proof}

\begin{thm}\label{secondcopula}
The copula (= cdf)  corresponding to the pdf \eqref{CopulaDensityOnMinusOneToOneSquare} is given by
\begin{equation*}
F(u,v)=\frac{u+v+1}{4}+\frac{u}{2\pi}
    \arcsin\left( \frac{v \sqrt{1-u^2}}{\sqrt{1 - u^2 v^2}}\right)
    +\frac{v}{2\pi}\arcsin\left(\frac{u \sqrt{1-v^2}}{\sqrt{1 - u^2 v^2}}\right),\ \ (u,v)\in C_2.
\label{NonLinearCopula1}
\end{equation*}
\end{thm}

\begin{figure}[h]
\centering
\includegraphics[scale=1]{DistributionPlotSec6-new.pdf}
\caption{Nonlinear transformation copula $F(u,v)$ in Theorem 6.3}
\label{fig6.3}
\end{figure}

\begin{proof}
Because $(U,V)$ is sign-change invariant and has uniform$[-1,1]$ marginals,
it follows from \eqref{F0} and \eqref{XYabsXY} in Lemma \ref{absvalues}
and from \eqref{UVdef} that for $(u,v)\in C_2$,
\begin{eqnarray*}
F(u,v)&=&\frac{u+v+1}{4}+\sigma(uv)\,P[0\le U\le |u|,\, 0\le V\le |v|]\label{UVabsUV}\\
&=&\frac{u+v+1}{4}+\frac{\sigma(uv)}{4}\,P[U^2\le u^2,\, V^2\le v^2]\label{UVabsUVsquared}\\
&=&\frac{u+v+1}{4}+\frac{\sigma(uv)}{4}\,P[(X,Y)\in E_1(u)\cap E_2(v)]\\
&=&\frac{u+v+1}{4}+\frac{\sigma(uv)}{4\pi}\mathrm{Area}(E_1(|u|)\cap E_2(|v|))\label{Buv}.
\end{eqnarray*}
The result now follows from Lemma \ref{AreaE1E2}.
\end{proof}
\medskip

\par\noindent
{\bf Problem:}
The construction \eqref{UVdef} extends readily to generate a copula
on $C_d$. For $d=3$, for example, let $(X,Y,Z)$ be uniformly distributed on the unit ball $B_3$ and define
\begin{equation}\label{UVW}
U:=\frac{X}{\sqrt{1-Y^2-Z^2}},\quad V:=\frac{Y}{\sqrt{1-X^2-Z^2}},\quad W:=\frac{Z}{\sqrt{1-X^2-Y^2}}.
\nonumber
\end{equation}
Then the marginal distributions of $U$, $V$,
and $W$ are each uniform$[-1,1]$ so the cdf $G(u,v,w)$ is a copula on $C_3$.
To find this copula one would need to determine
\begin{equation}\label{area3}
\mathrm{Volume}(E_1(u)\cap E_2(v)\cap E_3(w)),
\nonumber
\end{equation}
where now, for $0\le u,v,w\le 1$,
$E_1(u)$, $E_2(v)$, and $E_3(w)$ are the ellipsoids
\begin{eqnarray*}
E_1(u)&=&\left\{(x,y,z)\Bigm| \frac{x^2}{u^2}+y^2+z^2\le1 \right\},\\
E_2(v)&=&\left\{(x,y,z)\Bigm| x^2+\frac{y^2}{v^2}+z^2\le1 \right\},\\
E_3(w)&=&\left\{(x,y,z)\Bigm| x^2+y^2+\frac{y^2}{w^2}\le1 \right\}.
\end{eqnarray*}

\noindent{\bf Acknowledgement:}
We gratefully acknowledge several helpful suggestions by Ilya Vakser and Andrey Tovchigrechko.

\bigskip

\bibliographystyle{ims}


Abelian varieties as automorphism groups of smooth projective varieties

Abstract: We determine which complex abelian varieties can be realized as the automorphism group of a smooth projective variety.

\section{Introduction}

In this note we determine which complex abelian varieties $A$ can be realized as the automorphism group of a complex smooth projective variety.
{Given an abelian variety $A$, we denote by $\operatorname{Aut}_0(A)$ (respectively $\operatorname{Aut}(A)$)
the automorphism group of $A$ as an algebraic group (respectively as a projective variety).}
We prove that if {$\operatorname{Aut}_0(A)$} is infinite then $A$ can never be realized
{as the automorphism group of a smooth projective variety} (Theorem \ref{thm:InfiniteOrder}), while if {$\operatorname{Aut}_0(A)$} is finite there exists a smooth projective
variety $Y$ of dimension $2+\dim A$ such that $\Aut(Y)=A$ (Theorem \ref{thm:AutoY}).

\medskip

\textbf{Acknowledgments.} This work was motivated by a more general question posed by Michel Brion in Oberwolfach; the specific case of abelian varieties, which we analyse in this paper, was raised by Corrado De Concini.
We would like to thank Angelo Vistoli for some very useful correspondence regarding the comparison between étale and analytic cohomology.

\section{Abelian varieties with infinite automorphism group}
In this section we show that {no abelian variety with infinite
$\operatorname{Aut}_0(A)$ can be realized as the automorphism group of a smooth projective variety:}
\begin{theorem}\label{thm:InfiniteOrder}
Let $A$ be an abelian variety such that $\operatorname{Aut}_0(A)$, the automorphism group of $A$ as an algebraic group, is infinite.
Let $X$ be a smooth projective variety on which $A$ acts faithfully: then the automorphism group of $X$ is strictly larger than $A$.
\end{theorem}

The proof relies on the following result, due to Brion \cite[page 2]{MR2562620}:
\begin{theorem}\label{thm:Brion}
Let $X$ be a smooth projective variety on which an abelian variety $A$ acts faithfully.
There is a positive integer $n$ and a $A[n]$-invariant closed subscheme $Y$ of $X$ such that there is an $A$-equivariant isomorphism
\[
X \cong Y \times^{A[n]} A.
\]
\end{theorem}

\begin{proof}(of Theorem \ref{thm:InfiniteOrder})
Let $\iota : A \hookrightarrow \operatorname{Aut}(X)$ be the given action of $A$ on $X$ and write $X \cong Y \times^{A[n]} A$ as in Theorem \ref{thm:Brion}.
We can represent $X \cong Y \times^{A[n]} A$ more explicitly as the quotient
\[
X \cong \frac{Y \times A}{A[n]},
\]
where $t \in A[n]$ acts on $(y,a)$ as $t \cdot (y,a) = (\iota(t)(y),a-t)$. This quotient is well-behaved, because $A[n]$ is a finite group acting on $Y \times A$ with no fixed points. In particular, in order to give an (invertible) map $X \to X$ it is enough to give an (invertible) map $Y \times A \to Y \times A$ that is compatible with the action of $A[n]$.
Notice that, since $A[n]$ is finite and stable under the action of $\Aut_0(A)$,
there exists a nontrivial automorphism $\varphi\in \Aut_0(A) $ that acts trivially on $A[n]$.
We claim that the automorphism $\psi$ of $Y \times A$ given by $(y,a) \mapsto (y,\varphi(a))$
descends to an automorphism $\overline{\psi}$ of $X$; since $\overline{\psi}$ is not in the
image of $\iota$, this proves that $\operatorname{Aut}(X)$ is strictly larger than $\iota(A)$.
To see that $\overline{\psi}$ descends to $X$, it suffices to check that for every $t \in A[n]$ we have
$\psi(t \cdot (y,a)) = t \cdot \psi((y,a))$, that is,
\[
\psi((\iota(t)y,a-t)) = t \cdot (y,\varphi(a)) \Longleftrightarrow (\iota(t) (y),\varphi(a-t)) = (\iota(t) (y),\varphi(a)-t);
\]
this last equality holds since $\varphi$ is a group homomorphism and
$t$ is in $A[n]$, which $\varphi$ fixes pointwise.
\end{proof}

\section{Abelian varieties with finite automorphism group}

We will now prove that any abelian variety such that $\operatorname{Aut}_0(A)$ is finite can be realized as the automorphism group of a smooth projective variety $Y$. We first make some remaks on the structure of abelian varieties with finite automorphism group.

\begin{lemma}\label{lemma:Autfinito}
Let $A$ and $B$ two isogenous abelian varieties then $\Aut_0(A)$ is finite if and only if $\Aut_0(B)$ is finite.
\end{lemma}

\begin{proof}
Since being isogenous is a symmetric relation, it suffices to prove that if $A \to B$ is an isogeny and $\operatorname{Aut}_0(A)$ is infinite, then so is $\operatorname{Aut}_0(B)$.
Write $B \cong A/H$, where $H$ is a finite subgroup of $A$, and assume that $\Aut_0(A)$ is infinite.
Notice that every automorphism $\varphi$ of $A$ which leaves $H$ stable induces an automorphism $\bar \varphi$ of $B$, and that
$\bar\varphi$ is trivial if and only if $\varphi$ is trivial.
Let $n$ be the order of $H$; in particular, we have $H\subset A[n]$.
Any automorphism $\varphi$ of $A$ leaves $A[n]$ stable, so, since $\Aut_0(A)$ is infinite
and $A[n]$ is finite, the subgroup of automorphisms $\varphi$ which fix $A[n]$ pointwise is infinite.
Every such automorphism leaves $H$ stable, hence it descends to an automorphism of $B$, and since the map $\varphi \mapsto \bar{\varphi}$ is injective we deduce that $\Aut_0(B)$ is infinite.
\end{proof}

\begin{lemma}\label{lemma:UniquenessSubvarieties}
Let $A$ be an abelian variety such that $\operatorname{Aut}_0(A)$ is finite. Then any two simple abelian subvarieties $A_1, A_2$ of $A$ are isogenous
if and only if they coincide. Moreover, if $A_1$ is a simple abelian subvariety of $A$, then $\Aut_0(A_1)$ is finite.
\end{lemma}
\begin{proof}
Suppose by contradiction that we can find two distinct but {isogenous} simple abelian subvarieties $A_1, A_2$ of $A$.
By Poincaré's reducibility theorem, there is an abelian subvariety $C$ of $A$ such that the multiplication map $A_1 \times A_2 \times C \to A$
is an isogeny. Let $B$ be an abelian variety such that there exists isogenies $\varphi_i:B\longrightarrow A_i$ and define the isogeny
\[
\varphi: B^2 \times C \to A \quad\text{ by } \quad \varphi(b_1,b_2,c)=\varphi_1(b_1)+\varphi_2(b_2)+c.
\]
Now notice that $\psi(b_1,b_2,c)=(b_1,b_1+b_2,c)$ defines an automorphism of $B^2\times C$ of infinite order, and by the previous Lemma we
conclude that $\Aut_0(A)$ is also infinite, contradiction.
\begin{comment}
Let $H$ be the finite kernel of $\varphi$; in order to give an automorphism of $A$, it suffices to give an automorphism of $B^2 \times C$ that fixes $H$.
Let  $\psi$ be the endomorphism of $B^2\times C$ given by
$$
\psi(b_1,b_2,c)=(b_1,b_1+b_2,c)
$$
and notice that it is of infinite order.
Since $H$ is finite, it is contained in the torsion subgroup $(B^2 \times C)[n]$ for some $n>0$.
Since $(B^2 \times C)[n]$ is finite and $\psi$ stabilizes it as a set, a sufficiently large power $\psi^k$ of $\psi$ fixes
$(B^2 \times C)[n]$ pointwise and in particular fixes $H$. This implies that $\psi^k$ descends
to an automorphism of infinite order of $A$, contradicting our hypotheses.
Observe that any matrix $M =\begin{pmatrix}
\alpha & \beta \\ \gamma & \delta
\end{pmatrix} \in \operatorname{GL}_2(\mathbb{Z})$ induces an automorphism $\psi_M$ of $B^2 \times C$ by the formula
\[
\begin{array}{cccc}
\psi_M : & B^2 \times C & \to & B^2 \times C \\
& (b_1,b_2,c) & \mapsto & (\alpha b_1+ \beta b_2, \gamma b_1 + \delta b_2, c).
\end{array}
\]
Pick any $M \in \operatorname{SL}_2(\mathbb{Z})$ of infinite order. Since $H$ is finite, it is contained in the torsion subgroup $(B^2 \times C)[n]$ for some $n>0$.
Since $(B^2 \times C)[n]$ is finite and $\psi_M$ stabilizes it as a set, a sufficiently large power $\psi_M^k$ of $\psi_M$ fixes $(B^2 \times C)[n]$ pointwise. This implies that $\psi_M^k$ (which is of infinite order) descends to an automorphism of infinite order of $A$, contradicting our hypotheses.
\end{comment}
The proof that for any simple abelian subvariety $A_1$ of $A$ the group $\Aut_0(A_1)$ is finite is completely analogous.
\end{proof}

From now on we fix an abelian variety $A$ with finite automorphism group $\Aut_0(A)$. By the previous Lemma and Poincar\'e reducibility Theorem we know that there
exist uniquely
determined simple abelian subvarieties $A_1,\dots,A_h$ of $A$ such that the sum
$$
\sigma: A_1\times\cdots \times A_h\longrightarrow A  \qquad \sigma(a_1,\dots,a_h)=a_1+\dots+a_h
$$
is an isogeny. We denote by $\Sigma$ the finite kernel of this map and denote by $N$ its order. By Lemma \ref{lemma:UniquenessSubvarieties}, $A_i$ and $A_j$ are not isogenous
if $i\neq j$ and $\Aut_0(A_i)$ is finite for all $i$. Finally, notice that any abelian variety constructed in this way has finite automorphism group.

\subsection{Construction of the example}
Let $A$ be as above and choose a prime number $p\geq 7$ such that
\begin{itemize}
 \item[(*)] for $i=1,\dots,h$, for any subgroup $H$ of $A_i$ contained in $A[N]$, and for any nontrivial $\varphi\in \Aut_0(A_i/H)$,
       $p$ is larger than the order of $(A_i/H)^\varphi = \{ x \in A_i/H : \varphi(x)=x \}$.
\end{itemize}
Notice that if $\varphi$ is a nontrivial automorphism of a simple abelian variety then $\varphi$ has only finitely many fixed points, so
a prime number $p$ with this property exists.

Let $S/\mathbb{C}$ be a smooth hypersurface of degree $p$ in $\mathbb{P}^3$ with $\operatorname{Aut}(S) \cong \mathbb{Z}/p\mathbb{Z}$ and such
that every automorphism of $S$ acts on it without any fixed points;
an explicit example of such a hypersurface is given in Theorem \ref{thm:AutoS}.
Let $G = \operatorname{Aut}(S) \cong \mathbb{Z}/p\mathbb{Z}$ and set $X := S / G$.
{We now proceed to describe some basic properties of $X$ (§\ref{sect:PropertiesX}),
construct a certain smooth projective variety $Y$ of dimension $2+\dim A$ (§\ref{sect:ConstructionY}),
and prove that $Y$ has automorphism group isomorphic to $A$ (Theorem \ref{thm:AutoY} in §\ref{sect:AutY}).
}

\subsubsection{Properties of $X$}\label{sect:PropertiesX}

\begin{lemma}\label{lemma:XSmoothProjective}
$X$ is a smooth projective variety.
\end{lemma}
\begin{proof}
$X$ is smooth since $G$ acts on $S$ without fixed points, and is projective since any quotient of a projective variety by a finite group of automorphisms is projective (see \cite[Remarque on page 51]{MR0098097}).
\end{proof}
\begin{lemma}\label{lemma:TrivialAutomorphisms}
$X$ does not admit any nontrivial automorphisms.
\end{lemma}
\begin{proof}
Let $\varphi : X \to X$ be an automorphism. Composing with the natural projection $\pi : S \to X$, we obtain a map $\varphi \circ \pi : S \to X$ which, since $S$ is simply connected, lifts to a map $\tilde{\varphi} : S \to S$. Clearly $\tilde{\varphi}$ is algebraic, and it is easily seen to be a covering map, so it is an isomorphism since $S$ is connected and simply connected. It follows that $\tilde{\varphi} : S \to S$ is in $G$, hence (by passing to the quotient) it induces the identity on $X$. Since on the other hand $\tilde{\varphi}$ induces $\varphi$ on $X$, we get $\varphi = \operatorname{id}_X$ as claimed.
\end{proof}

\begin{lemma}\label{lemma:GeneralType}
$X$ has Kodaira dimension $2$.
\end{lemma}
\begin{proof}
Kodaira dimension is invariant under finite étale covers, so $\operatorname{kod}(X)=\operatorname{kod}(S)$. By adjunction, $K_S = O_{\mathbb{P}^3}(p-3-1)|_S$ is ample, so $\operatorname{kod}(S)=\operatorname{dim}(S)=2$.
\end{proof}

\begin{lemma}\label{lemma:MapsToA}
The Albanese variety of $X$ is trivial, therefore there are no non-constant maps from $X$ to any abelian variety.
\end{lemma}
\begin{proof}
Clearly $S$ is the universal cover of $X$, so $\pi_1(X)$ is isomorphic to $\operatorname{Aut}(S \to X) \cong \mathbb{Z}/p\mathbb{Z}$
and in particular is finite.
Since the Albanese variety of $X$ is dual to its Picard variety, one has
$\dim \operatorname{Alb}(X)=\operatorname{dim} \operatorname{H}^1(X,\mathcal{O}_X) =h^{1,0}(X)$; on the other hand, the fact that
$\pi_1(X)$ is finite implies that $H_1(X,\mathbb{Q})$ is trivial, so $h^{1,0}(X) \leq h^1(X)=\dim H^1(X,\mathbb{C})=0$, hence
$\operatorname{Alb}(X)$ is trivial as claimed.
\end{proof}

\subsubsection{A nontrivial {$A$-torsor} $Y \to X$}\label{sect:ConstructionY}

\begin{definition}\label{def:Y} Fix an isomorphism $\chi : G \to \mathbb{Z}/p\mathbb{Z}$ and a point $P$ such that
\begin{itemize}
 \item[(**)] $P$ is a $p$-torsion point of $A$ which is not contained
in any proper abelian subvariety of $A$.
\end{itemize}
The abelian subvarieties of $A$ are all of the form $A_{i_1}+\dots+A_{i_k}$, so a point with this property exists.
We let $\mathbb{Z}/p\mathbb{Z}$ act on the group generated by $P$ in the obvious way
(that is, for $n \in \mathbb{Z}$ the class of $n$ in $\mathbb{Z}/p\mathbb{Z}$ sends $P$ to $nP$).
We set $Y = (S \times A) / G$, where the action of $G$ on the product $S \times A$ is given by
\[
g \cdot (s,a) = (g \cdot s, a+\chi(g)P).
\]
\end{definition}
As in the proof of Lemma \ref{lemma:XSmoothProjective}, it is easy to see that $Y$ is a smooth projective variety;
moreover, $Y$ has a natural structure of principal space under $A$. Indeed
for each $b \in A$, the translation map
\[
\begin{array}{ccc}
S \times A & \to & S \times A \\
(s,a) & \mapsto & (s,a+b)
\end{array}
\]
commutes with the action of $G$, so it descends to an automorphism of $Y=(S \times A)/G$ that we denote by $y\mapsto b+y$ or by $\tau_b$.
This defines an action of $A$ on $Y$ which is free and transitive along the fibers of the map $Y\to X$.
Moreover, $Y \to X$ is an $A$-torsor in the analytic (and in fact even étale) topology: indeed, $S$ is an
étale covering of $X$, and the pullback of $Y$ to $S$ is trivial.

\begin{lemma}\label{lemma:YNonTrivial}
{The map $Y \to X$ does not admit a section (in the analytic topology).}
\end{lemma}
\begin{proof}
Notice that $Y \to X$ {admits a section} if and only if it is trivial as a torsor.
Indeed if $Y\to X$ has a section $s$ then the map $A\times X\to Y$ given by $(a,x)\mapsto a+s(x)$ is an isomorphism of {torsors}.
Let $\mathcal{A}$ be the sheaf of holomorphic functions on $X$ with values in $A$; {$A$-torsors on $X$ are classified by $H^1(X,\mathcal{A})$, where the cohomology is taken in the analytic category}. For any fixed $n>0$, consider the exact sequence of sheaves on $X$
\[
0 \to \mathcal{A}[n] \to \mathcal{A} \xrightarrow{[n]} \mathcal{A} \to 0
\]
and take cohomology to obtain the long exact sequence
\[
0 \to H^0(X,\mathcal{A}[n]) \to H^0(X,\mathcal{A}) \xrightarrow{[n]} H^0(X,\mathcal{A}) \to H^1(X,\mathcal{A}[n]) \to H^1(X,\mathcal{A}).
\]
{By Serre's GAGA principle,} all maps from $X$ to $A$ are algebraic, so by Lemma \ref{lemma:MapsToA} we have $H^0(X,\mathcal{A}) = A$, and $H^0(X,\mathcal{A}) \xrightarrow{[n]} H^0(X,\mathcal{A})$ is just $A \xrightarrow{[n]} A$, which is surjective. It follows in particular that the natural arrow
\begin{equation}\label{eq:injectivityGeneral}
H^1(X,\mathcal{A}[n]) \to H^1(X,\mathcal{A})
\end{equation}
is injective.
Consider $Z := (S \times \langle P \rangle) / G \hookrightarrow Y$, where $\langle P \rangle$ denotes the order $p$ subgroup of $A(\mathbb{C})$ generated by $P$.
By the injectivity of \eqref{eq:injectivityGeneral} (with $n=p$), proving that $Z$ is a nontrivial covering space of $X$ suffices to show that $Y \to X$ is a nontrivial {torsor}. But this is
clear, because the natural map $S \to S \times A \to (S \times \langle P \rangle) / G$ is injective and surjective, hence (since $S$ is compact) a homeomorphism. It follows
that $Z \cong S$ is a nontrivial cover of $X$ as desired.
\end{proof}

\subsection{Determination of $\operatorname{Aut}(Y)$}\label{sect:AutY}
In this section we show:

\begin{theorem}\label{thm:AutoY}
The automorphism group of $Y$ is isomorphic to $A$.
\end{theorem}

\subsubsection{Preliminaries on simple abelian varieties}
We shall need the following basic fact about simple abelian varieties.

\begin{lemma}\label{lemma:DeterminantOnCohomology}
Let $T$ be a projective complex torus. Let $A$ be the abelian variety obtained from $T$ by fixing an arbitrary origin; notice that $T$ is naturally a torsor under $A$. Finally let $\alpha$ be an automorphism of $T$ (as a projective variety) and assume that $A$ is simple. Then:
\begin{enumerate}
   \item if $\alpha$ is translation by a point of $A$, then the determinant of $(1-\alpha)_* : H_1(T,\mathbb{Q}) \to H_1(T,\mathbb{Q})$ is $0$;
   \item if $\alpha$ is \textit{not} translation by a point of $A$, then $\alpha$ has at least one fixed point and the determinant of
         $(1-\alpha)_* : H_1(T,\mathbb{Q}) \to H_1(T,\mathbb{Q})$ is the number of fixed points of $\alpha$.
\end{enumerate}
\end{lemma}
\begin{proof}
The statement of (1) is obvious, because translations induce the identity on $H_1(T,\mathbb{Q})$.
Assume now that $\alpha$ is not a translation and identify $T$ with $A$ by choosing a point $t_0 \in T$ as the origin. We prove first that $\alpha$ has at least one fixed point.
Letting $a=\alpha(t_0)-t_0$ we have $\alpha(t)=\varphi(t)+a$, where $\varphi\in \Aut_0(A)$ is different from the identity. Let $\psi=\varphi-\operatorname{id}_A :A\longrightarrow A$; it is an endomorphism of $A$, and since $A$ is simple and $\varphi$ is nontrivial the image of $\psi$
is $A$ itself. One checks that $b \in T$ is a fixed point of $\alpha$ if and only if $\psi(b)=-a$. As $\psi$ is surjective, such $b$ exist, and there are only finitely many of them because the set $\{ b : \psi(b)=-a\}$ is naturally a torsor under the finite group $\ker \psi$.
We can then choose the origin $t_0$ to be a fixed point of $\alpha$, in which case $\alpha$ belongs to $\operatorname{Aut}_0(A)$ and we have $\psi(t)=\alpha(t)-t$,
so that $A^\varphi$ is equal to the kernel of $\psi$ and its order is the degree of $\psi$. The lemma follows from the fact that
for a complex torus $H_n(\psi,\mathbb Q)=\det(\psi_*:H_1(T,\mathbb Q)\longrightarrow H_1(T,\mathbb Q))$.
\end{proof}

\subsubsection{Preliminaries on surfaces of Kodaira dimension 2}
We shall need the following consequence of \cite{MR2393263}.

\begin{lemma}\label{prop:ImagesOfAbelianVarieties}
Let $S$ be a surface of Kodaira dimension 2 and $A$ be an abelian variety.
The image of any morphism $f : A \to S$ is either a point or a (possibly singular) irreducible curve of geometric genus at most one.
\end{lemma}
\begin{proof}
Suppose by contradiction that $f$ is surjective. Then by \cite[Theorem 1.1]{MR2393263} the surface $S$ admits a finite étale cover which is a product of projective spaces and an abelian variety: since Kodaira dimension is invariant under finite étale covers, this contradicts the fact that the Kodaira dimension of $S$ is 2, because any such product has non-positive Kodaira dimension. So the image of $f$ can only be a single point or a curve, which is then automatically irreducible since $A$ is. Suppose that the image of $f$, call it $Z$, is a curve, and let $\tilde{Z}$ be its normalization. By the universal property of normalization, $f$
induces a (dominant, hence surjective) map from $A$ to $\tilde{Z}$; applying \cite[Theorem 1.1]{MR2393263} again we obtain that $\tilde{Z}$ is covered by either $\mathbb{P}^1$ or an elliptic curve, which proves the statement about the genus.
\end{proof}

\subsubsection{Proof of Theorem \ref{thm:AutoY}}
We already noticed that $A$ injects into $\operatorname{Aut}(Y)$.
For the other inclusion let $\varphi$ be an automorphism of $Y$. We prove first that $\varphi$ preserve the fibers of the map $\pi:Y\longrightarrow X$.
For each $x \in X$, let $Y_x$ be the fiber of $\pi$ over $x$ and let
\[
\varphi_x : Y_x \hookrightarrow Y \xrightarrow{\varphi} Y \xrightarrow{\pi} X.
\]
Suppose that for general $x$ the image of $Y_x$ is not reduced to a single point: then Lemma \ref{prop:ImagesOfAbelianVarieties} implies
that generically the image of $\varphi_x$ is a (possibly singular) curve of genus at most 1. By \cite[Proposition VII.2.1]{MR2030225},
a surface of Kodaira dimension 2 admits no algebraic system (of positive dimension) of effective divisors whose general member is a
(possibly singular) rational or elliptic curve. By Lemma \ref{lemma:GeneralType} we know that $X$ is a surface of Kodaira dimension 2, so it follows that $\varphi_x$ is
constant for all $x \in X$. In particular,
\[
Y \xrightarrow{\varphi} Y \to A\backslash Y = X
\]
descends to a map $\varphi_X : X \to X$, which is easily seen to be biregular (its inverse being $(\varphi^{-1})_X$), and hence an automorphism.
It follows from Lemma \ref{lemma:TrivialAutomorphisms} that $\varphi_X$ is the identity, which implies that the equality $\varphi (Y_x)=Y_x$
holds for all $x \in X$. Thus we see that for every $x \in X$ the automorphism $\varphi$ of $Y$ induces an automorphism $\varphi |_{Y_x}$ of $Y_x$.

\medskip

Thus, locally in the analytic topology, the automomorphism $\varphi$ can be described as follows.
For each $x\in X$ we can choose an open connected neighborhood $U \subset X$ of $x$ such that
$V=\pi^{-1}(U)$ can be identified with $U\times A$ {(as an $A$-torsor)} and $\varphi(u,a)=(u,\phi(u,a))$.
Let $r:U\longrightarrow A$ be defined by $r(u)=\phi(u,0)-0$. Then
$a\mapsto \phi(u,a)-r(u)$ is an automorphism of $A$ as an algebraic group,
and since $\Aut_0(A)$ is finite it must be equal to an automorphism $\phi$ independent of $u$.
Hence {$\varphi(u,a)=(u,\phi(a)+r(u))$} and $\psi=\phi-id_A$ is an endomorphism of $A$ as an algebraic group.
{Furthermore, since any two identifications of a fiber of $Y \to X$ with the trivial $A$-torsor differ only by a translation,
we see that the endomorphism $\psi$ thus obtained is independent of our choice of $U$ and of the local trivialization $\pi^{-1}(U) \cong U \times A$.
}

\medskip

We now prove the theorem by induction on $h$, the number of simple factors of $A$.
Assume first that $h=1$, so that $A$ is simple. For $x\in X$ we define
$$n(x) = \det \left( (1-\varphi|_{Y_x})_* \bigm\vert H_1(Y_x,\mathbb{Q}) \right);$$
it is a continuous function on $X$.
Since $X$ is connected and $\mathbb{Z}$ is discrete, it follows that $n(x)$ is actually constant:
let $n$ be the common value of the various $n(x)$. We show that $n=0$.
Suppose by contradiction that $n > 0$. Let $\tilde X=Y^\varphi$ and let $\tilde \pi$ be the restriction of $\pi$ to $\tilde X$.
We prove that $\tilde \pi$ is an $n$-to-1 covering of $X$. The fact that it is $n$-to-1 follows from Lemma \ref{lemma:DeterminantOnCohomology}.
The claim that it is a covering can be checked locally using the analytic topology: using the local description above we obtain
$V^\varphi=\{(u,a): \psi(a)=r(u)\}$, which is a covering of $U$.

If at least one of the connected components of $\tilde{X}$ is the trivial cover of $X$, then this gives a section of the projection map $\pi : Y \to X$,
contradicting Lemma \ref{lemma:YNonTrivial}.
Otherwise, take a connected subcover of $\tilde{X}$: this is a connected $m$-to-1 cover of $X$ for some $m \leq n$ which is smaller than $p$ by our assumption
(*) on $p$. This contradicts the fact that $\#\pi_1(X)=p$.

It follows that $n(x)=n=0$ for all $x$, hence by Lemma \ref{lemma:DeterminantOnCohomology} $\varphi|_{Y_x}$ is translation by a point $a(x) \in A$
(recall that $Y_x$ is naturally a torsor under $A$, so it makes sense to identify translations of $Y_x$ with elements of $A$). Now $x \mapsto a(x)$
gives a map $X \to A$, which is necessarily constant by Lemma \ref{lemma:MapsToA}, hence $\varphi$ is globally a translation by a point of $A$.

\medskip

We now prove the inductive step. Let $h>1$.
Since $\varphi$ preserves the fibers $Y_x$, composing with a translation by an element of $A$ we can assume that there exists
$y_0 \in Y$ such that $\varphi(y_0)=y_0$. We want to prove that in this case $\varphi$ is the identity.
\begin{comment}
Let $A'=A_1\backslash A$, and let $P'$ be the image of the point $P$ in $A'$ (see Definition \ref{def:Y}). If we define
$A_i'={(A_i\cap A_1)}\backslash A_i$, then $A_2',\dots, A_h'$ are the abelian factors of $A'$ and
$\Sigma'=(A_1\cap\Sigma)\backslash \Sigma$ is the kernel of the sum $\sigma':A_2'\times \dots\times A_h'\longrightarrow A'$.

{\color{red} Ho delle remore di natura formale; per esempio, strettamente parlando non ha senso intersecare $\Sigma$ con $A_1$, perché $\Sigma$ vive in $A_1 \times \cdots \times A_h$. Inoltre, per come capisco io l'operazione $\Sigma/(\Sigma \cap A_1)$, non penso venga il risultato giusto. Supponiamo che $A$ sia data da due curve ellittiche incollate lungo un sottogruppo ciclico (formalmente: siano $E_1, E_2$ due curve ellittiche (non isogene) e siano $P_1,P_2$ punti di $m$-torsione rispettivamente su $E_1, E_2$. Allora posso considerare il sottogruppo $H$ di $E_1 \times E_2$ dato da $\langle (P_1,P_2) \rangle$. Sia $A=(E_1 \times E_2)/H$). I fattori semplici di $A$ sono proprio $A_1=E_1, A_2=E_2$ (ognuna delle quali si immerge naturalmente in $A$), e $\Sigma$ è $H \subset E_1 \times E_2$. Ora $\Sigma \cap A_1$ (nel senso di $\Sigma \cap (A_1 \times \{0\})$) è banale, per cui $\Sigma/(\Sigma \cap A_1) = \Sigma$. Ma chiaramente $A'$ è semplice, per cui la "somma" è un isomorfismo, e $\Sigma'$ dovrebbe venire banale. Quello che avevi in mente (credo) è di prendere l'intersezione "in $A$", che penso sia la cosa giusta. Provo a scriverlo. \bigskip
\end{comment}
Let $\pi : A \to A':=A/A_1$ be the natural projection and set $A_i' := \pi(A_i)$ for $i=2,\ldots,h$. We let $P'=\pi(P)$ and write $\tilde{\pi} : A_1 \times \cdots \times A_h \to A_2' \times \cdots \times A_h'$ for the homomorphism
\[
\tilde{\pi}(a_1,\ldots,a_h) = (\pi(a_2), \ldots, \pi(a_h));
\]
finally, we set $\Sigma' := \tilde{\pi}(\Sigma)$. One then checks that the sum $\sigma' : A_2' \times \cdots \times A_h' \to A'$ is an isogeny with kernel $\Sigma'$.

Let $K = \ker(\Sigma \to \Sigma')$. For every $i=2,\ldots,h$, the intersection $A_1 \cap A_i$ embeds naturally into $K$,
so $N' \cdot \#(A_1 \cap A_i) \bigm\vert N' \cdot \#K = N$. It follows that every quotient of $A_i'=A_i/(A_1 \cap A_i)$ by a subgroup of $A_i'[N']$ is a quotient of $A_i$ by a subgroup of $A_i[N]$, so the analogue of condition (*) is satisfied by $A'$ and the prime $p$. It is immediate to check that (**) also holds for $A'$, $p$, and the point $P'$.
In particular, by induction, the automorphism group of $Y'=S\times ^G A'$ is equal to $A'$.

The projection map $S\times A\longrightarrow S\times A'$ is $G$-equivariant, so it induces a map $q:Y\longrightarrow Y'$ which we prove to be a
categorical quotient by the action of $A_1$. Indeed let $f:Y\longrightarrow Z$ be a $A_1$-invariant map.
It induces a $G\times A_1$-invariant map $f_1:S\times A\longrightarrow Z$
and therefore a map $f_2:S\times (A_1\times\dots\times A_h) \longrightarrow Z$ which is invariant by the action of {both $A_1$ and $\Sigma$ on the second factor.
Since the quotient of $A_1 \times \cdots \times A_h$ by the subgroup generated by $A_1$ and $\Sigma$ is $A'$,} the map $f_2$ induces a
regular map {$g_2:S\times A'\longrightarrow Z$ such that $f_2=g_2\circ(id_S\times \pi')$, where $\pi' := \pi \circ \sigma$ is the
natural map $A_1 \times \cdots \times A_h \to A'$. Since furthermore $f_2$} is $G$-invariant, {$g_2$} is also $G$-invariant,
hence it induces a map $g:Y'\longrightarrow Z$ such that $f=g\circ q$. Moreover, as $q$ is surjective, the map $g$ is unique.

We can now prove that $\varphi$ is the identity. For $a\in A$ denote by $\tau_a$ the translation by $a$ in $Y$.
Notice that for each $a$ and for each $x\in X$ there exists $\phi_x \in \Aut_0(A)$ such that
$$ \varphi\circ \tau_a \circ \varphi^{-1}=\tau_{\phi_x(a)}:Y_x\longrightarrow Y_x. $$
In particular, if $a\in A_1$, then $\phi_x(a)\in A_1$.
Being $Y'$ a categorical quotient of $Y$ by the action of $A_1$, we have that
$\varphi$ induces a map $\varphi':Y'\longrightarrow Y'$, which is an automorphism since $(\varphi^{-1})'$ is its inverse. Moreover, the image of $y_0$ in $Y'$ is
fixed by $\varphi'$, so $\varphi'$ is equal to the identity.

Hence $\varphi(y)-y \in A_1$ for all $y\in Y$. Arguing in the same way, but using $A_2$ instead of $A_1$, we obtain $\varphi(y)-y \in A_2$ for all $y$.
So $\varphi(y)-y\in A_1\cap A_2$ for all $y\in Y$, and since $A_1\cap A_2$ is finite and $\varphi(y_0)=y_0$ we obtain $\varphi(y)=y$ for all $y$. \qed

\subsection{A hypersurface in $\mathbb{P}^3$ with automorphism group $\mathbb{Z}/p\mathbb{Z}$}\label{sect:Hypersurface}In this section we explicitly construct, for every prime $p \geq 7$, an algebraic surface in $\mathbb{P}^3$ of degree $p$ whose automorphism group is cyclic of order $p$:
\begin{theorem}\label{thm:AutoS}
Let $p \geq 7$ be a prime number, and for $\lambda \in \mathbb{C}$ let $S_\lambda$ be the algebraic surface over $\mathbb{C}$ given by the zero locus in $\mathbb{P}^3$ of the homogeneous polynomial
\[
f_\lambda(x_1,x_2,x_3,x_4) := x_1^p + x_2^p + x_3^p + x_4^p + \lambda ( x_1^2 x_2^{p-4} x_3^2 + x_1^4 x_2^{p-6} x_4^2 ).
\]
The surface $S_\lambda$ is smooth for all but finitely many $\lambda \in \mathbb{C}$; if $\lambda \neq 0$, the automorphism group of $S_\lambda$ is cyclic of order $p$, generated by $[x_1 : x_2 : x_3 : x_4] \mapsto [x_1 : \zeta_p x_2 : \zeta_p^2 x_3 : \zeta_p^3 x_4]$, where $\zeta_p$ is a primitive $p$-th root of unity.
Moreover, each nontrivial element of $\operatorname{Aut}(S_\lambda)$ acts on $S_\lambda$ without any fixed points.
\end{theorem}

We start by noticing that for $\lambda=0$ the surface $S_0$ is smooth. Since being smooth is a Zariski-open condition in the defining polynomial, this shows that $S_\lambda$ is smooth away from a proper Zariski-closed subset of $\mathbb{C}$, that is, $S_\lambda$ is smooth for all but finitely many values of $\lambda$. From now on fix a nonzero value of $\lambda$ such that $S_\lambda$ is smooth, and to simplify the notation write $S$ for $S_\lambda$ and $f(x_1,x_2,x_3,x_4)$ for $f_\lambda(x_1,x_2,x_3,x_4)$.

By \cite[Theorem 2]{MR0168559} we know that all the automorphisms of $S$ are induced by
(linear) automorphisms of $\mathbb{P}^3$, so we only need to consider these.
Let $L : \mathbb{P}^3 \to \mathbb{P}^3$ be a linear transformation that satisfies
$L(S)=S$. We identify $L$ to the class $[M] \in \operatorname{PGL}_4(\mathbb{C})$ of a matrix $M=(M_{ij}) \in \operatorname{GL}_4(\mathbb{C})$. Furthermore, we let $e_1,\ldots,e_4$ be the canonical basis of $\mathbb{C}^4$ and denote by $\langle e_i \rangle$ the 1-dimensional $\mathbb{C}$-vector subspace of $\mathbb{C}^4$ generated by $e_i$.
We shall show Theorem \ref{thm:AutoS} in three steps: first we shall prove that $M$ either fixes or permutes the lines generated by $e_3$ and $e_4$; then we shall show that the same statement holds for the lines generated by $e_1$ and $e_2$; finally, we shall deduce from this that $M$ needs to be a diagonal matrix, at which point a direct computation concludes the proof. {This approach is inspired by \cite{MR2129679}.}

\subsubsection{Step 1: $M$ permutes $\langle e_3 \rangle$ and $\langle e_4 \rangle $}
The condition that $L(S)=S$ translates into the polynomial equality
\begin{equation}\label{eq:Automorphism}
f \circ M(x_1,\ldots,x_4) = \alpha f(x_1,\ldots,x_4)
\end{equation}
for some $\alpha \in \mathbb{C}^\times$.
Applying $\frac{\partial^2}{\partial x_i \partial x_j}$ to the two members of this equation and setting
\[
H_{ij}(x_1,\ldots,x_4) := \frac{\partial^2f}{\partial x_i \partial x_j}(x_1,\ldots,x_4)
\]
we find
\[
\sum_k \sum_m  M_{kj} M_{mi} H_{mk} (M(x_1,\ldots,x_4)) = \alpha H_{ij}(x_1,\ldots,x_4).
\]
Let $u, v$ be two vectors in $\mathbb{C}^4$. Multiplying the previous identity by $u_iv_j$ and summing over $i$ and $j$ we get
\begin{equation}\label{eq:SecondDerivatives}
\sum_{k,m} (Mv)_k (Mu)_{m} H_{mk}( M(x_1,\ldots,x_4)) = \alpha\sum_{i,j} H_{ij}(x_1,\ldots,x_4) u_i v_j.
\end{equation}
We now define a bilinear pairing
\[
\begin{array}{cccc}
\langle \cdot, \cdot \rangle : & \mathbb{C}^4 \times \mathbb{C}^4 & \to & \mathbb{C}[x_1,\ldots,x_4] \\
& (u,v) & \mapsto & \sum_{i,j} H_{i,j}(x_1,\ldots,x_4) u_i v_j,
\end{array}
\]
so that Equation \ref{eq:SecondDerivatives} reads
\[
\langle Mu, Mv \rangle \left(  M(x_1,\ldots,x_4) \right) = \alpha \langle u,v \rangle.
\]
In particular, since $M$ is invertible we obtain:
\begin{proposition}\label{prop:Orthogonality}
Let $u, v$ be vectors in $\mathbb{C}^4$. The equalities $\langle u,v \rangle=0$ and $\langle Mu, Mv \rangle=0$ are equivalent.
\end{proposition}

\begin{lemma}\label{lemma:Orthogonal}
Let $a, b \in \mathbb{C}^4$ be two nonzero vectors such that $\langle a,b \rangle=0$. Then there exist $\lambda, \mu \in \mathbb{C}^\times$ such that either $a=\lambda e_3, b=\mu e_4$, or $a = \lambda e_4, b = \mu e_3$ hold.
\end{lemma}

\begin{proof}
Write $a=(a_1,a_2,a_3,a_4)$ and $b=(b_1,b_2,b_3,b_4)$.
By direct inspection, one checks that, for $i=1,2,3,4$, the only second derivative of $f$ involving the monomial $x_i^{p-2}$ is $H_{ii}$. This immediately implies that $a_ib_i=0$ for $i=1,\ldots,4$, and by symmetry we can assume $a_1=0$.
The coefficients of the monomials $x_1 x_2^{p-5} x_3^2$, $x_1 x_2^{p-4} x_3$ and $x_1^3 x_2^{p-6} x_4$ in $\langle a,b \rangle$ are given by
$2 \lambda (p-4) (a_2 b_1 + a_1 b_2), 4 \lambda (a_3 b_1 + a_1 b_3)$ and $8 \lambda (a_4 b_1 + a_1 b_4)$ respectively, so under our assumptions $\langle a,b \rangle=0$, $\lambda \neq 0$ and $a_1=0$ we obtain $b_1a_2=b_1a_3=b_1a_4=0$. If we had $b_1 \neq 0$, this would imply $a=(0,0,0,0)$, contradicting our assumptions, so we must have $b_1=0$ as well. The situation is now again symmetric in $a,b$, so we might assume $a_2=0$. Arguing as before (but looking at the monomials $x_1^2 x_2^{p-5} x_3$ and $x_1^4 x_2^{p-7} x_4$) one finds $a_3b_2=a_4b_2=0$, so that $b_2=0$ as well. The conclusion now follows easily from the equalities $a_3b_3=a_4b_4=0$.
\end{proof}

\begin{corollary}\label{cor:Me1e2}
One of the following holds:
\begin{itemize}
\item $M \langle e_3 \rangle = \langle e_3 \rangle$ and $M \langle e_4 \rangle =\langle e_4 \rangle$;
\item $M \langle e_3 \rangle = \langle e_4 \rangle$ and $M \langle e_4 \rangle =\langle e_3 \rangle$.
\end{itemize}
\end{corollary}
\begin{proof}
Apply Proposition \ref{prop:Orthogonality} to $u=e_3$ and $v=e_4$: since $\langle e_3,  e_4\rangle=H_{34}=0$ we obtain
$
\langle M e_3, M e_4 \rangle = 0
$. The claim then follows from the previous lemma.
\end{proof}

\subsubsection{Step 2: $M$ permutes $\langle e_1 \rangle$ and $\langle e_2 \rangle$}

Arguing as in the previous section, it is easily seen that if we let $A : (\mathbb{C}^4)^p \to \mathbb{C}$ denote the multilinear form
\[
A: (u_1,\ldots,u_p) \mapsto \sum_{i_1,\ldots,i_p} \frac{\partial^p f}{\partial x_{i_p} \cdots \partial x_{i_p}}
(u_1)_{i_1} \cdots (u_p)_{i_p},
\
where $(u_i)_j$ is the $j$-th coordinate of $u_i$,
we have $A(Mu_1,\ldots,Mu_p) = \beta A(u_1,\ldots,u_p)$ for some $\beta \in \mathbb{C}^\times$;
notice that here we do not need to compose with $M$ on the left
hand side, because $p$-th derivatives of $f$ are just scalars. Suppose that $M \langle e_3 \rangle = \langle e_3 \rangle$ and $M \langle e_4 \rangle= \langle e_4 \rangle$; the case $M \langle e_3 \rangle = \langle e_4 \rangle$ and $M \langle e_4 \rangle = \langle e_3 \rangle$ is completely analogous. Rescaling $M$ if necessary (which we can do, since we are only interested in its projective class) we can assume $Me_3=e_3$.
Choosing $u_1=\cdots=u_{p-1}=e_3$ and $u_p=e_1$ we have
\[
\beta A(e_3,\ldots,e_3,e_1) = \beta \frac{\partial^p f}{\partial x_3^{p-1} \partial x_1}=0,
\]
from which we deduce
\[
0 = A(Me_3,\ldots,Me_3,Me_1) = A(e_3,\ldots,e_3,Me_1)= \sum_{i_p} \frac{\partial^p f}{\partial x_3^{p-1} \partial x_{i_p}} (Me_1)_{i_p};
\]
since the only nonvanishing partial derivative of the form $\frac{\partial^p f}{\partial x_3^{p-1} \partial x_{i_p}}$ is $\frac{\partial^p f}{\partial x_3^p}$, this implies $M_{31}=0$. Similary, the choice $(e_4,\ldots,e_4,e_1)$ shows $M_{41}=0$, while the choices $(e_3,\ldots,e_3,e_2)$ and $(e_4,\ldots,e_4,e_2)$ give $M_{32}=M_{42}=0$.
It follows that $M$ sends the 2-plane $\{x_3=x_4=0\}$ to itself; in particular, $M$ induces an automorphism of the finite set of points in $\mathbb{P}^3$ defined by the equations
\[
f(x_1,x_2,x_3,x_4)=0, \quad x_3=x_4=0 \quad \Longleftrightarrow \quad x_3=x_4=0, \quad x_1^p+x_2^p=0.
\]
From this it is immediate to deduce:
\begin{corollary}\label{cor:Me3e4}
One of the following holds:
\begin{itemize}
\item $M \langle e_1\rangle = \langle e_1 \rangle$ and $M \langle e_2 \rangle = \langle e_2 \rangle$;
\item $M \langle e_1 \rangle = \langle e_2 \rangle$ and $M \langle e_2 \rangle = \langle e_1 \rangle$.
\end{itemize}
\end{corollary}

\subsubsection{Step 3: determination of $\operatorname{Aut}(S)$}
Corollaries \ref{cor:Me1e2} and \ref{cor:Me3e4}
tell us that $M$ either fixes or permutes the lines $\langle e_1 \rangle, \langle e_2 \rangle$, and that the same holds for the lines $\langle e_3 \rangle, \langle e_4 \rangle$.
One checks easily that if $M$ exchanges $\langle e_1 \rangle$ with $ \langle e_2 \rangle$, and/or it exchanges $\langle e_3 \rangle$ with $\langle e_4 \rangle$, then $f \circ M$ is not a scalar multiple of $f$, so that $M$ needs to be a diagonal matrix. Normalize $M$ so that $M_{11}=1$ and write $M=\operatorname{diag}(1,\mu_2,\mu_3,\mu_4)$: replacing in Equation \eqref{eq:Automorphism} and comparing the coefficients of $x_1^p$ on the two sides we find $\alpha=1$. Comparing the coefficients of $x_i^p$ for $i=2,3,4$ we then obtain $\mu_i^p=1$ for $i=2,3,4$, so that $\mu_2,\mu_3,\mu_4$ are $p$-th roots of unity. It is now immediate to check that the only automorphisms of $S$ are represented by the powers of the (order $p$) matrix
$
\begin{pmatrix}
1 \\ & \zeta_p \\ & & \zeta_p^2 \\ & & & \zeta_p^3
\end{pmatrix},
$
where $\zeta_p$ is a primitive $p$-th root of unity. The fixed points (in $\mathbb{P}^3$) for the action of this matrix (or any of its powers, with the exception of the identity) are $[1:0:0:0], [0:1:0:0], [0:0:1:0], [0:0:0:1]$, none of which lies on the hypersurface $f(x_1,x_2,x_3,x_4)=0$.
This concludes the proof of Theorem \ref{thm:AutoS}.

\bibliographystyle{alpha}