# A PAIR CORRELATION PROBLEM, AND COUNTING LATTICE POINTS WITH THE ZETA FUNCTION

CHRISTOPH AISTLEITNER, DANIEL EL-BAZ AND MARC MUNSCH

*Dedicated to the memory of Professor Aleksandar Ivić (1949–2020)*

<figure><img src="image_1.png" /></figure>

Check for
updates

**Abstract.** The pair correlation is a localized statistic for sequences in the unit interval. Pseudo-random behavior with respect to this statistic is called Poissonian behavior. The metric theory of pair correlations of sequences of the form $(a_n\alpha)_{n\ge1}$ has been pioneered by Rudnick, Sarnak and Zaharescu. Here $\alpha$ is a real parameter, and $(a_n)_{n\ge1}$ is an integer sequence, often of arithmetic origin. Recently, a general framework was developed which gives criteria for Poissonian pair correlation of such sequences for almost every real number $\alpha$, in terms of the additive energy of the integer sequence $(a_n)_{n\ge1}$. In the present paper we develop a similar framework for the case when $(a_n)_{n\ge1}$ is a sequence of reals rather than integers, thereby pursuing a line of research which was recently initiated by Rudnick and Technau. As an application of our method, we prove that for every real number $\theta > 1$, the sequence $(n^\theta\alpha)_{n\ge1}$ has Poissonian pair correlation for almost all $\alpha \in \mathbb{R}$.

## 1 Introduction and statement of results

A sequence $(y_n)_{n\ge1}$ of real numbers is called uniformly distributed (or equidistributed) modulo one if for all intervals $A \subset [0, 1)$ the asymptotic equality

$$ \lim_{N \to \infty} \frac{1}{N} \sum_{n=1}^{N} \mathbf{1}_{A}(y_n) = \lambda(A) \quad (1) $$

holds. Here $\mathbf{1}_A$ is the indicator function of $A$, extended periodically with period 1, and $\lambda$ denotes Lebesgue measure. Uniform distribution theory has a long history, going back to the seminal paper of Hermann Weyl [Wey16]. For general background, see [DT97, KN74]. Uniform distribution of a sequence can be seen as a pseudorandomness property, in the sense that a sequence $(Y_n)_{n\ge1}$ of independent, identically distributed random variables having uniform distribution on $[0, 1)$ satisfies

<footnote>*Keywords and phrases:* Pair correlation, Riemann zeta function, Lattice points, Diophantine inequality</footnote>

<footnote>*Mathematics Subject Classification:* Primary 11K06, 11J83, 11M06; Secondary 11B05, 11J25, 11J71</footnote>
(1) almost surely as a consequence of the Glivenko–Cantelli theorem; thus a deterministic sequence $(y_n)_{n \ge 1}$ which is uniformly distributed mod 1 exhibits the same behavior as a typical realization of a random sequence.

A sequence $(y_n)_{n \ge 1}$ is said to have Poissonian pair correlation if for all real numbers $s \ge 0$,

$$ \lim_{N \to \infty} \frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} \mathbf{1}_{[-s/N, s/N]}(y_n - y_m) = 2s. $$

This notion is motivated by questions from theoretical physics, and plays a key role in the Berry–Tabor conjecture; see [Mar01] for more information. Just like equidistribution, Poissonian pair correlation can be seen as a pseudo-randomness property, since a random sequence $(Y_n)_{n \ge 1}$ as above almost surely has Poissonian pair correlation. However, clearly the two properties are of a rather different nature. While equidistribution is a “large-scale” statistic (where the test interval always remains the same), the pair correlation is a highly localized statistic (where the size of the test interval shrinks in proportion with $N$). Note that the two properties are not independent: it is known that a sequence having Poissonian pair correlation necessarily must be equidistributed [ALP18, GL17, Mar20], whereas the opposite implication is generally false. An illustrative example is the sequence $(n\alpha)_{n \ge 1}$, which is equidistributed if and only if $\alpha \notin \mathbb{Q}$, but which fails to have Poissonian pair correlation for any $\alpha$ (see [LS20] for a more general result along those lines).

The theory of uniform distribution modulo one can be said to be relatively well understood (at least in the one-dimensional case). Many specific sequences are known which are uniformly distributed mod one. In contrast, only very few specific results are known in the pair correlation setting. A notable exception is the sequence $(\sqrt{n})_{n \in \mathbb{Z}_{\ge 1} \setminus \square}$, which is known to have Poissonian pair correlation [EMV15]. The sequence $(n^2\alpha)_{n \ge 1}$ is conjectured to have Poissonian pair correlation under mild Diophantine assumptions on $\alpha$, but only partial results are known in this direction [Hea10, MY18, RSZ01, Tru10]. Lacking specific examples, it is natural to turn to a metric theory instead. Let $(a_n)_{n \ge 1}$ be a sequence of distinct integers, let $\alpha \in \mathbb{R}$, and consider sequences of the form $(a_n\alpha)_{n \ge 1}$. The metric theory of such sequences with respect to equidistribution is very simple: for every such $(a_n)_n$, the sequence $(a_n\alpha)_n$ is uniformly distributed mod 1 for almost all $\alpha$ [Wey16]. The situation with respect to pair correlation is much more delicate. Pioneering work in this area was carried out by Rudnick, Sarnak and Zaharescu [RS98, RZ99]. As noted above, $(n\alpha)_{n \ge 1}$ does not have Poissonian pair correlation for any $\alpha$. However, for any polynomial $p \in \mathbb{Z}[X]$ of degree at least 2, the pair correlation of $(p(n)\alpha)_n$ is Poissonian for almost all $\alpha$. For related results, see for example [BPT01, CLZ15, RZ02].

Recently, a simple criterion was established in [ALL17] which allows to decide whether the sequence $(a_n\alpha)_n$ has Poissonian pair correlation for almost all $\alpha$ for many naturally arising integer sequences $(a_n)_n$. Let $E_N$ denote the number of
solutions $(n_1, n_2, n_3, n_4)$ of the equation

$$a_{n_1} - a_{n_2} + a_{n_3} - a_{n_4} = 0, \qquad (2)$$

subject to $1 \le n_1, n_2, n_3, n_4 \le N$. This quantity is called the *additive energy* in the additive combinatorics literature (see [Gry13, TV06]). Note that $N^2 \le E_N \le N^3$ for every $(a_n)_n$ and every $N$. The criterion is as follows. If a sequence $(a_n)_n$ satisfies $E_N \ll N^{3-\epsilon}$ for some $\epsilon > 0$, then $(a_n\alpha)_n$ has Poissonian pair correlation for almost all $\alpha$. If in contrast $E_N \gg N^3$, then the conclusion fails to be true. For further refinements of this criterion, and for remaining open problems, see [ALT19, BCGW18, BW20, LS20]. We emphasize that all that was written in this paragraph requires $(a_n)_n$ to be a sequence of integers.

Very little is known in the metric theory of pair correlation of sequences $(x_n\alpha)_n$ when $(x_n)_n$ is a sequence of reals rather than integers. One step in this general direction is [CLZ15], where $(x_n)_n$ is allowed to take rational values and the results obtained depend on the size of the denominators of these rationals. A general result was obtained recently in [RT20], where the authors gave a criterion formulated in terms of the number of solutions of a certain Diophantine inequality. The criterion is as follows: for a sequence $(x_n)_n$, assume that there exist some $\epsilon > 0$ and $\delta > 0$ such that the number of integer solutions $(n_1, n_2, n_3, n_4, j_1, j_2)$ of the equation

$$|j_1(x_{n_1} - x_{n_2}) - j_2(x_{n_3} - x_{n_4})| < N^{\epsilon}, \qquad (3)$$

subject to $1 \le |j_1|, |j_2| \le N^{1+\epsilon}$, $1 \le n_1, n_2, n_3, n_4 \le N$, $n_1 \ne n_2$, $n_3 \ne n_4$, is of order $\ll N^{4-\delta}$, then $(x_n\alpha)_n$ has Poissonian pair correlation for almost all $\alpha$. It is verified in [RT20] that this condition is satisfied for lacunary sequences. A condition in the spirit of (3) arises very naturally when studying this sort of problem (cf. also [RZ99]); we will encounter a variant of this condition in Equation (16) below. In particular, it is very natural that in the integer case one has to count solutions of Diophantine equations, while in the real-number setting one has to count solutions of Diophantine inequalities. The problem with (3) is that it is in general rather difficult to verify whether this condition is satisfied for a given sequence or not, with issues being caused in particular by the presence of the coefficients $j_1$ and $j_2$. The purpose of the present paper is to give a simplified criterion, in the spirit of the criterion of [ALL17] which was specified in terms of the number of solutions of the equation (2).

**Theorem 1.** Let $(x_n)_{n \ge 1}$ be a sequence of positive real numbers for which there exists a constant $c > 0$ such that $x_{n+1} - x_n \ge c$, $n \ge 1$. Let $E_N^*$ denote the number of solutions $(n_1, n_2, n_3, n_4)$ of the inequality

$$|x_{n_1} - x_{n_2} + x_{n_3} - x_{n_4}| < 1, \qquad (4)$$

subject to $n_i \le N$, $i = 1, 2, 3, 4$. Assume that there exists some $\delta > 0$ such that $E_N^* \ll N^{183/76-\delta}$ as $N \to \infty$. Then the sequence $(x_n\alpha)_{n \ge 1}$ has Poissonian pair correlation for almost all $\alpha \in \mathbb{R}$.
The exponent $183/76 \approx 2.408$ in the conclusion of the theorem comes from a bound for the 178/13th moment of the Riemann zeta function on the critical line due to Ivić [Ivi03], building on earlier work of Heath-Brown [Hea78]. Conditionally under the Lindelöf hypothesis, our bound for $E_N^*$ can be relaxed to $E_N^* \ll N^{3-\epsilon}$ for any $\epsilon > 0$, which would be in accordance with the results known for the integer case.

