\section{Introduction}

We consider a quantum particle at time $t=0$ with wave function $\psi_0$ in a region $\Omega\subset\mathbb{R}^3$ of physical space, surrounded by detectors placed along the boundary $\partial \R$, and ask for the probability distribution $\mu$ of $Z=(T,{\boldsymbol{X}})$, where $T\geq 0$ is the time at which a detector clicks, and ${\boldsymbol{X}}$ is the location on $\partial \R$ where the particle gets detected; if no detector ever clicks, we write $Z=\infty$. As we have discussed in \cite{detect-rule}, quantum mechanics in principle makes a prediction for $\mu$, assuming that the wave function of all detectors were exactly known and we could solve the Schr\"odinger equation for a macroscopic number of particles. In practice, however, it is desirable to have a simple mathematical rule for computing $\mu$ from $\psi_0$ for an \emph{ideal detector}. For comparison, consider the standard rule that if we make a quantum measurement of position at time $t$ on a particle with wave function $\varphi_t$ then the outcome has $|\varphi_t|^2$ distribution; also this rule is certainly highly idealized and requires no detailed information about the nature or state of the measurement apparatus. A particularly convincing rule for $\mu$ from $\psi_0$ for ideal detectors along a surface $\partial \R$ in the non-relativistic case is the \emph{absorbing boundary rule} \cite{detect-rule}, based on evolving $\psi_0$ according to the Schr\"odinger equation with an ``absorbing'' boundary condition (ABC) on $\partial \R$ first considered by Werner \cite{Wer87},
\begin{equation}\label{ABC}
\boldsymbol{n}(\boldsymbol{x}) \cdot \nabla \psi(\boldsymbol{x}) = i\kappa \psi(\boldsymbol{x})
\end{equation}
at all $\boldsymbol{x}\in\partial \R$, where $\boldsymbol{n}(\boldsymbol{x})$ is the outward unit normal vector on $\partial \R$ and $\kappa>0$ a constant (the detector's wave number of sensitivity).
The ABC entails that the probability current across $\partial \R$ always points outward. The rule can be expressed by saying that $Z=(T,{\boldsymbol{X}})$ is the time and place where the random Bohmian trajectory starting out with the $|\psi_0|^2$ distribution, and guided by a wave function evolving according to the Schr\"odinger equation with ABC, hits $\partial \R$. (We note that this trajectory is different from what it would be in the absence of detectors.) Other rules for the detection time distribution have been proposed in the literature; see, e.g., \cite{ML00,MSE08,MRC09} for an overview.

In this paper, we describe a relativistic analog of the absorbing boundary rule using the Dirac equation  instead of the non-relativistic Schr\"odinger equation. A key element is again a boundary condition on the detecting surface, an absorbing boundary condition for the Dirac equation (ABCD). Like in the non-relativistic case, the Hamiltonian with an ABCD is not self-adjoint and the time-evolution not unitary. That is because it is part of the setup that, upon detection, the particle gets absorbed (or removed from consideration); since $|\psi_t(\boldsymbol{x})|^2\, d^3\boldsymbol{x}$ represents the probability that the particle is located at time $t$ in the volume $d^3\boldsymbol{x}$, we have that $\|\psi_t\|^2 = \int_{\Omega} d^3\boldsymbol{x} \, |\psi_t(\boldsymbol{x})|^2$ equals the probability that the particle has not been detected up to time $t$. And again, the time and place of detection can be expressed as the time and place where the Bohmian trajectory hits the boundary. Among the family of suitable boundary conditions, there is one simplest choice for every Lorentz frame (detector rest frame); we propose this choice as part of the definition of \emph{ideal} detector. We note that the absorbing boundary rule (in either its non-relativistic or the present relativistic form) is inequivalent to the detection time observables considered by Werner in \cite{Wer86}.

While in the non-relativistic case, the boundary condition involves the parameter $\kappa$ (such that $\hbar^2\kappa^2/2m$ represents the particle energy at which the detector is maximally efficient), the  condition proposed here (the ``ideal'' ABCD) does not involve such a parameter, although it does involve two other parameters, representing the choice of the detector rest frame, which do not show up in the non-relativistic case. In fact, the ideal ABCD does not possess a non-relativistic limit. For this reason, we also consider another family of absorbing boundary conditions for the Dirac equation that we call \emph{semi-ideal ABCDs} which do possess a non-relativistic limit, given by the boundary condition
\begin{equation}\label{ABC2}
{\boldsymbol{\sigma}}\cdot \nabla \phi(\boldsymbol{x}) = i\kappa\, \boldsymbol{n}(\boldsymbol{x})\cdot {\boldsymbol{\sigma}}\, \phi(\boldsymbol{x})\,,
\end{equation}
where $\phi$ is a $\mathbb{C}^2$-valued non-relativistic wave function governed by the Pauli equation, and ${\boldsymbol{\sigma}}=(\sigma_1,\sigma_2,\sigma_3)$ are the three Pauli matrices. Also \eqref{ABC2} is a non-relativistic absorbing boundary condition, similar but not equivalent to \eqref{ABC}.

We also describe in this paper variants of the ABCD rule for moving detectors, for curved space-time, and for several detectable particles; similar extensions of the non-relativistic rule to moving detectors and several particles are described in \cite{detect-several}. For a discussion of another boundary condition for the Dirac equation that leads to a self-adjoint Hamiltonian, see \cite{FR15}; for a general discussion of boundary conditions for Dirac operators, see \cite{BB13}; for a discussion of detection time on a lattice, see \cite{DTT15}; an uncertainty relation between detection time and energy in the non-relativistic case is derived in \cite{uncertainty}.

It is often assumed that only positive-energy states (i.e., wave functions from the subspace of Hilbert space $L^2(\mathbb{R}^3,\mathbb{C}^4)$ associated with the positive half of the spectrum of the free Dirac Hamiltonian) are physical. For our discussion we need to drop this assumption, as the initial wave function $\psi_0$ needs to be concentrated inside $\Omega$, while positive-energy states cannot be so concentrated but have nonzero tails over all of $\mathbb{R}^3$.

The remainder of this paper is organized as follows. In Section~\ref{sec:ruleflat} we state and discuss the proposed rule for the Dirac equation in flat space-time for ideal detectors at rest relative to some fixed Lorentz frame. In Section~\ref{sec:rulecurved}, we provide a version of the rule for curved space-time and moving ideal detectors. In Section~\ref{sec:several}, we discuss the case of several particles. In Section~\ref{sec:semi}, we discuss another model that we call \emph{semi-ideal detectors} and derive that its non-relativistic limit agrees with the non-relativistic absorbing boundary rule based on \eqref{ABC}.

\section{Absorbing Boundary Rule for a Single Dirac Particle in Minkowski Space-Time}
\label{sec:ruleflat}

We begin with the simplest case, of a single particle in flat space-time and detectors along a surface in 3-dimensional space that is at rest in some Lorentz frame, which we take to be the rest frame of each detector in the continuous family of detectors placed along $\partial \R$.

\subsection{Statement of the Proposed Rule}

Let $\psi:[0,\infty)\times \Omega \to \mathbb{C}^4$ be the unique solution of the Dirac equation
\begin{equation}\label{Dirac}
ic\hbar \gamma^\mu \partial_\mu \psi = mc^2 \psi
\end{equation}
with initial condition
\begin{equation}\label{ic}
\psi(0,\boldsymbol{x}) = \psi_0(\boldsymbol{x}) \quad \text{for all }\boldsymbol{x}\in\Omega,
\end{equation}
and boundary condition
\begin{equation}\label{abcd1}
\boldsymbol{n}(\boldsymbol{x}) \cdot {\boldsymbol{\alpha}}\: \psi(t,\boldsymbol{x}) = \psi(t,\boldsymbol{x})
\quad \text{for all }t> 0,\boldsymbol{x}\in\partial \R.
\end{equation}
We call \eqref{abcd1} an absorbing boundary condition for the Dirac equation (ABCD).
As we explain elsewhere \cite{detect-thm}, the Hille-Yosida theorem implies that this initial value-boundary value problem \eqref{Dirac}--\eqref{abcd1} has a unique solution (in the appropriate sense) for every $\psi_0\in L^2(\Omega,\mathbb{C}^4)$.

Let $j^\mu$ be the 4-vector field usually called the probability current of $\psi$,
\begin{equation}\label{jdef}
j^\mu = \overline{\psi}\gamma^\mu \psi
\end{equation}
or, equivalently,
\begin{equation}
j=(j^0,{\boldsymbol{j}})=(|\psi|^2, \psi^\dagger{\boldsymbol{\alpha}}\psi)\,.
\end{equation}
Suppose $\|\psi_0\|^2=\int_\Omega d^3\boldsymbol{x}\, |\psi_0(\boldsymbol{x})|^2 =1$, where $|\psi|^2$ means $\sum_{s=1}^4 |\psi_s|^2$. Then, the proposed rule asserts, the distribution $\mu$ of $Z$ satisfies
\begin{align}\label{probnjR1}
  \mu \Bigl( t_1 \leq T<t_2, {\boldsymbol{X}} \in B \Bigr)
  &= \int\limits_{t_1}^{t_2} dt \int\limits_{B} d^2\boldsymbol{x} \; \boldsymbol{n}(\boldsymbol{x}) \cdot
  {\boldsymbol{j}}^{\psi_t}(\boldsymbol{x})\\
  &= \int\limits_{t_1}^{t_2} dt \int\limits_{B} d^2\boldsymbol{x}\; |\psi_t(\boldsymbol{x})|^2 \label{probpsiR}
\end{align}
for any $0\leq t_1<t_2$ and any set $B\subseteq \partial \R$. Note that the formulas \eqref{probnjR1} and \eqref{probpsiR} are equivalent by virtue of the ABCD \eqref{abcd1}. Furthermore,
\begin{equation}\label{Zinfty}
  \mu (Z=\infty) =
  1-\int\limits_{0}^{\infty} dt \int\limits_{\partial \R} d^2\boldsymbol{x} \; \boldsymbol{n}(\boldsymbol{x}) \cdot
  {\boldsymbol{j}}^{\psi_t}(\boldsymbol{x})\,.
\end{equation}
This completes the statement of the rule.

Equivalently, the rule can be expressed in terms of the Bohmian trajectory $X^\mu(s)$ \cite{Bohm53}, i.e., the integral curve of $j^\mu$,
\begin{equation}
\frac{dX^\mu}{ds}\propto j^\mu(X(s))
\end{equation}
with arbitrary curve parameter $s$ (which can, but does not have to, be taken to be coordinate time $t$ or proper time $\tau$, the latter as long as the integral curve is timelike) and initial point
\begin{equation}
X(0)=\bigl(0,{\boldsymbol{X}}(0)\bigr)\,,
\end{equation}
where ${\boldsymbol{X}}(0)\in\Omega$ is random with $|\psi_0|^2$ distribution. The Bohmian trajectory is guided by $\psi$ evolving according to \eqref{Dirac} and \eqref{abcd1}. The space-time point $Z=(T,{\boldsymbol{X}})$ of detection is where the Bohmian trajectory hits $[0,\infty)\times \partial \R$ (and ends), and $Z=\infty$ if the Bohmian trajectory never hits the boundary.

Should the experiment be terminated at time $t>0$ without any detection having occurred, then the collapsed wave function, which becomes the initial wave function at time $t$ for any subsequent experiment, is $\psi_t/\|\psi_t\|$.

\subsection{Discussion}

Here are some comments on the mathematics involved. The boundary condition \eqref{abcd1} means that the spinor $\psi(t,\boldsymbol{x})\in\mathbb{C}^4$ is an eigenvector of the matrix $\alpha_n(\boldsymbol{x})=\boldsymbol{n}(\boldsymbol{x})\cdot {\boldsymbol{\alpha}}$, the Dirac alpha matrix associated with the normal direction, with eigenvalue $+1$. It is known that, for every unit vector $\boldsymbol{e}\in\mathbb{R}^3$, $\boldsymbol{e}\cdot{\boldsymbol{\alpha}}$ has eigenvalues $\pm1$, each with multiplicity 2. Thus, \eqref{abcd1} constrains $\psi(t,\boldsymbol{x})$ to lie in a particular 2-dimensional subspace of $\mathbb{C}^4$. It is common \cite{BB13} that boundary conditions on the Dirac equation are of this form, requiring $\psi$ at each boundary point to lie in a particular 2-dimensional subspace of 4-dimensional spin space.

The ABCD \eqref{abcd1} implies that the Bohmian particle, when it reaches the boundary, moves at the speed of light in the outward direction normal to the boundary. So the Bohmian particle can cross $\partial \R$ only outward.

Since it also follows from the ABCD \eqref{abcd1} that
$\boldsymbol{n}(\boldsymbol{x})\cdot {\boldsymbol{j}}^{\psi_t}(\boldsymbol{x})=|\psi(t,\boldsymbol{x})|^2$,
the density of $\mu$ relative to $dt\, d^2\boldsymbol{x}$ is non-negative. Since the Dirac equation \eqref{Dirac} implies that
\begin{equation}\label{continuity}
\partial_\mu j^\mu =0\,,
\end{equation}
we obtain by integrating \eqref{continuity} over $\Omega$ and $t\in[0,\infty)$ and applying the divergence theorem that
\begin{equation}\label{eq1}
\int\limits_0^\infty dt \int\limits_{\partial \R} d^2\boldsymbol{x} \; \boldsymbol{n}(\boldsymbol{x}) \cdot {\boldsymbol{j}}^\psi(\boldsymbol{x},t)
= \int\limits_{\Omega} d^3\boldsymbol{x}\, |\psi_0(\boldsymbol{x})|^2 - \lim_{t\to\infty} \int\limits_{\Omega} d^3\boldsymbol{x}\, |\psi(\boldsymbol{x},t)|^2 \,,
\end{equation}
so, using that $\|\psi_0\|=1$,
\begin{equation}
\mu(Z=\infty) = \lim_{t\to\infty} \int_{\Omega} d^3\boldsymbol{x}\, |\psi(\boldsymbol{x},t)|^2 \geq 0\,.
\end{equation}
Thus, $\mu$ is a probability measure. It can also be expressed as
\begin{equation}
\mu(\cdot) = \scp{\psi_0}{E(\cdot)|\psi_0}
\end{equation}
with the POVM $E$ on $[0,\infty)\times \partial \R \cup \{\infty\}$ given by
\begin{align}
E(dt \times d^2\boldsymbol{x}) &=  W_t^\dagger |\boldsymbol{x}\rangle\langle\boldsymbol{x}| W_t \,dt\, d^2\boldsymbol{x} \label{Edef1}\\
E(\{\infty\})&= \lim_{t\to\infty} W_t^\dagger W_t\label{Edef2}
\end{align}
with $W_t$ the time evolution operator on $L^2(\Omega,\mathbb{C}^4)$,
\begin{equation}\label{Wdef}
\psi_t=W_t\psi_0\,,
\end{equation}
according to \eqref{Dirac} and \eqref{abcd1}. As explained in detail in \cite{detect-thm}, the operators $(W_t)_{t\geq 0}$ form a semigroup, $W_sW_t=W_{s+t}$ and $W_0=I$, and are in general not unitary but contractions, $\|W_t\psi\|\leq \|\psi\|$. In fact, $\|W_t\psi_0\|^2=\mu(T>t \text{ or }Z=\infty)$. The time evolution is of the form $W_t=\exp(-iHt/\hbar)$, but the Hamiltonian $H$, defined by \eqref{Dirac} and \eqref{abcd1}, is not self-adjoint.

\section{Rule for Curved Space-Time and Moving Detectors}
\label{sec:rulecurved}

The version formulated in this section is the natural generalization of the one in Section~\ref{sec:ruleflat} to curved space-time and an arbitrary space-time shape of the detecting surface.

\subsection{Setup}

Let $(\mathscr{M},g)$ be a Lorentzian manifold, suppose it is globally hyperbolic, and let $\mathscr{S}$ be a vector bundle of Dirac spin spaces over $(\mathscr{M},g)$ \cite{PR84}. Thus, for every $x\in \mathscr{M}$, $\mathscr{S}_x$ is a 4-dimensional complex vector space, the spin space at $x$, equipped with an indefinite inner product $\overline{\phi}\psi$ and Dirac matrices $\gamma^\mu$; moreover, the bundle is equipped with a covariant derivative $\nabla_{\!\mu}$; the inner product and the Dirac matrices are parallel relative to $\nabla_{\!\mu}$. With any spacelike hypersurface $\Sigma$ in $\mathscr{M}$ there is associated a Hilbert space $\mathscr{H}_\Sigma=L^2(\Sigma,\mathscr{S})$ of cross-sections $\psi:\Sigma\to\mathscr{S}$ with the positive definite inner product
\begin{equation}
\scp{\phi}{\psi}_{\Sigma} = \int_{\Sigma} d^3x \, \overline{\phi(x)}\, \gamma^\mu(x) \,n^\Sigma_\mu(x)\, \psi(x)\,,
\end{equation}
where $d^3x$ is the volume of a hypersurface element defined by the Riemannian 3-metric on $\Sigma$, and $n^\Sigma_\mu(x)$ is the future-pointing unit normal vector to $\Sigma$ at $x$. As usual, the norm is $\|\psi\|^2_\Sigma=\scp{\psi}{\psi}_\Sigma$. The Dirac equation in curved space-time reads
\begin{equation}\label{curvedDirac}
ic\hbar \gamma^\mu(x) \nabla_{\!\mu} \psi(x) = mc^2 \psi(x)
\end{equation}
and defines, for any two Cauchy hypersurfaces $\Sigma,\Sigma'$, a unitary isomorphism $\mathscr{H}_{\Sigma}\to \mathscr{H}_{\Sigma'}$. The probability current 4-vector field is, in analogy to \eqref{jdef}, given by
\begin{equation}\label{jdef2}
j^\mu = \overline{\psi} \, \gamma^\mu \, \psi
\end{equation}
and satisfies, by virtue of the Dirac equation \eqref{curvedDirac}, the continuity equation
\begin{equation}\label{continuity2}
\nabla_{\!\mu} j^\mu =0
\end{equation}
analogous to \eqref{continuity}.