Theorem 1 applies, for example, to all sequences of the form $x_n = p(n)$, $n \ge 1$, where $p$ is a quadratic polynomial with real coefficients. For such a sequence $(x_n)_n$ we have $E_N^* \ll N^{2+\epsilon}$ for any $\epsilon > 0$ by Lemma 5.2 of [BKW10]. Theorem 1 also applies to $x_n = p(n)$ for every polynomial $p \in \mathbb{R}[X]$ of degree $d \ge 3$, under the additional assumption that the coefficient of $x^{d-1}$ is rational<fnref n="1" />; the required bound for $E_N^*$ then follows, after eliminating this coefficient, from Lemma 7 below (with the choice of $\theta = d$ and $\gamma = N^{d-2}$). The extra assumption on the second coefficient is most likely redundant, but we have not been able to establish the necessary bound for $E_N^*$ without it. A famous open conjecture in additive combinatorics asserts that $E_N \ll N^{2+\epsilon}$ for all convex sequences $(x_n)_n$, which would provide many further applications of our theorem; however, unfortunately the best current bound in this direction (Shkredov's 32/13 $\approx 2.46$ from [Shk13]) is just beyond the range of applicability of our theorem.

Bounding the number of solutions of (4) is necessary to control the variance of the pair correlation function. When carefully reading the Proof of Theorem 1 it becomes visible that not all solutions of (4) contribute equally to the variance, but that rather a 4-tuple $(n_1, n_2, n_3, n_4)$ with $x_{n_1} - x_{n_2} + x_{n_3} - x_{n_4} = \gamma$ for some $\gamma \in (-1, 1)$ has a stronger effect on the variance the smaller the absolute value of $\gamma$ is. This suggests to consider the quantity $E_{N,\gamma}^*$, which is defined as the number of solutions $(n_1, n_2, n_3, n_4)$ of the inequality

$$|x_{n_1} - x_{n_2} + x_{n_3} - x_{n_4}| < \gamma,$$

for $\gamma \in (0, 1]$ and subject to $n_i \le N$, $i = 1, 2, 3, 4$. Very informally speaking, one might expect that $E_{N,\gamma}^*$ scales as $E_{N,\gamma}^* \approx \gamma E_N^*$ for a “randomly behaved” real sequence $(x_n)_n$, except for the contribution of the trivial solutions $n_1 = n_2$ and $n_3 = n_4$ which always is of order $N^2$. The following theorem states that being able to control $E_{N,\gamma}^*$ as a function of $\gamma$ indeed allows us to deduce metric pair correlations in some cases where the condition on the additive energy in Theorem 1 fails to hold.

**Theorem 2.** Let $(x_n)_{n \ge 1}$ be a sequence of positive real numbers for which there exists $c > 0$ such that $x_{n+1} - x_n \ge c$, $n \ge 1$. Assume that there exists some $\delta > 0$ such that for all $\eta > 0$ we have

$$E_{N,\gamma}^* \ll_{\eta,\delta} N^{2+\eta} + \gamma N^{3-\delta} \quad (5)$$

as $N \to \infty$, uniformly for $\gamma \in (0, 1]$. Then the sequence $(x_n\alpha)_{n \ge 1}$ has Poissonian pair correlation for almost all $\alpha \in \mathbb{R}$.

<footnote n="1">¹ This could be relaxed to assuming some Diophantine condition on this coefficient.</footnote>
**Theorem 2** is tailor-made for an application to the sequence $x_n = n^\theta$. For that sequence, Equation (5) holds provided $\theta > 1$ as a consequence of Lemma 7 below. Since that particular problem was a key driving force for writing the present paper, we formulate this conclusion as a theorem rather than just as a corollary.

**Theorem 3.** For every real number $\theta > 1$ and almost every $\alpha \in \mathbb{R}$, the sequence $(n^\theta \alpha)_{n \ge 1}$ has Poissonian pair correlation.

As noted above, the conclusion of Theorem 3 is not true when $\theta = 1$. It seems plausible that the conclusion of the theorem is valid again for $0 < \theta < 1$. However, this cannot be proved with the methods used in the present paper, which break down in the case of a sequence $(x_n)_n$ whose order of growth is only linear or even slower. We will address this aspect at the very end of the paper, where we also formulate some further open problems.

In conclusion we note that Technau and Yesha recently obtained a result which is somewhat similar to our Theorem 3, but which is “metric” in the exponent rather than in a multiplicative parameter. More precisely, they showed that $(n^\theta)_n$ has Poissonian pair correlation for almost all $\theta > 7$. Their paper also contains similar results on higher correlations, which require a larger value of $\theta$. From a technical perspective, their problem is rather different from ours. For details see their paper [TY20].

## 2 Preliminaries

As in the introduction, let $\mathbf{1}_{[-s/N, s/N]}(x)$ denote the indicator function of the interval $[-s/N, s/N]$, extended with period 1. That is,

$$ \mathbf{1}_{[-s/N, s/N]}(x) = \begin{cases} 1 & \text{if } x - \langle x \rangle \in [-s/N, s/N], \\ 0 & \text{otherwise,} \end{cases} $$

where $\langle x \rangle$ denotes the nearest integer to $x$. We wish to show that under the assumptions of Theorem 1 we have for almost all $\alpha \in \mathbb{R}$

$$ \frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} \mathbf{1}_{[-s/N, s/N]}(x_m \alpha - x_n \alpha) \to 2s \quad (6) $$

as $N \to \infty$ for all $s \ge 0$. It is well-known that for any $s$ and $N$, and for any positive integer $K$ there exist trigonometric polynomials $f_{K,s,N}^+(x)$ and $f_{K,s,N}^-(x)$ of degree at most $K$ such that

$$ f_{K,s,N}^-(x) \le \mathbf{1}_{[-s/N, s/N]}(x) \le f_{K,s,N}^+(x) \quad (7) $$

for all $x$, and such that

$$ \int_0^1 f_{K,s,N}^{\pm}(x) dx = 2s/N \pm \frac{1}{K+1}. \quad (8) $$
Furthermore, the jth Fourier coefficient $c_j$ of $f_{K,s,N}^{-}$ satisfies

$$|c_j| \le \min \left( \frac{2s}{N}, \frac{1}{\pi|j|} \right) + \frac{1}{K+1} \qquad (9)$$

for all j, and an analogous bound holds for the Fourier coefficients of $f_{K,s,N}^{+}$. These trigonometric polynomials are called Selberg polynomials, and their construction is described in detail in Chapter 1 of [Mon94].

Instead of establishing the required convergence relation (6) for indicator functions, we will rather work with the trigonometric polynomials $f_{K,s,N}^{+}$ and $f_{K,s,N}^{-}$ instead, which is technically more convenient. More precisely, in order to obtain (6) it suffices to prove the following. For every fixed positive integer $r$, and for every fixed real number $s \ge 0$, we have

$$\frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} f_{rN, s, N}^{+} (x_m \alpha - x_n \alpha) \sim N \int_{0}^{1} f_{rN, s, N}^{+} (x) dx \qquad (10)$$

as $N \to \infty$, for almost all $\alpha \in \mathbb{R}$, and the same is true when $f^{+}$ is replaced by $f^{-}$. The desired result for indicator functions then follows from (7) and (8) and letting $r \to \infty$ (see [RS98, RT20] for more details).

To establish (10) we prove that the “expected value” (with respect to $\alpha$) of the left-hand side is asymptotic to the right-hand side, and that the “variance” of the left-hand side of (10) is not too large. An application of Chebyshev’s inequality together with the Borel–Cantelli lemma then gives the desired result. As usual in such problems, controlling the expectation is easier than controlling the variance.

We will obtain the required bound for the expectation in Section 3, and the bound for the variance in Sections 4 and 5. In Section 6 we conclude the Proof of Theorem 1. Section 7 contains all of the necessary modifications for the Proof of Theorem 2, and in Section 8 we show that the sequence $(n^{\theta})_{n}$ indeed allows an application of Theorem 2. Finally, in Section 9 we discuss limitations of our method, and outline open problems and directions for future research.

### 3 Proof of Theorem 1: Controlling expectations

Throughout the argument we assume that a positive integer $r$ and a positive real $s$ are fixed. We write $f_{N}$ for the function $f_{rN, s, N}^{+}$, as defined in the previous section (or for the function $f_{rN, s, N}^{-}$—both cases work in exactly the same way). We want to control the “expected value” with respect to $\alpha$ of the left-hand side of (10) as $N \to \infty$. In the case when $(x_n)_{n \ge 1}$ is an integer sequence everything is periodic with period 1, and it is appropriate to integrate over $\alpha \in [0, 1]$ with respect to the Lebesgue measure. In our case, when $(x_n)_{n \ge 1}$ is a sequence of reals, we do not have such periodicity. We thus have to integrate over all $\alpha \in \mathbb{R}$ with respect to an appropriate measure $\mu$, which is absolutely continuous with respect to the
Lebesgue measure, so that a $\mu$-almost everywhere conclusion implies a Lebesgue-almost everywhere conclusion. A good choice for the measure $\mu$ is the measure whose density with respect to the Lebesgue measure is given by

$$d\mu(x) = \frac{2(\sin(x/2))^2}{\pi x^2} dx. \qquad (11)$$

The Fourier transform of $x \mapsto \frac{2(\sin(x/2))^2}{\pi x^2}$ is a non-negative real function which is supported on the interval $(-1, 1)$, and which is uniformly bounded by $1/\sqrt{2\pi}$. Note that the measure $\mu$ is normalized such that $\mu(\mathbb{R}) = 1$.

Expanding the function $f_N(x)$ into a Fourier series

$$\sum_{j \in \mathbb{Z}} c_j e^{2\pi i j x},$$

by construction we have $c_j = 0$ when $|j| > rN$, and $|c_j| \le 2s/N + 1/(rN) \ll N^{-1}$ for all $j$ (recall that $r$ and $s$ are assumed to be fixed). Moreover we have $c_0 = \int_0^1 f_N(x)dx$. Using the fact that the Fourier transform of the measure $\mu$ is supported on $(-1, 1)$ and uniformly bounded, we obtain

$$\left| \int_{\mathbb{R}} \frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} f_N(x_m \alpha - x_n \alpha) \, d\mu(\alpha) - N \int_0^1 f_N(x) dx \right| \\ \le \underbrace{\left(N - \frac{N(N-1)}{N}\right)}_{=1} \underbrace{\int_0^1 f_N(x) dx}_{\ll N^{-1}} \\ + \frac{1}{N} \sum_{1 \le |j| \le rN} |c_j| \left| \int_{\mathbb{R}} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} e^{2\pi i j(x_m \alpha - x_n \alpha)} \, d\mu(\alpha) \right| \\ \ll N^{-1} + N^{-2} \underbrace{\sum_{1 \le |j| \le rN} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} \mathbf{1}(|j(x_m - x_n)| < 1)}_{\ll N \text{ due to the growth assumption on } (x_n)_{n \ge 1}} \\ \ll N^{-1}, \qquad (12)$$

where we estimated $|c_j|$ using (9). Thus we have

$$\int_{\mathbb{R}} \frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} f_N(x_m \alpha - x_n \alpha) \, d\mu(\alpha) = N \int_0^1 f_N(x) dx + O(1/N), \qquad (13)$$

as desired.

Controlling the variances is more difficult, and will be done in the next two sections.
# 4 Proof of Theorem 1: Controlling variances

We keep the setup as in Section 3 above, that is, we assume that $r$ and $s$ are fixed, and we write $f_N$ for either $f_{rN,s,N}^+$ or $f_{rN,s,N}^-$. Furthermore, we write $h_N$ for the centered version of $f_N$, that is, for the function

$$h_N(x) = f_N(x) - \int_0^1 f(x) \, dx = \sum_{\substack{j \in \mathbb{Z}, \\ j \neq 0}} c_j e^{2\pi i jx}. \quad (14)$$

We wish to estimate the “variance” of our localized counting function, or more precisely the quantity

$$\mathrm{Var}(h_N, \mu) := \int_{\mathbb{R}} \left( \frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} h_N(x_m \alpha - x_n \alpha) \right)^2 d\mu(\alpha). \quad (15)$$

The following bound on $\mathrm{Var}(h_N, \mu)$ is the crucial ingredient in our Proof of Theorem 1.

**LEMMA 1.** For every $\varepsilon > 0$ we have, as $N \to \infty$,

$$\mathrm{Var}(h_N, \mu) \ll \max \left( N^{-\varepsilon/8} + N^{-183/89+3\varepsilon} (E_N^*)^{76/89}, E_N^* N^{-2.49+4\varepsilon} \right).$$

For the convenience of the reader, we note at this point that our assumption that there is some $\delta > 0$ such that $E_N^* \ll N^{183/76-\delta}$ ensures together with Lemma 1 that there is some $\delta' > 0$ such that $\mathrm{Var}(h_N, \mu) \ll N^{-\delta'}$, which is sufficient to deduce Theorem 1 (see Section 6 for details). We also note that conditionally under the Lindelöf hypothesis the bound which follows from our method is $\mathrm{Var}(h_N, \mu) \ll N^{-3+\varepsilon} E_N^*$.

# 5 Proof of Lemma 1: Lattice point counting via the Riemann zeta function

## 5.1 A first reduction.
Squaring out in (15) and using again the properties of the Fourier transform of the measure $\mu$, we can bound $\mathrm{Var}(h_N, \mu)$ by

$$\int_{\mathbb{R}} \frac{1}{N^2} \sum_{\substack{1 \le n_1, n_2, n_3, n_4 \le N, \\ n_1 \ne n_2, n_3 \ne n_4}} \sum_{\substack{j_1, j_2 \in \mathbb{Z}, j_1, j_2 \ne 0 \\ |j_1|, |j_2| \le rN}} \underbrace{|c_{j_1} c_{j_2}|}_{\ll N^{-2}} e^{2\pi i \alpha (j_1(x_{n_1} - x_{n_2}) - j_2(x_{n_3} - x_{n_4}))} d\mu(\alpha) \\ \ll \frac{1}{N^4} \sum_{\substack{1 \le n_1, n_2, n_3, n_4 \le N, \\ n_1 > n_2, n_3 > n_4}} \sum_{1 \le j_1, j_2 \le rN} \mathbf{1}(|j_1(x_{n_1} - x_{n_2}) - j_2(x_{n_3} - x_{n_4})| < 1),$$

thereby essentially arriving at (3). For technical reasons, in this paper we prefer to localize the variables $j_1, j_2$ into dyadic regions and thus apply the Cauchy-Schwarz
inequality to (15). To simplify later formulas we also replace the differences $x_{n_1}-x_{n_2}$ and $x_{n_3}-x_{n_4}$ by their respective absolute values using the parity of $h_N$. Then, writing $U$ for the smallest integer for which $2^U \ge rN$, we can bound $\text{Var}(h_N, \mu)$ by

$$ 
\begin{aligned}
& \int_{\mathbb{R}} \left( \frac{1}{N} \sum_{u=1}^{U} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} \sum_{2^{u-1} \le |j| < 2^u} c_j e^{2\pi i j |x_m - x_n| \alpha} \right)^2 d\mu(\alpha) \\
& \ll \frac{1}{N^2} \int_{\mathbb{R}} \left( \sum_{k=1}^{U} 1 \right) \left( \sum_{u=1}^{U} \left| \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} \sum_{2^{u-1} \le |j| < 2^u} c_j e^{2\pi i j |x_m - x_n| \alpha} \right| \right)^2 d\mu(\alpha) \\
& \ll \frac{\log N}{N^4} \sum_{u=1}^{U} \sum_{\substack{1 \le n_1, n_2, n_3, n_4 \le N, \\ n_1 \ne n_2, n_3 \ne n_4}} \sum_{2^{u-1} \le j_1, j_2 < 2^u} \mathbf{1} \left( |j_1| |x_{n_1} - x_{n_2}| - |j_2| |x_{n_3} - x_{n_4}| < 1 \right).
\end{aligned}
\quad (16) $$

Thus we have reduced the problem of estimating the variance to a problem of bounding the number of solutions of a Diophantine inequality.

**5.2 Counting solutions by using the Riemann zeta function.** We will relate the counting problem in Equation (16) to the problem of bounding a twisted moment of the Riemann zeta function. Before we return to the proof, we point out the difference between the real-number case (in this paper) and the corresponding results for the case of $(x_n)_{n \ge 1}$ being an integer sequence. In the integer case, the problem of estimating the variance of the pair correlation function can be reduced to counting solutions of $j_1(x_{n_1} - x_{n_2}) = j_2(x_{n_3} - x_{n_4})$. Note that this is in accordance with the situation in the present paper, where we count solutions to $|j_1(x_{n_1} - x_{n_2}) - j_2(x_{n_3} - x_{n_4})| < 1$, with the difference that in the integer case “< 1” implies “= 0”. The number of solutions of the counting problem in the integer case is essentially governed by what is called a GCD sum. It is known that such sums have a connection with the Riemann zeta function (see [Ais16, Hil09]), and strong estimates for such sums were obtained in [ABS15, BS17, BT19]. Our approach below is motivated by a beautiful argument of Lewko and Radziwiłł [LR17], who showed how the relevant GCD sum can be estimated in terms of a twisted moment of a random model of the Riemann zeta function on the critical line.<fnref n="2" /> The randomization was crucial in their argument for different reasons, one being that the required distributional estimates for extreme values of the actual Riemann zeta function are not known unconditionally. Their argument relied crucially on the fundamental theorem of arithmetic, and thus on the fact that they were dealing with integer sequences. In the real-number case the situation is much more delicate. We will relate

<footnote n="2">² See also [BMT21] for links between twisted moments of character sums and GCD sums. Furthermore, see [Shk20] for a very recent paper of Shkredov, where he applies GCD sums and methods from [LR17] to give upper bounds for the maximal length of arithmetic progressions contained in sets with small product set.</footnote>
our counting problem to a convolution formula for the Riemann zeta function. The kernel will be chosen for its good properties with respect to the Fourier transform (positivity and localized support) which allow to overcount without substantial loss. To summarize, in our argument below we will use of a combination of ideas from [Ais16, BS17, BS18, BT19] and [LR17].

Let $(x_n)_{n \ge 1}$ be the sequence from the statement of Theorem 1. Let $M = N^2 - N$, and let $\{z_1, \dots, z_M\}$ be the multi-set of all absolute differences $\{|x_m - x_n| : 1 \le m, n \le N, m \ne n\}$, meaning that we allow repetitions in the definition. For a given positive integer $u$ with $2^u \le 2rN$, we wish to estimate

$$ \sum_{1 \le m, n \le M} \sum_{2^{u-1} \le j_1, j_2 < 2^u} \mathbf{1}(|j_1 z_m - j_2 z_n| < 1). \quad (17) $$

We write $\zeta(\sigma + it)$ for the Riemann zeta function. We also write $\Phi(t) = e^{-t^2/2}$, and note that this function has a positive Fourier transform given by $\Phi = \sqrt{2\pi}\Phi$. Throughout the proof $\varepsilon > 0$ is a small constant, and we take it for granted that $N$ is “large”.

Our argument proceeds by splitting into two cases depending on the size of $\min\{z_m, z_n\}$. We first treat the case when $z_m, z_n$ are both at least of size $N^{1.01}$. We then treat the case when one of $z_m$ or $z_n$ is “small”, which because of our dyadic splitting essentially amounts to saying that both variables are small.

• Case 1: Counting solutions for $z_m, z_n \ge N^{1.01}$.

Let $u$ be given such that $2^{u-1} \le j_1, j_2 \le 2^u$. Set $T = 2^u N^{1+\varepsilon/2}$. For any integer $k \ge N^{1.01}$, we set

$$ b_k = \sum_{m=1}^{M} \mathbf{1}(z_m \in [k, k+1)), \quad (18) $$

while for $k < N^{1.01}$ we set $b_k = 0$. Clearly we have

$$ \sum_{k=1}^{\infty} b_k \le M = N^2 - N. $$

We split the interval $[1, \infty)$ into a disjoint union of intervals $I_h$ for $h \ge 0$, where

$$ I_h = \left[ \left\lfloor \left( 1 + \frac{1}{T} \right)^h \right\rfloor, \left\lceil \left( 1 + \frac{1}{T} \right)^{h+1} \right\rceil \right), $$

and we set

$$ a_h = \left( \sum_{k \in I_h} b_k^2 \right)^{1/2}, \quad h \ge 0. \quad (19) $$

Finally, we define a function

$$ P(t) = \sum_{h=0}^{\infty} a_h \left(1 + \frac{1}{T}\right)^{iht}. \quad (20) $$

The following four lemmas correspond to the key steps in our Case 1 analysis.
LEMMA 2 (Controlling the square-integral of $P$ in terms of the additive energy).
We have

$$ \int_{\mathbb{R}} |P(t)|^2 \Phi(t/T) dt \ll T E_N^*. \quad (21) $$

LEMMA 3 (Counting solutions of Diophantine inequalities in terms of $a_h$). We have

$$ \sum_{2^{u-1} \le j_1, j_2 < 2^u} \sum_{1 \le m, n \le M} 1(|j_1 z_m - j_2 z_n| < 1) \ll \sum_{2^{u-1} \le j_1, j_2 \le 2^u} \sum_{\substack{h_1, h_2 \ge 0, \\ |(1+\frac{1}{T})^{h_1-h_2} - \frac{j_2}{j_1}| \le \frac{4}{T}}} a_{h_1} a_{h_2}. $$

The next step, in Lemma 4 below, is to relate the sum on the right-hand side of the equation above to a complex integral, where the term $\sum_{j_1, j_2 \ge 1} (j_1 j_2)^{-1/2} (j_1 / j_2)^{it}$ can be roughly interpreted as $|\zeta(1/2 + it)|^2$. However, instead of simply using a truncated expression or an approximate functional equation for $\zeta(1/2 + it)$, we will rather use of a convolution formula to improve the analysis near $t = 0$. To do so, we introduce the function $K$ defined by

$$ K(u) := \frac{\sin^2((1 + \varepsilon/4)u \log N)}{\pi u^2(1 + \varepsilon/4)(\log N)}, $$

whose Fourier transform is given by

$$ \hat{K}(\xi) = \max \left( 1 - \frac{|\xi|}{2(1 + \varepsilon/4) \log N}, 0 \right). $$

A similar idea was also fruitfully used in a paper of Bondarenko and Seip [BS18]. The function $K$ is chosen in such a way that we have $\hat{K}(\log j_1 j_2) \gg 1 - \frac{2(\log r N)}{2(1+\varepsilon/4)\log N} \gg 1$ (where we suppress the dependence on the constants $\varepsilon$ and $r$).

LEMMA 4 (Counting solutions of the Diophantine inequality by complex integration). We have

$$ \sum_{2^{u-1} \le j_1, j_2 < 2^u} \sum_{\substack{h_1 \ge 0, h_2 \ge 0, \\ |(1+\frac{1}{T})^{h_1-h_2} - \frac{j_2}{j_1}| \le \frac{4}{T}}} a_{h_1} a_{h_2} \\ \ll \frac{2^u}{T} \int_{\mathbb{R}} \sum_{j_1, j_2 \ge 1} \frac{\hat{K}(\log j_1 j_2)}{(j_1 j_2)^{1/2}} \left(\frac{j_1}{j_2}\right)^{it} |P(t)|^2 \Phi(t/T) dt. $$

LEMMA 5 (Estimating the complex integral). We have

$$ \frac{2^u}{T} \int_{\mathbb{R}} \sum_{j_1, j_2 \ge 1} \frac{\hat{K}(\log j_1 j_2)}{(j_1 j_2)^{1/2}} \left(\frac{j_1}{j_2}\right)^{it} |P(t)|^2 \Phi(t/T) dt \ll N^{4-\varepsilon/4} + N^{173/89+2\varepsilon} (E_N^*)^{76/89}. \quad (22) $$

We now prove these four lemmas.
*Proof of Lemma 2.* We have