Now let $\Sigma_0$ be a spacelike hypersurface in $\mathscr{M}$, and let $S$ be a subset of the future of $\Sigma_0$ with (piecewise smooth) boundary $\partial S$, see Figure~\ref{fig:S}. The detectors are placed along $\tilde{\partial} S=\partial S\setminus \Sigma_0$.
We are given an initial 1-particle wave function $\psi_0$, which is a cross-section of $\mathscr{S}$ defined on $\Sigma_0\cap \partial S$, and the rule defines, in terms of $\psi_0$, the probability distribution of $Z\in \tilde{\partial} S \cup\{\infty\}$, where $Z$ is either the space-time point of detection or $\infty$ if the particle is never detected.

\begin{figure}[h]
\begin{center}
\includegraphics[width=.4 \textwidth]{detect-dirac-fig1.pdf}
\end{center}
\caption{Example of a space-time region $S$ in the future of a spacelike hypersurface $\Sigma_0$, shown here in 1+1 dimensions. This example happens to have finite volume, although $S$ is not required to.}
\label{fig:S}
\end{figure}

Note that this setup includes, also in flat space-time, the possibility of a moving detecting surface, as opposed to the resting detecting surface (of the form $[0,\infty)\times \partial \R$ in some Lorentz frame) considered in Section~\ref{sec:ruleflat}. This possibility can be thought of either as moving detectors or as detectors that get switched on at a certain time; the latter scenario makes it clear that the detecting surface can very well have spacelike regions. The detecting surface $\tilde{\partial} S$ can be subdivided into the region $\partial_{\mathrm{ti}} S$ where it is timelike, the region $\partial_{\mathrm{psp}} S$ where it is spacelike or null and lies in the past of $S$, and the region $\partial_{\mathrm{fsp}} S$ where it is spacelike or null and lies in the future of $S$; see Firgure~\ref{fig:dS}. The particle may be detected in $\partial_{\mathrm{ti}} S$ or $\partial_{\mathrm{fsp}} S$, but it cannot reach $\partial_{\mathrm{psp}} S$ along a timelike curve from its starting region $\Sigma_0\cap \partial S$. (There may be parts of $\partial_{\mathrm{ti}} S$ and $\partial_{\mathrm{fsp}} S$ that cannot be reached either.)

\begin{figure}[h]
\begin{center}
\includegraphics[width=.4 \textwidth]{detect-dirac-fig2.pdf}
\end{center}
\caption{Parts of the boundary of $S$. Thick curves: $\partial_{\mathrm{ti}} S$; dashed: $\partial_{\mathrm{fsp}} S$; dotted: $\partial_{\mathrm{psp}} S$. Lightlike directions are drawn at $45^\circ$.
The dash--dotted line (which is at $45^\circ$) means that $\psi$ vanishes to the right of it because $\psi$
cannot reach that region within $S$ from $\Sigma_0$ propagating no faster than light (assuming that $\psi$ vanishes on $\partial_{\mathrm{psp}} S$).}
\label{fig:dS}
\end{figure}

\subsection{Detector Frame}

As mentioned already, the rule requires that we specify the rest frame of every detector in $\partial_{\mathrm{ti}} S$. This will be encoded in a future-pointing timelike unit vector field $u^\mu$ defined on and tangent to $\partial_{\mathrm{ti}} S$. Note that there are two degrees of freedom for the choice of $u^\mu$ for every $x\in\partial_{\mathrm{ti}} S$, as the tangent space $T_x\partial S$ is 3-dimensional and the future-pointing timelike unit vectors in there form a 2-dimensional set.

To get an understanding of why the choice of a vector $u^\mu$ is necessary, let us consider the setting of Section~\ref{sec:ruleflat}, in which there is a Lorentz frame in which the set $\Omega$ is time independent. Specifically, consider $\Omega=\{x^1<0\}\subset \mathbb{R}^3$, so the detecting surface is plane $\partial \R=\{x^1=0\}$, corresponding to the timelike hyperplane $\{x^1=0\}$ in space-time. Apart from the Lorentz frame $\Lambda$ in which Equations \eqref{Dirac}--\eqref{abcd1} hold, let us consider another Lorentz frame $\hat\Lambda$ with coordinates $\hat{x}^\mu$, arising from $\Lambda$ by a boost in the $x^2$ direction,
\begin{align}
\hat{x}^0 &= x^0\cosh \xi+x^2\sinh \xi  & \hat{x}^1&=x^1\nonumber\\
\hat{x}^2 &= x^0\sinh \xi +x^2\cosh \xi & \hat{x}^3 &=x^3\,.
\end{align}
Since the space-time set $\{x^1<0\}$ is the same as $\{\hat{x}^1<0\}$, the timelike hyperplane where the detectors are located, $\{\hat{x}^1=0\}$, looks static also in the $\hat{x}$ coordinates. However, the boundary condition \eqref{abcd1}, which holds in the $\Lambda$ coordinates, does not hold in the $\hat\Lambda$ coordinates. This can be seen, e.g., from the Bohmian trajectories: While \eqref{abcd1} entails that everywhere on the boundary ($x^1=0$), the coordinate velocity ${\boldsymbol{v}}$ defined by
\begin{equation}
v^i = \frac{dX^i}{dX^0}\,, \quad i=1,2,3,
\end{equation}
is pointing in the $x^1$ direction,
\begin{equation}
{\boldsymbol{v}}=(1,0,0)\,,
\end{equation}
the coordinate velocity in $\hat\Lambda$ must be
\begin{equation}
\hat{\boldsymbol{v}}=\biggl( \frac{1}{\cosh\xi}, \tanh\xi,0 \biggr)
\end{equation}
(as can be seen easily by transforming the lightlike wordline $X(s)=(s,s,0,0)$ to $\hat{X}(s) = (s\cosh\xi,s,s\sinh\xi,0)$). Since $\hat{\boldsymbol{v}}$ does not point in the $\hat{x}^1$ direction, \eqref{abcd1} cannot hold in $\hat\Lambda$. That is why we need to specify the Lorentz frame in which \eqref{abcd1} holds or, equivalently, in which the coordinate velocity on the boundary is perpendicular to the boundary.

So let a future-timelike tangent vector field $u^\mu$ on $\partial_{\mathrm{ti}} S$ be given.

\subsection{Statement of the Rule}

Let $\psi$ be the cross-section of $\mathscr{S}$ defined on $S$ that is the unique solution of the Dirac equation \eqref{curvedDirac}
with the initial condition
\begin{equation}\label{iccurvedD}
\psi(x) = \psi_0(x) \quad \text{for all }x\in\Sigma_0\cap \partial S\,,
\end{equation}
the boundary condition
\begin{equation}\label{abcd2}
n^{\partial S}_\mu(x)\, \gamma^\mu(x)\, \psi(x)
=u_\mu(x)\, \gamma^\mu(x)\, \psi(x)
\quad  \text{for all }x\in\partial_{\mathrm{ti}} S,
\end{equation}
where $n^{\partial S}(x)$ is the (spacelike) unit normal vector to $\partial S$ at $x\in\partial_{\mathrm{ti}} S$,
and the further condition
\begin{equation}\label{psp}
\psi(x)=0 \quad \text{for all }x\in \partial_{\mathrm{psp}} S\,.
\end{equation}
The last condition plays mathematically the role of a further part of the initial condition; for our purposes it states that the detecting surface does not emit particles, so the one particle we consider has to start in the region $\Sigma_0\cap \partial S$, not in $\partial_{\mathrm{psp}} S$.
As a shorthand for \eqref{abcd2}, we also write
\begin{equation}\label{abcd3}
n\!\!\!/ \, \psi = u\!\!\!/\, \psi \quad \text{on }\partial_{\mathrm{ti}} S,
\end{equation}
where $n$ is for $n^{\partial S}$, and the slash $/$ denotes the contraction with $\gamma^\mu$.

Let $Z$ be the space-time point on $\partial S$ where the particle gets detected, if it gets detected, and $Z=\infty$ it the particle never gets detected. We take for granted that $\psi_0$ vanishes outside $\partial S$ and $\|\psi_0\|_{\Sigma_0}=1$ and assume for simplicity that the lightlike portion of $\partial_{\mathrm{fsp}} S$ has dimension $<3$ (see Remarks~\ref{rem:curvedBohm} and \ref{rem:forms} below for a formulation that does not require this assumption).

Our rule asserts that $Z$ has probability density
\begin{equation}
n^{\partial S}_\mu j^\mu = \overline{\psi}\, n\!\!\!/\: \psi
\end{equation}
on $\tilde{\partial} S$ relative to the volume measure on $\tilde{\partial} S$ defined by the 3-metric on $\tilde{\partial} S$. That is, the distribution $\mu$ of $Z$ is given by
\begin{align}
\mu(Z\in d^3x) &= \overline{\psi(x)}\, n\!\!\!/(x)\, \psi(x)\, d^3x \label{curvedmudef1}\\
\mu(Z=\infty) &= 1-\int_{\tilde{\partial} S} d^3x\, \overline{\psi(x)}\, n\!\!\!/(x)\, \psi(x)\,.\label{curvedmudef2}
\end{align}
This completes the statement of the rule. Its relativistic invariance is manifest.

\subsection{Remarks}
\label{sec:rem}

\begin{enumerate}

\item \label{rem:curvedBohm} Equivalently, the rule can be expressed in terms of the Bohmian trajectory $X^\mu(s)$, guided by $\psi$ evolving according to \eqref{curvedDirac}, \eqref{abcd2}, and \eqref{psp}, and starting at a random point in $\Sigma_0\cap \partial S$ with probability density given by $|\psi_0|^2$ (i.e., by $\overline{\psi_0} \, n\!\!\!/\,^{\Sigma_0} \, \psi_0$). Namely, $Z$ is the random point at which the Bohmian trajectory hits $\tilde{\partial} S$ if it hits $\tilde{\partial} S$, and $Z=\infty$ otherwise. This formulation also applies when the lightlike portion of $\partial_{\mathrm{fsp}} S$ has dimension 3, a case in which \eqref{curvedmudef1} does not apply because $d^3x=0$.