$$
\begin{align*}
\int_{\mathbb{R}} |P(t)|^2 \Phi(t/T) dt &= \int_{\mathbb{R}} \sum_{h_1, h_2 \ge 0} a_{h_1} a_{h_2} \left(1 + \frac{1}{T}\right)^{(h_1-h_2)it} \Phi(t/T) dt \\
&= T \sum_{h_1, h_2 \ge 0} a_{h_1} a_{h_2} \int_{\mathbb{R}} \exp\left(\left(\log\left(1 + \frac{1}{T}\right)\right) T(h_1 - h_2)iy\right) \Phi(y) dy \\
&= T \sum_{h_1, h_2 \ge 0} a_{h_1} a_{h_2} \hat{\Phi}\left(\left(\log\left(1 + \frac{1}{T}\right)\right) T(h_1 - h_2)\right) \\
&\ll T \sum_{h_1, h_2 \ge 0} a_{h_1} a_{h_2} \hat{\Phi}\left(\frac{h_1 - h_2}{2}\right) \tag{23}
\end{align*}
$$

$$
\begin{align*}
&\ll T \sum_{h=0}^{\infty} a_h^2 \\
&\ll T \sum_{k=1}^{\infty} b_k^2 \ll T \sum_{\substack{m,n=1 \\ |z_m-z_n|<1}}^{M} 1 = T E_N^*.
\end{align*}
\tag{24}
$$

Here we used that $T \log(1 + 1/T) \ge 1/2$ for sufficiently large $N$ (note that large $N$ implies large $T$), and the Cauchy-Schwarz inequality together with the rapid decay of $\hat{\Phi}$ to pass from (23) to (24). We will further comment on the construction of $P(t)$ at the very end of our Case 1 analysis. $\square$

*Proof of Lemma 3.* Let $j_1$ and $j_2$ be fixed, and assume without loss of generality that $j_1 \ge j_2$. Let $k \ge N^{1.01}$ be an integer in $I_{h_1}$, and assume that $z_m \in [k, k+1)$. Then the inequality $|j_1 z_m - j_2 z_n| < 1$ is only possible when

$$
\left| \left[ \frac{j_1 k}{j_2} \right] - z_n \right| < 4 \qquad (25)
$$

(recall that $j_1/j_2 \le 2$ because $j_1, j_2$ are located in the same dyadic interval). We write $\ell(k) = \lfloor j_1 k / j_2 \rfloor$. Recall that $j_1/j_2 \ge 1$ by assumption, so the mapping $k \mapsto \ell(k)$ is injective. Thus we have

$$
\begin{align*}
& \sum_{z_m \in I_{h_1}, z_n \in I_{h_2}} 1(|j_1 z_m - j_2 z_n| < 1) \\
& \ll \sum_{k \in I_{h_1}} \sum_{z_m \in [k, k+1)} \sum_{\substack{z_n \in I_{h_2}, \\ |\ell(k) - z_n| < 4}} 1 \\
& \ll \sum_{k \in I_{h_1}} \sum_{z_m \in [k, k+1)} \sum_{\substack{z_n \in I_{h_2}, \\ z_n \in [\ell(k) + v, \ell(k) + v+1]}} 1 \\
& \ll \sum_{-4 \le v \le 3} \sum_{k \in I_{h_1} \text{ such that } \ell(k) + v \in I_{h_2}} b_k b_{\ell(k)+v}
\end{align*}
$$
$$ \begin{aligned} & \ll \left( \sum_{k \in I_{h_1}} b_k^2 \right)^{1/2} \left( \sum_{\ell \in I_{h_2}} b_\ell^2 \right)^{1/2} \\ & \ll a_{h_1} a_{h_2} \end{aligned} \qquad (26) $$

by Cauchy–Schwarz.

When $j_1$ and $j_2$ are fixed, there can only be solutions of $|j_1z_m - j_2z_n| < 1$ with $z_m \in I_{h_1}$ and $z_n \in I_{h_2}$ for particular pairs $(h_1, h_2)$. Assume that $z_m \in I_{h_1}$ and $z_n \in I_{h_2}$ such that $|j_1z_m - j_2z_n| < 1$. Recall that $j_1 \ge j_2$ by assumption, so we have $\left| \frac{z_m}{z_n} - \frac{j_2}{j_1} \right| < \frac{1}{j_1z_n}$ and consequently $\frac{z_m}{z_n} \le \frac{j_2}{j_1} + \frac{1}{j_1z_n} \le 2$. Since $z_m \in I_{h_1}$ and $z_n \in I_{h_2}$, the quotient $z_m/z_n$ is somewhere between $(1+1/T)^{h_1-h_2-1}$ and $(1+1/T)^{h_1-h_2+1}$, so that

$$ \frac{z_m}{z_n} \le \underbrace{(1+1/T)^{h_1-h_2}}_{\le 2(1+1/T) \le 3, \text{ since } z_m/z_n \le 2} (1+1/T) \le (1+1/T)^{h_1-h_2} + \frac{3}{T}. \qquad (27) $$

Similarly

$$ \frac{z_m}{z_n} \ge (1 + 1/T)^{h_1 - h_2} - \frac{3}{T}. \qquad (28) $$

Since $j_1 \ge 2^{u-1}$ and $z_n \ge N^{1.01}$ by assumption, we have $\left| \frac{z_m}{z_n} - \frac{j_2}{j_1} \right| \le \frac{1}{2^{u-1}N^{1.01}} \le \frac{1}{T}$, where the last inequality follows from our choice of $T$. Overall, together with (27) and (28) this shows that the inequality $|j_1z_m - j_2z_n| < 1$ for $z_m \in I_{h_1}$ and $z_n \in I_{h_2}$ is only possible when

$$ \left| \left( 1 + \frac{1}{T} \right)^{h_1 - h_2} - \frac{j_2}{j_1} \right| \le \frac{4}{T}. $$

Note that, for fixed $j_1, j_2$, this is an inequality which only depends on $h_1, h_2$ and not on $z_m, z_n$ anymore. Thus in combination with (26) we obtain

$$ \sum_{1 \le m, n \le M} \mathbf{1}(|j_1 z_m - j_2 z_n| < 1) \ll \sum_{\substack{h_1, h_2 \ge 0, \\ |(1+\frac{1}{T})^{h_1-h_2} - \frac{j_2}{j_1}| \le \frac{4}{T}}} a_{h_1} a_{h_2} $$

for all fixed $j_1$ and $j_2$. When summing over $j_1$ and $j_2$, we obtain the conclusion of Lemma 3. □
*Proof of Lemma 4.* By the properties of $\hat{K}$ and $\hat{\Phi}$ we have

$$ 
\begin{aligned}
& \sum_{2^{u-1} \le j_1, j_2 < 2^u} \frac{1}{(j_1 j_2)^{1/2}} \sum_{\substack{h_1 \ge 0, h_2 \ge 0, \\ \left|\left(1+\frac{1}{T}\right)^{h_1-h_2} \frac{j_2}{j_1}\right| \le \frac{4}{T}}} a_{h_1} a_{h_2} \\
& \ll \sum_{j_1, j_2 \ge 1} \frac{\hat{K}(\log j_1 j_2)}{(j_1 j_2)^{1/2}} \sum_{h_1, h_2 \ge 0} a_{h_1} a_{h_2} \hat{\Phi}\left(T \log\left(\frac{j_1}{j_2}(1+1/T)^{h_1-h_2}\right)\right) \\
& \ll \frac{1}{T} \int_{\mathbb{R}} \sum_{j_1, j_2 \ge 1} \frac{\hat{K}(\log j_1 j_2)}{(j_1 j_2)^{1/2}} \left(\frac{j_1}{j_2}\right)^{it} |P(t)|^2 \Phi(t/T) dt. \quad (29)
\end{aligned}
$$

Note that we crucially used the fact that in all three lines of the displayed equation above, all terms in the summations are non-negative, because $\hat{K}, \hat{\Phi}$ and $\hat{\Phi}$ are all non-negative. Thus, using that $2^u \ll (j_1j_2)^{1/2} \ll 2^u$, we have

$$ 
\begin{aligned}
& \sum_{2^{u-1} \le j_1, j_2 \le 2^u} \sum_{\substack{h_1, h_2 \ge 0, \\ \left|\left(1+\frac{1}{T}\right)^{h_1-h_2} \frac{j_2}{j_1}\right| \le \frac{4}{T}}} a_{h_1} a_{h_2} \\
& \ll \frac{2^u}{T} \int_{\mathbb{R}} \sum_{j_1, j_2 \ge 1} \frac{\hat{K}(\log j_1 j_2)}{(j_1 j_2)^{1/2}} \left(\frac{j_1}{j_2}\right)^{it} |P(t)|^2 \Phi(t/T) dt,
\end{aligned}
$$

as claimed.

□

In order to prove Lemma 5, we need the following technical tool, which is Lemma 5.3 of [BT19].

LEMMA 6. Let $\sigma \in (-\infty, 1)$ and let $F$ be a holomorphic function in the strip $y = \Im z \in [\sigma - 2, 0]$, such that

$$ \sup_{\sigma-2 \le y \le 0} |F(x + iy)| \ll \frac{1}{x^2 + 1}. \quad (30) $$

Then for all $s = \sigma + it \in \mathbb{C}$, $t \ne 0$, we have

$$ 
\begin{aligned}
& \sum_{k, \ell \ge 1} \frac{\hat{F}(\log k\ell)}{k^s \ell^s} \\
& = \int_{\mathbb{R}} \zeta(s + iu) \overline{\zeta(s - iu)} F(u) du \\
& \quad + 2\pi\zeta(1 - 2it) F(is - i) + 2\pi\zeta(1 + 2it) F(is + i).
\end{aligned}
$$

The Proof of this lemma (Lemma 6) is only briefly sketched in [BT19]. However, a detailed proof of a similar lemma is given in [BS18, Lemma 1]; for the proof of our lemma one can exactly follow the argument given there, just using the function $f: z \mapsto \zeta(z + it)\zeta(z - it)K(i\sigma - iz)$ instead of the one considered there.
*Proof of Lemma 5.* For simplicity of writing we define

$$G(t) = \sum_{j_1, j_2 \ge 1} \frac{\hat{K}(\log j_1 j_2)}{(j_1 j_2)^{1/2}} \left(\frac{j_1}{j_2}\right)^{it}.$$

We note that $K$ can be extended analytically and satisfies assumption (30). Furthermore

$$|K(t - i/2)| \ll N^{1+\epsilon/4}/(t+1)^2, \quad |K(-t - i/2)| \ll N^{1+\epsilon/4}/(t+1)^2. \quad (31)$$

We first note that we have the pointwise bound

$$|P(t)|^2 \le |P(0)|^2 = \left(\sum_{h \ge 0} a_h\right)^2 \ll \left(\sum_{n=1}^{\infty} b_n\right)^2 \ll N^4, \quad t \in \mathbb{R}. \quad (32)$$

This allows us to see that

$$\int_{-1}^{1} G(t) |P(t)|^2 \Phi(t/T) dt \ll N^5$$

and thus we can restrict our integration domain to $|t| \ge 1$ below. By Lemma 6 we have

$$\int_{|t| \ge 1} G(t) |P(t)|^2 \Phi(t/T) dt = \text{Int}_1 + \text{Int}_2 + \text{Int}_3,$$

where

$$\text{Int}_1 = \int_{|t| \ge 1} |P(t)|^2 \Phi(t/T) \int_{\mathbb{R}} \zeta(1/2 + it + iu) \zeta(1/2 - it + iu) K(u) du dt,$$

$$\text{Int}_2 = 2\pi \int_{|t| \ge 1} \zeta(1 - 2it) K(-t - i/2) |P(t)|^2 \Phi(t/T) dt,$$

$$\text{Int}_3 = 2\pi \int_{|t| \ge 1} \zeta(1 + 2it) K(t - i/2) |P(t)|^2 \Phi(t/T) dt.$$

Using (31) together with the easy estimate $|\zeta(1 + it)| \ll \log t$, we obtain

$$\begin{aligned} \text{Int}_2 &\ll N^{1+\epsilon/4} N^4 \int_{t \ge 1} \frac{\log t}{t^2} dt \\ &\ll N^{5+\epsilon/4}. \end{aligned}$$

Exactly the same estimate holds for Int$_3$.

The classical convexity bound $|\zeta(1/2 + it)| \ll |t|^{1/4}$ gives $|\zeta(1/2 + it + iu)\zeta(1/2 - it + iu)| \ll (|t| + |u|)^{1/2} \ll |t|^{1/2} + |u|^{1/2}$. Hence we can bound the contribution of the domain $|u| \ge T$ to Int$_1$ by
$$ \int_{\mathbb{R}} |t|^{1/2} |P(t)|^2 \Phi(t/T) \underbrace{\int_{|u| \ge T} K(u) \, du \, dt}_{\ll T^{-1}} + \int_{\mathbb{R}} |P(t)|^2 \Phi(t/T) \underbrace{\int_{|u| \ge T} |u|^{1/2} K(u) \, du \, dt}_{\ll T^{-1/2}} \ll \int_{\mathbb{R}} \left( \frac{|t|^{1/2}}{T} + \frac{1}{T^{1/2}} \right) |P(t)|^2 \Phi(t/T) dt \ll T^{1/2} E_N^*, $$

where we used $K(u) \ll u^{-2}$ and Lemma 2 together with the quick decay of $\Phi$. Let $A = \frac{178}{13}$. Then by Ivić's theorem [Ivi03, Theorem 8.3] we have

$$ \int_0^T |\zeta(1/2 + it)|^A dt \ll T^{2+\frac{3(A-12)}{22}+\varepsilon} = T^{29/13+\varepsilon}. \quad (33) $$

The contribution to $\text{Int}_1$ of small $u$ is

$$ \int_{|u| \le T} K(u) \left( \int_{\mathbb{R}} |\zeta(1/2 + it + iu)||\zeta(1/2 - it + iu)||P(t)|^2 \Phi(t/T) \, dt \right) \, du. $$

To estimate the term inside the brackets, we use Hölder's inequality with parameters $1/A + 1/A + 1/B = 1$, so that $B = \frac{A}{A-2} = \frac{89}{76}$, and write $|P(t)|^2 = |P(t)|^{2-2/B} |P(t)|^{2/B}$. By Lemma 2 together with (32) and (33) we deduce

$$ \begin{aligned} & \int_{\mathbb{R}} |\zeta(1/2 + it + iu)||\zeta(1/2 - it + iu)||P(t)|^2 \Phi(t/T) \, dt \\ & \ll \left( \int_{\mathbb{R}} |\zeta(1/2 + it + iu)|^A \Phi(t/T) \, dt \right)^{1/A} \left( \int_{\mathbb{R}} |\zeta(1/2 - it + iu)|^A \Phi(t/T) \, dt \right)^{1/A} \times \\ & \qquad \times |P(0)|^{2(1-1/B)} \left( \int_{\mathbb{R}} |P(t)|^2 \Phi(t/T) \, dt \right)^{1/B} \\ & \ll (T^{29/13+\varepsilon})^{2/A} N^{4(1-1/B)} T^{1/B} (E_N^*)^{1/B} \\ & = (T^{29/13+\varepsilon})^{2/A} N^{8/A} T^{(A-2)/A} (E_N^*)^{(A-2)/A}. \end{aligned} \quad (34) $$

Integrating over $u$ we deduce that

$$ \begin{aligned} \text{Int}_1 & \ll (T^{29/13+\varepsilon})^{2/A} N^{8/A} T^{(A-2)/A} (E_N^*)^{(A-2)/A} \\ & = T^{105/89+\varepsilon} N^{52/89} (E_N^*)^{76/89}. \end{aligned} $$
Using Lemma 3 and (29) and inserting our bounds for $\mathrm{Int}_1$, $\mathrm{Int}_2$ and $\mathrm{Int}_3$ we obtain

$$
\begin{align*}
& \sum_{1 \le m, n \le M} \sum_{2^{u-1} \le j_1, j_2 < 2^u} \mathbf{1}(|j_1 z_m - j_2 z_n| < 1) \\
& \ll \frac{2^u}{T} \int_{\mathbb{R}} G(t) |P(t)|^2 \Phi(t/T) dt \\
& \ll \frac{2^u}{T} \left( N^{5+\varepsilon/4} + T^{105/89+\varepsilon} N^{52/89} (E_N^*)^{76/89} + T^{1/2} E_N^* \right) \\
& \ll \frac{2^u}{T} \left( N^{5+\varepsilon/4} + T^{105/89+\varepsilon} N^{52/89} (E_N^*)^{76/89} \right),
\end{align*}
$$

where we used that the term $T^{1/2}E_N^*$ is dominated by the other two summands.
Substituting $T = 2^u N^{1+\varepsilon/2}$ and using that $2^u \ll N$, we finally get the upper bound

$$
\begin{align*}
& \sum_{1 \le m, n \le M} \sum_{2^{u-1} \le j_1, j_2 < 2^u} \mathbf{1}(|j_1 z_m - j_2 z_n| < 1) \\
& \ll \frac{1}{N^{1+\varepsilon/2}} (N^{5+\varepsilon/4} + N^{52/89} N^{2 \times 105/89 + 2\varepsilon} (E_N^*)^{76/89}) \\
& = N^{4-\varepsilon/4} + N^{173/89 + 2\varepsilon} (E_N^*)^{76/89}.
\end{align*}
$$

This establishes Lemma 5.

<figure><img src="image_2.png" /></figure>

We note here that conditionally under the Lindelöf hypothesis we could estimate
the integral in line (34) much more efficiently, by using a pointwise bound for the
zeta function and estimating the remaining integral with Lemma 2.

Combining Lemmas 3, 4 and 5, we have shown in our Case 1 analysis that the
contribution of pairs $z_m, z_n$ with $z_m, z_n \ge N^{1.01}$ to the counting problem (17) is

$$
\sum_{1 \le m, n \le M, z_m, z_n \ge N^{1.01}} \sum_{2^{u-1} \le j_1, j_2 < 2^u} \mathbf{1}(|j_1 z_m - j_2 z_n| < 1) \ll N^{4-\varepsilon/4} + N^{173/89+2\varepsilon} (E_N^*)^{76/89}. \tag{35}
$$

Before we move on to Case 2, we make some further comments on our argument
above. Intuitively, it would seem more natural to work with $Q(t) := \sum_m z_m^{it}$ rather
than with the more complicated function $P(t)$. However, from a technical point of
view the key problem in the whole argument is to be able to choose an appropriate
value of $T$ which balances the contribution to our integrals of those values of $t$
for which $|t|$ is “large” (this gets worse when $T$ is larger, since the bound for the
zeta function grows polynomially in $T$) against the contribution coming from those
$t$ for which $|t|$ is “small” (this contribution can only be compensated in the final
estimate when $T$ is sufficiently large). When working directly with $Q$, the size of $T$
would need to depend on the size of the $z_m$ in order to be able to control $\int |Q|^2$.
The “orthogonalization” procedure leading to our definition of $P(t)$ gives us more
freedom in our choice of $T$. The whole problem discussed in this paragraph occurs
only in the real-number setting, in contrast to the integer setting.
• Case 2: Counting solutions for $z_m, z_n$ with $\min\{z_m, z_n\} < N^{1.01}$.

First we consider the contribution to (17) of those $z_m$ and $z_n$ for which $\max\{z_m, z_n\} < 4N^{1/4}$. We assumed that $x_{n+1} - x_n \ge c > 0$, so $z_n \ge c$ for all $n$. Furthermore, we deduce that among $z_1, \dots, z_M$ there are at most $\ll N^{5/4}$ many elements which are smaller than $4N^{1/4}$ (we suppress the dependence of the implied constant on $c$). Note that whenever $j_1$ and $z_m, z_n$ are fixed, there are at most $\ll 1$ many possible choices for $j_2$ such that $|j_1z_m - j_2z_n| < 1$, again since $z_n \ge c$. Thus the total contribution of pairs $z_m, z_n$ with $\max\{z_m, z_n\} < 4N^{1/4}$ to our counting problem is at most

$$ \sum_{\substack{1 \le m, n \le M \\ \max\{z_m, z_n\} < 4N^{1/4}}} \sum_{2^{u-1} \le j_1, j_2 < 2^u} \mathbf{1}(|j_1z_m - j_2z_n| < 1) \ll (N^{5/4})^2 2^u \ll N^{7/2}. $$

Now consider the case when $\max\{z_m, z_n\} \ge 4N^{1/4}$. Recall that we have localized $j_1, j_2$ into a dyadic interval in the counting problem. This implies a similar localization for $z_m$ and $z_n$, since $j_1/j_2 \in [1/2, 2]$ and $|j_1z_m - j_2z_n| < 1$ are only possible if we have $z_m/z_n \in [1/4, 4]$, given the fact that $\max\{z_m, z_n\} \ge 4N^{1/4}$.

Thus we can restrict ourselves in the counting problem (17) to the case when $z_m \in [4N^\beta, 8N^\beta)$ for some $\beta \ge 1/4$, and when consequently $z_n$ needs to be in $[N^\beta, 32N^\beta)$. Note that there are $\ll \log N$ many intervals of this form necessary to cover the whole relevant range $[N^{1/4}, N^{1.01}]$, and clearly we only need to consider $1/4 \le \beta \le 1.01$. We count more solutions if we relax the condition to $z_m, z_n \in [N^\beta, 32N^\beta)$. Thus, let us consider

$$ \sum_{\substack{1 \le m, n \le M, \\ z_m, z_n \in [N^\beta, 32N^\beta]}} \sum_{2^{u-1} \le j_1, j_2 < 2^u} \mathbf{1}(|j_1z_m - j_2z_n| < 1) \quad (36) $$

for some $\beta \in [1/4, 1.01]$. We set up everything as in Case 1, but now we define $T = 2^u N^\beta$. We define the $b_k$'s as before but restricting ourselves to those $z_m$ contained in $[N^\beta, 32N^\beta)$. That is, we set

$$ b_k = \sum_{\substack{1 \le m \le M, \\ z_m \in [N^\beta, 32N^\beta]}} \mathbf{1}(z_m \in [k, k+1)). $$

Note that previously we had $\sum_k b_k = M \le N^2$, whereas now we have a stronger bound. Applying the Cauchy-Schwarz inequality we obtain

$$ \sum_k b_k \ll \sqrt{E_N^*} N^{\beta/2}. \quad (37) $$