\item It seems very plausible that the Dirac equation \eqref{curvedDirac} on $S$ together with initial and boundary conditions \eqref{iccurvedD}--\eqref{psp} possesses a unique solution for every $\psi_0\in L^2(\Sigma_0\cap \partial S,\mathscr{S})$; it would be of interest to have a rigorous proof. It then follows from \eqref{continuity2} that for every spacelike hypersurface $\Sigma$ in the future of $\Sigma_0$, the restriction $\psi_{\Sigma}$ of $\psi$ to $\Sigma$ is an element of $L^2(\Sigma\cap S,\mathscr{S})$ whose norm squared $\|\psi_\Sigma\|^2_\Sigma$
equals the probability that either $Z=\infty$ or $Z$ lies in the future of $\Sigma$. The time evolution operator $W_{\Sigma_1}^{\Sigma_2}$ defined for $\Sigma_2$ in the future of $\Sigma_1$ (in the future of $\Sigma_0$, where the ``future of $\Sigma$'' includes $\Sigma$ itself) by $\psi_{\Sigma_2}=W_{\Sigma_1}^{\Sigma_2}\psi_{\Sigma_1}$ is a contraction $\mathscr{H}_{\Sigma_1}\to\mathscr{H}_{\Sigma_2}$.

\item External field. In the presence of an external electromagnetic field, the $\nabla_{\!\mu}$ in the Dirac equation \eqref{curvedDirac} must be replaced by $\nabla_{\!\mu}+i\tfrac{e}{c}A_\mu(x)$, where $e$ is the particle's charge. The boundary condition does not change.

\item Let us take a closer look at the boundary condition \eqref{abcd2}. For every $x\in \partial_{\mathrm{ti}} S$, there is a local Lorentz frame (i.e., an orthonormal basis of $T_x\mathscr{M}$) for which $u=(1,0,0,0)$ and $n=(0,1,0,0)$; that is because $u^\mu n_\mu=0$ as a consequence of the fact that $u$ is tangent to $\partial S$ and $n$ is normal to it. In this frame, the boundary condition becomes $\gamma^1\psi=\gamma^0\psi$ or, what is equivalent by virtue of the relation $\alpha_1=(\gamma^0)^{-1}\gamma^1$,
\begin{equation}
\alpha_1 \psi=\psi\,,
\end{equation}
which is the ABCD \eqref{abcd1} considered in Section~\ref{sec:ruleflat}. It follows that the system of equations \eqref{Dirac}--\eqref{abcd1} in Section~\ref{sec:ruleflat} is a special case of the system \eqref{curvedDirac}, \eqref{iccurvedD}--\eqref{psp} in this section.

\item\label{rem:L} Subbundle. It also follows that the ABCD $n\!\!\!/\psi=u\!\!\!/\psi$ is equivalent to requiring that
\begin{equation}
\psi(x) \in \mathscr{L}_x
\end{equation}
for all $x\in\partial_{\mathrm{ti}} S$, where $\mathscr{L}_x$ is a subspace (of complex dimension 2) of the spin space $\mathscr{S}_x$, namely the eigenspace with eigenvalue $+1$ of the matrix $\alpha_1$ in the local Lorentz frame with $u=(1,0,0,0)$ and $n=(0,1,0,0)$. Put differently, $\mathscr{L}_x$ is the kernel (nullspace) of the operator $n\!\!\!/ - u\!\!\!/:\mathscr{S}_x \to \mathscr{S}_x$. Together, the $\mathscr{L}_x$ form a rank-2 subbundle $\mathscr{L}$ of $\mathscr{S}$ over $\partial_{\mathrm{ti}} S$.

\item \label{rem:forms} Differential forms. The distribution of $Z$ over $\tilde{\partial} S$ can also be expressed as a differential 3-form $\omega$ on $\tilde{\partial} S$. Two advantages of the formalism of differential forms in this context are its independence of any additional structure of space-time such as a Lorentzian metric or Galilean structure and, relatedly, that it has no difficulties with lightlike surfaces in Lorentzian manifolds, for which the concept of 3-volume becomes degenerate.

To say that $\omega$ expresses the probability distribution $\mu$ of $Z$ means that
\begin{equation}\label{muomega}
\mu(Z\in B)= \int_B \omega
\end{equation}
for any subset $B\subseteq\tilde{\partial} S$ (assumed to have outward orientation).
In our case, the differential form $\omega$ is given by
\begin{equation}\label{omegaj}
\omega_{\lambda\mu\nu} = j^\sigma \varepsilon_{\sigma\lambda\mu\nu}
\end{equation}
or $\omega = \iota_j \varepsilon$, where $\varepsilon$ is the 4-form that represents space-time volume (given by the Levi-Civita symbol in any Galilean/Lorentzian frame, or $\varepsilon = dt\wedge dx^1 \wedge dx^2 \wedge dx^3$), and $\iota_j$ means inserting the 4-vector field $j$ into the first slot of a differential form. Equation \eqref{muomega} with \eqref{omegaj} is a reformulation of \eqref{curvedmudef1} that remains valid when the lightlike portion of $\partial_{\mathrm{fsp}} S$ has dimension 3.

In the non-relativistic (Galilean) case \cite{detect-several}, with $\Sigma_0$ a $t=\mathrm{const.}$ hypersurface and possibly moving detectors, \eqref{muomega} and \eqref{omegaj} are still valid with $j=(|\psi|^2, {\boldsymbol{j}}^{\psi})$ and ${\boldsymbol{j}}^\psi = (\hbar/m) \Im(\psi^* \nabla \psi)$; that is, \eqref{muomega} agrees with Equation (13) of \cite{detect-several}.

In general, for Bohm-like world lines that are integral curves of a current 4-vector field $j$, the integral $\int_B\omega$ with $\omega$ as in \eqref{omegaj} and $B$ an oriented piece of hypersurface has the following meaning \cite{Tum01}. The integral equals the expected number of signed crossings of the random world line through $B$, where a crossing against the orientation of $B$ is counted negative. In the present case, if $B$ is any piece of $\tilde{\partial} S$, the particle world line can cross $B$ only outward and at most once, so $\int_B \omega$ equals the probability of crossing $B$, in agreement with \eqref{muomega}.

\item\label{rem:muprob} We now verify that $\mu$ as defined in \eqref{curvedmudef1} and \eqref{curvedmudef2} is a probability measure. Since $\mu(\tilde{\partial} S \cup \{\infty\})=1$ by construction, it only remains to verify that $\mu$ is non-negative. The density of $\mu$ relative to $d^3x$, $\overline{\psi} \,n\!\!\!/\,\psi$, is non-negative because it is equal, by virtue of $n\!\!\!/\;\psi=u\!\!\!/\:\psi$, to $\overline{\psi} \,u\!\!\!/\,\psi$, and this is equal to $|\psi|^2$ in any Lorentz frame with $u=(1,0,0,0)$.