We define $(a_h)_{h \ge 0}$ and $P(t)$ as in Case 1; see (19) and (20). In the present case the inequality $|j_1z_m - j_2z_n| < 1$ is only possible when $\left|\frac{z_m}{z_n} - \frac{j_2}{j_1}\right| \le \frac{1}{2^{u-1}N^\beta}$. By construction, $2^{u-1}N^\beta$ becomes large in comparison with $T$, and we can continue to
argue as in Case 1. Note that now, as a consequence of (37), we have $|P(0)|^2 \ll E_N^* N^\beta$ instead of $|P(0)|^2 \ll N^4$ as in Case 1. Proceeding as in Case 1 we obtain

$$ \mathrm{Int}_2, \mathrm{Int}_3 \ll N^{1+\beta+\varepsilon/4} E_N^* + T^{1/2} E_N^*. $$

As for $\mathrm{Int}_1$, we now obtain, again writing $A = \frac{178}{13}$,

$$ \begin{aligned} \mathrm{Int}_1 &\ll (T^{29/13+\varepsilon})^{2/A} (E_N^* N^\beta)^{2/A} (E_N^*)^{1-2/A} T^{1-2/A} \\ &= E_N^*(T^{29/13+\varepsilon})^{2/A} (N^\beta)^{2/A} T^{1-2/A}. \end{aligned} $$

Substituting $T = 2^u N^\beta \ll N^{1+\beta}$, we conclude that (36) is bounded by

$$ \begin{aligned} &\ll \frac{2^u}{T} (E_N^* N^{1+\beta+\varepsilon/4} + E_N^* T^{105/89+\varepsilon} N^{13\beta/89}) \\ &\ll E_N^* N^{1+\varepsilon/4} + E_N^* N^{105/89+29\beta/89+(1+\beta)\varepsilon}. \end{aligned} $$

Recall that we only need to consider $\beta \le 1.01$, and that there are at most $\ll \log N$ many different values of $\beta$ to consider for Case 2. Hence it follows that, for every fixed value of $u$, (36) is bounded by

$$ \ll E_N^* N^{1.51+3\varepsilon}. \qquad (38) $$