We now show that $\mu(\{\infty\})\geq 0$, or, equivalently, $\mu(\tilde{\partial} S)\leq 1$. Consider first the case that $S$ is compact (as in Figure~\ref{fig:S}). We integrate the equation $\nabla_{\!\mu}j^\mu=0$ over $S$. The divergence theorem in Lorentzian manifolds yields that
\begin{equation}\label{divthm}
\int_S d^4x\, \nabla_{\!\mu}j^\mu = \int_{\partial S} d^3x\, n^{\partial S}_\mu \, j^\mu
\end{equation}
with outward orientation (in deviation from our previous convention that $n^{\partial S}$ points to the future on spacelike parts of $\partial S$).\footnote{Equivalently, this step can be expressed (and may be more transparent) in terms of differential forms, noting that $(\nabla_{\!\mu}j^\mu)\varepsilon=d\omega$ with $d$ the exterior derivative (because $\nabla_{\!\rho}\varepsilon_{\sigma\lambda\mu\nu}=0$). The left-hand side of \eqref{divthm} equals $\int_S d\omega$, which by Stokes' theorem for differential forms equals $\int_{\partial S} \omega$, which equals the right-hand side of \eqref{divthm}.} The left-hand side vanishes, and the right-hand side equals (paying attention to the orientation)
\begin{equation}
-\|\psi_0\|^2+\int_{\tilde{\partial} S}d^3x\, n^{\partial S}_\mu \, j^\mu\,.
\end{equation}
Since $j$ vanishes on $\partial_{\mathrm{psp}} S$ by \eqref{psp}, we obtain that
\begin{equation}\label{intboutiboufsp}
\int_{\partial_{\mathrm{ti}} S\cup\partial_{\mathrm{fsp}} S}d^3x\, n^{\partial S}_\mu \, j^\mu =1\,.
\end{equation}
Since now all remaining spacelike portions are future boundaries of $S$, the outward orientation agrees with the future orientation, so we can return to our previous convention of taking $n^{\partial S}$ future-pointing when timelike (or lightlike), and the derivation of $\mu(\tilde{\partial} S)\leq 1$ is complete. Its version with differential forms also applies when the lightlike parts of $\partial S$ have dimension 3.

Now suppose that $S$ is not compact; for example, we may think of Minkowski space-time with $\Omega$ a 3-dimensional half space in some Lorentz frame. Then the previous reasoning still applies if we interpret $\partial S$ as including ``boundaries at infinity,'' such as parts of $\mathscr{I}^+$. (Alternatively, we may consider a limit in which a sequence of compact $S$'s approaches the desired non-compact $S$.) Since we defined $\partial_{\mathrm{ti}} S\cup\partial_{\mathrm{fsp}} S$ to contain only finite space-time points, not ideal points at infinity, some parts of the full future boundary are missing in \eqref{intboutiboufsp}, with the consequence that ``$=1$'' must be replaced by ``$\leq 1$'' in \eqref{intboutiboufsp}.

\item\label{rem:POVM} POVM. As in the flat case, the distribution of $Z$ is given by a POVM $E(\cdot)$ on $\tilde{\partial} S$ acting on $\mathscr{H}_{\Sigma_0}$,
\begin{equation}
\mu(\cdot) = \scp{\psi_0}{E(\cdot)|\psi_0}\,.
\end{equation}
A perhaps useful explicit specification of this POVM can be given in terms of the operators $J_{\mathrm{ti}}$ and $J_{\mathrm{fsp}}$ that map $\psi_0$ to the restriction of $\psi$ to $\partial_{\mathrm{ti}} S$ and $\partial_{\mathrm{fsp}} S$, respectively. Although $\partial_{\mathrm{ti}} S$ is not spacelike, it is naturally associated with a Hilbert space $\mathscr{H}_{\mathrm{ti}}=\mathscr{H}_{\partial_{\mathrm{ti}} S,u}$ (= ``$L^2(\partial_{\mathrm{ti}} S,\mathscr{L})$''), consisting of cross-sections $\phi:\partial_{\mathrm{ti}} S \to \mathscr{L}$ (with $\mathscr{L}$ as defined in Remark~\ref{rem:L}) and equipped with the inner product
\begin{equation}
\scp{\phi}{\chi}_{\partial_{\mathrm{ti}} S} = \int_{\partial_{\mathrm{ti}} S} d^3x \, \overline{\phi(x)} \, n\!\!\!/(x) \, \chi(x)\,.
\end{equation}
While this inner product would not be positive definite for cross-section of $\mathscr{S}$, it is for cross-sections of $\mathscr{L}$ because for them $n\!\!\!/\, \chi = u\!\!\!/\, \chi$, and $\overline{\chi} \, u\!\!\!/\, \chi\geq 0$ as pointed out already at the beginning of Remark~\ref{rem:muprob}. Both $\mathscr{H}_{\mathrm{ti}}$ and $\mathscr{H}_{\mathrm{fsp}}=\mathscr{H}_{\partial_{\mathrm{fsp}} S}$, being spaces of functions, are equipped with natural PVMs $P_{\mathrm{ti}}$ and $P_{\mathrm{fsp}}$, where $P_{\#}(B)$ is the multiplication with the characteristic function of $B$. Both $J_{\mathrm{ti}}:\mathscr{H}_{\Sigma_0}\to \mathscr{H}_{\mathrm{ti}}$ and $J_{\mathrm{fsp}}:\mathscr{H}_{\Sigma_0}\to\mathscr{H}_{\mathrm{fsp}}$ are contractions, so $J:= J_{\mathrm{ti}}\oplus J_{\mathrm{fsp}}$ is a contraction $\mathscr{H}_{\Sigma_0}\to\mathscr{H}_{\mathrm{ti}}\oplus\mathscr{H}_{\mathrm{fsp}}$. We find that
\begin{equation}
E(B) = J^*\Bigl[P_{\mathrm{ti}}(B\cap \partial_{\mathrm{ti}} S) \oplus P_{\mathrm{fsp}}(B\cap \partial_{\mathrm{fsp}} S)\Bigr] J
\end{equation}
for $B\subseteq \partial_{\mathrm{ti}} S\cup \partial_{\mathrm{fsp}} S$, and we can see that this is a positive operator. Since $E$ vanishes on $\partial_{\mathrm{psp}} S$, it only remains to specify $E(\{\infty\})$, which is
\begin{equation}
E(\{\infty\}) = I-J^*J\,.
\end{equation}

\item Conditional Distribution.
Given that no detection has occurred between $\Sigma_0$ and another spacelike hypersurface $\Sigma_1$ in the future of $\Sigma_0$, the conditional distribution of the detection space-time point $Z$ is $\scp{\tilde\psi_1}{E_1(\cdot)|\tilde\psi_1}_{\Sigma_1}$, where $\tilde\psi_1=\psi_1/\|\psi_1\|$ with $\psi_1=W_{\Sigma_0}^{\Sigma_1}\psi_0$, and $E_1(\cdot)$ is the POVM constructed from $\Sigma_1$ and $S\cap \mathrm{future}(\Sigma_1)$ in the same way as $E(\cdot)$ is constructed from $\Sigma_0$ and $S$. This follows from two facts: First, that the conditional distribution of a random variable $Z$ in a probability space $(X,\mu)$, given that $Z\in C\subset X$, is $\mu(Z\in B|Z\in C) = \mu(B\cap C)/\mu(C)$. And second, that $\psi_0$ and $\psi_1$ actually lead to the same wave function $\psi$, up to a global factor $\|\psi_1\|$, on $S\cap \mathrm{future}(\Sigma_1)$.
\end{enumerate}

\section{Semi-Ideal Detectors}
\label{sec:semi}

The Dirac equation with the ideal ABCD \eqref{abcd1} or, more generally, \eqref{abcd2}, does not possess a non-relativistic limit, as is obvious already from the fact that the ideal ABCD forces the particle to cross the detecting surface $\partial \R$ at the speed of light. However, \eqref{abcd1} is not the only absorbing boundary condition for the Dirac equation, and another type of these conditions, which we call the \emph{semi-ideal ABCDs}, do possess a non-relativistic limit. They are of the form
\begin{equation}\label{abcd4}
\bigl(\boldsymbol{n}(\boldsymbol{x})\cdot {\boldsymbol{\alpha}} + \theta \beta\bigr) \psi = \sqrt{1+\theta^2}\,\psi
\end{equation}
with $\theta\in \mathbb{R}$ a constant, or, in a more general formulation that also applies to curved space-time and moving detectors,
\begin{equation}\label{abcd5}
(n\!\!\!/ + \theta) \psi = \sqrt{1+\theta^2} \, u\!\!\!/\, \psi\,.
\end{equation}
Obviously, the ideal ABCD \eqref{abcd1} or \eqref{abcd2} is included for $\theta=0$.

\subsection{Discussion of the Semi-Ideal ABCD}

It is well known that the eigenvalues of the matrix ${\boldsymbol{v}}\cdot {\boldsymbol{\alpha}}+ \theta \beta$ (for arbitrary ${\boldsymbol{v}}\in\mathbb{R}^3$ and $\theta\in\mathbb{R}$) are $\pm\sqrt{|{\boldsymbol{v}}|^2+\theta^2}$, each with multiplicity two. This is often formulated as the fact that the free Dirac Hamiltonian in Fourier form, which is multiplication by the matrix
\begin{equation}
H({\boldsymbol{k}}) = c\hbar{\boldsymbol{k}} \cdot {\boldsymbol{\alpha}} + mc^2 \beta\,,
\end{equation}
has eigenvalues
\begin{equation}
E({\boldsymbol{k}})=\pm \sqrt{c^2\hbar^2|{\boldsymbol{k}}|^2 + m^2c^4}\,,
\end{equation}
each with multiplicity two. In our case with ${\boldsymbol{v}}=\boldsymbol{n}(\boldsymbol{x})$, we obtain that $\boldsymbol{n}(\boldsymbol{x})\cdot {\boldsymbol{\alpha}} +\theta\beta$ has eigenvalues $\pm \sqrt{1+\theta^2}$, and the ABCD \eqref{abcd4} amounts to saying that $\psi(\boldsymbol{x})$ for $\boldsymbol{x}\in\partial \R$ has to lie in the 2-dimensional eigenspace $\mathscr{L}_{\boldsymbol{x}}$ with the positive eigenvalue.

In order to show that \eqref{abcd4} implies that the current points outward at every point of the boundary $\partial \R$, we use the Dirac basis in spin space (a.k.a.\ standard representation), in which
\begin{equation}
\gamma^0=\beta = \begin{pmatrix} I_2&0\\0&-I_2\end{pmatrix}\,,
\quad \gamma^i=\gamma^0\alpha^i = \begin{pmatrix} 0&\sigma_i\\ -\sigma_i &0\end{pmatrix}\,,
\quad \alpha^i = \begin{pmatrix} 0&\sigma_i\\ \sigma_i&0 \end{pmatrix}
\end{equation}
with $I_2$ the $2\times 2$ identity matrix, and write
\begin{equation}
\psi = \begin{pmatrix} \psi_1 \\ \psi_2 \end{pmatrix}
\end{equation}
with 2-spinors $\psi_1,\psi_2$. We obtain that, at any $\boldsymbol{x}\in\partial \R$,
\begin{align}
\boldsymbol{n} \cdot {\boldsymbol{j}}
&= \psi^\dagger (\boldsymbol{n}\cdot{\boldsymbol{\alpha}}) \psi\\
& =  \psi^\dagger (-\theta\beta + \sqrt{1+\theta^2}) \psi\\
&= \Bigl(\sqrt{1+\theta^2}-\theta\Bigr) \psi_1^\dagger \psi_1
+ \Bigl(\sqrt{1+\theta^2}+\theta \Bigr)\psi_2^\dagger \psi_2 \geq 0
\end{align}
since $\sqrt{1+\theta^2}\geq |\theta|$. This is what we claimed.

We note further that the current ${\boldsymbol{j}}$ is, in fact, everywhere orthogonal to the boundary. To see this, one can derive that eigenvectors $\psi$ of $\boldsymbol{n}\cdot {\boldsymbol{\alpha}}+\theta\beta$ with eigenvalue $\sqrt{1+\theta^2}$ satisfy
\begin{equation}\label{abcd6}
\psi_2 = \Bigl(\sqrt{1+\theta^2}-\theta\Bigr)(\boldsymbol{n}\cdot{\boldsymbol{\sigma}})\psi_1\,.
\end{equation}
Thus, for any vector ${\boldsymbol{v}}$ orthogonal to $\boldsymbol{n}$,
\begin{align}
{\boldsymbol{v}} \cdot {\boldsymbol{j}}
&= \psi^\dagger ({\boldsymbol{v}}\cdot{\boldsymbol{\alpha}}) \psi\\
&= \psi_1^\dagger ({\boldsymbol{v}}\cdot {\boldsymbol{\sigma}})\psi_2 + \psi_2^\dagger ({\boldsymbol{v}}\cdot {\boldsymbol{\sigma}}) \psi_1\\
&= \Bigl(\sqrt{1+\theta^2}-\theta\Bigr) \psi_1^\dagger\Bigl[ ({\boldsymbol{v}}\cdot {\boldsymbol{\sigma}})(\boldsymbol{n}\cdot{\boldsymbol{\sigma}}) +  (\boldsymbol{n}\cdot{\boldsymbol{\sigma}})({\boldsymbol{v}}\cdot {\boldsymbol{\sigma}}) \Bigr] \psi_1 \\
&=0
\end{align}
because, as is well known,
\begin{equation}
(\boldsymbol{a}\cdot {\boldsymbol{\sigma}})(\boldsymbol{b}\cdot{\boldsymbol{\sigma}}) = (\boldsymbol{a}\cdot \boldsymbol{b}) I_2 + i (\boldsymbol{a} \times \boldsymbol{b})\cdot {\boldsymbol{\sigma}}
\end{equation}
for any $\boldsymbol{a},\boldsymbol{b}\in\mathbb{R}^3$ and $\times$ the cross product in $\mathbb{R}^3$.

Also with the semi-ideal ABCD \eqref{abcd4}, the Dirac Hamiltonian generates, for a time-independent region $\Omega$, a contraction semigroup $W_t=e^{-iHt/\hbar}$ with non-self-adjoint $H$ \cite{detect-thm}.

\subsection{Non-Relativistic Limit}

In the non-relativistic limit $c\to\infty$ of the Dirac equation, one considers wave functions superposed of energy eigenstates with energies infinitesimally above $mc^2$. Such wave functions obey, in the Dirac basis again, \cite{Dirac}
\begin{equation}\label{psi21}
\psi_2 \approx -\frac{i\hbar}{2mc} {\boldsymbol{\sigma}}\cdot \nabla \psi_1\,,
\end{equation}
so in particular (since $c$ is large) $\psi_2\ll \psi_1$. In the limit, the 2-spinor wave function
\begin{equation}
\phi=\psi_1
\end{equation}
satisfies the Pauli equation, which coincides with the usual Schr\"odinger equation in the absence of external fields.
The current has $j^0=|\phi|^2$ and spacelike components \cite[p.~216]{BH93}
\begin{equation}\label{jPauli1}
{\boldsymbol{j}} = \frac{\hbar}{m} \Im[\phi^\dagger\nabla \phi] + \frac{\hbar}{2m} \nabla\times [\phi^\dagger{\boldsymbol{\sigma}} \phi]\,,
\end{equation}
or, equivalently,
\begin{equation}\label{jPauli2}
{\boldsymbol{j}} = -i\frac{\hbar}{2m} \Bigl[ \phi^\dagger {\boldsymbol{\sigma}}({\boldsymbol{\sigma}}\cdot \nabla)\phi - ({\boldsymbol{\sigma}}\cdot \nabla \phi)^\dagger {\boldsymbol{\sigma}} \phi\Bigr]\,.
\end{equation}

The semi-ideal ABCD \eqref{abcd4}, or equivalently \eqref{abcd6}, becomes, if we express $\psi_2$ using \eqref{psi21} and write $\phi$ for $\psi_1$, the non-relativistic ABC \eqref{ABC2}, which we repeat here for convenience:
\begin{equation}\label{ABC3}
{\boldsymbol{\sigma}}\cdot \nabla \phi(\boldsymbol{x}) = i\kappa \,\boldsymbol{n}(\boldsymbol{x})\cdot{\boldsymbol{\sigma}}\, \phi(\boldsymbol{x})\,,
\end{equation}
with
\begin{equation}\label{kappatheta}
\kappa = \frac{2mc}{\hbar} \Bigl( \sqrt{1+\theta^2}-\theta \Bigr)\,.
\end{equation}
Since the inverse relation to \eqref{kappatheta} is
\begin{equation}\label{thetakappa}
\theta = \frac{r^{-1}-r}{2} \quad \text{with } r= \frac{\hbar\kappa}{2mc}\,,
\end{equation}
a non-relativistic limit of the ABCD \eqref{abcd4} exists if we let $\theta\to \infty$ as $c\to\infty$ in such a way that $0<\lim(\theta/c)<\infty$; then,
\begin{equation}
\lim \frac{\theta}{c} = \frac{m}{\hbar\kappa}\,.
\end{equation}

We note that also \eqref{ABC3} leads to an outward-pointing current, as inserting it into \eqref{jPauli2}, contracting with $\boldsymbol{n}$, and using $(\boldsymbol{n}\cdot{\boldsymbol{\sigma}})^2=I_2$ yields
\b
\boldsymbol{n}\cdot{\boldsymbol{j}}
= \frac{\hbar\kappa}{m}  \phi^\dagger \phi\geq 0\,.
\e
We note further that also \eqref{ABC3} leads to a contraction semigroup.

\section{Rule for Several Particles}
\label{sec:several}

In this section, we outline the combination of the version of the absorbing boundary rule developed here for the Dirac equation with that developed in \cite{detect-several} for $n>1$ particles. For simplicity, we suppose that all particles have the same detecting surface $\tilde{\partial} S$. For now, we also suppose that they do not interact (and will make a remark about the interacting case later).

\subsection{Statement of the Rule}
\label{sec:rule-several}

Let $\mathscr{F}$ be any spacelike foliation of the space-time $(\mathscr{M},g)$ containing the initial hypersurface $\Sigma_0$, and choose a coordinate system in which $\{x^0=\mathrm{const.}\}$ are the leaves of $\mathscr{F}$, which we will denote $\Sigma_t$ with $t=x^0$,and the initial hypersurface is $\{x^0=0\}$. For $t\geq0$, the wave function $\psi_t(x_1,\ldots,x_n)\in \mathscr{S}_{x_1}\otimes \cdots \otimes \mathscr{S}_{x_n}$ is defined on the set $(\Sigma_t\cap S)^n$. Note that the union $\cup_{t\geq0} (\Sigma_t\cap S)^n$ forms the set $\mathscr{C}\subset S^n$ of configurations that are simultaneous relative to $\mathscr{F}$, and all $\psi_t$ together define a function $\psi$ on $\mathscr{C}$.

The Dirac equation for a single particle naturally defines a time evolution for the $n$-particle wave function $\psi$, which can be formulated as follows. Consider the $n$ multi-time Dirac equations \cite{Blo34,PT13a}
\begin{equation}\label{multiDirac}
ic\hbar\gamma^\mu_j(x_j) \nabla_{\!j\mu} \phi(x_1,\ldots,x_n) = m_j c^2 \phi(x_1,\ldots,x_n)
\end{equation}
for $j=1,\ldots,n$ and functions $\phi$ on $S^n$ with $\phi(x_1,\ldots,x_n)\in\mathscr{S}_{x_1}\otimes \cdots\otimes \mathscr{S}_{x_n}$, where $\gamma^\mu_j(x_j)$ acts on $\mathscr{S}_{x_j}$ and $\nabla_{\!j\mu}$ denotes the (covariant) derivative relative to $x^\mu_j$. They are supplemented by the boundary conditions
\begin{equation}\label{abcd7}
n\!\!\!/(x_j)\, \phi(x_1,\ldots,x_n) = u\!\!\!/(x_j) \,\phi(x_1,\ldots,x_n)
\quad\text{when }x_j \in \partial_{\mathrm{ti}} S
\end{equation}
and quasi-initial conditions
\begin{equation}\label{psp2}
\phi(x_1,\ldots,x_n) = 0
\quad \text{when }x_j \in \partial_{\mathrm{psp}} S\,.
\end{equation}
As explained in \cite{PT13a}, multi-time equations on $S^n$ define a time evolution on the set $\mathscr{C}$ of simultaneous configurations; this is the time evolution of $\psi$. (We will come back to $\phi$ later.)
That is, $\psi$ satisfies \eqref{multiDirac}-summed-over-$j$,  \eqref{abcd7}, and \eqref{psp2}, and is uniquely determined by these equations from initial data $\psi_0$ on $(\Sigma_0\cap S)^n$. As before, we take for granted that $\|\psi_0\|=1$.

Let $T^1$ be the first time (according to $\mathscr{F}$) at which a particle is registered by a detector, $Z^1\in \Sigma_{T^1}\cap \tilde{\partial} S$ the space-time point where it got registered, and $I^1\in \{1,\ldots,n\}$ the label of the registered particle. Our proposed rule asserts, in analogy to that in \cite{detect-several}, that their joint probability distribution is
\begin{equation}\label{probnpsi}
\mathrm{Prob}\bigl( I^1=j, Z^1\in d^3x_j \bigr) = d^3x_j \!\!\!\!\! \int\limits_{(\Sigma_{x_j^0}\cap S )^{n-1}} \!\!\!\!\!
d^3x_1 \cdots \widehat{d^3x_j}\cdots d^3x_n\: \overline{\psi} \, n\!\!\!/(x_j)\, \psi
\end{equation}
where $d^3x_j$ is a hypersurface element of $\tilde{\partial} S$ and the hat $\widehat{\quad}$ denotes omission. In the event of a detection of particle $j$ at time $t$, the wave function collapses according to
\begin{equation}\label{collapse1}
\psi_t'(x_1,\ldots,x_n) = \mathscr{N} \, \psi_t(x_1,\ldots,x_{j-1},Z^1,x_{j+1},\ldots,x_n) \, \delta^3(x_j-Z^1)
\end{equation}
with $\mathscr{N}$ a normalization factor and $\delta^3$ the 3-dimensional Dirac delta function (relative to the 3-metric on $\Sigma_t$). If we remove particle $j$ from consideration after detection, we may equivalently proceed with the following collapsed wave function of $n-1$ particles:
\begin{equation}\label{collapse2}
\psi_t'(x_1,\ldots,\widehat{x_j}, \ldots, x_n) = \mathscr{N}\, \psi_t(x_1,\ldots,x_{j-1},Z^1,x_{j+1},\ldots,x_n)\,.
\end{equation}
Now proceed in the time evolution along $\mathscr{F}$ as above but with $n-1$ particles and starting out on $\Sigma_t$ with $\psi_t'$. Note that $\psi'$, while it has only $n-1$ space-time arguments, still has $n$ spin indices, i.e., $4^n$ components, or values in $\mathscr{S}_{x_1}\otimes \cdots \otimes \mathscr{S}_{x_n}$ including $\mathscr{S}_{x_j}$. However, since the future time evolution does not act on $\mathscr{S}_{x_j}$, the equations given above still equally apply, even if $\psi'$ carries a further index. This completes the statement of the rule.

Sorting the detection events by particle label, we write $Z_j$ (with lower index, where our use of upper and lower indices is unrelated to that in relativity for contravariant and covariant components of a 4-vector, as this index does not label space-time components) for $Z^k$ if $I^k=j$ (and $Z_j=\infty$ if all $I^k\neq j$); that is, $Z_j$ is the space-time point where particle $j$ was detected. The outcome of the experiment is thus the $n$-tuple $(Z_1,\ldots,Z_n) \in (\tilde{\partial} S\cup\{\infty\})^n$. Importantly, as we will show below, its probability distribution is independent of the choice of the foliation $\mathscr{F}$.

\subsection{Discussion}

We begin with the brief remark that the Bohmian trajectories that occur in this scenario are those for which $\mathscr{F}$ is the preferred foliation of space-time (``time foliation''): First, already the Bohmian motion of a free $n$-particle system depends on $\mathscr{F}$ \cite{HBD}, and second, the collapses must be taken instantaneous relative to $\mathscr{F}$ because they represent the decoherence due to the entanglement with the detectors, which the Bohmian particles at spacelike separation from a detection event $Z_j$ will feel instantaneously relative to $\mathscr{F}$. Note that the trajectories depend on $\mathscr{F}$ while the distribution of the $Z_j$ does not.

Now let us leave aside the Bohmian trajectories and turn to the multi-time wave function $\phi$. As discussed in \cite{Blo34,PT13a}, the multi-time equations \eqref{multiDirac}, here with \eqref{abcd7} and \eqref{psp2}, actually possess a unique solution $\phi$ on $S^n$. While multi-time equations can be inconsistent, the equations here are easily seen to be consistent because of the absence of interaction; and while multi-time wave functions are usually defined only on spacelike configurations, in this case it is defined on all of $S^n$, again because of absence of interaction. If desired, the solution can be expressed in coordinates as
\begin{equation}\label{phiW}
\phi(x_1,\ldots,x_n) = (W_{\Sigma_0}^{\Sigma_{x_1^0}}\otimes \cdots \otimes W_{\Sigma_0}^{\Sigma_{x_n^0}} \psi_0)(\boldsymbol{x}_1,\ldots,\boldsymbol{x}_n)\,.
\end{equation}
However, both the equations \eqref{multiDirac}--\eqref{psp2} and the solution $\phi$ are independent of any choice of coordinates or foliation. We note that the $\psi_t$ of Section~\ref{sec:rule-several} is just the restriction of $\phi$ to $\Sigma_t^n$, and that the post-collapse wave function $\psi'$ of $n-1$ particles, evolved to later times from its initial value given in \eqref{collapse2}, is also a restriction of $\phi$, except for the normalization factor $\mathscr{N}$. Namely, plug $Z^1$ into $x_j$ and let the other $x_k$ lie in the future of $\Sigma_t$. It follows from \eqref{probnpsi} that
\begin{equation}\label{probnphi}
\mathrm{Prob}\Bigl(Z_1\in d^3x_1,\ldots,Z_n\in d^3x_n\Bigr) = \overline{\phi} \: n\!\!\!/(x_1)\cdots n\!\!\!/(x_n)\: \phi \: d^3x_1 \cdots d^3x_n
\end{equation}
with $\phi=\phi(x_1,\ldots,x_n)$. That is, the use of multi-time wave functions allows us to compute the probability distribution in a direct way without recourse to collapse.

From here, it also follows easily that the joint distribution of $Z_1,\ldots,Z_n$ is given by a POVM $E_n(\cdot)$ on $(\tilde{\partial} S\cup\{\infty\})^n$ acting on $\mathscr{H}_{\Sigma_0}^{\otimes n}$,
\begin{equation}\label{probnE}
\mathrm{Prob}\Bigl( (Z_1,\ldots,Z_n)\in B \Bigr) =\scp{\psi_0}{E_n(B)|\psi_0}\,,
\end{equation}
where $E_n(\cdot)$ is the product POVM on $n$ copies of the 1-particle POVM $E(\cdot)$ introduced in Remark~\ref{rem:POVM} in Section~\ref{sec:rem} (and in \eqref{Edef1}, \eqref{Edef2} for the flat and stationary case). The defining property of a product POVM is
\begin{equation}\label{EnE}
E_n(B_1\times\cdots\times B_n) = E(B_1)\otimes \cdots \otimes E(B_n)\,.
\end{equation}
(The existence and uniqueness of the product POVM follows from Corollary 7 in Section~4.4 of \cite{DGTZ05}.) Note that while $E_n$ is a product measure, the probability distribution \eqref{probnE} is not unless $\psi_0$ is a product function; in other words, $Z_1,\ldots,Z_n$ are not independent unless the $n$ particles are disentangled.

We now turn to the case with interaction. It is not easy to write down an explicit example of such a dynamics because because instantaneous interaction by means of potentials would break the relativistic invariance, so the interaction should be implemented by the creation and annihilation of other particles, which lies beyond the scope of this paper. However, we make a few remarks on what we expect to happen in this case. While the joint distribution of $Z_1,\ldots,Z_n$ will still be given by a POVM $E_n(\cdot)$ as in \eqref{probnE}, the POVM will no longer be of product form \eqref{EnE}. While a multi-time wave function is still well defined \cite{PT13b}, it is no longer of the form \eqref{phiW}, and it is defined only on spacelike configurations. As a consequence, \eqref{probnphi} will be valid only when the $x_j$ are mutually spacelike (while $n(x_j)$ may be either spacelike or timelike). When $x_j$ and $x_k$ are timelike separated, then also the calculation with $\phi$ needs to take collapses into account. The formulation of the rule given in Section~\ref{sec:rule-several} in terms of a single-time wave function $\psi_t$ is still correct, except that its time evolution is different from the non-interacting one expressed by \eqref{multiDirac}. For setting up a model with interaction, it makes a difference whether detected particles get absorbed or are allowed to further interact with the undetected particles; also, whether the particles emitted by the $n$ Dirac particles get detected (or absorbed, reflected, or transmitted) on $\tilde{\partial} S$. It would be of interest to study such a model.

\medskip

\noindent \textit{Acknowledgments.}
I thank Julian Schmidt and Stefan Teufel for helpful discussions.


\section{Introduction}

The familiar van der Corput $k$-th derivative estimate for exponential
sums (Titchmarsh \cite[Theorems 5.9, 5.11, \& 5.13]{Titch}, for
example), may be stated as follows. Let $k\ge 2$ be an integer, and
suppose that $f(x):[0,N]\to\R$ has continuous
derivatives of order up to $k$ on $(0,N)$.
Suppose further that $0<\lambda_k\le f^{(k)}(x)\le A\lambda_k$ on
$(0,N)$. Then
\beql{vdc}
\sum_{n\le N}e(f(n))\ll A^{2^{2-k}}N\lambda_k^{1/(2^k-2)}
+N^{1-2^{2-k}}\lambda_k^{-1/(2^k-2)},
\eeq
with an implied constant independent of $k$. One usually
chooses $k$ so that the first term dominates, and one often has
$A^{2^{2-k}}\ll 1$, so that the bound is merely
$O(N\lambda_k^{1/(2^k-2)})$. Clearly one can only get a non-trivial
bound when $\lk<1$.  A typical application is the series of
estimates
\[\zeta(\sigma+it)\ll t^{1/(2^k-2)}\log t,\;\;\;
\left(\sigma=1-\frac{k}{2^k-2},\; t\ge 2\right)\]
for $k=2,3,\ldots$. Again the implied constant is independent of $k$.

One can improve on the standard $k$-th derivative bound somewhat.  Thus
Robert and Sargos \cite{RS} show roughly that if $k=4$ then
\[\sum_{n\le N}e(f(n))\ll_{\ep} N^{\ep}(N\lambda_4^{1/13}+\lambda_4^{-7/13}),\]
for any $\ep>0$.  In the corresponding version of (\ref{vdc}) one
would have a term $N\lambda_4^{1/14}$ in place of
$N\lambda_4^{1/13}$. Similarly for $k=8$ and 9, Sargos \cite[Theorems
3 \& 4]{Sarg} gives bounds
\[\sum_{n\le N}e(f(n))\ll_{\ep}
N^{\ep}(N\lambda_8^{1/204}+\lambda_8^{-95/204}),\]
and
\[\sum_{n\le N}e(f(n))\ll_{\ep}
N^{\ep}(N\lambda_9^{7/2640}+\lambda_9^{-1001/2640}),\]
respectively.  Here the exponents $1/204$ and $7/2640$ should be
compared with the values $1/254$ and $1/510$ produced by (\ref{vdc}).
\bigskip

There are quite different approaches to exponential sums, using
estimates for the Vinogradov mean value integral
\beql{Jd}
J_{s,l}(P)=\int_0^1\ldots\int_0^1\left|\sum_{n\le P}e(\alpha_1
  n+\ldots +\alpha_l n^l)\right|^{2s}d\bal,
\eeq
see Vinogradov \cite{VB}, \cite{Vk}, and Korobov \cite{Kv}, amongst others.
The first of these methods is described by Titchmarsh \cite[Chapter 6]{Titch}
for example. The Vinogradov-Korobov machinery has been used by Ford
\cite[Theorem 2]{ford} to show that
\beql{f1}
\sum_{N<n\le 2N}n^{-it}\ll N^{1-1/134k^2}
\eeq
for $N^k\ge t\ge 2$.  (Ford's result is somewhat more precise, and
more general.)  One may think of this as corresponding very roughly to
a bound of the form (\ref{vdc}) with first term
$N\lambda_k^{1/134k^2}$.

A slightly refined version of
the original method of Vinogradov \cite{VB} coupled with
new estimates for
the Vinogradov mean value integral,
leads to distinctly stronger
bounds.  For example, Wooley \cite[Theorem 1.2]{W3} gives
\[J_{s,l}(P)\ll_{\ep,l} P^{2s-l(l+1)/2+\ep}\;\;\;(s\ge l(l-1),\]
and Robert \cite[Theorem 10]{Rob} used this to show that if $k\ge 4$ then
\[\sum_{n\le N}e(f(n))\ll_{A,k,\ep}
N^{1+\ep}(\lk^{1/2(k-1)(k-2)}+N^{-1/2(k-1)(k-2)})\]
for $N\ge\lk^{-(k-1)/(2k-3)}$.  This is a remarkable improvement on the
  classical $k$-th derivative estimate.  The exponent of $\lk$ is
  better than $1/(2^k-2)$ for all $k\ge 4$, and decreases
  quadratically rather than exponentially.

The purpose of this paper is to further refine
the original method of Vinogradov \cite{VB} and to input the very
recent optimal bounds for
the Vinogradov mean value integral, due to Wooley \cite{W3} (for $l=3$), and to
Bourgain, Demeter and Guth \cite{BDG} (for $l\ge 4$). These theorems
show that
\beql{bdg}
J_{s,l}(P)\ll_{\ep,l} P^{2s-l(l+1)/2+\ep}\;\;\;(s\ge\tfrac12 l(l+1),\;
l\ge 1),
\eeq
the cases $l=1$ and $l=2$ being elementary. The range for $s$ is
optimal, and it
is this feature that represents the dramatic culmination of many
previous works over the past 80 years. Unfortunately neither result
gives an explicit dependence on $l$ and $s$, nor gives an explicit form
for the factor $P^{\ep}$. Results prior to the advent of Wooley's
efficient congruencing method had required $s$ to be
larger, but had given an explicit dependence on $l$.  Thus for
example, Ford
\cite[Theorem 3]{ford} implies in particular that
\[J_{s,l}(P)\ll l^{13l^3/4} P^{2s-l(l+1)/2+l^2/1000}\;\;\;
(s\ge\tfrac52 l^2,\; l\ge 129),\]
in which one has an additional term $l^2/1000$ in the exponent, and
more restrictive conditions on $l$ and $s$. An important application
of bounds for Weyl sums is to the zero-free region for $\zeta(s)$, as
described by Ford.  However for this it is crucial to have a suitable
dependence on the parameter $l$, so that the new result of Bourgain,
Demeter and Guth is not applicable.

Our first result gives a new $k$-th derivative estimate
\begin{theorem}\label{esV}
Let $k\ge 3$ be an integer, and
suppose that $f(x):[0,N]\to\R$ has continuous
derivatives of order up to $k$ on $(0,N)$.  Suppose further that
\[0<\lambda_k\le f^{(k)}(x)\le A\lambda_k,\;\;\; x\in (0,N). \]
Then
\[
\sum_{n\le N}e(f(n)) \ll_{A,k,\ep}
N^{1+\ep}(\lk^{1/k(k-1)}+N^{-1/k(k-1)}+N^{-2/k(k-1)}\lk^{-2/k^2(k-1)}).
\]
\end{theorem}

If one thinks of $N\lk^{1/k(k-1)}$ as being the leading term here,
then one needs to compare the exponent $1/k(k-1)$ with the
corresponding exponent $1/(2^k-2)$ in (\ref{vdc}).  These agree for
$k=3$, but for larger values of $k$ the new exponent tends to zero far
more slowly than the old one.  It may perhaps be something of a
surprise that an analysis via Vinogradov's mean value integral
reproduces the same term $N\lambda_3^{1/6}$ as in the classical
third-derivative estimate.

We should emphasize that the stength of Theorem \ref{esV} comes almost
entirely from the new bound (\ref{bdg}). One could have injected
(\ref{bdg}) into the method of Robert \cite{Rob}, to produce an
estimate with the same terms $\lk^{1/k(k-1)}+N^{-1/k(k-1)}$ as in
Theorem \ref{esV}, but valid only for $N\ge \lk^{-(k-1)/(2k-3)}$.
Our result, incorporating a slightly better way of using the
Vinogradov mean value, gives the terms $\lk^{1/k(k-1)}+N^{-1/k(k-1)}$
in the substantially longer range $N\ge \lk^{-2/k}$.  However for our
application to Theorems \ref{EP}--\ref{rich} below, Robert's range
would have been very nearly sufficient.

The secondary terms in the bound given by Theorem \ref{esV}
are somewhat awkward. The
classical estimate (\ref{vdc}) leads easily to an exponent pair,
\[\left(\frac{1}{2^k-2}\,,\,\frac{2^k-k-1}{2^k-2}\right)\]
in which the term $N^{1-2^{2-k}}\lk^{-1/(2^k-2)}$ has no effect.
However the situation with Theorem \ref{esV} is more complicated. None
the less we are able to produce a series of new exponent pairs.

Before stating the
result we remind the reader of the necessary background. Let $s$ and
$c$ be positive constants, and let $\mathcal{F}(s,c)$ be the set of
quadruples $(N,I,f,y)$ where $y\ge N^s$ are positive real numbers, $I$
is a subinterval of $(N,2N]$, and $f$ is an infinitely differentiable
function on $I$, with
\[\left|f^{(n+1)}(x)-\frac{d^n}{dx^n}(yx^{-s})\right|\le c
\left|\frac{d^n}{dx^n}(yx^{-s})\right|\]
for $x\in I$, for all $n\ge 0$.  We then say that $(p,q)$ is an
exponent pair, if $p$ and $q$ lie in the range
$0\le p\le\tfrac12\le q\le 1$, and for each $s$ there
is a corresponding $c=c(p,q,s)>0$ such that
\[\sum_{n\in I}e(f(n))\ll_{p,q,s} (yN^{-s})^p N^q,\]
uniformly for all quadruples $(N,I,f,y)\in\mathcal{F}(s,c)$.

We then have the following.
\begin{theorem}\label{EP}
For any integer $k\ge 3$ and any real $\ep>0$ there is an exponent
pair given by
\beql{pqp}
p=\frac{2}{(k-1)^2(k+2)},
\eeq
and
\beql{pqq}
q=\frac{k^3+k^2-5k+2}{k(k-1)(k+2)}+\ep
=1-\frac{3k-2}{k(k-1)(k+2)}+\ep.
\eeq
\end{theorem}

In fact we are able to handle a much weaker condition on $f$. Let
$\mathcal{A}=(a_k)_3^\infty$ and $\mathcal{B}=(b_k)_3^\infty$ be
sequences of  positive real numbers, and let
$\mathcal{G}(\mathcal{A},\mathcal{B})$ be the set of
quadruples $(N,I,g,T)$ where $T\ge N$ are positive real numbers, $I$
is a subinterval of $(N,2N]$, and $g$ is an infinitely differentiable
function on $I$, with
\[a_kTN^{-k}\le\left|g^{(k)}(x)\right|\le b_k TN^{-k}\]
for $x\in I$, and for all $k\ge 3$. We then have the following.
\begin{theorem}\label{EP0}
For any integer $k\ge 3$ and any real $\ep>0$, let $p$ and $q$ be given
by (\ref{pqp}) and (\ref{pqq}). Then
\[\sum_{n\in I}e(g(n))\ll_{k,\ep,\mathcal{A},\mathcal{B}} (TN^{-1})^p N^q,\]
uniformly for $(N,I,g,T)\in\mathcal{G}(\mathcal{A},\mathcal{B})$.
\end{theorem}
If $(N,I,g,y)\in\mathcal{F}(s,\tfrac14)$, then
$(N,I,g,yN^{1-s})\in\mathcal{G}(\mathcal{A},\mathcal{B})$ with
\[a_k=\frac{3\times 2^{1-2-k}}{4s(s+1)\ldots(s+k-2)},\;\;\;
b_k=\frac{5}{4s(s+1)\ldots(s+k-2)}.\]
The sequences $\mathcal{A}$ and $\mathcal{B}$ depend only on $s$, and
we immediately see that Theorem \ref{EP} follows from Theorem \ref{EP0}.

We next present a slightly weaker version of Theorem \ref{EP0}, which is
somewhat more immediately intelligible. It will be convenient to write
$T=N^{\tau}$.
\begin{theorem}\label{EP1}
Let sequences $\mathcal{A}$ and $\mathcal{B}$, and a real number
$\ep>0$ be given, then
\[\sum_{n\in I}e(g(n))\ll_{\ep,\mathcal{A},\mathcal{B}} N^{1-49/(80\tau^2)+\ep},\]
uniformly for quadruples
$(N,I,g,T)\in\mathcal{G}(\mathcal{A},\mathcal{B})$ with $N\le T^{1/2}$.
\end{theorem}
The constant $49/80$ arises from the use of an exponent pair
\[(\tfrac{1}{20}\,,\,\tfrac{33}{40})=A^2BA^2B(0,1)\]
when $\tau=\tfrac72$. One could improve the constant slightly by employing a
better exponent pair. As will be clear from the proof, the constant
$\tfrac{49}{80}$ may be replaced by $1-\delta$ for any small
$\delta>0$, if we restrict to sufficiently large values
$\tau\ge\tau(\delta)$.

As an example of Theorem \ref{EP1}, if $t\ge 2$ we find that
\beql{f2}
\sum_{n\in I}n^{-it}\ll_{\ep} N^{1-49/80\tau^2+\ep},
\eeq
for $\tau=(\log t)/(\log N)\ge 2$. This should be compared with
(\ref{f1}). Using (\ref{f2}) we produce the following result.
\begin{theorem}\label{rich}
Let $\kappa=\tfrac{8}{63}\sqrt{15}=0.4918\ldots$.  Then for any fixed
$\ep>0$ we have
\beql{rich1}
\zeta(\sigma+it)\ll_{\ep} t^{\kappa(1-\sigma)^{3/2}+\ep}
\eeq
uniformly for $t\ge 1$ and $\tfrac12\le\sigma\le 1$. Moreover we have
\beql{rich2}
\zeta(\sigma+it)\ll_{\ep} t^{\tfrac12 (1-\sigma)^{3/2}+\ep}
\eeq
uniformly for $t\ge 1$ and $0\le\sigma\le 1$.
\end{theorem}

One sees from the proof that $\kappa$ may be reduced to
$2/\sqrt{27}+\delta=0.3849\ldots$ for any small $\delta>0$, if we
restrict $\sigma$ to a suitably small range
$\sigma(\delta)\le\sigma\le 1$.
The corresponding result in the work of Ford \cite[Theorem 1]{ford}
states that
\[|\zeta(\sigma+it)|\le 76.2 t^{4.45(1-\sigma)^{3/2}}(\log t)^{2/3}\]
for $t\ge 3$ and $\tfrac12\le\sigma\le 1$.  Thus we have reduced the
constant $4.45$ to $0.4918\ldots$. Unfortunately our result yields no useful
information when $\sigma$ tends to 1, which is a critical situation in
many applications.  Moreover we do not have the explicit order constant
that Ford finds.

As Ford explains, there are a number of interesting corollaries, for which we
merely have to replace the constant $B=4.45$ by $B=0.492$ in the
arguments given in \cite[Pages 566 and 567]{ford}.
We can feed our bound into the zero-density theorem of Montgomery
\cite[Theorem 12.3]{mont} (with $1-\alpha=4.93(1-\sigma)$ as used by Ford
\cite[Page 566]{ford}) to give the following.
\begin{corollary}
We have
\[N(\sigma,T)\ll_{\ep}T^{6.42(1-\sigma)^{3/2}+\ep}\]
for $\tfrac{9}{10}\le\sigma\le 1$.
\end{corollary}

For moments of the Riemann Zeta-function we have:
\begin{corollary}
For any positive integer $k$ one has
\[\int_{0}^T|\zeta(\sigma+it)|^{2k}dt\sim
T\sum_{1}^\infty d_k(n)^2n^{-2\sigma},\]
as $t\to\infty$, for any fixed $\sigma\ge 1-0.534 k^{-2/3}$.
\end{corollary}

For the generalized divisor problem we have:
\begin{corollary}
For any positive integer $k$ the error term $\Delta_k(x)$ in the
generalized divisor problem satisfies
\[\Delta(x)\ll_k x^{1-0.849 k^{-2/3}}.\]
\end{corollary}

In Section \ref{sec2} we will reduce the proof of Theorem \ref{esV}
to a two-variable
counting problem involving fractional parts of the derivatives $f^{(j)}(n)$.
Section \ref{secN} shows how this counting problem is tackled, and
finally Section \ref{sPT} completes the proof of our theorems.

{\bf Acknowledgement.}  This work was supported by EPSRC grant number
EP/K021132X/1

\section{Initial Steps}\label{sec2}

Our goal in the first stage of the proof is to estimate the sum
\[\Sigma=\sum_{n\le N}e(f(n))\]
in terms of $J_{s,l}(P)$, together with, a counting function involving
the fractional parts of numbers of the form $f^{(j)}(n)/j!$.
\begin{lemma}\label{l1}
Let $k\ge 2$ be an integer, and
suppose that $f(x):[0,N]\to\R$ has continuous
derivatives of order up to $k$ on $(0,N)$.  Suppose further that
\[0<\lambda_k\le f^{(k)}(x)\le A\lambda_k,\;\;\; x\in (0,N), \]
and that $A\lk\le \tfrac14$.  Then
\[\Sigma\ll H+k^2 N^{1-1/s}\mathcal{N}^{1/2s}
\left\{H^{-2s+k(k-1)/2}J_{s,k-1}(H)\right\}^{1/2s},\]
where $H=[(A\lk)^{-1/k}]$ and
\[\mathcal{N}=\#\left\{m,n\le N:
\left|\left|\frac{f^{(j)}(m)}{j!}-\frac{f^{(j)}(n)}{j!}\right|\right|
\le 2H^{-j}\mbox{ for } 1\le j\le k-1\right\}.\]
\end{lemma}

If $J_{s,k-1}(H)\ll_{\ep,k} H^{2s-k(k-1)/2+\ep}$ as in (\ref{bdg})
the estimate in the lemma reduces to
\beql{basic}
\Sigma\ll_{\ep,k} H+N^{1-1/s+\ep}\mathcal{N}^{1/2s}.
\eeq
Here we would want to choose $s$ to be as small as possible, and since
we are taking $l=k-1$ this means that we will have $s=k(k-1)/2$.

The lemma is clearly trivial if $H\ge N$, and we may therefore suppose
for the proof that $H\le N$.
For any positive integer $H\le N$ we will have
\[H\Sigma=\sum_{h\le H}\sum_{-h<n\le N-h}e(f(n+h))=
\sum_{h\le H}\sum_{1\le n\le N-H}e(f(n+h))+O(H^2),\]
so that
\beql{Sigma1}
\Sigma=H^{-1}\sum_{n\le N-H}\;\sum_{h\le H}e(f(n+h))+O(H).
\eeq

We proceed to approximate $f(n+h)$ by the polynomial
\[f_n(h):=f(n)+f'(n)h+\ldots+\frac{f^{(k-1)}(n)}{(k-1)!}h^{k-1}.\]
To do this we set $g_n(x)=f(n+x)-f_n(x)$ and use summation by parts to
obtain the bound
\[\sum_{h\le H}e(f(n+h))\ll |S_n(H)|+\int_0^H |S_n(x)g_n'(x)|dx,\]
where we have written
\[S_n(x)=\sum_{h\le x}e(f_n(h))\]
for convenience.

If $0\le x\le H$ we may use Taylor's Theorem with Lagrange's form of
the remainder to show that
\[f'(n+x)=f_n'(x)+\frac{f^{(k)}(\xi)}{k!}x^{k-1}\]
for some $\xi\in(n,n+x)\subseteq(0,N)$. It follows that
\[g_n'(x)\ll A\lk H^{k-1}\]
on $[0,H]$. With the choice $H=[(A\lk)^{-1/k}]$
we find that
\[\sum_{h\le H}e(f(n+h))\ll |S_n(H)|+H^{-1}\int_0^H |S_n(x)|dx.\]
The bound (\ref{Sigma1}) now yields
\[\Sigma\ll H+H^{-1}\sum_{n\le N-H}|S_n(H)|+
H^{-2}\int_0^H \left\{\sum_{n\le N-H} |S_n(x)|\right\}dx.\]
It then follows that there is a positive integer $H_0\le H$ such that
\beql{b2}
\Sigma\ll H+H^{-1}\sum_{n\le N-H}|S_n(H_0)|.
\eeq

Now suppose that $\bal\in [0,1]^{k-1}$ and
\beql{bal}
||f^{(j)}(n)/j!-\alpha_j||\le H^{-j}\mbox{ for } 1\le j\le k-1
\eeq
where
\[||\theta||=\min_{n\in\Z}|\theta-n|\]
as usual. We proceed to replace $f_n(h)$ by
\[f(h;\bal)=\alpha_1 h+\ldots+\alpha_{k-1}h^{k-1}\]
as follows. Firstly we remove the constant term $f(n)$ from
$f_n(h)$. This has no effect on $|S_n(H_0)|$. Next, we replace each
coefficient $f^{(j)}(n)/j!$ by $c_j$, say, with
$f^{(j)}(n)/j!-c_j\in\Z$, so that $|c_j-\alpha_j|\le H^{-j}$, and
denote the resulting polynomial by $f^*_n(h)$. If we write
\[S^*_n(H_0)=\sum_{h\le H_0}e(f^*_n(h))\]
then clearly $|S_n(H_0)|=|S^*_n(H_0)|$.  Moreover
\[\frac{d}{dx}\left(f(x;\bal)-f^*_n(x)\right)\ll k^2\max_{j\le
  k-1}|c_j-\alpha_j|H^{j-1} \ll k^2H^{-1}.\]
It therefore follows on summing by parts that
\[S^*_n(H_0)\ll |S(H_0;\bal)|+k^2 H^{-1}\int_0^{H_0}|S(x;\bal)|dx,\]
where we have set
\[S(x;\bal)=\sum_{h\le x}e(f(h;\bal)).\]
We may therefore conclude that
\[S^*_n(H_0)\ll 2^{-k}H^{k(k-1)/2}\left\{\int_{\bal}|S(H_0;\bal)|
d\bal+k^2 H^{-1}\int_0^{H_0}\int_{\bal}|S(x;\bal)|
d\bal dx\right\},\]
where the integral over $\bal$ is for vectors in $[0,1]^{k-1}$
satisfying (\ref{bal}).

For each $\bal\in [0,1]^{k-1}$ we now define
\[\nu(\bal)=\#\{n\le N-H: ||f^{(j)}(n)/j!-\alpha_j||\le H^{-j}\mbox{
  for } 1\le j\le k-1\}.\]
We then find that
\beql{b1}
\sum_{n\le N-H}|S_n(H_0)|\ll 2^{-k}H^{k(k-1)/2}\left\{I(H_0)+
k^2 H^{-1}\int_0^{H_0}I(x)dx\right\},
\eeq
with
\[I(x)=\int_0^1\ldots\int_0^1|S(x;\bal)|\nu(\bal)
d\bal.\]
We easily see that
\[\int_0^1\ldots\int_0^1\nu(\bal)d\bal=
2^{k-1}H^{-k(k-1)/2}(N-H),\]
and that
\[\int_0^1\ldots\int_0^1\nu(\bal)^2d\bal\le
2^{k-1}H^{-k(k-1)/2}\mathcal{N},\]
where $\mathcal{N}$ is defined in Lemma \ref{l1}.
Moreover
\[\int_0^1\ldots\int_0^1|S(x;\bal)|^{2s}d\bal=J_{s,k-1}(x)\]
in the notation of (\ref{Jd}). Since $J_{s,k-1}(P)$ is non-decreasing
in $P$ this last integral may be bounded by $J_{s,k-1}(H)$.

Hence, by H\"{o}lder's inequality, for any positive integer $s$ we have
\[I(x)\ll
2^kH^{-k(k-1)/2}N^{1-1/s}\mathcal{N}^{1/2s}
\left\{H^{k(k-1)/2}J_{s,k-1}(H)\right\}^{1/2s}.\]
Thus (\ref{b1}) yields
\[\sum_{n\le N-H}|S_n(H_0)|\ll k^2 N^{1-1/s}\mathcal{N}^{1/2s}
\left\{H^{k(k-1)/2}J_{s,k-1}(H)\right\}^{1/2s}\]
and (\ref{b2}) gives us
\[\Sigma\ll H+k^2 N^{1-1/s}\mathcal{N}^{1/2s}
\left\{H^{-2s+k(k-1)/2}J_{s,k-1}(H)\right\}^{1/2s}\]
as required.

\section{The counting function $\mathcal{N}$}\label{secN}

Naturally our next task is to bound $\mathcal{N}$. The original
approach taken by Vinogradov, as described in Titchmarsh \cite[Chapter
6]{Titch}, merely used an $L^{\infty}$ bound for $\nu(\bal)$. One
discards all the information on $f^{(j)}(n)/j!$ for $j\le k-2$ and
uses only the case $j=k-1$. One then employs a standard procedure given by
the following trivial variant of \cite[Lemma 6.11]{Titch}, for
example.
\begin{lemma}\label{space}
Let $N$ be a positive integer, and suppose that $g(x):[0,N]\to\R$ has a
continuous derivative on $(0,N)$.  Suppose further that
\[0<\mu\le g'(x)\le A_0\mu,\;\;\; x\in(0,N).\]
Then
\[\#\{n\le N: ||g(n)||\le \theta\}\ll (1+A_0\mu N)(1+\mu^{-1}\theta).\]
\end{lemma}
We fix $m$ and take
\[g(x)=\frac{f^{(k-1)}(x)-f^{(k-1)}(m)}{(k-1)!}\]
and $\mu=\lk/(k-1)!$, $A_0=A$. This leads to a bound
\[\mathcal{N}\ll (k-1)!N(1+AN\lk)(1+H^{1-k}\lk^{-1}).\]
Under the assumption $A\lk\le \tfrac14$ in Lemma \ref{l1} we have
\[H^{1-k}\lk^{-1}\asymp (A\lk)^{1-1/k}\lk^{-1}=A(A\lk)^{-1/k}\ge A\ge 1,\]
whence our bound produces
\beql{addd}
\mathcal{N}\ll A^2(k-1)!N\lk^{-1/k}(1+N\lk).
\eeq
If one inserts this into (\ref{basic}) with $s=k(k-1)/2$
one gets an estimate
\begin{eqnarray*}
\Sigma&\ll_{\ep,k}&
(A\lk)^{-1/k}+N^{1-1/s+\ep}\{A^2N\lk^{-1/k}(1+N\lk)\}^{1/2s}\\
&\ll_{\ep,k}&
AN^{\ep}\{\lk^{-1/k}+N^{1-1/k(k-1)}\lk^{-1/k^2(k-1)}+N\lk^{1/k^2}\}.
\end{eqnarray*}
In fact the first term can be dropped, giving
\beql{vv}
\Sigma\ll_{\ep,k}AN^{\ep}\{N\lk^{1/k^2}+N^{1-1/k(k-1)}\lk^{-1/k^2(k-1)}\}.
\eeq
To see this we note that we have
\[\Sigma\ll N^{1-1/k(k-1)}\lk^{-1/k^2(k-1)}\]
trivially
unless
\[N^{1-1/k(k-1)}\lk^{-1/k^2(k-1)}\le N.\]
In this latter case however one sees that
\[\lk^{-1/k}\le N^{1-1/k(k-1)}\lk^{-1/k^2(k-1)}.\]

We may therefore regard (\ref{vv}) as being the result that
Vinogradov's method achieves, given the results of Wooley \cite{W3}
and Bourgain, Demeter and Guth \cite{BDG}.  It is already a remarkable
improvement on (\ref{vdc}), replacing the critical exponent
$1/(2^k-2)$ by $1/k^2$.  Thus, in appropriate circumstances, we get
an improvement as soon as $k\ge 5$.  Our goal in this section is to
make the following small further sharpening in the estimation of
$\mathcal{N}$.
\begin{lemma}\label{N2}
When $k\ge 3$ we have
\[\mathcal{N}\ll
\big((k-1)!A\big)^4(N+\lk N^2+\lk^{-2/k})\log N.\]
\end{lemma}
Apart from the term $\lk^{-2/k}$, which is insignificant in
applications, this represents an improvement of (\ref{addd}) by a
factor $\ll_{A,k}\lk^{1/k}$.

On the one hand our proof will use the fact that $\mathcal{N}$ is a counting
function of two variables $m$ and $n$.  On the other we shall use information
about both $f^{(k-1)}$ and $f^{(k-2)}$.  The reader may find it
slightly surprising in the light of this that our bound depends on
$\lk$ only, and not on estimates for other derivatives $f^{(j)}$. The
introduction of $\mathcal{N}$, and our procedure for estimating it,
are the only really new aspects to this paper.
\bigskip

We begin our analysis by assuming that $k\ge 3$ and noting that
$\mathcal{N}$ is at most
\[\mathcal{N}_1=\#\left\{m,n\le N:
\left|\left|\frac{f^{(j)}(m)}{j!}-\frac{f^{(j)}(n)}{j!}\right|\right|
\le 2H^{-j}\mbox{ for } j=k-2, k-1\right\}.\]
We proceed to show that it suffices to consider pairs $m,n$ of integers that
are relatively close. It will be convenient to write $B=4H^{2-k}$ and
$C=4H^{1-k}$ and to set
\[g_1(x)=\frac{f^{(k-2)}(x)}{(k-2)!},\;\;\; g_2(x)=\frac{f^{(k-1)}(x)}{(k-1)!}.\]
We also define the doubly-periodic function
\[\phi(x,y)=\max\left(1-B^{-1}||x||,0\right)\max\left(1-C^{-1}||y||,0\right),\]
so that
\[\mathcal{N}_1\ll \sum_{m,n\le N}
\phi\big(g_1(m)-g_1(n)\,,\, g_2(m)-g_2(n)\big).\]
The function $\phi(x,y)$ has an absolutely convergent Fourier series
\[\phi(x,y)=\sum_{r,s\in\Z}c_{r,s}e(rx+sy)\]
with non-negative coefficients
\[c_{r,s}=BC\left(\frac{\sin(\pi rB)\sin(\pi sC)}{\pi^2
    rsBC}\right)^2.\]
Thus
\begin{eqnarray*}
\mathcal{N}_1&\ll& \sum_{r,s\in\Z}c_{r,s}\sum_{m,n\le N}
e\big(r(g_1(m)-g_1(n))+s(g_2(m)-g_2(n))\big)\\
&=& \sum_{r,s\in\Z}c_{r,s}\left|\sum_{n\le N}
e\big(rg_1(n)+sg_2(n)\big)\right|^2.
\end{eqnarray*}
Let $K$ be a positive integer parameter, to be chosen later.
We proceed to partition the range $(0,N]$ into $K$ intervals
$I_i=(a_i,b_i]$ for $i\le K$, having integer endpoints, and length
$b_i-a_i\le 1+N/K$.  An application of Cauchy's inequality then yields
\begin{eqnarray*}
\mathcal{N}_1 &\ll & K\sum_{i\le K}\sum_{r,s\in\Z}c_{r,s}
\left|\sum_{n\in I_i}e\big(rg_1(n)+sg_2(n)\big)\right|^2\\
&=&K\sum_{i\le K}\sum_{r,s\in\Z}c_{r,s}\sum_{m,n\in I_i}
e\big(r(g_1(m)-g_1(n))+s(g_2(m)-g_2(n))\big)\\
&=&K\sum_{i\le K}\sum_{m,n\in I_i}
\phi\big(g_1(m)-g_1(n)\,,\, g_2(m)-g_2(n)\big)\\
&\le & K\twosum{m,n\le N}{|m-n|\le 1+N/K}
\phi\big(g_1(m)-g_1(n)\,,\, g_2(m)-g_2(n)\big).
\end{eqnarray*}
We may therefore conclude that $\mathcal{N}_1\ll K\mathcal{N}_2$,
where $\mathcal{N}_2$ counts pairs of integers $m,n\le N$ with
$|m-n|\le 1+N/K$ for which
\[\left|\left|\frac{f^{(j)}(m)}{j!}-\frac{f^{(j)}(n)}{j!}\right|\right|
\le 4H^{-j}\mbox{ for } j=k-2, k-1.\]

If $|m-n|\le 1+N/K$ we will have
\[\left|\frac{f^{(k-1)}(m)}{(k-1)!}-\frac{f^{(k-1)}(n)}{(k-1)!}\right|\le
\frac{|m-n|}{(k-1)!}\sup |f^{(k)}|\le A\lk(1+N/K),\]
by the mean-value theorem. We will choose
\[K=1+[4A\lk N],\]
so that
\[\left|\frac{f^{(k-1)}(m)}{(k-1)!}-\frac{f^{(k-1)}(n)}{(k-1)!}\right|\le
\frac{1}{2},\]
in view of our assumption that $A\lk\le\tfrac14$. Thus if
\[\left|\left|\frac{f^{(k-1)}(m)}{(k-1)!}-
\frac{f^{(k-1)}(n)}{(k-1)!}\right|\right|\le 4H^{1-k}\]
we must have
\[\left|\frac{f^{(k-1)}(m)}{(k-1)!}-\frac{f^{(k-1)}(n)}{(k-1)!}\right|
\le 4H^{1-k}.\]
However the mean-value theorem also tells us that
\[\left|\frac{f^{(k-1)}(m)}{(k-1)!}-\frac{f^{(k-1)}(n)}{(k-1)!}\right|\ge
\frac{|m-n|}{(k-1)!}\inf |f^{(k)}|\ge \lk\frac{|m-n|}{(k-1)!}.\]
We therefore conclude that
\[|m-n|\le \frac{4(k-1)!}{\lk H^{k-1}}\]
for any pair $m,n$ counted by $\mathcal{N}_2$.

There are $N$ pairs $m=n$ counted by $\mathcal{N}_2$. We consider
the remaining pairs with $m>n$, the alternative case producing the
same estimates by symmetry.  Then $m=n+d$ with $1\le d\le D$, where
\[D=\min\left(N\,,\,\left[\frac{4(k-1)!}{\lk H^{k-1}}\right]\right).\]
For each available value of $d$ we estimate the number of
corresponding integers $n$ via Lemma \ref{space}, taking
\[g(x)=\frac{f^{(k-2)}(x+d)-f^{(k-2)}(x)}{(k-2)!}.\]
Then
\[g'(x)=\frac{f^{(k-1)}(x+d)-f^{(k-1)}(x)}{(k-2)!},\]
so that
\[d\frac{\lk}{(k-2)!}\le d\frac{\inf |f^{(k)}|}{(k-2)!}\le g'(x)\le
d\frac{\sup |f^{(k)}|}{(k-2)!}\le d\frac{A\lk}{(k-2)!},\]
by the mean-value theorem. We therefore apply the lemma with $\mu=
\lk d/(k-2)!$ and $A_0=A$. This shows that each $d\ge 1$ contributes
\begin{eqnarray*}
&\ll& (k-2)!(1+AN\lk d)(1+H^{2-k}\lk^{-1}d^{-1})\\
&\ll& (k-2)!(1+AN\lk D)(D+H^{2-k}\lk^{-1})d^{-1}\\
&\ll& \big((k-1)!\big)^3A(1+NH^{1-k})H^{2-k}\lk^{-1}d^{-1}\\
&\ll& \big((k-1)!A\big)^3(1+N\lk^{1-1/k})\lk^{-2/k}d^{-1}.
\end{eqnarray*}
Summing for $d\le D$ we therefore find that
\[\mathcal{N}_2\ll N+\big((k-1)!A\big)^3(1+N\lk^{1-1/k})\lk^{-2/k}\log D.\]
Since $k\ge 3$, $\lk\le 1$ and $D\le N$ this simplifies to give
\[\mathcal{N}_2\ll \big((k-1)!A\big)^3(N+\lk^{-2/k})\log N,\]
whence
\[\mathcal{N}\le\mathcal{N}_1\ll K\mathcal{N}_2\ll (1+A\lk N)
\big((k-1)!A\big)^3(N+\lk^{-2/k})\log N\]
\[\ll \big((k-1)!A\big)^4(N+\lk N^2+\lk^{-2/k}+N\lk^{1-2/k})\log N.\]
Since $k\ge 3$ and $\lk\le 1$ we have $N\lk^{1-2/k}\le N$, and
Lemma \ref{N2} follows.

\section{Proof of the Theorems}\label{sPT}

If we insert Lemma \ref{N2} into Lemma \ref{l1}, and use the bound
(\ref{bdg}) with the choices $l=k-1, s=k(k-1)/2$, we see that
\[\Sigma\ll_{A,k,\ep}
N^{\ep}(\lk^{-1/k}+N^{1-1/k(k-1)}+N\lk^{1/k(k-1)}
+N^{1-2/k(k-1)}\lk^{-2/k^2(k-1)}).\]
The term $\lk^{-1/k}$ may be omitted, since the resulting bound
\[\Sigma\ll_{A,k,\ep}N^{\ep}(N^{1-1/k(k-1)}+N\lk^{1/k(k-1)}
+N^{1-2/k(k-1)}\lk^{-2/k^2(k-1)})\]
holds trivially when $N\le N^{1-2/k(k-1)}\lk^{-2/k^2(k-1)}$, while
\[\lk^{-1/k}\le N^{1-2/k(k-1)}\lk^{-2/k^2(k-1)}\]
when $N\ge N^{1-2/k(k-1)}\lk^{-2/k^2(k-1)}$.
This suffices for Theorem \ref{esV}.

We turn next to Theorem \ref{EP0}.  Suppose that
$(N,I,g,T)\in\mathcal{G}(\mathcal{A},\mathcal{B})$, and let $I$ have end points
$N_0$ and $N_0+N_1$, so that $N_1\le N$. We apply Theorem \ref{esV}
to the function $f(x)=g(N_0+x)$, taking $\lambda_k=a_kTN^{-k}$ and
$A=b_k/a_k$. (Since $f^{(k)}$ is differentiable it is continuous, and
hence it cannot change sign if $|f^{(k)}(x)|\ge a_kTN^{-k}>0$. Taking
complex conjugates of our sum if necessary we may therefore assume
that $f^{(k)}(x)$ is positive on $I$.) It follows that if $k\ge 3$ then
\begin{eqnarray*}
\lefteqn{\sum_{n\in I}e(g(n))}\\
&\ll_{\ep,k,\mathcal{A},\mathcal{B}}&
N^{1+\ep}(\lambda_k^{1/k(k-1)}+N^{-1/k(k-1)}+N^{-2/k(k-1)}\lambda_k^{-2/k^2(k-1)})\\
&\ll_{\ep,k,\mathcal{A},\mathcal{B}}&
N^{1+\ep}(N^{-1/(k-1)}T^{1/k(k-1)}+N^{-1/k(k-1)}+T^{-2/k^2(k-1)}).
\end{eqnarray*}

We use the above bound for
\[\frac{(k-1)^2+1}{k}\le\tau<\frac{k^2+1}{k+1}\]
where we define $\tau$ by  $T=N^{\tau}$.  For this range of $\tau$ we
find that
\begin{eqnarray}\label{mid}
\lefteqn{\max\left(\frac{\tau-k}{k(k-1)}\,,\,\frac{-1}{k(k-1)}\,,\,
\frac{-2\tau}{k^2(k-1)}\right)}\nonumber\\
&=&\left\{\begin{array}{cc}
-1/k(k-1), & ((k-1)^2+1)/k \le \tau\le k-1,\\
(\tau-k)/k(k-1), & k-1\le\tau<(k^2+1)/(k+1),\end{array}\right.\\
&\le & A_k\tau+B_k,\nonumber
\end{eqnarray}
where the coefficients $A_k$ and $B_k$ are chosen so that
\[A_k\frac{(k-1)^2+1}{k}+B_k=\frac{-1}{k(k-1)}\]
and
\[A_k\frac{k^2+1}{k+1}+B_k=\frac{(k^2+1)/(k+1)-k}{k(k-1)}=\frac{-1}{k(k+1)}.\]
One then calculates that
\[A_k=\frac{2}{(k-1)^2(k+2)}\;\;\;\mbox{and}\;\;\;
B_k=-\frac{3k^2-3k+2}{k(k-1)^2(k+2)}.\]
If we now define $\phi(\tau):[2,\infty)\to\R$ by taking
$\phi(\tau)=A_k\tau+B_k$ on
\[\left[\frac{(k-1)^2+1}{k}\,,\,\frac{k^2+1}{k+1}\right)\]
for each integer $k\ge 3$, we conclude that
\beql{nd}
\sum_{n\in I}e(g(n))\ll_{\ep,\tau_0,\mathcal{A},\mathcal{B}}
N^{1+\phi(\tau)+\ep},
\eeq
uniformly for $2\le\tau\le\tau_0$. The function $\phi$ is continuous,
and since the coefficients $A_k$ are monotonic decreasing $\phi$ is
also convex. It follows that $\phi(\tau)\le A_k\tau+B_k$ for any
$\tau\in[2,\infty)$ and any $k\ge 3$. Thus
\[\sum_{n\in I}e(g(n))\ll_{\ep,\tau_0,\mathcal{A},\mathcal{B}}
N^{1+B_k+\ep}T^{A_k}=(TN^{-1})^p N^q,\]
with $p,q$ given by (\ref{pqp}) and (\ref{pqq}).
As before, this is uniform in any
finite range $2\le\tau\le\tau_0$.  However if we set $\tau_0=1+(1-q)/p$
then $\tau_0$ will depend on $\ep$ and $k$ alone.  Moreover, if
$\tau\ge\tau_0$ then we trivially have
\[\sum_{n\in I}e(g(n))\ll N\le (TN^{-1})^pN^q.\]
Finally, if $\tau\le 2$ we use the well known exponent pair
$(\tfrac16,\tfrac23)$ to show that
\[\sum_{n\in I}e(g(n))\ll T^{1/6}N^{1/2}.\]
When $k\ge 3$ one easily verifies that $q\ge p+1/2$ and $p+q\ge 5/6$
for the values (\ref{pqp}) and (\ref{pqq}), whence $T^{1/6}N^{1/2}\le
T^pN^{q-p}$ for $N\ge T^{1/2}$. It then follows that
\[\sum_{n\in I}e(g(n))\ll T^{1/6}N^{1/2} \le T^pN^{q-p}\]
for the remaining range $1\le\tau\le 2$.  This completes the proof of
Theorem \ref{EP0}.
\bigskip

We move now to the proof of Theorem \ref{EP1}. Let
$\tau_0=\sqrt{49/80\ep^2}$.  Then if $\tau\ge\tau_0$ we will trivially
have
\[\sum_{n\in I}e(g(n))\ll N\le N^{1-49/80\tau^2+\ep}.\]
When $\tau\le\tau_0$ we begin by handling
the range $\tfrac{13}{3}\le\tau\le\tau_0$, for which we claim that
$\phi(\tau)\le -49/80\tau^2$. This will clearly suffice, in view of the
estimate (\ref{nd}). Since $\phi(\tau)$ is piecewise linear, while the
function $-49/80\tau^2$ is convex, it suffices to verify that
$\phi(\tau)\le -49/80\tau^2$ at each of the points
$\tau=(k^2+1)/(k+1)$, for $k\ge 5$. This condition is equivalent to
\[\frac{(k^2+1)^2}{k(k+1)^3}\ge\frac{49}{80}.\]
However the fraction on the right is increasing for $k\ge 5$, and
takes the value $169/270>49/80$ at $k=5$.

When $\tfrac72\le\tau\le\tfrac{13}{3}$ we will use the bounds
\[\sum_{n\in I}e(g(n))\ll_{\ep} N^{1-1/20+\ep},\;\;\;(\tfrac72\le\tau\le
4)\]
and
\[\sum_{n\in I}e(g(n))\ll_{\ep} N^{1-(5-\tau)/20+\ep},\;\;\;(4\le\tau\le
\tfrac{13}{3})\]
which come from the case $k=5$ of (\ref{mid}). Note that the first of
these is valid in the longer range $\tfrac{17}{5}\le\tau\le 4$, but we
shall only use it when $\tfrac72\le\tau\le 4$. We therefore need to
verify that $-1/20\le-49/80\tau^2$ for $\tfrac72\le\tau\le 4$ and that
$-(5-\tau)/20\le-49/80\tau^2$ for $4\le\tau\le \tfrac{13}{3}$. This
is routine, but we observe that we have equality at $\tau=\tfrac72$.

We next consider the case in which $\tfrac{59}{22}\le\tau\le\tfrac72$,
for which we use the bound
\[\sum_{n\in I}e(g(n))\ll (T/N)^{1/20}N^{33/40}=N^{1+(2\tau-9)/40}\]
corresponding to the exponent pair
$(\tfrac{1}{20},\tfrac{33}{40})$. (This pair is $A^2BA^2B(0,1)$ in the
usual notation, see Titchmarsh \cite[\S 5.20]{Titch}, for example.)
Again, it is routine to check that
\[\frac{2\tau-9}{40}\le-\frac{49}{80\tau^2},\;\;\;(\frac{59}{22}\le
\tau\le\frac{7}{2}).\]

Finally we examine the range $2\le\tau\le\tfrac{59}{22}$, and here we
use the bound
\[\sum_{n\in I}e(g(n))\ll (T/N)^{1/9}N^{13/18}=N^{1+(2\tau-7)/18}\]
corresponding to the exponent pair
$(\tfrac{1}{9},\tfrac{13}{18})$. (This pair is $ABA^2B(0,1)$ in the
usual notation, see Titchmarsh \cite[\S 5.20]{Titch}, for example.)
Another routine check shows that
\[\frac{2\tau-7}{18}\le-\frac{49}{80\tau^2},\;\;\;(2\le
\tau\le\frac{59}{22}),\]
thereby completing the proof of Theorem \ref{EP1}.
\bigskip

We turn now to Theorem \ref{rich}. If $\tau\ge 2$ we may use
(\ref{f2}) along with a
partial summation to obtain
\[\sum_{n\in J}n^{-\sigma-it}\ll_{\ep} N^{1-49/80\tau^2-\sigma+\ep}\le
t^{(1-\sigma)\tau^{-1}-\tfrac{49}{80}\tau^{-3}+\ep/2}\]
for any $\sigma\in[\tfrac12,1]$, and for any interval $J\subseteq(N,2N]$.
As a function of
$\tau\in(0,\infty)$ the exponent of $t$ is maximal at
\[\tau=\sqrt{\frac{147}{80(1-\sigma)}},\]
whence
\[\sum_{n\in J}n^{-\sigma-it}\ll_{\ep}t^{\kappa(1-\sigma)^{3/2}+\ep/2}.\]
Using a dyadic subdivision of $(0,N]$ we therefore have
\[\sum_{n\le N}n^{-\sigma-it}\ll_{\ep}t^{\kappa(1-\sigma)^{3/2}+3\ep/4}\]
for any $N\le t^{1/2}$. A further summation by parts then shows that
\[\sum_{n\le M}n^{-1+\sigma-it}
\ll_{\ep}M^{2\sigma-1}t^{\kappa(1-\sigma)^{3/2}+\ep}
\ll_{\ep}t^{\sigma-\tfrac12+\kappa(1-\sigma)^{3/2}+\ep}\]
for any $M\le t^{1/2}$. The required bound (\ref{rich1}) then follows
from the approximate functional equation for $\zeta(s)$.

The bound (\ref{rich2}) follows from (\ref{rich1}) when
$\tfrac12\le\sigma\le 1$, since $\kappa<\tfrac12$. For the remaining
range we use the functional equation, which shows that
\[\zeta(\sigma+it)\ll t^{\tfrac12-\sigma}|\zeta(1-\sigma+it)|\ll_{\ep}
t^{\tfrac12-\sigma+\tfrac12 \sigma^{3/2}+\ep}.\]
However one can readily verify that
\[\frac{1}{2}-\sigma+\frac{\sigma^{3/2}}{2}\le\frac{(1-\sigma)^{3/2}}{2}\]
for $0\le\sigma\le\tfrac12$, which completes the proof of Theorem
\ref{rich}.