We could have used a different value of $A$ here (for instance $A = 12$, namely Heath-Brown's bound on the twelfth moment from [Hea78]) to arrive at (38). However, we kept the same parameters as in Case 1 to simplify the writing.

Finally, inserting in (16) resp. (17) the bound (35) from Case 1 together with the bound (38) yields

$$ \mathrm{Var}(h_N, \mu) \ll \max \left( N^{-\varepsilon/8} + N^{-183/89+3\varepsilon} (E_N^*)^{76/89}, E_N^* N^{-2.49+4\varepsilon} \right), $$

which concludes the Proof of Lemma 1.

# 6 Proof of Theorem 1: conclusion of the proof

The crucial ingredient in the Proof of Theorem 1 is the variance bound from Lemma 1. We record that we have, for every sufficiently small $\varepsilon > 0$,

$$ \mathrm{Var}(h_N, \mu) \ll \max \left( N^{-\varepsilon/8} + N^{-183/89+3\varepsilon} (E_N^*)^{76/89}, E_N^* N^{-2.49+4\varepsilon} \right). $$

Inserting $E_N^* \ll N^{183/76-\delta}$ shows that for every sufficiently small $\delta' > 0$ we have

$$ \mathrm{Var}(h_N, \mu) \ll N^{-\delta'}. \qquad (39) $$

Everything else now follows from a standard procedure. To be a bit more specific, convergence in (10) can be established using the estimate for the expectations in Section 3, and using the variance bound (39) together with Chebyshev's inequality and the Borel-Cantelli lemma. From that we get a convergence result for almost all
$\alpha \in \mathbb{R}$, for fixed values of $r$ and $s$. One notes that there are only countably many possible values of $r$, and that by continuity/monotonicity it is sufficient to consider countably many values of $s$. Since a countable union of sets of measure zero has measure zero as well, almost all $\alpha \in \mathbb{R}$ have the property that (10) holds for all $r$ and all $s$, as desired. We refer the reader to [ALL17] or [RT20], where this argument is carried out in full detail. It applies without any modifications to the situation in the present paper.

## 7 Proof of Theorem 2

We set up the same machinery as in the Proof of Theorem 1. Controlling the expectations, as in Section 3 above, is unproblematic. The crucial part is again the variance estimate. As in Section 4, we are led to the counting problem

$$ \sum_{1 \le m, n \le M} \sum_{2^{u-1} \le j_1, j_2 < 2^u} 1(|j_1 z_m - j_2 z_n| < 1) $$

where $\{z_1, \dots, z_M\}$ is the multi-set of all the absolute differences $\{|x_m - x_n| : 1 \le m, n \le N, m \neq n\}$. As above, $u$ is a positive integer with $2^u \le 2rN$, and $M = N^2 - N$.

As in the general argument before, we can easily dispose of the contribution of those $z_m, z_n$ for which $\max\{z_m, z_n\} < 4N^{1/4}$. Thus again we can localize $z_m$ and $z_n$, and restrict ourselves to counting

$$ \sum_{\substack{1 \le m, n \le M, \\ z_m, z_n \in [N^\beta, 32N^\beta)}} \sum_{2^{u-1} \le j_1, j_2 < 2^u} 1(|j_1 z_m - j_2 z_n| < 1) \quad (40) $$

for some $1/4 \le \beta \le 1.01$, where we understand that the final range for $\beta = 1.01$ extends over all of $[N^\beta, \infty)$ rather than only $[N^\beta, 32N^\beta)$. Note that in total at most $\ll \log N$ many different values of $\beta$ need to be considered. Let $\beta \in [1/4, 1.01]$ be fixed. Let $u$ in (40) be fixed. For integers $k \ge 0$ we define

$$ b_k = \sum_{\substack{1 \le m \le M, \\ z_m \in [N^\beta, 32N^\beta)}} 1\left(z_m \in \left[\frac{k}{2^u}, \frac{k+1}{2^u}\right)\right), \quad (41) $$

where again in case $\beta = 1.01$ the range $[N^\beta, 32N^\beta)$ is understood to be replaced by $[N^\beta, \infty)$. Note the difference in comparison with (18). There we collected all $z_m$ in a range of the form $[k, k+1)$, since we could only control the number of solutions of the specific inequality (4), which has “< 1” on the right-hand side. In contrast we can now control the number of solutions on a finer scale, and can accordingly set shorter ranges for the grouping of the $z_m$ (where $\gamma = 2^{-u}$).

Let $\varepsilon > 0$ be a small constant (chosen depending on the size of $\eta$ and $\delta$ in the statement of the theorem). Set $T = 2^u N^{\min\{\beta-\varepsilon, 1+\varepsilon\}}$. Unlike the argument in the
general case in Section 4, we do not explicitly distinguish between Case 1 and Case 2, but have implicitly included this distinction into the way that $T$ is defined. As in Section 4, we split the interval $[1, \infty)$ into a disjoint union $\bigcup_{h=0}^{\infty} I_h$, where

$$I_h = \left[ \left\lfloor \left(1 + \frac{1}{T}\right)^h \right\rfloor, \left\lfloor \left(1 + \frac{1}{T}\right)^{h+1} \right\rfloor \right),$$

and set

$$a_h = \left( \sum_{k: k/2^u \in I_h} b_k^2 \right)^{1/2}, \quad h \ge 0,$$

as well as

$$P(t) = \sum_{h=0}^{\infty} a_h \left(1 + \frac{1}{T}\right)^{iht}.$$

Then by construction we again have

$$\int_{\mathbb{R}} |P(t)|^2 \Phi(t/T) dt \ll T \sum_{h=0}^{\infty} a_h^2$$

as during the Proof of Lemma 2 in Section 4, but now we can continue to estimate this by (5) and obtain

$$T \sum_{h \ge 0} a_h^2 \ll T \sum_{k \ge 0} b_k^2 \ll T E_{N, 2^{-u}}^* \ll T (N^{2+\eta} + 2^{-u}N^{3-\delta}). \quad (42)$$

Now we establish the necessary upper bound on $|P(0)|$. Trivially we always have $|P(0)| \ll N^2$. For small values of $\beta$ we obtain a better estimate. Note that

$$
\begin{align*}
|P(0)| &= \sum_h a_h \le \sum_k b_k = \# \{m : z_m \in [N^\beta, 32N^\beta)\} \\
&\le \sum_{a=[N^\beta]}^{[32N^\beta]} \# \{m : z_m \in [a, a+1)\} \\
&\ll N^{\beta/2} \left( \sum_{a=[N^\beta]}^{[32N^\beta]} \left( \# \{m : z_m \in [a, a+1)\} \right)^2 \right)^{1/2} \tag{43} \\
&\ll N^{\beta/2} \sqrt{E_{N,1}^*} \\
&\ll N^{\beta/2} N^{(3-\delta)/2} = N^{3/2+\beta/2-\delta/2}, \tag{44}
\end{align*}
$$

where for (43) we used Cauchy-Schwarz and for (44) our assumption (5) in the statement of Theorem 2. So overall we have $|P(0)| \ll N^{\min\{2,3/2+\beta/2-\delta/2\}}$.
As in the Proof of Lemma 3 in Section 4, we assume without loss of generality that $j_1 \ge j_2$. Recall that $j_1, j_2 \ge 2^{u-1}$ by assumption. Assume that $z_m \in [k/2^u, (k+1)/2^u)$ for some $k$. Then the inequality $|j_1z_m - j_2z_n| \le 1$ is only possible when

$$ \left| \frac{\left[ \frac{j_1 k}{j_2} \right]}{2^u} - z_n \right| \le \frac{5}{2^u}, $$

which is a version of (25) that is adapted to the construction in (41). Arguing as in the lines leading to (26), this again gives

$$ \sum_{z_m \in I_{h_1}, z_n \in I_{h_2}} \mathbf{1}(|j_1 z_m - j_2 z_n| < 1) \ll a_{h_1} a_{h_2}, $$

which perfectly resembles (26) but where now the $b_k$ and $a_h$ are defined in a different way according to (41). Note that $T$ is chosen in such a way that $j_1z_m$ and $j_2z_n$ exceed $T$; indeed, by assumption $j_1, j_2 \ge 2^{u-1}$ and $z_m, z_n \ge N^\beta$, while $T \le 2^u N^{\beta-\epsilon}$ by definition. Thus we can continue the argument as in Section 5. It turns out that in this setting in order to bound $\mathrm{Int}_1$ it is sufficient to use $|\zeta(1/2 + it)| \ll |t|^{1/6}$ (which is essentially the Weyl–Hardy–Littlewood bound), rather than the more elaborate argument relying on estimates for moments of the Riemann zeta function. We obtain

$$ \begin{aligned} & \sum_{\substack{1 \le m, n \le M, \\ z_m, z_n \in [N^\beta, 32N^\beta)}} \sum_{2^{u-1} \le j_1, j_2 < 2^u} \mathbf{1}(|j_1 z_m - j_2 z_n| < 1) \\ & \ll \frac{2^u}{T} (|P(0)|^2 N^{1+\varepsilon/4} + T^{1+1/3} E_{N, 2^{-u}}^*) \\ & \ll \frac{2^u N^{\min\{4,3+\beta-\delta\}} N^{1+\varepsilon/4}}{T} + (2^u)^{4/3} N^{(1+\varepsilon)/3} N^{2+\eta} + (2^u)^{4/3} N^{(1+\varepsilon)/3} 2^{-u} N^{3-\delta} \\ & \ll N^{1+\varepsilon/4+\min\{4,3+\beta-\delta\}-\min\{\beta-\varepsilon, 1+\varepsilon\}} + N^{11/3+\varepsilon+\eta} + N^{11/3+\varepsilon-\delta} \\ & \ll N^{4-\varepsilon/2} \end{aligned} $$

if $\varepsilon$ was chosen sufficiently small (with respect to $\delta$). Here we used (5) as well as $T \ll 2^u N^{1+\varepsilon}$, $2^u \ll N$. Noting that we need to consider at most $\ll \log N$ different values of $\beta$, this gives the necessary variance estimate. The remaining part of the Proof of Theorem 3 can be carried out exactly as in the Proof of Theorem 1. We remark that any subconvex bound for the Riemann zeta function would be sufficient to derive the same conclusion.

## 8 Proof of Theorem 3

We assume that $\theta > 1$ is fixed, and consider the sequence $(x_n)$ defined by $x_n = n^\theta$, $n \ge 1$. Note that with this definition we have $x_{n+1} - x_n \ge 1$ for all $n \ge 1$, so the assumption $x_{n+1} - x_n \ge c$ of Theorem 2 is satisfied in this case with $c = 1$. The
following lemma of Robert and Sargos shows that the necessary bound on $E_{N,\gamma}^*$ also is satisfied for this sequence.<fnref n="3" />

LEMMA 7 ([RS06, Theorem 2]). Let $\theta \neq 0, 1$ be a fixed real number. For any $\gamma > 0$ and $B \ge 2$, let $\mathcal{N}(B, \gamma)$ denote the number of 4-tuples $(n_1, n_2, n_3, n_4) \in \{B+1, B+2, \dots, 2B\}^4$ for which

$$|n_1^\theta - n_2^\theta + n_3^\theta - n_4^\theta| \le \gamma. \qquad (45)$$

Then for every $\varepsilon > 0$,

$$\mathcal{N}(B, \gamma) \ll_{\varepsilon} B^{2+\varepsilon} + \gamma B^{4-\theta+\varepsilon}.$$

The restriction to a dyadic range for $(n_1, n_2, n_3, n_4)$ in the statement of the lemma does not actually play a role. This is easily seen by interpreting the number of solutions of the inequality as an $L^4$-norm. Indeed, generalizing the definition in (11) and setting

$$d\mu_{2\gamma}(x) = \frac{(\sin(\gamma x))^2}{\pi\gamma x^2} dx,$$

we have a measure whose Fourier transform is a (normalized) tent function on $[-2\gamma, 2\gamma]$. Let $E_{N,\gamma}^*$ denote the number of solutions of (45), subject to $(n_1, n_2, n_3, n_4) \in \{1, \dots, N\}^4$. Assume for simplicity of writing that $N$ is a power of 2, i.e. $N = 2^L$ for some $L \ge 1$. Then applying Hölder's inequality we obtain

$$
\begin{align*}
E_{N,\gamma}^* & \ll \int_{\mathbb{R}} \left( \sum_{\ell=0}^{L} \sum_{2^{\ell-1} < n \le 2^\ell} e^{2\pi i n^\theta x} \right)^4 d\mu_{2\gamma}(x) \\
& \ll \int_{\mathbb{R}} \left( \sum_{\ell=0}^{L} 1 \right)^3 \sum_{\ell=0}^{L} \left| \sum_{2^{\ell-1} < n \le 2^\ell} e^{2\pi i n^\theta x} \right|^4 d\mu_{2\gamma}(x) \\
& \ll (\log N)^3 \sum_{\ell=0}^{L} \int_{\mathbb{R}} \left| \sum_{2^{\ell-1} < n \le 2^\ell} e^{2\pi i n^\theta x} \right|^4 d\mu_{2\gamma}(x) \\
& \ll (\log N)^3 \sum_{\ell=0}^{L} E_{2^{\ell-1}, 2\gamma} \\
& \ll_{\varepsilon} N^{2+\varepsilon} + \gamma N^{4-\theta+\varepsilon}, \tag{46}
\end{align*}
$$

<footnote n="3">³ We thank Niclas Technau for pointing out to us that the estimate in Lemma 7 is also contained as a special case in a general result in a very recent paper of Huang [Hua20]. Huang’s result gives improved error terms, but for our application this does not play a role. However, the generality of Huang’s results could allow further applications of our method in the spirit of our Theorem 3.</footnote>
which is obtained by interpreting the integrals in line (46) in terms of solutions of the Diophantine inequality (45), and applying Lemma 7 with parameters $2\gamma$ and $B = 2^{\ell-1}$. Thus we have

$$E_{N,\gamma}^* \ll_{\epsilon} N^{2+\epsilon} + \gamma N^{4-\theta+\epsilon} \quad (47)$$

for any $\epsilon > 0$. Consequently all assumptions of Theorem 2 are satisfied, and we can conclude that $(n^\theta \alpha)_n$ has Poissonian pair correlation for almost all $\alpha$.

## 9 Closing remarks

As remarked in the introduction, our method breaks down completely when the growth order of the sequence $(x_n)_{n \ge 1}$ is only linear or even slower. Not only does the “lattice point counting with the zeta function” argument from Section 5 fail to work in this situation, but there is a much more fundamental reason why the whole approach based on calculating first and second moments (expectations and variances, as in Sections 3 and 4) fails to work in this setup. To give a brief sketch of what causes the problem, assume that $(x_n)_n$ is a sequence of reals such that $x_n \le n$, $n \ge 1$. Assume that we want to bound the variance in analogy with (15), so say we want to show that

$$\int_{\mathbb{R}} \left( \frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} \mathbf{1}_{[-1/N, 1/N]}(x_m \alpha - x_n \alpha) \right)^2 d\mu(\alpha) \quad (48)$$

tends to zero as $N \to \infty$ (where we write the original indicator function instead of its approximation by a trigonometric polynomial, and where for simplicity of writing we set $s = 1$). By our assumption on the growth of $(x_n)_n$, all differences $x_m - x_n$ appearing in the sum above are uniformly bounded by $N$. Thus we have $\mathbf{1}_{[-1/N, 1/N]}(x_m \alpha - x_n \alpha) = 1$ throughout the range $\alpha \in [-1/N^2, 1/N^2]$, for all $m, n \le N$. Consequently

$$\begin{align*}
& \int_{\mathbb{R}} \left( \frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} \mathbf{1}_{[-1/N, 1/N]}(x_m \alpha - x_n \alpha) \right)^2 d\mu(\alpha) \\
\ge & \int_{-1/N^2}^{1/N^2} \left( \frac{1}{N} \sum_{\substack{1 \le m, n \le N, \\ m \ne n}} 1 \right)^2 d\mu(\alpha) \\
\gg & 1.
\end{align*}$$
Thus the variance fails to tend to zero for a slowly growing $(x_n)_n$, due to the fact that the contribution of small values of $\alpha$ to the variance integral is too large.<fnref n="4" /> The argument used in Section 5 fails to work in a similar way for slowly growing $(x_n)_n$, since the error terms coming from the contribution to the integrals of values of $t$ near zero become too large.<fnref n="5" /> Consequently, it seems that for establishing Poissonian pair correlation of $(x_n\alpha)_n$ for almost all $\alpha$ for slowly growing $(x_n)_n$ some genuine new ideas are necessary. Note that we cannot simply remove all values of $\alpha$ near zero from the variance integral (48) by replacing $\mu$ with some other measure which vanishes for small $\alpha$, since such a measure would fail to have non-negative Fourier transform (thereby causing major problems in other places). Note also that all these problems with slowly growing sequences $(x_n)_n$ are a novel aspect which only shows up in the real-number setup—in contrast, when $(a_n)_n$ is an integer sequence which grows at most linearly, then $(a_n\alpha)_n$ is known to fail to have Poissonian pair correlation for any $\alpha$, because the additive energy of $(a_1, \dots, a_N)$ necessarily is of maximal possible order (cf. [LS20]).

We emphasize that the fact that our method fails to work in the case of slowly growing sequences $(x_n)_n$ should not be understood as indicating that in such a case $(x_n\alpha)_n$ should necessarily fail to have Poissonian pair correlation for almost all $\alpha$. Quite on the contrary, there are good reasons to expect that also for slowly growing $(x_n\alpha)_n$ one should in “generic” situations obtain Poissonian pair correlation for almost all $\alpha$. It seems that the property of having Poissonian pair correlation for $(x_n\alpha)_n$ for almost all $\alpha$ can only be prevented by a certain (“small-scale”) combinatorial obstruction, in such a way that the case of integer sequences $(x_n)_n$ with slowly growing $(x_n)_n$ can be seen as a degenerate situation exhibiting exactly this type of combinatorial obstruction (coming from the fact that in the integer setup everything which is smaller than one in absolute value necessarily equals zero). We believe that these are very interesting phenomena, and we propose the following open problems.

**Open Problem 1.** Let $\theta \in (0, 1)$. Show that $(n^\theta \alpha)_{n \ge 1}$ has Poissonian pair correlation for almost all $\alpha$. Note that Lemma 7 is still valid for this range of $\theta$.

<footnote n="4">⁴ A similar argument appears at the end of [RS98], where it is used to show that the $L^2$ approach fails to work in the case of the triple correlation of $(n^2\alpha)_n$; cf. also [TW20].</footnote>

<footnote n="5">⁵ It might be difficult to spot at a quick glance, so we briefly comment on where the speed of growth of $(x_n)_n$ was used in our argument in Sections 5 and 7. There is a term $|P(0)|^2 N^{1+\epsilon/4}$ coming from the contribution of values of $t$ near the origin to the integral. This term is divided by $T$ at the end of the calculation, so we cannot take $T$ too small since we need $N|P(0)|^2 N^{1+\epsilon/4}/N^4 T \to 0$. On the other hand, we cannot take $T$ too large, since we need $T \ll 2^u z_n$ to be able to detect the solutions of our Diophantine inequality. To balance everything out, we need to be able to assure that there are not too many small values of $z_n$ (i.e., not too many differences $x_m - x_n$ which are “small”). In our Proof of Theorem 1 our assumption on the order of the additive energy takes care of this: it is easy to see that an upper bound on $E_N^*$ implies an upper bound on the number of “small” differences $x_m - x_n$, which is what we used in Case 2 of Section 5. In the setting of Theorem 2 a similar argument based on the energy assumption allowed us to control the number of small differences $x_m - x_n$; the relevant equations there are (37) and (44).</footnote>
**Open Problem 2.** Let $x_n = n + \log n$. Show that $(x_n\alpha)_n$ has Poissonian pair correlation for almost all $\alpha$. We note that it is possible to establish a variant of Lemma 7 for this setting (with exponent 3 in place of $4 - \theta$).

**Open Problem 3.** Let $x_n = n \log n$, $n \ge 1$. Show that $(x_n\alpha)_n$ has Poissonian pair correlation for almost all $\alpha$.

Clearly the exponent $183/76 - \delta$ in the statement of Theorem 1 is not optimal, and most likely it can be improved to $3 - \delta$ (which is the case conditionally under the Lindelöf hypothesis). It seems to us that the method of Bloom and Walker [BW20], which led to a quantitative improvement of the results of [ALL17], cannot be used here. Their method relied on sum-product estimates, which, roughly speaking, leads to an integrand $|P(t)|^2$ being replaced by $|P(t)|^4$. In the case of integer sequences (when working with the random model of the zeta function) one has perfect orthogonality, so that $\int |P|^4$ can be efficiently bounded. In our setting the situation is quite different – we have constructed our function $P(t)$ in such a way that the diagonal contribution dominates when calculating $\int |P|^2$, but we do not have orthogonality for $\int |P|^4$ and cannot efficiently bound this integral.

**Open Problem 4.** Show that Theorem 1 remains valid under the weaker assumption $E_N^* \ll N^{3-\delta}$ for some $\delta > 0$. Show that this can be further relaxed to assuming $E_{N,\gamma}^* \ll \gamma N^{4-\delta}$, for all $\gamma$ in a range from roughly $1/N$ to $1$. It might even be the case that only values of $\gamma$ near a critical size of roughly $1/N$ are relevant. Note that if the condition $E_{N,\gamma}^* \ll \gamma N^{4-\delta}$ uniformly for $\gamma \in [1/N, 1]$ truly is the “right” condition, then this would give a unified picture for the real-sequence case as well as for the integer-sequence case. Indeed, in the latter case clearly $E_{N,\gamma}^* = E_N$ for all $\gamma < 1$ and thus the condition would reduce to $E_N \ll N^{3-\delta}$, in accordance with the criterion stated after (2).

As noted, in the case of an integer sequence $(x_n)_n$ it is known that $(x_n\alpha)_n$ cannot have Poissonian pair correlation for almost all $\alpha$ when $E_N \gg N^3$. It would be interesting to obtain an analogous result in the case of real sequences.

**Open Problem 5.** Show that unlike in the integer case, it is possible for an increasing sequence $(x_n)_{n \ge 1}$ of reals that $E_N^* \gg N^3$ and that $(x_n\alpha)_n$ has Poissonian pair correlation for almost all $\alpha$ (compare Open Problems 1 and 2 above, where $E_N^* \gg N^3$). Establish a criterion (stated for example in terms of $E_{N,\gamma}^*$) which ensures that $(x_n\alpha)_n$ does not have Poissonian pair correlation for almost all $\alpha$. A candidate for such a criterion is that $E_{N,\gamma}^* \gg \gamma N^4$ for some $\gamma = \gamma(N)$ for infinitely many $N$, where maybe one also has to assume that these values of $\gamma$ are of size $\gamma \approx 1/N$.
Acknowledgments

CA is supported by the Austrian Science Fund (FWF), Projects F-5512, I-3466, I-4945 and Y-901. DE is supported by FWF Projects F-5512 and Y-901. MM is supported by FWF Project P-33043. We thank Winston Heap, Olivier Robert, Zeev Rudnick, Ilya Shkredov, Igor Shparlinski, Athanasios Sourmelidis and Niclas Technau for discussions and comments.

**Funding** Open access funding provided by Graz University of Technology.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

[Ais16] C. AISTLEITNER. Lower bounds for the maximum of the Riemann zeta function along vertical lines. *Math. Ann.*, (1-2)365 (2016),473–496

[ABS15] C. AISTLEITNER, I. BERKES, and K. SEIP. GCD sums from Poisson integrals and systems of dilated functions. *J. Eur. Math. Soc. (JEMS)*, (6)17 (2015), 1517–1546

[ALP18] C. AISTLEITNER, T. LACHMANN, and F. PAUSINGER. Pair correlations and equidistribution. *J. Number Theory*, 182(2018), 206–220.

[ALT19] C. AISTLEITNER, T. LACHMANN, and N. TECHNAU. There is no Khintchine threshold for metric pair correlations. *Mathematika*, (4)65 (2019),929–949

[ALL17] C. AISTLEITNER, G. LARCHER, and M. LEWKO. Additive energy and the Hausdorff dimension of the exceptional set in metric pair correlation problems. *Israel J. Math.*, (1)222 (2017),463–485

[BPT01] I. BERKES, W. PHILIPP, and R. TICHY. Pair correlations and *U*-statistics for independent and weakly dependent random variables. *Illinois J. Math.*, (2)45 (2001), 559–580

[BCGW18] T. F. BLOOM, S. CHOW, A. GAFNI, and A. WALKER. Additive energy and the metric Poissonian property. *Mathematika*, (3)64 (2018),679–700

[BW20] T. F. BLOOM and A. WALKER. GCD sums and sum-product estimates. *Israel J. Math.*, (1)235 (2020),1–11
[BS17] A. BONDARENKO and K. SEIP. Large greatest common divisor sums and extreme values of the Riemann zeta function. *Duke Math. J.*, (9)166 (2017),1685–1701

[BS18] A. BONDARENKO and K. SEIP. Extreme values of the Riemann zeta function and its argument. *Math. Ann.*, 372(3-4):999–1015, 2018.

[BKW10] J. BRÜDERN, K. KAWADA, and T. D. WOOLEY. Additive representation in thin sequences, VIII: Diophantine inequalities in review. In *Number theory*, volume 6 of *Ser. Number Theory Appl.*, World Sci. Publ., Hackensack, NJ, (2010) pages 20–79.

[CLZ15] S. CHAUBEY, M. LANIUS, and A. ZAHARESCU. Pair correlation of fractional parts derived from rational valued sequences. *J. Number Theory*, 151 (2015), 147–158

[BT19] R. de la BRETÈCHE and G. TENENBAUM. Sommes de Gál et applications. *Proc. Lond. Math. Soc. (3)*, (1)119 (2019), 104–134

[BMT21] R. de la BRETÈCHE, M. MUNSCH, and G. TENENBAUM. Small Gál sums and applications. *Journal of the London Mathematical Society*, (1)103 (2021), 336–352

[DT97] M. DRMOTA and R. F. TICHY. *Sequences, discrepancies and applications*, volume 1651 of *Lecture Notes in Mathematics*. Springer-Verlag, Berlin, 1997.

[EMV15] D. EL-BAZ, J. MARKLOF, and I. VINOGRADOV. The two-point correlation function of the fractional parts of $\sqrt{n}$ is Poisson. *Proc. Amer. Math. Soc.*, (7)143 (2015) 2815–2828,

[GL17] S. GREPSTAD and G. LARCHER. On pair correlation and discrepancy. *Arch. Math. (Basel)*, (2)109 (2017), 143–149

[Gry13] D. J. GRYNKIEWICZ. *Structural additive theory*, volume 30 of *Developments in Mathematics*. Springer, Cham, (2013).

[Hea78] D. R. HEATH-BROWN. The twelfth power moment of the Riemann-function. *Quart. J. Math. Oxford Ser. (2)*, (116)29 (1978), 443–462

[Hea10] D. R. HEATH-BROWN. Pair correlation for fractional parts of $\alpha n^2$. *Math. Proc. Cambridge Philos. Soc.*, (3)148 (2010), 385–407

[Hil09] T. HILBERDINK. An arithmetical mapping and applications to $\Omega$-results for the Riemann zeta function. *Acta Arith.*, (139)139 (2009), 341–367

[Hua20] J.-J. HUANG. The density of rational points near hypersurfaces. *Duke Math. J.*, (11)169 (2020), 2045–2077

[Ivi03] A. IVIĆ. *The Riemann zeta-function*. Dover Publications, Inc., Mineola, NY, 2003. Theory and applications, Reprint of the 1985 original [Wiley, New York; MR0792089 (87d:11062)].

[KN74] L. KUIPERS and H. NIEDERREITER. *Uniform distribution of sequences*. Wiley-Interscience [John Wiley & Sons], New York-London-Sydney, (1974).

[LS20] G. LARCHER and W. STOCKINGER. Pair correlation of sequences ($\{a_n\}_{n \in \mathbb{N}}$ with maximal additive energy. *Math. Proc. Cambridge Philos. Soc.*, (2)168 (2020), 287–293

[LS20] G. LARCHER and W. STOCKINGER. Pair correlation of sequences with maximal additive energy. *Math. Proc. Cambridge Philos. Soc.*, (2)168 (2020), 287–293

[LR17] M. LEWKO and M. RADZIWILŁ. Refinements of Gál’s theorem and applications. *Adv. Math.*, 305 (2017), 280–297

[Mar01] J. MARKLOF. The Berry-Tabor conjecture. In *European Congress of Mathematics, Vol. II (Barcelona, 2000)*, volume 202 of *Progr. Math.*, Birkhäuser, Basel
(2001) pages 421–427.

[Mar20] J. MARKLOF. Pair correlation and equidistribution on manifolds. *Monatsh. Math.*, (2)191 (2020), 279–294

[MY18] J. MARKLOF and N. YESHA. Pair correlation for quadratic polynomials mod 1. *Compos. Math.*, (5)154 (2018), 960–983

[Mon94] H. L. MONTGOMERY. Ten lectures on the interface between analytic number theory and harmonic analysis, volume 84 of CBMS Regional Conference Series in Mathematics. Published for the Conference Board of the Mathematical Sciences, Washington, DC; by the American Mathematical Society, Providence, RI, (1994).

[RS06] O. ROBERT and P. SARGOS. Three-dimensional exponential sums with monomials. *J. Reine Angew. Math.*, 591 (2006), 1–20

[RS98] Z. RUDNICK and P. SARNAK. The pair correlation function of fractional parts of polynomials. *Comm. Math. Phys.*, (1)194 (1998), 61–70

[RSZ01] Z. RUDNICK, P. SARNAK, and A. ZAHARESCU. The distribution of spacings between the fractional parts of $n^2\alpha$. *Invent. Math.*, (1)145 (2001), 37–57

[RT20] Z. RUDNICK and N. TECHNAU. The metric theory of the pair correlation function of real-valued lacunary sequences. *Illinois J. Math.*, (2020). Advance publication.

[RZ99] Z. RUDNICK and A. ZAHARESCU. A metric result on the pair correlation of fractional parts of sequences. *Acta Arith.*, (3)89 (1999), 283–293

[RZ02] Z. RUDNICK and A. ZAHARESCU. The distribution of spacings between fractional parts of lacunary sequences. *Forum Math.*, (5)14 (2002), 691–712

[Shk20] I. SHKREDOV. On some applications of GCD sums to arithmetic combinatorics, 2020. Preprint, arXiv:2010.03020.

[Shk13] I. D. SHKREDOV. Some new results on higher energies. *Trans. Moscow Math. Soc.*, (2013), pages 31–63.

[TV06] T. TAO and V. VU. *Additive combinatorics*, volume 105 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, (2006).

[TW20] N. TECHNAU and A. WALKER. On the triple correlations of fractional parts of $n^2\alpha$, 2020. Preprint, arXiv:2005.01490.

[TY20] N. TECHNAU and N. YESHA. On the correlations of $n^\alpha$ mod 1, 2020. Preprint, arXiv:2006.16629.

[Tru10] J. L. TRUELSEN. Divisor Problems and the Pair Correlation for the Fractional Parts of $n^2\alpha$. *International Mathematics Research Notices*, (16)2010 (2010) 3144–3183, 01.

[Wey16] H. WEYL. Über die Gleichverteilung von Zahlen mod. Eins. *Math. Ann.*, (3)77 (1916), 313–352

C. AISTLEITNER · D. EL-BAZ · M. MUNSCH

Graz University of Technology, Steyrergasse 30, 8010 Graz, Austria.

daniel.elbaz.88@gmail.com

C. AISTLEITNER

aistleitner@math.tugraz.at
M. Munsch
munsch@math.tugraz.at

Received: October 13, 2020

Revised: February 12, 2021

Accepted: February 22, 2021