Averages of arithmetic functions over polynomials in many variables

Kevin Destagnol Laboratoire de mathématiques d’Orsay
Université Paris Saclay
Orsay
France [email protected] and Efthymios Sofos Mathematics Department
Glasgow University
Glasgow
G12 8QQ
UK [email protected]

Abstract.

We estimate the average of any arithmetic function $k$ over the values of any smooth polynomial in many variables provided only that $k$ has a distribution in arithmetic progressions of fixed modulus. We give several applications of this result including the analytic Hasse principle for an intersection of two cubics in at least 21 variables and asymptotics for the number of integer solutions of a non-algebraic variety.

2020 Mathematics Subject Classification:

11N37, 11P55.

1. Introduction

Numerous challenges within the realm of Diophantine geometry can be rephrased in terms of computing averages of arithmetic functions, denoted as $k:\mathbb{N}\to\mathbb{C}$ , over the values attained by integer polynomials. Examples include, among others, the Hasse principle for conic bundles [9], Manin’s conjecture for Châtelet surfaces [3], [4], [6] and [19] or prime and squarefree values of polynomials (see [24], [23] or [20] for example).

To estimate the average of $k$ even over linear univariate polynomials one must plainly assume that $k$ has an average over every fixed arithmetic progression $a\left(\textnormal{mod}\ q\right)$ . We may thus assume that

\sum_{\begin{subarray}{c}m\in\mathbb{N}\cap[1,x]\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}k(m)\approx\rho(a,q)% \int_{1}^{x}\omega(t)\mathrm{d}t

for some $\rho(a,q)\in\mathbb{C}$ and a function $\omega$ of class $\mathscr{C}^{1}$ . We may then use the Cramér–Granville model to predict the average of $k$ over the values of any polynomial $f$ as done, for example in [20, Appendix] when $k$ is the indicator function of the primes. The statement analogous to [20, Eq. (A.1)] is that the conditional expectation of $k(m)$ given that $m\in[1,x]$ lies in the progression $a\left(\textnormal{mod}\ q\right)$ is $\approx q\rho(a,q)\int_{1}^{x}\omega(t)\mathrm{d}t$ . Letting $W$ be an integer that is divisible by all primes $p$ below some $z\to\infty$ the arguments in [20, pg. 32] lead to the heuristic

	$\displaystyle\sum_{\begin{subarray}{c}\mathbf{t}\in(\mathbb{Z}\cap[-P,P])^{n}% \\ f(\mathbf{t})>0\end{subarray}}k(f(\mathbf{t}))$	$\displaystyle\approx\sum_{\mathbf{a}\in(\mathbb{Z}/W\mathbb{Z})^{n}}\sum_{% \begin{subarray}{c}\\ \mathbf{t}\in(\mathbb{Z}\cap[-P,P])^{n}\\ \mathbf{t}\equiv\mathbf{a}\left(\textnormal{mod}\ W\right)\end{subarray}}W\rho% (f(\mathbf{t}),W)\omega(f(\mathbf{t}))$
		$\displaystyle\approx\int_{[-P,P]^{n}}\omega(f(\mathbf{y}))\mathrm{d}\mathbf{y}% \ \ \ \sum_{\mathbf{a}\in(\mathbb{Z}/W\mathbb{Z})^{n}}\frac{\rho(f(\mathbf{a})% ,W)}{W^{n-1}}.$

We will prove this heuristic for all polynomials $f$ in sufficiently many variables. The cases with a small number of variables are harder; in this case if one is allowed to assume that the function is non-negative and has a weak multiplicative property then work of Nair–Tenenbaum [31] or [5] provide sharp upper bounds.

1.1. The general result

Let $\omega:[1,\infty)\to\mathbb{C}$ be $\mathscr{C}^{1}$ and for any $q\in\mathbb{N}$ and $a\in\mathbb{Z}/q\mathbb{Z}$ let $\rho(a,q)$ be a number in $\mathbb{C}$ . Assume we are given an arithmetic function $k:\mathbb{N}\to\mathbb{C}$ and define

{E_{k,\omega,\rho}(x,q):=\sup_{y\in\mathbb{R}\cap[1,x]}\max_{a\in\mathbb{Z}/q% \mathbb{Z}}\left|\sum_{\begin{subarray}{c}1\leqslant m\leqslant y\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}k(m)-\rho(a,q)\int_{1}% ^{y}\omega(t)\mathrm{d}t\right|.}

(1.1)

We shall henceforth denote $E_{k,\omega,\rho}(x,q)$ by $E(x,q)$ to simplify notation. The function $k$ is equidistributed in progressions modulo $q$ exactly when $E(x,q)$ has smaller order of magnitude than $\int_{1}^{x}\omega$ .

Example 1.1.

Fix any $A>0$ . When $k$ is the indicator function of the primes we choose $\omega(m)=1/\log(m+1)$ and $\rho(a,q)=\mathds{1}(\gcd(a,q)=1)/\varphi(q)$ . The Siegel–Walfisz theorem [28, Eq.(5.77)] is equivalent to

E(x,q)=O\left(\frac{\int_{1}^{x}\omega(t)\mathrm{d}t}{(\log x)^{A}}\right)

with an implied constant independent of $q$ .

Definition 1.1.

Let $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ be an integer polynomial. Let $\mathscr{B}\subset[-1,1]^{n}$ be of the form $\mathscr{B}=\prod_{j=1}^{n}[a_{j},a^{\prime}_{j}]$ for some $a_{j}$ with $\max_{j}|a^{\prime}_{j}-a_{j}|\leqslant 1$ and define

{b=2\max\{|f(\mathbf{t})|:\mathbf{t}\in\mathscr{B}\}.}

(1.2)

Definition 1.2.

For any $z\geqslant 2$ , assume we have a function $m_{p}(z):\{\textrm{primes}\}\to\mathbb{Z}_{\geqslant 1}$ and let

W_{z}:=\prod_{p\leqslant z}p^{m_{p}(z)}.

We shall assume that $m_{p}(z)$ is suitably large so that $\widetilde{\varepsilon}(z)\to 0$ as $z\to\infty$ , where

\widetilde{\varepsilon}(z):=\sum_{p\leqslant z}\frac{1}{p^{1+m_{p}(z)}}.

The number $\widetilde{\varepsilon}(z)$ is an upper bound for the probability that a random integer is divisible by a high power of at least one small prime. We abbreviate $m_{p}(z)$ by $m_{p}$ .

Our main technical tool expresses the error term as a function of $E(x,W_{z})$ .

Theorem 1.1.

Assume that the positive integers $d,n\in\mathbb{N}$ satisfy $n>2^{d}(d-1)$ and let $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ be a smooth form of degree $d$ . Let $s\in\{-1,1\}$ and let $\mathscr{B}\subset\mathbb{R}^{n}$ be as in Definition 1.1. For any $z\geqslant 2$ and $m_{p}$ as in Definition 1.2, let $W_{z}$ be as in Definition 1.2 and assume that $\widetilde{\varepsilon}(z)\to 0$ . Then there exist $c,\delta>0$ that only depend on $f$ such that for any function $k:\mathbb{N}\to\mathbb{C}$ and all large enough $P$ we have

	$\displaystyle\frac{1}{P^{n}}\sum_{\begin{subarray}{c}\mathbf{t}\in\mathbb{Z}^{% n}\cap P\mathscr{B}\\ sf(\mathbf{t})>0\end{subarray}}k(sf(\mathbf{t}))$	$\displaystyle=\Bigg{(}\int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}% \\ sf(\mathbf{t})>P^{-d}\end{subarray}}\omega(P^{d}sf(\mathbf{t}))\mathrm{d}% \mathbf{t}\Bigg{)}\sum_{\mathbf{t}\in(\mathbb{Z}/W_{z}\mathbb{Z})^{n}}\frac{% \rho(sf(\mathbf{t}),W_{z})}{W_{z}^{n-1}}$
		$\displaystyle+O\left(\frac{\\|k\\|_{1}}{P^{d}}(P^{-\delta}+\widetilde{% \varepsilon}(z)+z^{-c})+\frac{E(bP^{d},W_{z})W_{z}}{P^{d}}\right),$

where the implied constant is independent of $k$ . Here $b$ is defined in (1.2) and

\|k\|_{1}:=\sum_{\nu\in\mathbb{N}\cap[1,bP^{d}]}|k(\nu)|.

If the densities $\rho(a,W_{z})$ have a multiplicative structure then the sum over $\mathbf{t}$ is essentially a truncated Euler product over the primes $p\leqslant z$ . The main idea of the proof is to not go through the usual route of major and minor arcs of the circle method but instead approximate the values of $f(\mathbf{t})$ by a Cramér–Granville model. Let $\mathrm{e}(z)=\mathrm{e}^{2\pi iz}$ for $z\in\mathbb{C}$ . We shall show in Lemmas 2.6-2.7-2.8 that

{\int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}\\ sf(\mathbf{t})>P^{-d}\end{subarray}}\omega(P^{d}sf(\mathbf{t}))\mathrm{d}% \mathbf{t}=\int_{\gamma\in\mathbb{R}}\left(\int_{\mathscr{B}}\mathrm{e}(\gamma f% (\mathbf{t})\mathrm{d}\mathbf{t})\right)\int_{[P^{-d},b]}\omega(P^{d}\mu)% \mathrm{e}(-s\gamma\mu)\mathrm{d}\mu\mathrm{d}\gamma.}

(1.3)

This expression is useful in applications as it has $\omega$ and $f$ separated.

Remark.

The number of variables can be reduced by half if one has better error terms regarding the distribution of $k$ on arithmetic progressions. This was done in [20] for the indicator function of primes and square-free integers.

Remark.

A convenient feature of Theorem 1.1 is that it allows the user to make choices for $z$ and the exponents $m_{p}(z)$ ; this is useful in situations where equidistribution for $k$ is easier for moduli of specific factorisation. For example, in certain Diophantine applications it is much easier to prove equidistribution for square-full moduli. Furthermore, $z$ is allowed to go to infinity arbitrarily slow; this means that one only needs to prove equidistribution modulo very small moduli.

1.2. A more accessible version

Theorem 1.1 makes no assumptions on the arithmetic function $k$ . We give a version that is easier to use if $k$ has certain equidistribution properties:

Theorem 1.2.

Assume that the positive integers $d,n\in\mathbb{N}$ satisfy $n>2^{d}(d-1)$ and let $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ be a smooth form of degree $d$ . Let $s\in\{-1,1\}$ , let $\mathscr{B}\subset\mathbb{R}^{n}$ be as in Definition 1.1 and let $b$ be given by (1.2). Assume that $k:\mathbb{N}\to\mathbb{C}$ is any function for which

•

$\displaystyle\sup_{q\in\mathbb{N}}\sum_{a=1}^{q}|\rho(a,q)|<\infty$ or $\rho(a,q)\in\mathbb{R}\cap[0,\infty)$ for all $q\in\mathbb{N},a\in\mathbb{Z}/q\mathbb{Z}$ ;
•

$\displaystyle\int_{1}^{x}\omega(t)\mathrm{d}t$ is non-zero for all large enough $x$ ;
•

For each fixed $q\in\mathbb{N}$ we have $\displaystyle\lim_{x\to\infty}\frac{E(x,q)}{\displaystyle\int_{1}^{x}\omega(t)% \mathrm{d}t}=0$ ;
•

We have $\displaystyle\liminf_{P\to\infty}P^{d}\Bigg{|}\int\limits_{\begin{subarray}{c}% \mathbf{t}\in\mathscr{B}:sf(\mathbf{t})>P^{-d}\end{subarray}}\omega(P^{d}sf(% \mathbf{t}))\mathrm{d}\mathbf{t}\Bigg{|}\cdot\Bigg{|}\int_{1}^{bP^{d}}\omega(t% )\mathrm{d}t\Bigg{|}^{-1}\neq 0$ ;
•

We have $\displaystyle\limsup_{x\to\infty}\frac{\displaystyle\sum_{m\leqslant x}|k(m)|}% {\displaystyle\left|\int_{1}^{x}\omega(t)\mathrm{d}t\right|}\neq\infty$ .

Define for any $z>1$ the integer $T_{z}=\prod_{p\leqslant z}p^{t_{p}(z)}$ , where $t_{p}(z)$ are arbitrary positive integers such that for each fixed prime $p$ one has $\lim_{z\to\infty}t_{p}(z)=\infty$ . Then the limit

{\sigma(f):=\lim_{z\to\infty}T_{z}^{-n+1}\sum_{\mathbf{t}\in(\mathbb{Z}/T_{z}% \mathbb{Z})^{n}}\rho(sf(\mathbf{t}),T_{z})}

(1.4)

exists and is independent of the choice of $t_{p}(z)$ . Furthermore,

\lim_{P\to\infty}\frac{1}{P^{n}\displaystyle\int\limits_{\begin{subarray}{c}% \mathbf{t}\in\mathscr{B}:sf(\mathbf{t})>P^{-d}\end{subarray}}\omega(P^{d}sf(% \mathbf{t}))\mathrm{d}\mathbf{t}}\sum_{\begin{subarray}{c}\mathbf{t}\in\mathbb% {Z}^{n}\cap P\mathscr{B}\\ sf(\mathbf{t})>0\end{subarray}}k(sf(\mathbf{t}))=\displaystyle\sigma(f).

A useful feature of Theorem 1.2 is that we only assume equidistribution in arithmetic progressions of fixed modulus and without an explicit error term.

1.3. Analytic Hasse principle for intersections

We say that a family of systems of polynomial equations with integer coefficients satisfies the analytic Hasse principle if the number of integer solutions in an expanding box converges to a constant, strictly positive, multiple of the analogous real and $p$ -adic densities. We will prove this for a specific intersection of two degree $d$ polynomials. Our intersection is given by

{a^{2}+b^{2}=f(x_{1},\ldots,x_{n}),c^{2}+d^{2}=1+f(x_{1},\ldots,x_{n}),}

(1.5)

where $f$ is smooth homogeneous of degree $d$ in $n>2^{d}(d-1)$ variables. There is recent analytic work in this area by Rydin-Myerson [36] who studied intersections of $R$ forms. For intersections of two cubic forms this was recently improved by Northey–Vishe [32] by relaxing the assumption on the number of variables so as to work provided that there are at least $39$ . For systems of two diagonal cubic forms work of Brüdern and Wooley [11] proves asymptotics when one has at least $13$ variables. The three simplest cases of our next result are

•

two quadratic equations in at least $9$ variables,
•

two cubic equations in at least $21$ variables,
•

two quartic equations in at least $53$ variables.

For any box $\mathscr{B}\subset[-1,1]^{n}$ and any $P\geqslant 1$ let

N_{f}(P):=\sum_{\begin{subarray}{c}\mathbf{x}\in\mathbb{Z}^{n}\cap P\mathscr{B% }\\ f(\mathbf{x})>0\end{subarray}}r(f(\mathbf{x}))r(1+f(\mathbf{x})),

where $r$ is the number of representations as the sum of two integer squares.

Theorem 1.3.

Assume that $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ is a smooth homogeneous degree $d$ polynomial in $n>2^{d}(d-1)$ variables. For any box $\mathscr{B}\subset[-1,1]^{n}$ inside of which $f$ only takes non-negative values we have

\lim_{P\to\infty}\frac{N_{f}(P)}{P^{n}}=\pi^{2}\mathrm{vol}\left(\{\mathbf{t}% \in\mathscr{B}\,:\,f(\mathbf{t})>0\}\right)\prod_{p}\sigma_{p}(f),

where $\sigma_{p}(f)$ is given by

\lim_{k\to\infty}p^{-k(n+2)}\sharp\left\{(\mathbf{t},x_{1},x_{2},y_{1},y_{2})% \in(\mathbb{Z}/p^{k}\mathbb{Z})^{n+4}:x_{1}^{2}+x_{2}^{2}=f(\mathbf{t}),y_{1}^% {2}+y_{2}^{2}=1+f(\mathbf{t})\right\}.

In particular, the analytic Hasse Principle holds for (1.5) as soon as $f$ assumes at least one strictly positive value.

The proof relies on proving a shifted convolution result (Theorem 3.1) and then feeding it into Theorem 1.1. The proof of Theorems 3.1 and 1.3 are respectively in §3.4, 3.5 and §3.6.

1.4. Subsets of integers

During the $2024$ conference on rational points in BIRS Chennai Pieropan [7] asked the following (see [7]): assume we are given a set $\mathscr{A}\subset\mathbb{Z}$ whose density we can estimate as

\sharp\{a\in\mathscr{A}:|a|\leqslant T\}=f(T)+O(g(T)),

where $f(T)$ has the form $c_{0}T^{c}(\log T)^{e}$ for some real constants $c_{0},c$ and $e$ and where

\displaystyle\lim_{T\to\infty}g(T)/(T^{c}(\log T)^{e})=0.

For an arbitrary polynomial $F\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ define

N(T;\mathscr{A},F)=\sharp\left\{\mathbf{x}\in\mathbb{Z}^{n}\cap[-T,T]^{n}:F(% \mathbf{x})\in\mathscr{A}\right\}.

The question is under what conditions on $\mathscr{A},f,g,n$ and $\deg(F)$ can one prove asymptotics for $N(T;\mathscr{A},F)$ as $T\to\infty$ ?

Our result states that the condition $n>2^{\deg(F)}(\deg(F)-1)$ is sufficient, as long as one has the additional assumption that we can count asymptotically the elements of $\mathscr{A}$ on any fixed arithmetic progression $r\left(\textnormal{mod}\ q\right)$ . To see why this is necessary, just consider the polynomial $F(x)=r+qx$ . For convenience of notation we will assume that $\mathscr{A}$ consists of strictly positive integers; in applications this does not cause any problems as one can split in cases according to the sign.

Theorem 1.4.

Assume that the positive integers $d,n\in\mathbb{N}$ satisfy $n>2^{d}(d-1)$ and let $F\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ be a smooth form of degree $d$ . Let $\mathscr{B}\subset\mathbb{R}^{n}$ be as in Definition 1.1 and let $b$ be given by (1.2). Assume that $\mathscr{A}\subset\mathbb{Z}_{\geqslant 0}$ is a non-empty set for which there exists $\omega:[1,\infty)\to(0,\infty)$ in $\mathscr{C}^{1}$ and for all $q\in\mathbb{N}$ and $r\in\mathbb{Z}/q\mathbb{Z}$ a real number $\rho(q,r)\geqslant 0$ such that we have

\lim_{x\to\infty}\sup_{y\in\mathbb{R}\cap[1,x]}\max_{a\in\mathbb{Z}/q\mathbb{Z% }}\left|\frac{\sharp\{a\in\mathscr{A}:a\leqslant y,a\equiv r\left(\textnormal{% mod}\ q\right)\}}{\displaystyle\int_{1}^{y}\omega(t)\mathrm{d}t}-\rho(r,q)% \right|=0.

Assume, in addition, that

\liminf_{P\to\infty}P^{d}\Bigg{|}\int\limits_{\begin{subarray}{c}\mathbf{t}\in% \mathscr{B}:|F(\mathbf{t})|>P^{-d}\end{subarray}}\omega(P^{d}F(\mathbf{t}))% \mathrm{d}\mathbf{t}\Bigg{|}\cdot\Bigg{|}\int_{1}^{bP^{d}}\omega(t)\mathrm{d}t% \Bigg{|}^{-1}\neq 0.

\sigma(F):=\lim_{z\to\infty}T_{z}^{-n+1}\sum_{\mathbf{t}\in(\mathbb{Z}/T_{z}% \mathbb{Z})^{n}}\rho(F(\mathbf{t}),T_{z})

exists and is independent of the choice of $t_{p}(z)$ . Furthermore,

\lim_{P\to\infty}\frac{\sharp\{\mathbf{x}\in\mathbb{Z}^{n}\cap P\mathscr{B}:F(% \mathbf{x})\in\mathscr{A}\}}{P^{n}\displaystyle\int\limits_{\begin{subarray}{c% }\mathbf{t}\in\mathscr{B}:F(\mathbf{t})>P^{-d}\end{subarray}}\omega(P^{d}F(% \mathbf{t}))\mathrm{d}\mathbf{t}}=\displaystyle\sigma(F).

The proof is given in §3.7 and is a direct application of Theorem 1.2.

1.5. Points on non-algebraic varieties

Exponential Diophantine equations form a considerable area of research in number theory with many results directed at proving finiteness of solutions, see for example the work of Bugeaud–Mignotte–Siksek [12] and the book of Shorey and Tijdeman [38]. Here we change perspective and count asymptotically the number of integer solutions of

2^{x_{0}}=f(x_{1},\ldots,x_{n}),

where $f$ is an integer polynomial in many variables.

Theorem 1.5.

For every smooth homogeneous $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ of degree $d$ with $n>2^{d}(d-1)$ and for every box $\mathscr{B}\subset\mathbb{R}^{n}$ as in Definition 1.1 inside which $f$ assumes at least one strictly positive value, we have

\lim_{P\to\infty}\frac{\sharp\left\{\mathbf{x}\in\mathbb{Z}^{n}\cap P\mathscr{% B},x_{0}\in\mathbb{N}\cap\left[1,\log(bP^{d})/\log 2\right]:f(\mathbf{x})=2^{x% _{0}}\right\}}{P^{n-d}\log P}=\frac{d\sigma_{\infty}(f)}{\log 2}\sigma_{2}(f)% \prod_{\begin{subarray}{c}p>2\end{subarray}}\sigma_{p}^{\prime}(f),

where $\sigma_{\infty}(f),\sigma_{2}(f)$ are the usual Hardy–Littlewood real and $2$ -adic densities associated to $f=0$ , $b$ is as in (1.1) and for an odd prime $p$ we have

\sigma_{p}^{\prime}(f):=\lim_{m\to\infty}\frac{\sharp\left\{\mathbf{y}\in(% \mathbb{Z}/p^{m}\mathbb{Z})^{n},h\in\mathbb{Z}/\varphi(p^{m})\mathbb{Z}:f(% \mathbf{y})=2^{h}\right\}}{(p-1)p^{mn}}.

The proof is given in §3.10 by applying Theorem 1.2.

1.6. Chowla’s conjecture

Let $\mu$ denote the Möbius function. Chowla’s conjecture states that for any $f\in\mathbb{Z}[x]$ for which $f\neq c_{1}g^{2}$ for $c_{1}\in\mathbb{Q}$ and $g\in\mathbb{Z}[x]$ one has

\frac{1}{P}\sum_{\begin{subarray}{c}m\in\mathbb{N}\cap[1,P]\\ f(m)\neq 0\end{subarray}}\mu(|f(\mathbf{m})|)\to 0

as $P\to\infty$ . It has been the focus of intensive investigation recently, see the work of Matomäki–Radziwiłł–Tao [29] and Tao–Teräväinen [39], for example. We prove the analogue of the conjecture for all polynomials in a sufficiently large number of variables.

Theorem 1.6.

Chowla’s conjecture holds for all smooth homogeneous $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ of degree $d$ when $n>2^{d}(d-1)$ . In particular, for every box $\mathscr{B}\subset\mathbb{R}^{n}$ as in Definition 1.1 we have

\lim_{P\to\infty}\frac{1}{P^{n}}\sum_{\begin{subarray}{c}\mathbf{x}\in\mathbb{% Z}^{n}\cap P\mathscr{B}\\ f(\mathbf{x})\neq 0\end{subarray}}\mu(|f(\mathbf{x})|)=0.

The proof is given in §3.1.

1.7. Sums of Fourier coefficients of modular forms over values of multivariable polynomials

Asymptotics for averages of coefficients of modular forms over values of irreducible quadratic polynomials in one variable have been studied by Blomer [2], Templier [40] and Templier–Tsimerman [41]. Bounds in the case of higher degree polynomials were later given by Chiriac and Yang [14].

Let $f$ be a holomorphic cusp form of weight $k$ for $\textrm{SL}_{2}(\mathbb{Z})$ with Fourier expansion

f(z)=\sum_{m=1}^{\infty}\lambda_{f}(m)m^{(k-1)/2}\mathrm{e}(mz)

for $\Im(z)>0$ . Let $f$ be a normalized Hecke eigenform so that $\lambda_{f}(1)=1$ . We have the bound $\lambda_{f}(n)\ll\tau(n)$ by the work the Ramanujan–Petersson conjecture, proved by Deligne [17] when $k\geqslant 2$ and by Deligne–Serre[18] when $k=1$ .

Theorem 1.7.

For all smooth homogeneous $F\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ of degree $d$ with $n>2^{d}(d-1)$ and all boxes $\mathscr{B}\subset\mathbb{R}^{n}$ as in Definition 1.1 we have

\lim_{P\to\infty}\frac{1}{P^{n}}\sum_{\begin{subarray}{c}\mathbf{x}\in\mathbb{% Z}^{n}\cap P\mathscr{B}\\ F(\mathbf{x})\neq 0\end{subarray}}\lambda_{f}(|F(\mathbf{x})|)=0.

The proof is given in §3.2.

1.8. Divisor function over values of polynomials

Let us continue with estimating the average of the divisor function $\tau$ over polynomials. In the case of one variable polynomials, asymptotics are known only for polynomials of degree $1$ and $2$ (see Hooley [26]), while, for binary forms asymptotics are known for degree $\leqslant 4$ by work of Daniel [15]. For a prime $p$ and positive integers $m,\ell$ we define $\tau_{p^{m}}(\ell):=1+\min\{v_{p}(\ell),m\}$ . This is a local model for the divisor function $\tau$ since $\tau(\ell)=\prod_{p\mid\ell}\lim_{m\to\infty}\tau_{p^{m}}(\ell)$ .

Theorem 1.8.

For every smooth homogeneous $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ of degree $d$ with $n>2^{d}(d-1)$ and for every box $\mathscr{B}\subset\mathbb{R}^{n}$ as in Definition 1.1 we have

\lim_{P\to\infty}\frac{1}{P^{n}\log P}\sum_{\begin{subarray}{c}\mathbf{x}\in% \mathbb{Z}^{n}\cap P\mathscr{B}\\ f(\mathbf{x})\neq 0\end{subarray}}\tau(|f(\mathbf{x})|)=d\mathrm{vol}(\mathscr% {B})\prod_{p}\left(1-\frac{1}{p}\right)\lim_{m\to\infty}\frac{1}{p^{nm}}\sum_{% \mathbf{t}\in(\mathbb{Z}/p^{m})^{n}}\tau_{p^{m}}(f(\mathbf{t})).

The proof is given in §3.3 as an application of an equidistribution result of Pongsriiam–Vaughan [35] and Theorem 1.2.

1.9. Rational points on a hyperelliptic hypersurface

For an integer polynomial $f$ of degree at least $5$ the hyperelliptic curve $y^{2}=f(x)$ has finitely many rational points by Faltings’ theorem [22]. On the contrary when $f$ has many variables one expects infinitely many rational points; here we shall give an asymptotic estimate.

Theorem 1.9.

Fix an integer $k\geqslant 2$ . For every smooth homogeneous $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ of degree $d$ with $n>2^{d}(d-1)$ and for every box $\mathscr{B}\subset\mathbb{R}^{n}$ as in Definition 1.1 inside which $f$ assumes at least one strictly positive value, we have

	$\displaystyle\lim_{P\to\infty}$	$\displaystyle\frac{\sharp\left\{\mathbf{x}\in\mathbb{Z}^{n}\cap P\mathscr{B}:f% (\mathbf{x})\textrm{ is a }k\textrm{-th power of a positive integer}\right\}}{% P^{n-d(1-1/k)}}$
		$\displaystyle=\prod_{p}\lim_{m\to\infty}\frac{\sharp\{y,\mathbf{t}\in(\mathbb{% Z}/p^{m}\mathbb{Z})^{n+1}:y^{k}=f(\mathbf{t})\}}{p^{mn}}\lim_{\delta\to 0_{+}}% \int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}\\ f(\mathbf{t})>\delta\end{subarray}}\frac{\mathrm{d}\mathbf{t}}{kf(\mathbf{t})^% {1-1/k}}.$

The proof is given in §3.8 as an application of Theorem 1.2.

1.10. $m$ -full numbers represented by polynomials

Let $m\in\mathbb{N}$ , $m\geqslant 2$ . An integer is called $m$ -full if it is divided by the $m$ -th power of every of its prime divisors. Values of polynomials that are $m$ -full have been studied by various authors; for example, Pasten [33] proved that an integer polynomial should take very few $m$ -full values conditionally on Vojta’s conjecture and have applications to the study of Campana points (see for example [34]). The picture for multivariable polynomials is rather different; indeed here one expects infinitely many solutions. We provide asymptotics.

Theorem 1.10.

Fix an integer $m\geqslant 2$ . For every smooth homogeneous $f\in\mathbb{Z}[x_{1},\ldots,x_{n}]$ of degree $d$ with $n>2^{d}(d-1)$ and for every box $\mathscr{B}\subset\mathbb{R}^{n}$ as in Definition 1.1 inside which $f$ assumes at least one strictly positive value, we have

\lim_{P\to\infty}\frac{\sharp\left\{\mathbf{x}\in\mathbb{Z}^{n}\cap P\mathscr{% B}:f(\mathbf{x})\textrm{ is }m\textrm{-full and positive}\right\}}{P^{n-d(1-1/% m)}}=\sigma(f)\lim_{\delta\to 0_{+}}\int\limits_{\begin{subarray}{c}\mathbf{t}% \in\mathscr{B}\\ f(\mathbf{t})>\delta\end{subarray}}\frac{\mathrm{d}\mathbf{t}}{mf(\mathbf{t})^% {1-1/m}},

where $\sigma(f)$ is given by

\sum_{\begin{subarray}{c}k_{2},\ldots,k_{m}\in\mathbb{N}\end{subarray}}\frac{% \mu(k_{2}\cdots k_{m})^{2}}{k_{2}^{1+1/m}k_{3}^{1+2/m}\cdots k_{m}^{2-1/m}}% \prod_{p}\lim_{m\to\infty}\frac{\sharp\{y,\mathbf{t}\in(\mathbb{Z}/p^{m}% \mathbb{Z})^{n+1}:y^{m}k_{2}^{m+1}\cdots k_{m}^{2m-1}=f(\mathbf{t})\}}{p^{mn}}.

The proof is given in §3.9 as an application of Theorem 1.2.

Acknowledgements.

We would like to thank Ping Xi for asking questions during a seminar in Xi’an Jiaotong University that started activities that led to Theorem 1.1. Parts of this investigation too place when the second named author visited the university Paris-Saclay during April $2022$ , the generous hospitality and support of which is greatly appreciated.

2. Proof of Theorems 1.1 and 1.2

The Hardy–Littlewood singular series for the problem of representing an integer $\nu$ by the polynomial $f$ is given by

\mathfrak{S}(\nu):=\prod_{\begin{subarray}{c}p=2\\ p\textrm{ prime}\end{subarray}}^{\infty}\lim_{m\to\infty}\frac{\sharp\{\mathbf% {t}\in(\mathbb{Z}/p^{m}\mathbb{Z})^{n}:f(\mathbf{t})\equiv\nu\left(\textnormal% {mod}\ p^{m}\right)\}}{p^{m(n-1)}}.

We start by approximating it by its multiplicative model

\mathfrak{S}^{\flat}(\nu):=\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/W_{z}\mathbb% {Z})^{n}:f(\mathbf{t})\equiv\nu\left(\textnormal{mod}\ p^{m}\right)\}}{W_{z}^{% n-1}}.

Lemma 2.1.

Let $f$ be as in Theorem 1.1 and assume that $\widetilde{\varepsilon}(z)\to 0$ . There exists $c=c(f)>0$ such that for all $\nu\in\mathbb{Z}$ we have

\mathfrak{S}(\nu)-\mathfrak{S}^{\flat}(\nu)\ll\widetilde{\varepsilon}(z)+z^{-c},

where the implied constant is independent of $\nu$ .

Proof.

By [1, §7] we have

\lim_{m\to\infty}\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/p^{m}\mathbb{Z})^{n}:f% (\mathbf{t})=\nu\}}{p^{m(n-1)}}=1+\sum_{r\in\mathbb{N}}p^{-rn}\sum_{\begin{% subarray}{c}a\in\mathbb{Z}/p^{r}\mathbb{Z}\\ p\nmid a\end{subarray}}S_{a,p^{r}}\mathrm{e}(-a\nu/p^{r}),

where

{S_{a,q}:=\sum_{\mathbf{t}\in(\mathbb{Z}/q\mathbb{Z})^{n}}\mathrm{e}\left(% \frac{a}{q}f(\mathbf{t})\right).}

(2.1)

The assumption $n>2^{d}(d-1)$ shows that the constant

{c:=\frac{n2^{-d}}{(d-1)}-1}

(2.2)

is strictly positive. Using [1, pg. 256, Eq. (19)] with $\varepsilon=c$ we see that

{|S_{a,p^{r}}|\ll p^{r\left(n-c-2\right)}}

(2.3)

with an implied constant that depends at most on $f$ . We infer that for all $m\geqslant 1$ one has

{\left|\sum_{r\geqslant m}p^{-rn}\sum_{\begin{subarray}{c}a\in\mathbb{Z}/p^{r}% \mathbb{Z}\\ p\nmid a\end{subarray}}S_{a,p^{r}}\mathrm{e}(-a\nu/p^{r})\right|\ll\sum_{r% \geqslant m}p^{-r(1+c)}\ll p^{-m(1+c)},}

(2.4)

where the implied constant depends at most on $f$ . In the special case that $m=1$ this yields

{\left|\lim_{m\to\infty}\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/p^{m}\mathbb{Z}% )^{n}:f(\mathbf{t})\equiv\nu\left(\textnormal{mod}\ p^{m}\right)\}}{p^{m(n-1)}% }-1\right|\ll p^{-(1+c)}}

(2.5)

uniformly in $\nu$ . Hence, for any $z>1$ we have

\log\prod_{p>z}\lim_{m\to\infty}\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/p^{m}% \mathbb{Z})^{n}:f(\mathbf{t})\equiv\nu\left(\textnormal{mod}\ p^{m}\right)\}}{% p^{m(n-1)}}=\sum_{p>z}\log\left(1+O(p^{-1-c})\right)\ll\sum_{p>z}p^{-1-c}\ll z% ^{-c}

and therefore

\prod_{p>z}\lim_{m\to\infty}\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/p^{m}% \mathbb{Z})^{n}:f(\mathbf{t})\equiv\nu\left(\textnormal{mod}\ p^{m}\right)\}}{% p^{m(n-1)}}=1+O\left(z^{-c}\right).

We have $\mathfrak{S}(\nu)=O(1)$ independently of $\nu$ due to [1, Corollary, pg.256]. Thus,

\mathfrak{S}(\nu)=\left(\prod_{p\leqslant z}\lim_{m\to\infty}\frac{\sharp\{% \mathbf{t}\in(\mathbb{Z}/p^{m}\mathbb{Z})^{n}:f(\mathbf{t})\equiv\nu\left(% \textnormal{mod}\ p^{m}\right)\}}{p^{m(n-1)}}\right)+O(z^{-c})

with an implied constant independent of $\nu$ . By (2.4) we see that for any $m_{p}\geqslant 1$ ,

\lim_{m\to\infty}\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/p^{m}\mathbb{Z})^{n}:f% (\mathbf{t})\equiv\nu\left(\textnormal{mod}\ p^{m}\right)\}}{p^{m(n-1)}}-\sum_% {0\leqslant r\leqslant m_{p}}p^{-rn}\!\!\!\sum_{\begin{subarray}{c}a\in\mathbb% {Z}/p^{r}\mathbb{Z}\\ p\nmid a\end{subarray}}S_{a,p^{r}}\mathrm{e}(-a\nu/p^{r})\ll p^{-(1+c)(1+m_{p}% )}.

This means that

\prod_{p\leqslant z}\lim_{m\to\infty}\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/p^% {m}\mathbb{Z})^{n}:f(\mathbf{t})\equiv\nu\left(\textnormal{mod}\ p^{m}\right)% \}}{p^{m(n-1)}}=V\prod_{p\leqslant z}\left(\sum_{0\leqslant r\leqslant m_{p}}p% ^{-rn}\sum_{\begin{subarray}{c}a\in\mathbb{Z}/p^{r}\mathbb{Z}\\ p\nmid a\end{subarray}}S_{a,p^{r}}\mathrm{e}(-a\nu/p^{r})\right),

for some quantity $V$ that satisfies

V=\prod_{p\leqslant z}\left(1+O\left(p^{-(1+c)(1+m_{p})}\right)\right)=\exp% \left(O\left(\sum_{p\leqslant z}p^{-(1+c)(1+m_{p})}\right)\right)=1+O(% \widetilde{\varepsilon}(z))

due to the assumption $\widetilde{\varepsilon}(z)\to 0$ . Recalling that $W_{z}=\prod_{p\leqslant z}p^{m_{p}},$ the last equation in [1, pg. 259] gives

{\sum_{0\leqslant r\leqslant m_{p}}p^{-rn}\sum_{\begin{subarray}{c}a\in\mathbb% {Z}/p^{r}\mathbb{Z}\\ p\nmid a\end{subarray}}S_{a,p^{r}}\mathrm{e}(-a\nu/p^{r})=\frac{\sharp\{% \mathbf{t}\in(\mathbb{Z}/p^{m_{p}})^{n}:f(\mathbf{t})\equiv\nu\left(% \textnormal{mod}\ p^{m}\right)\}}{p^{(n-1)m_{p}}},}

(2.6)

which concludes the proof via the Chinese remainder theorem. ∎

The Hardy–Littlewood singular integral for representing $\nu$ by $f$ is

J(\nu):=\int_{\mathbb{R}}I(\mathscr{B};\gamma)\mathrm{e}(-\gamma\nu)\mathrm{d}% \gamma,\ \ \ \textrm{ where }\ \ \ I(\mathscr{B};\gamma):=\int_{\mathscr{B}}% \mathrm{e}\left(\gamma f(\mathbf{t})\right)\mathrm{d}\mathbf{t}.

We are now in position to start averaging the arithmetic function over the polynomial values.

Lemma 2.2.

Let $f$ be as in Theorem 1.1 and assume that $\widetilde{\varepsilon}(z)\to 0$ . There exists $c=c(f)>0$ such that for all large $P$ one has

	$\displaystyle\sum_{\begin{subarray}{c}\mathbf{t}\in\mathbb{Z}^{n}\cap P% \mathscr{B}\\ sf(\mathbf{t})>0\end{subarray}}k(sf(\mathbf{t}))=$	$\displaystyle\frac{P^{n-d}}{W_{z}^{n-1}}\sum_{\mathbf{t}\in(\mathbb{Z}/W_{z})^% {n}}\int_{\mathbb{R}}I(\mathscr{B};\gamma)\left\{\sum_{\begin{subarray}{c}\nu% \in\mathbb{N}\cap[1,bP^{d}]\\ \nu\equiv sf(\mathbf{t})\left(\textnormal{mod}\ W_{z}\right)\end{subarray}}k(% \nu)\mathrm{e}(-\gamma s\nu P^{-d})\right\}\mathrm{d}\gamma$
	$\displaystyle+$	$\displaystyle O\left(P^{n-d}\left(P^{-\delta}+\widetilde{\varepsilon}(z)+z^{-c% }\right)\\|k\\|_{1}\right).$

Proof.

For $P$ suitable large and all $\mathbf{t}\in P\mathscr{B}$ we have $|f(\mathbf{t})|\leqslant bP^{d}$ , hence,

{\sum_{\begin{subarray}{c}\mathbf{t}\in\mathbb{Z}^{n}\cap P\mathscr{B}\\ sf(\mathbf{t})>0\end{subarray}}k(sf(\mathbf{t}))=\sum_{\nu\in\mathbb{N}\cap[1,% bP^{d}]}k(\nu)\sharp\{\mathbf{t}\in\mathbb{Z}^{n}\cap P\mathscr{B}:f(\mathbf{t% })=s\nu\}.}

(2.7)

By [1, Lem.5.5] there is $\delta=\delta(f)>0$ such that

\sharp\{\mathbf{t}\in\mathbb{Z}^{n}\cap P\mathscr{B}:f(\mathbf{t})=s\nu\}=P^{n% -d}\mathfrak{S}(s\nu)J(s\nu P^{-d})+O(P^{n-d-\delta}),

where the implied constant is independent of $s$ and $\nu$ . By Lemma 2.1 we can write the right-hand side as

P^{n-d}\mathfrak{S}^{\flat}(s\nu)\int_{\mathbb{R}}I(\mathscr{B};\gamma)\mathrm% {e}(-\gamma s\nu P^{-d})\mathrm{d}\gamma+O\left(P^{n-d}\left(P^{-\delta}+% \widetilde{\varepsilon}(z)+z^{-c}\right)\right),

where we used that $J(s\nu P^{-d})$ is bounded only in terms of $f$ . Substituting into (2.7) we obtain the claimed error term. In the ensuing main term we can interchange the order of summation and integration owing to the fact that $\int_{\mathbb{R}}|I(\mathscr{B};\gamma)|\mathrm{d}\gamma<\infty$ by [1, Lem. 5.2]. This gives

P^{n-d}\int_{\mathbb{R}}I(\mathscr{B};\gamma)\left\{\sum_{\nu\in\mathbb{N}\cap% [1,bP^{d}]}k(\nu)\mathfrak{S}^{\flat}(s\nu)\mathrm{e}(-\gamma s\nu P^{-d})% \right\}\mathrm{d}\gamma.

By definition of the function $\mathfrak{S}^{\flat}(s\nu)$ , we obtain

P^{n-d}W_{z}^{-n+1}\int_{\mathbb{R}}I(\mathscr{B};\gamma)\sum_{\mathbf{t}\in(% \mathbb{Z}/W_{z}\mathbb{Z})^{n}}\left\{\sum_{\begin{subarray}{c}\nu\in\mathbb{% N}\cap[1,bP^{d}]\\ \nu\equiv sf(\mathbf{t})\left(\textnormal{mod}\ W_{z}\right)\end{subarray}}k(% \nu)\mathrm{e}(-\gamma s\nu P^{-d})\right\}\mathrm{d}\gamma,

as claimed. ∎

The presence of $I(\mathscr{B};\gamma)$ in the main term shows that only small $\gamma$ matter. Therefore, a simple approach to deal with those will suffice.

Lemma 2.3.

For any $x\geqslant 1$ , $a,q\in\mathbb{Z}$ with $q\neq 0$ and any $\mu\in\mathbb{R}$ we have

\sum_{\begin{subarray}{c}1\leqslant\nu\leqslant x\\ \nu\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}k(\nu)\mathrm{e}(\mu% \nu)=\rho(a,q)\int_{1}^{x}\omega(t)\mathrm{e}(\mu t)\mathrm{d}t+O\big{(}(1+|% \mu|x)E(x,q)\big{)},

where the implied constant is absolute.

Proof.

By partial summation we obtain

A(x)\mathrm{e}(\mu x)-2\pi i\mu\int_{1}^{x}A(t)\mathrm{e}(\mu t)\mathrm{d}t,\ % \ \ \textrm{ where }\ \ \ A(t):=\sum_{\begin{subarray}{c}1\leqslant\nu% \leqslant t\\ \nu\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}k(\nu).

Alluding to (1.1) we then obtain

\rho(a,q)\left\{\mathrm{e}(\mu x)\int_{1}^{x}\omega(t)\mathrm{d}t-2\pi i\mu% \int_{1}^{x}\int_{1}^{t}\omega(u)\mathrm{e}(\mu t)\mathrm{d}u\mathrm{d}t\right% \}+O\left((1+|\mu|x)E(x,q)\right)

and infer that the difference inside the brackets equals $\displaystyle\int_{1}^{x}\omega(t)\mathrm{e}(\mu t)\mathrm{d}t.$ ∎

Define

H:=\int_{\mathbb{R}}I(\mathscr{B};\gamma)\bigg{(}\int_{P^{-d}}^{b}\omega(P^{d}% y)\mathrm{e}(-\gamma sy)\mathrm{d}y\bigg{)}\mathrm{d}\gamma.

Lemma 2.4.

Keep the setting of Lemma 2.2. There exists $c=c(f)>0$ such that for all large $P$ one has

	$\displaystyle\sum_{\begin{subarray}{c}\mathbf{t}\in\mathbb{Z}^{n}\cap P% \mathscr{B}\\ sf(\mathbf{t})>0\end{subarray}}$	$\displaystyle k(sf(\mathbf{t}))-\frac{P^{n}H}{W_{z}^{n-1}}\sum_{\mathbf{t}\in(% \mathbb{Z}/W_{z})^{n}}\rho(sf(\mathbf{t}),W_{z})$
	$\displaystyle\ll$	$\displaystyle P^{n-d}W_{z}E(bP^{d},W_{z})+P^{n-d}\left(P^{-\delta}+\widetilde{% \varepsilon}(z)+z^{-c}\right)\\|k\\|_{1}.$

Proof.

We inject Lemma 2.3 into Lemma 2.2. The contribution of the error term is

\ll\frac{P^{n-d}}{W_{z}^{n-1}}\sum_{\mathbf{t}\in(\mathbb{Z}/W_{z})^{n}}\int_{% \mathbb{R}}|I(\mathscr{B};\gamma)|\left(1+\left|\gamma sP^{-d}\right|P^{d}% \right)E(bP^{d},W_{z})\mathrm{d}\gamma.

By [1, Lem. 5.2] with $\varepsilon$ being the positive constant $c$ in (2.2) we obtain

{|I(\mathscr{B};\gamma)|\ll\min\{1,|\gamma|^{-c-2}\},}

(2.8)

hence, $\displaystyle\int_{\mathbb{R}}|I(\mathscr{B};\gamma)|(1+|\gamma|)\mathrm{d}% \gamma<\infty$ . This gives the estimate

\ll E(bP^{d},W_{z})\frac{P^{n-d}}{W_{z}^{n-1}}\sum_{\mathbf{t}\in(\mathbb{Z}/W% _{z})^{n}}\int_{\mathbb{R}}I(\mathscr{B};\gamma)(1+|\gamma|)\mathrm{d}\gamma% \ll E(bP^{d},W_{z})P^{n-d}W_{z}.

∎

To deal with $H$ we generalise the Fourier analysis approach in [20]. For $k\in\mathbb{N}$ let $\varphi_{k}:\mathbb{R}\to\mathbb{R}$ be $\varphi_{k}(\gamma)=\pi^{-1/2}k\exp(-k^{2}\gamma^{2}),$ whose Fourier transform is $\widehat{\varphi}_{k}(\gamma)=\exp(-\pi^{2}k^{-2}\gamma^{2})$ . Define

H_{k}=\int_{\gamma\in\mathbb{R}}\widehat{\varphi}_{k}(\gamma)I(\mathscr{B};% \gamma)\int_{P^{-d}}^{b}\omega(P^{d}y)\mathrm{e}(-\gamma sy)\mathrm{d}y\mathrm% {d}\gamma.

Lemma 2.5.

For fixed $P\geqslant 1$ we have $H_{k}\to H$ as $k\to+\infty$ .

Proof.

By (2.8) and $|\widehat{\varphi}_{k}(\gamma)|=O(1)$ one sees that

\lim_{k\to\infty}\int_{|\gamma|>\log k}(\widehat{\varphi}_{k}(\gamma)-1)I(% \mathscr{B};\gamma)\int_{P^{-d}}^{b}\omega(P^{d}y)\mathrm{e}(-\gamma sy)% \mathrm{d}y\mathrm{d}\gamma=0.

In the region $|\gamma|\leqslant\log k$ one has

\widehat{\varphi}_{k}(\gamma)=\exp(-\pi^{2}k^{-2}\gamma^{2})=1+O(k^{-2}\gamma^% {2})=1+O\left(\frac{\log^{2}k}{k^{2}}\right),

with an implied constant that is independent of $\gamma$ . Hence,

\int_{|\gamma|>\log k}(\widehat{\varphi}_{k}(\gamma)-1)I(\mathscr{B};\gamma)% \int_{P^{-d}}^{b}\omega(P^{d}y)\mathrm{e}(-\gamma sy)\mathrm{d}y\mathrm{d}% \gamma\ll\frac{\log^{2}k}{k^{2}}\int_{\mathbb{R}}|I(\mathscr{B};\gamma)|% \mathrm{d}\gamma\ll\frac{\log^{2}k}{k^{2}},

which concludes the proof. ∎

For $k\in\mathbb{N}$ and $\mathbf{t}\in\mathscr{B}$ let

F_{k}(\mathbf{t}):=\int\limits_{\begin{subarray}{c}y\in[P^{-d},b]\\ |f(\mathbf{t})-sy|\leqslant 1/\sqrt{k}\end{subarray}}\omega(P^{d}y)\varphi_{k}% (f(\mathbf{t})-sy)\mathrm{d}y.

Lemma 2.6.

We have

H=\int_{\mathscr{B}}\lim_{k\to\infty}F_{k}(\mathbf{t})\mathrm{d}t.

Proof.

By the definition of $I(\mathscr{B};\boldsymbol{\gamma})$ we can write

H_{k}=\int_{\mathbf{t}\in\mathscr{B}}\int_{y\in[P^{-d},b]}\omega(P^{d}y)\int_{% \gamma\in\mathbb{R}}\widehat{\varphi}_{k}(\gamma)\mathrm{e}(\gamma\left\{f(% \mathbf{t})-sy\right\})\mathrm{d}\gamma\mathrm{d}y\mathrm{d}\mathbf{t},

where we swapped the order of integration owing to $\int_{\mathbb{R}}|\widehat{\varphi}_{k}(\boldsymbol{\gamma})|\mathrm{d}\gamma<\infty$ . By Fourier inversion

H_{k}=\int_{\mathbf{t}\in\mathscr{B}}\int_{y\in[P^{-d},b]}\omega(P^{d}y)% \varphi_{k}(f(\mathbf{t})-sy)\mathrm{d}y\mathrm{d}\mathbf{t}.

For $|y|>1/\sqrt{k}$ we see that

{\varphi_{k}(y)=\pi^{-1/2}k\exp(-k^{2}y^{2})\ll k\exp(-k).}

(2.9)

Hence, by Lemma 2.5 we get

H=\lim_{k\to\infty}\int_{\mathscr{B}}F_{k}(\mathbf{t})\mathrm{d}t.

To conclude the proof it will suffice to use the dominated convergence theorem. Its assumptions hold since $F_{k}(\mathbf{t})$ is bounded independently of $k$ and $\mathbf{t}$ . Indeed,

\left|F_{k}(\mathbf{t})\right|\ll\left(\max_{y\in[P^{-d},b]}|\omega(P^{d}y)|% \right)\int_{\begin{subarray}{c}\mathbb{R}\end{subarray}}\varphi_{k}(f(\mathbf% {t})-sy)\mathrm{d}y

and note that a linear change of variables makes the integral equal to $\displaystyle\int_{\mathbb{R}}\varphi_{k}(z)\mathrm{d}z=1$ . ∎

Lemma 2.7.

We have

H=\int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}\\ sf(\mathbf{t})>P^{-d}\end{subarray}}\lim_{k\to\infty}F_{k}(\mathbf{t})\mathrm{% d}t.

Proof.

The $\mathbf{t}$ with $sf(\mathbf{t})=P^{-d}$ have measure $0$ , hence, can be ignored. It suffices to show that if $sf(\mathbf{t})<P^{-d}$ then $\lim_{k}F_{k}(\mathbf{t})=0$ . Indeed, for such $\mathbf{t}$ we shall see that, for $k$ big enough, the range of integration in the definition of $F_{k}$ is empty. Indeed, if there is some $y$ in $F_{k}$ then

y=y-sf(\mathbf{t})+sf(\mathbf{t})\leqslant|y-sf(\mathbf{t})|+sf(\mathbf{t})% \leqslant\frac{1}{\sqrt{k}}+sf(\mathbf{t}),

hence, as $k\to\infty$ , we get $y<P^{-d}$ . This contradicts the definiton of $F_{k}$ , which includes the condition $y\in[P^{-d},b]$ . ∎

Lemma 2.8.

Let $\mathbf{t}\in\mathscr{B}$ . If $sf(\mathbf{t})>P^{-d}$ then $\displaystyle\lim_{k\to\infty}F_{k}(\mathbf{t})=\omega(sP^{d}f(\mathbf{t}))$ .

Proof.

For each $y$ such that $|y-sf(\mathbf{t})|\leqslant 1/\sqrt{k}$ we have

-\frac{1}{\sqrt{k}}+P^{-d}\leqslant-|y-sf(\mathbf{t})|+P^{-d}<y-sf(\mathbf{t})% +sf(\mathbf{t})\leqslant|y-sf(\mathbf{t})|+b\leqslant\frac{1}{\sqrt{k}}+b,

thus, as $k\to\infty$ , we can ignore the condition $y\in[P^{-d},b]$ in the definition of $F_{k}$ . Hence,

\lim_{k\to\infty}F_{k}(\mathbf{t})=\lim_{k\to\infty}\int_{|f(\mathbf{t})-sy|% \leqslant 1/\sqrt{k}}\omega(P^{d}y)\varphi_{k}(f(\mathbf{t})-sy)\mathrm{d}y.

Since $\omega$ is continuous and $\varphi_{k}$ has exponential decay, we infer that

\lim_{k\to\infty}F_{k}(\mathbf{t})=\omega(P^{d}sf(\mathbf{t}))\lim_{k\to\infty% }\int_{|f(\mathbf{t})-sy|\leqslant 1/\sqrt{k}}\varphi_{k}(f(\mathbf{t})-sy)% \mathrm{d}y.

It now suffices to prove that the last limit over $k\to\infty$ equals $1$ . By (2.9) we can dispose of the condition $|f(\mathbf{t})-sy|\leqslant 1/\sqrt{k}$ as $k\to+\infty$ and we can then use the fact $\displaystyle\int_{\mathbb{R}}\varphi_{k}(\mu)\mathrm{d}\mu=1.$ ∎

Bringing together Lemmas 2.6, 2.7 and 2.8 yields

H=\int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}\\ sf(\mathbf{t})>P^{-d}\end{subarray}}\omega(P^{d}sf(\mathbf{t}))\mathrm{d}% \mathbf{t}.

Injecting this into the main term present in Lemma 2.4 concludes the proof of Theorem 1.1.

2.1. Proof of Theorem 1.2

We first prove that if $\rho(a,q)\geqslant 0$ for all $a,q$ then $\displaystyle\sup_{q\in\mathbb{N}}\sum_{a=1}^{q}|\rho(a,q)|<\infty$ . Thus, to complete the proof it suffices to work under the latter assumption only. We write

\sum_{m\leqslant x}k(m)=\sum_{a\in\mathbb{Z}/q\mathbb{Z}}\sum_{\begin{subarray% }{c}m\leqslant x\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}k(m)

and divide through by $\displaystyle\int_{1}^{x}\omega(t)\mathrm{d}t$ . Letting $x\to\infty$ and using the assumption

\lim_{x\to+\infty}\frac{E(x,q)}{\displaystyle\int_{1}^{x}\omega(t)\mathrm{d}t}=0

we deduce that for all $a$ and $g$ we have

\displaystyle\sum_{a\in\mathbb{Z}/q\mathbb{Z}}\rho(a,q)=\lim_{x\to+\infty}% \frac{\displaystyle\sum_{m\leqslant x}k(m)}{\displaystyle\int_{1}^{x}\omega(t)% \mathrm{d}t}.

(2.10)

Thus, if $\rho(a,q)\geqslant 0$ for all $a,q$ then $\displaystyle\sum_{a=1}^{q}|\rho(a,q)|$ is uniformly bounded for all $q$ .

Next we prove that the limit in (1.4) exists. Fixing the value of $f(\mathbf{t})$ we see that

T_{z}^{-n+1}\sum_{\mathbf{t}\in(\mathbb{Z}/T_{z}\mathbb{Z})^{n}}\rho(sf(% \mathbf{t}),T_{z})=\sum_{\nu=1}^{T_{z}}\rho(\nu,T_{z})\prod_{p\leqslant z}% \frac{\sharp\left\{\mathbf{t}\in(\mathbb{Z}/p^{t_{p}(z)}\mathbb{Z})^{n}:sf(% \mathbf{t})\equiv\nu\left(\textnormal{mod}\ p^{t_{p}(z)}\right)\right\}}{p^{t_% {p}(z)(n-1)}},

which, by Lemma 2.1 equals

\sum_{\nu=1}^{T_{z}}\rho(\nu,T_{z})\left(\mathfrak{S}(\nu)+o_{z\to\infty}(1)% \right).

By assumption, $\displaystyle\sum_{\nu=1}^{T_{z}}|\rho(\nu,T_{z})|$ is bounded independently of $z$ , hence, we obtain

\sum_{\nu=1}^{T_{z}}\rho(\nu,T_{z})\mathfrak{S}(\nu)+o_{z\to\infty}(1)=\sum_{q% =1}^{\infty}q^{-n}\sum_{\begin{subarray}{c}a\in\mathbb{Z}/q\mathbb{Z}\\ \gcd(a,q)=1\end{subarray}}S_{a,q}\sum_{\nu=1}^{T_{z}}\rho(\nu,T_{z})\mathrm{e}% (-a\nu/q)+o_{z\to\infty}(1),

where we used the equation before [1, Lemma 5.5] to express $\mathfrak{S}$ via $S_{a,q}$ and the bound (2.4) to change order of summation. The new series over $q$ converges absolutely due to (2.4) and

\left|\sum_{\nu=1}^{T_{z}}\rho(\nu,T_{z})\mathrm{e}(-a\nu/q)\right|\leqslant% \sum_{\nu=1}^{T_{z}}|\rho(\nu,T_{z})|<\infty.

This proves that the limit defining $\sigma(f)$ exists and is independent of the choice of $t_{p}(z)$ due to 2.10.

The remaining claims in Theorem 1.2 will be deduced from Theorem 1.1 with the following choice for $m_{p}(z)$ . For any prime $p$ , define $m_{p}(z)$ as the largest integer satisfying $p^{m_{p}(z)}\leqslant z$ and let $W_{z}=\displaystyle\prod_{p\leqslant z}p^{m_{p}(z)}$ . For a fixed integer $k\geqslant 1$ note that the primes $p$ for which $m_{p}(z)=k$ are exactly those in the interval $(z^{1/(k+1)},z^{1/k}]$ . Summing over all integers $m$ in that interval gives

	$\displaystyle\widetilde{\varepsilon}(z)$	$\displaystyle\leqslant\sum_{1\leqslant k\leqslant(\log z)/(\log 2)}\sum_{z^{1/% (k+1)}<m\leqslant z^{1/k}}m^{-k-1}\leqslant\sum_{1\leqslant k\leqslant(\log z)% /(\log 2)}\int_{z^{1/(k+1)}}^{\infty}\frac{\mathrm{d}t}{t^{k+1}}$
		$\displaystyle\ll\sum_{1\leqslant k\leqslant(\log z)/(\log 2)}\frac{z^{-k/(k+1)% }}{k}\ll\sum_{k\leqslant(\log z)/(\log 2)}\frac{1}{\sqrt{z}}\ll\frac{\log z}{% \sqrt{z}},$

hence, $\widetilde{\varepsilon}(z)=o(1)$ . Since $m_{p}(z)=[\log z/\log p]$ we infer that $m_{p}(z)\to\infty$ as $z\to\infty$ and $p$ is fixed. Fix any $\varepsilon>0$ . Then there exists $z$ such that

\left|\displaystyle\sigma(f)-W_{z}^{-n+1}\sum_{\mathbf{t}\in(\mathbb{Z}/W_{z}% \mathbb{Z})^{n}}\rho(sf(\mathbf{t}),W_{z})\right|<\varepsilon\ \ \textrm{ and % }\ \ |\widetilde{\varepsilon}(z)+z^{-c}|<\varepsilon.

Hence, by Theorem 1.1 we obtain

	$\displaystyle\Bigg{\|}\frac{1}{P^{n}\displaystyle\int\omega(P^{d}sf(\mathbf{t})% )\mathrm{d}\mathbf{t}}\sum_{\begin{subarray}{c}\mathbf{t}\in\mathbb{Z}^{n}\cap P% \mathscr{B}\\ sf(\mathbf{t})>0\end{subarray}}k(sf(\mathbf{t}))-\sigma(f)\Bigg{\|}$
	$\displaystyle\ll\varepsilon+\frac{\\|k\\|_{1}}{P^{d}\displaystyle\int\omega(P^{d% }sf(\mathbf{t}))\mathrm{d}\mathbf{t}}(P^{-\delta}+\varepsilon)+W_{z}\frac{E(bP% ^{d},W_{z})}{P^{d}\displaystyle\int\omega(P^{d}sf(\mathbf{t}))\mathrm{d}% \mathbf{t}},$

where the integral is over $\mathbf{t}\in\mathscr{B}$ with $sf(\mathbf{t})>P^{-d}$ and the implied constant depends only on $f$ and $k$ . By the last two assumptions in the theorem the error term becomes

\ll\varepsilon+P^{-\delta}+\varepsilon+W_{z}\frac{E(bP^{d},W_{z})}{P^{d}% \displaystyle\int_{1}^{bP^{d}}\omega(t)\mathrm{d}t}.

Using the third assumption shows that for $P$ big enough, we have

\Bigg{|}\frac{1}{P^{n}\displaystyle\int\omega(P^{d}sf(\mathbf{t}))\mathrm{d}% \mathbf{t}}\sum_{\begin{subarray}{c}\mathbf{t}\in\mathbb{Z}^{n}\cap P\mathscr{% B}\\ sf(\mathbf{t})>0\end{subarray}}k(sf(\mathbf{t}))-\sigma(f)\Bigg{|}\ll\varepsilon,

where the implied constant depends at most on $f$ and $k$ . Taking arbitrary small $\varepsilon$ concludes the proof.

3. Applications

3.1. The proof of Theorem 1.6

For fixed $B>0$ Davenport [16, Lemma 6] proved that

\sum_{\begin{subarray}{c}1\leqslant m\leqslant x\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}\mu(m)\ll\frac{x}{(% \log x)^{B}}

for all $x\geqslant 1$ , $q\in\mathbb{N}$ and $a\in\mathbb{Z}/q\mathbb{Z}$ , where the implied constant is independent of $x,a$ and $q$ . Letting $\rho(a,q)=0$ and $\omega(t)=1$ we see that $E(x,q)\ll x(\log x)^{-B}$ . One can now easily verify all assumptions of Theorem 1.2, which suffices for the proof.

3.2. The proof of Theorem 1.7

Recall the Wilton-type bound

\sup_{\alpha\in\mathbb{R}}\left|\sum_{m\leqslant x}\lambda_{f}(m)\mathrm{e}(% \alpha m)\right|\ll_{f}\sqrt{x}\log x

that is proved in [27, Theorem 5.3]. Using additive characters modulo $q$ we obtain

\sup_{a\left(\textnormal{mod}\ q\right)}\left|\sum_{\begin{subarray}{c}m% \leqslant x\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}\lambda_{f}(m)\right|% \ll_{f}\sqrt{x}\log x.

The proof now follows similarly as in §3.1

3.3. The proof of Theorem 1.8

By [35, Theorem 1.1] one has for all $x\geqslant 1$ and $a,q\in\mathbb{N}$ and fixed $\varepsilon>0$ that

\sum_{\begin{subarray}{c}m\leqslant x\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}\tau(m)=\frac{x}{q}% \sum_{r\mid q}\frac{c_{r}(a)}{r}\left(\log\frac{x}{r^{2}}+2\gamma-1\right)+O_{% \varepsilon}\left((x^{1/3}+q^{1/2})x^{\varepsilon}\right),

where $\gamma$ is Euler’s constant, $c_{r}(a)$ is Ramanujan’s sum and the implied constant depends at most on $\varepsilon$ . We apply Theorem 1.2, hence, we do not need to record the dependence of the error term on $a,q$ . The right hand side is $\rho(a,q)x\log x+O(x)$ , where $q\rho(a,q)=\sum_{r\mid q}c_{r}(a)/r$ and the implied constant depends at most on $q$ and $a$ . Dividing through by $x\log x$ and letting $x\to\infty$ shows that $\rho(a,q)\geqslant 0$ . Letting $\omega(t)=\log t$ one has

E(x,q)=\sum_{\begin{subarray}{c}m\leqslant x\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}\tau(m)-\rho(a,q)\int_% {1}^{x}\omega(t)\mathrm{d}t=O(x).

It is now easy to verify all remaining assumptions of Theorem 1.2. The real density in the main term is

\int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}\\ sf(\mathbf{t})>P^{-d}\end{subarray}}\log\left(P^{d}sf(\mathbf{t})\right)% \mathrm{d}\mathbf{t}=d(\log P)\mathrm{vol}(\mathbf{t}\in\mathscr{B}:sf(\mathbf% {t})>P^{-d})+O(1).

This is asymptotic to $d(\log P)\mathrm{vol}(\mathbf{t}\in\mathscr{B}:sf(\mathbf{t})>0)$ since $\mathrm{vol}(\mathbf{t}\in\mathscr{B}:|f(\mathbf{t})|\ll P^{-d})\ll P^{-d}$ by [8, Lemma 1.19].

It remains to study $\sigma(f)$ . Since $c_{r}(a)$ is multiplicative with respect to $r$ we write $\sigma(f)$ as

\prod_{p=2}^{\infty}\lim_{m\to\infty}p^{-(n-1)m}\sum_{\mathbf{t}\in(\mathbb{Z}% /p^{m}\mathbb{Z})^{n}}\rho(sf(\mathbf{t}),p^{m})=\prod_{p=2}^{\infty}\lim_{m% \to\infty}p^{-nm}\sum_{\mathbf{t}\in(\mathbb{Z}/p^{m}\mathbb{Z})^{n}}\sum_{r=0% }^{m}p^{-r}c_{p^{r}}(sf(\mathbf{t})).

Using the fact that $c_{p^{r}}(a)=c_{p^{r}}(b)$ when $a\equiv b\left(\textnormal{mod}\ p^{r}\right)$ we obtain

\sum_{r=0}^{m}p^{-r(n+1)}\sum_{\mathbf{x}\in(\mathbb{Z}/p^{r}\mathbb{Z})^{n}}c% _{p^{r}}(sf(\mathbf{x}))=\sum_{r=0}^{m}p^{-r(n+1)}\sum_{\begin{subarray}{c}a% \in\mathbb{Z}/p^{r}\mathbb{Z}\\ p\nmid a\end{subarray}}S_{a,p^{r}}.

(3.1)

By (2.6) we have the following for each $m\geqslant 1$ ,

\sum_{0\leqslant r\leqslant m}p^{-rn}\sum_{\begin{subarray}{c}a\in\mathbb{Z}/p% ^{r}\mathbb{Z}\\ p\nmid a\end{subarray}}S_{a,p^{r}}=\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/p^{m% })^{n}:f(\mathbf{t})=0\}}{p^{(n-1)m}}=:\Gamma_{m}

from which it immediately follows that for all $m\geqslant 1$ one has

p^{-mn}\sum_{\begin{subarray}{c}a\in\mathbb{Z}/p^{m}\mathbb{Z}\\ p\nmid a\end{subarray}}S_{a,p^{m}}=\Gamma_{m}-\Gamma_{m-1}.

Therefore, (3.1) is equal to

		$\displaystyle 1+\sum_{r=1}^{m}p^{-r}(\Gamma_{r}-\Gamma_{r-1})=\frac{\Gamma_{m}% }{p^{m+1}}+\left(1-\frac{1}{p}\right)\sum_{k=0}^{m}\frac{\Gamma_{k}}{p^{k}}$
	$\displaystyle=$	$\displaystyle\frac{\Gamma_{m}}{p^{m+1}}+\left(1-\frac{1}{p}\right)\sum_{k=0}^{% m}\frac{\sharp\{\mathbf{t}\in(\mathbb{Z}/p^{k})^{n}:f(\mathbf{t})=0\}}{p^{nk}}$
	$\displaystyle=$	$\displaystyle\frac{\Gamma_{m}}{p^{m+1}}+\left(1-\frac{1}{p}\right)\frac{1}{p^{% nm}}\sum_{\mathbf{t}\in(\mathbb{Z}/p^{m})^{n}}\sum_{k=0}^{m}\mathds{1}(p^{k}% \mid f(\mathbf{t})).$

Recalling the definition of $\tau_{p^{m}}(n)$ and taking the limit as $m\to\infty$ , this becomes

\left(1-\frac{1}{p}\right)\lim_{m\to\infty}\frac{1}{p^{nm}}\sum_{\mathbf{t}\in% (\mathbb{Z}/p^{m})^{n}}\tau_{p^{m}}(f(\mathbf{t})),

which concludes the proof of Theorem 1.8.

3.4. A shifted convolution problem in arithmetic progressions

Let $\chi$ be the non-principal Dirichlet character modulo $4$ so that Dirichlet convolution $1\ast\chi(n)$ equals $4$ times the number of representations of $n$ as a sum of two integer squares for any integer $n$ . For $a,q\in\mathbb{N}$ we shall estimate

{S(x;q,a):=\sum_{\begin{subarray}{c}m\leqslant x\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}r(m)r(m+1)}

(3.2)

under the conditions

4\mid q,\ \ v_{p}(a)\leqslant v_{p}(q)-1\ \forall p\mid q.

(3.3)

The case $q=1$ can be treated using a variety of methods, for example, the first proof of Estermann [21] uses Kloosterman sums while the proofs in [13] and [28, Corollary 15.12, part (2)] use spectral methods. We were unable to locate a reference that provides an explicit error term for general arithmetic progressions and we therefore give a proof here.

Theorem 3.1.

For all $a,q\in\mathbb{N},x\geqslant 1$ satisfying $q\leqslant x^{1/2}$ and (3.3) we have

S(x;q,a)=\pi^{2}x\frac{\eta_{q}(a)\eta_{q}(a+1)}{q^{3}}\mathds{1}(a\equiv 0,1% \left(\textnormal{mod}\ 4\right))\prod_{p\nmid q}\left(1-\frac{1}{p^{2}}\right% )+O\left(x^{5/6}(\log x)^{20}q^{3/2}\tau(q)^{4}\right),

where the implied constant is absolute and for any $b$ integer

\eta_{q}(b):=\sharp\left\{\mathbf{y}\in(\mathbb{Z}/q\mathbb{Z})^{2}:y_{1}^{2}+% y_{2}^{2}\equiv b\left(\textnormal{mod}\ q\right)\right\}.

(3.4)

As a first step we recall that $r(m)$ vanishes for $m\equiv 3\left(\textnormal{mod}\ 4\right)$ , thus,

{S(x;q,a)=S_{0}(x;q,a)+S_{1}(x;q,a),}

(3.5)

where

S_{i}(x;q,a):=\sum_{\begin{subarray}{c}m\leqslant x,m\equiv i\left(\textnormal% {mod}\ 4\right)\\ m\equiv a\left(\textnormal{mod}\ q\right)\end{subarray}}r(m)r(m+1).

Indeed, either $m$ is odd and for $r(m)r(m+1)$ not to vanish, $m\equiv 1\left(\textnormal{mod}\ 4\right)$ or $m+1$ is odd and for $r(m)r(m+1)$ not to vanish, $m\equiv 0\left(\textnormal{mod}\ 4\right)$ . By Dirichlet’s hyperbola trick when $k$ is odd we have

\frac{1}{4}r(k)=\sum_{de=k}\chi(d)=\sum_{\begin{subarray}{c}d\mid k\\ d\leqslant\sqrt{x}\end{subarray}}\chi(d)+\chi(k)\sum_{\begin{subarray}{c}e\mid k% \\ e<k/\sqrt{x}\end{subarray}}\chi(e).

This shows that

{S_{1}(x;q,a)=4S_{1}^{-}(x;q,a)+4S_{1}^{+}(x;q,a),}

(3.6)

where

S_{1}^{-}(x;q,a):=\sum_{d\leqslant\sqrt{x}}\chi(d)\sum_{\begin{subarray}{c}m% \leqslant x,m\equiv 1\left(\textnormal{mod}\ 4\right)\\ m\equiv a\left(\textnormal{mod}\ q\right),d\mid m\end{subarray}}r(m+1)

and

S_{1}^{+}(x;q,a):=\sum_{e\leqslant\sqrt{x}}\chi(e)\sum_{\begin{subarray}{c}m% \in(e\sqrt{x},x],m\equiv 1\left(\textnormal{mod}\ 4\right)\\ m\equiv a\left(\textnormal{mod}\ q\right),e\mid m\end{subarray}}r(m+1).

Since $4\mid q$ by (3.3) we infer that the three congruences in the sums $S_{1}^{-},S_{1}^{+}$ are soluble if and only if $4\mid a-1$ and $\gcd(d,q)\mid a.$ Denoting the least common multiple by $[\cdot,\cdot]$ we see that under the aforementioned conditions there exists a unique $t\in\mathbb{Z}/[d,q]\mathbb{Z}$ such that

m\equiv 1\left(\textnormal{mod}\ 4\right),m\equiv a\left(\textnormal{mod}\ q% \right),m\equiv 0\left(\textnormal{mod}\ d\right)\iff m\equiv t\left(% \textnormal{mod}\ [d,q]\right).

(3.7)

Hence,

S_{1}^{-}(x;q,a)=\mathds{1}(a\equiv 1\left(\textnormal{mod}\ 4\right))\sum_{% \begin{subarray}{c}d\leqslant\sqrt{x}\\ \gcd(d,q)\mid a\end{subarray}}\chi(d)\sum_{\begin{subarray}{c}m\leqslant x\\ m\equiv t\left(\textnormal{mod}\ [d,q]\right)\end{subarray}}r(m+1).

Lemma 3.1.

For all $a,q$ satisfying (3.3) and all $x\geqslant 1$ we have

S_{1}^{-}(x;q,a)=\mathds{1}(a\equiv 1\left(\textnormal{mod}\ 4\right))\pi x% \mathscr{M}++O\left(q\tau(q)^{4}x^{5/6}(\log x)^{19}\right),

where the implied constant is absolute and

\mathscr{M}:=\sum_{\begin{subarray}{c}d\leqslant\sqrt{x}\\ \gcd(d,q)\mid a\end{subarray}}\chi(d)\frac{\sharp\left\{\mathbf{y}\in(\mathbb{% Z}/[d,q]\mathbb{Z})^{2}:y_{1}^{2}+y_{2}^{2}\equiv a+1\left(\textnormal{mod}\ q% \right),y_{1}^{2}+y_{2}^{2}\equiv 1\left(\textnormal{mod}\ d\right)\right\}}{[% d,q]^{2}}.

Proof.

Applying [42, Theorem, page 262] gives the error term

\ll\sum_{\begin{subarray}{c}d\leqslant\sqrt{x}\end{subarray}}\left(\left([d,q]% ^{1/2}+x^{1/3}\right)\gcd(a,[d,q])^{1/2}\tau([d,q])^{4}(\log x)^{4}\right).

By (3.3) we have $\gcd(a,[d,q])\leqslant q$ . Indeed, if a prime $p\mid a,d$ but $p\nmid q$ , since $d\mid m$ , we deduce that $p\mid q$ which is a contradiction. Hence, the error term is

\ll q^{1/2}\tau(q)^{4}(\log x)^{4}\sum_{\begin{subarray}{c}d\leqslant\sqrt{x}% \end{subarray}}\left((dq)^{1/2}+x^{1/3}\right)\tau(d)^{4}\ll q\tau(q)^{4}x^{5/% 6}(\log x)^{19}.

The main term supplied by [42, Theorem, page 262] equals

\mathds{1}(a\equiv 1\left(\textnormal{mod}\ 4\right))\pi x\sum_{\begin{% subarray}{c}d\leqslant\sqrt{x}\\ \gcd(d,q)\mid a\end{subarray}}\chi(d)\frac{\eta_{t+1}([d,q])}{[d,q]^{2}},

which takes the required shape in light of (3.7) and the fact that the congruence modulo 4 is implied by the one modulo $q$ provided that $a\equiv 1\left(\textnormal{mod}\ 4\right).$ ∎

One can similarly prove the estimate

{S_{1}^{+}(x;q,a)=\mathds{1}(a\equiv 1\left(\textnormal{mod}\ 4\right))\pi% \left(x\mathscr{M}-E\sqrt{x}\right)+O\left(q\tau(q)^{4}x^{5/6}(\log x)^{19}% \right),}

(3.8)

where the implied constant is absolute and

E:=\sum_{\begin{subarray}{c}e\leqslant\sqrt{x}\\ \gcd(e,q)\mid a\end{subarray}}\chi(e)e\frac{\sharp\left\{\mathbf{y}\in(\mathbb% {Z}/[e,q]\mathbb{Z})^{2}:y_{1}^{2}+y_{2}^{2}\equiv a+1\left(\textnormal{mod}\ % q\right),y_{1}^{2}+y_{2}^{2}\equiv 1\left(\textnormal{mod}\ d\right)\right\}}{% [e,q]^{2}}.

Taking into account the oscillation of $\chi(e)$ , we next show that $E$ makes a negligible contribution.

Lemma 3.2.

For all $a,q$ satisfying (3.3) and all $x\geqslant 1$ we have

E\ll q^{3/2}(\log q)\tau(q)(\log x),

where the implied constant is absolute.

Proof.

Write $e=e_{0}e_{1}$ where $e_{1}$ is coprime to $q$ and each prime factor of $e_{0}$ divides $q$ . Then the condition $\gcd(e,q)\mid a$ becomes $\gcd(e_{0},q)\mid a$ and by (3.3) we infer that $e_{0}\mid q$ . Thus,

E=\sum_{\begin{subarray}{c}e_{0}\leqslant\sqrt{x}\\ e_{0}\mid\gcd(a,q)\end{subarray}}\chi(e_{0})e_{0}\frac{\sharp\left\{\mathbf{y}% \in(\mathbb{Z}/q\mathbb{Z})^{2}:y_{1}^{2}+y_{2}^{2}\equiv a+1\left(\textnormal% {mod}\ q\right)\right\}}{q^{2}}E_{1}(\sqrt{x}/e_{0}),

where

E_{1}(T):=\sum_{\begin{subarray}{c}e_{1}\leqslant T\\ \gcd(e_{1},q)=1\end{subarray}}\chi(e_{1})\frac{\sharp\left\{\mathbf{y}\in(% \mathbb{Z}/e_{1}\mathbb{Z})^{2}:y_{1}^{2}+y_{2}^{2}\equiv 1\left(\textnormal{% mod}\ e_{1}\right)\right\}}{e_{1}}.

By [25, pages 27–28] we infer that $E_{1}(T)$ equals

\sum_{\begin{subarray}{c}e_{1}\leqslant T\\ \gcd(e_{1},q)=1\end{subarray}}\chi(e_{1})\prod_{p\mid e_{1}}\left(1-\frac{\chi% (p)}{p}\right)=\sum_{\begin{subarray}{c}k\leqslant T\\ \gcd(k,q)=1\end{subarray}}\frac{\mu(k)\chi(k)}{k}\sum_{\begin{subarray}{c}e_{1% }\leqslant T,k\mid e_{1}\\ \gcd(e_{1},q)=1\end{subarray}}\chi(e_{1})=\sum_{\begin{subarray}{c}k\leqslant T% \\ \gcd(k,q)=1\end{subarray}}\frac{\mu(k)}{k}\sum_{\begin{subarray}{c}t\leqslant T% /k\\ \gcd(t,q)=1\end{subarray}}\chi(t).

We may view $\mathds{1}(\gcd(t,q)=1)\chi(t)$ as a character modulo $q$ (remember here that $4\mid q$ ) in $t$ . Hence, the Pólya–Vinogradov estimate gives the bound $E_{1}(T)\ll\sqrt{q}(\log q)(\log T)$ . Hence,

E\ll\sqrt{q}(\log q)(\log x)\sum_{e_{0}\mid\gcd(a,q)}e_{0}\ll q^{3/2}(\log q)(% \log x)\tau(q)

since $\gcd(a,q)\mid q$ . ∎

Lemma 3.3.

Assume that $q\in\mathbb{N}$ is divisible by $4$ , that $a\equiv 1\left(\textnormal{mod}\ 4\right)$ and that $a\in\mathbb{Z}/q\mathbb{Z}$ is such that $v_{p}(a)<v_{p}(q)$ for all $p\mid q$ . Then

\sum_{d_{0}\mid\gcd(a,q)}\chi(d_{0})=\frac{1}{2}\frac{\eta_{q}(a)}{q}\prod_{p% \mid q}\left(1-\frac{\chi(p)}{p}\right)^{-1},

where $\eta_{q}(a)$ is defined in (3.4).

Proof.

The sum over $d_{0}$ can be written as

\prod_{\begin{subarray}{c}p\equiv 1\left(\textnormal{mod}\ 4\right)\\ p\mid\gcd(a,q)\end{subarray}}(1+v_{p}(a))\prod_{\begin{subarray}{c}p\equiv 3% \left(\textnormal{mod}\ 4\right)\\ p\mid\gcd(a,q)\end{subarray}}\frac{1+(-1)^{v_{p}(a)}}{2}.

Now we can replace the condition $p\mid\gcd(a,q)$ by $p\mid q$ in both products because if $v_{p}(a)=0$ then the analogous terms in the products equal $1$ . By [25, pages 27–28] and the fact that $v_{p}(a)<v_{p}(q)$ we then obtain

\prod_{\begin{subarray}{c}p\equiv 1\left(\textnormal{mod}\ 4\right)\\ p\mid q\end{subarray}}\frac{\eta_{p^{v_{p}(q)}}(a)}{p^{v_{p}(q)}(1-1/p)}\prod_% {\begin{subarray}{c}p\equiv 3\left(\textnormal{mod}\ 4\right)\\ p\mid q\end{subarray}}\frac{\eta_{p^{v_{p}(q)}}(a)}{p^{v_{p}(q)}(1+1/p)}.

Lastly, by [25, Equation (8.4)] and our assumptions $4\mid q$ and $a\equiv 1\left(\textnormal{mod}\ 4\right)$ this becomes

\frac{1}{2}\prod_{\begin{subarray}{c}p\mid q\end{subarray}}\frac{\eta_{p^{v_{p% }(q)}}(a)}{p^{v_{p}(q)}(1-\chi(p)/p)},

which completes the proof. ∎

Lemma 3.4.

For all $a,q$ satisfying (3.3) and all $q\geqslant x^{1/2},$ we have

\mathscr{M}=\frac{\pi}{8}\frac{\eta_{q}(a)\eta_{q}(a+1)}{q}\prod_{p\nmid q}% \left(1-\frac{1}{p^{2}}\right)+O\left(\tau(q)q^{3/2}\frac{(\log x)^{2}}{\sqrt{% x}}\right),

where the implied constant is absolute.

Proof.

Writing $d=d_{0}d_{1}$ as in the proof of Lemma 3.2 we see that $\mathscr{M}$ equals

\frac{\sharp\left\{\mathbf{y}\in(\mathbb{Z}/q\mathbb{Z})^{2}:y_{1}^{2}+y_{2}^{% 2}\equiv a+1\left(\textnormal{mod}\ q\right)\right\}}{q^{2}}\sum_{\begin{% subarray}{c}d_{0}\mid\gcd(a,q)\end{subarray}}\chi(d_{0})\sum_{\begin{subarray}% {c}d_{1}\leqslant\sqrt{x}/d_{0}\\ \gcd(d_{1},q)=1\end{subarray}}\frac{\chi(d_{1})}{d_{1}}\prod_{p\mid d_{1}}% \left(1-\frac{\chi(p)}{p}\right).

We used [25, pages 27–28] and the fact that $y_{1}^{2}+y_{2}^{2}\equiv 1\left(\textnormal{mod}\ d_{0}\right)$ is implied by $y_{1}^{2}+y_{2}^{2}\equiv a+1\left(\textnormal{mod}\ q\right)$ due to $d_{0}\mid\gcd(a,q)$ . Note that the condition $d_{0}\leqslant\sqrt{x}$ is implied by $d_{0}\mid q$ owing to the assumption $q\geqslant x^{1/2}$ . The sum over $d_{1}$ equals

\sum_{\begin{subarray}{c}k\leqslant\sqrt{x}/d_{0}\\ \gcd(k,q)=1\end{subarray}}\frac{\mu(k)}{k^{2}}\sum_{\begin{subarray}{c}t% \leqslant\sqrt{x}/(kd_{0})\\ \gcd(t,q)=1\end{subarray}}\frac{\chi(t)}{t}.

The sum over $t$ can be seen as converging to the value of an $L$ -function at $1$ with Dirichlet character $t\mapsto\chi(t)\mathds{1}(\gcd(t,q)=1)$ . Hence, by [37, Lemma 16] for example, it equals

\sum_{\begin{subarray}{c}t=1\\ \gcd(t,q)=1\end{subarray}}^{\infty}\frac{\chi(t)}{t}+O\left(q^{1/2}\log(q)% \frac{kd_{0}}{\sqrt{x}}\right),

with an absolute implied constant. Thus, the main term becomes

	$\displaystyle\sum_{\begin{subarray}{c}k\leqslant\sqrt{x}/d_{0}\\ \gcd(k,q)=1\end{subarray}}\frac{\mu(k)}{k^{2}}\left(\sum_{\begin{subarray}{c}t% =1\\ \gcd(t,q)=1\end{subarray}}^{\infty}\frac{\chi(t)}{t}+O\left(q^{1/2}\log(q)% \frac{kd_{0}}{\sqrt{x}}\right)\right)=$	$\displaystyle\prod_{p\nmid q}\left(1-\frac{1}{p^{2}}\right)\sum_{\begin{% subarray}{c}t=1\\ \gcd(t,q)=1\end{subarray}}^{\infty}\frac{\chi(t)}{t}$
		$\displaystyle+O\left(q^{1/2}\log(q)\frac{d_{0}}{\sqrt{x}}\log x\right),$

where we used the standard bound $L(1,\psi)\ll\log r$ for a non-principal character $\psi$ modulo $r$ . Finally, we use Lemma 3.3 to deal with the sum over $d_{0}$ and the calculation

\prod_{p\mid q}\left(1-\frac{\chi(p)}{p}\right)^{-1}\sum_{\begin{subarray}{c}t% =1\\ \gcd(t,q)=1\end{subarray}}^{\infty}\frac{\chi(t)}{t}=\sum_{t=1}^{\infty}\frac{% \chi(t)}{t}=\frac{\pi}{4}.

∎

3.5. The proof of Theorem 3.1

Injecting Lemma 3.2 into (3.8) gives an asymptotic for $S_{1}^{+}$ in terms of $\mathscr{M}$ . By Lemma 3.1 we have a similar asymptotic for $S_{1}^{-}$ . Putting these estimates into (3.6) and then using Lemma 3.4 to estimate $\mathscr{M}$ yields an asymptotic for $S_{1}$ . A similar argument works in an identical manner for $S_{0}$ and that concludes the proof of Theorem 3.1.

3.6. The proof of Theorem 1.3

We shall need the following lemma.

Lemma 3.5.

For every prime $p$ , any integer $m\geqslant 1$ and any $x\geqslant 1$ with $p^{m}\leqslant x^{1/4}$ we have

\sum_{\begin{subarray}{c}|n|\leqslant x\\ p^{m}\mid n\end{subarray}}r(n)r(n+1)\ll\frac{m}{p^{m}}x,

where the implied constant is absolute.

Proof.

We have

\sum_{\begin{subarray}{c}|n|\leqslant x\\ p^{m}\mid n\end{subarray}}r(n)r(n+1)\ll\sum_{\begin{subarray}{c}|n|\leqslant x% \\ p^{m}\mid n\end{subarray}}r(n(n+1))

where we used the coprimality of $n,1+n$ to bound $r(n)r(1+n)\ll r(n(1+n))$ . Note that we must have $n\equiv 0,1\left(\textnormal{mod}\ 4\right)$ in order to have $r(n)r(n+1)\neq 0.$ The result would follow from the work of Nair [30], however, its application is prohibited by the fact that the polynomial $n(n+1)$ has $2$ as a fixed divisor. It suffices to work with the multiplicative function $r^{\prime}=r/4$ . Then the sum in the lemma is

\ll\sum_{\begin{subarray}{c}|t|\leqslant x/p^{m}\end{subarray}}r^{\prime}(p^{m% }t)r^{\prime}(p^{m}t+1)\leqslant(m+1)\sum_{\begin{subarray}{c}|t|\leqslant x/p% ^{m}\end{subarray}}r^{\prime}(t(p^{m}t+1)),

since $r^{\prime}(ab)\leqslant\tau(a)r^{\prime}(b)$ holds for all integers $a,b$ . Using the bound $r^{\prime}(s)\leqslant\tau(s)\ll s^{1/10}$ we see that the contribution of the terms for which $2^{k}\mid t$ for some $k$ satisfying $2^{k}>x^{1/4}$ is

\ll(m+1)\sum_{\begin{subarray}{c}|t|\leqslant x/p^{m}\\ 2^{k}\mid t\end{subarray}}x^{1/5}\ll mx^{1/5}\left(\frac{x}{p^{m}2^{k}}+1% \right)\ll mx^{19/20}p^{-m}+mx^{1/5}.

The assumption $p^{m}\leqslant x^{1/4}$ implies that $x^{1/5}\leqslant xp^{-m}$ , hence, the bound is satisfactory.

When $k\geqslant 1$ and $2^{k}\leqslant x^{1/4}$ we can write, in the case $t=2^{k}t^{\prime}$ with $t^{\prime}$ odd to obtain

\ll m\sum_{1\leqslant k\leqslant(\log x)/(4\log 2)}k\sum_{\begin{subarray}{c}|% t^{\prime}|\leqslant x/(p^{m}2^{k})\\ 2\nmid t^{\prime}\end{subarray}}r^{\prime}(t^{\prime}(p^{m}2^{k}t^{\prime}+1))

where we used $r^{\prime}(2^{k}s)\ll kr^{\prime}(s)$ . If $t^{\prime}\equiv 3\left(\textnormal{mod}\ 4\right)$ then $r^{\prime}(t^{\prime})=0$ , hence, we may assume that $t^{\prime}\equiv 1\left(\textnormal{mod}\ 4\right)$ . Then we can write $t^{\prime}=4s+1$ , thus, we get the bound

\ll m\sum_{0\leqslant k\leqslant(\log x)/(4\log 2)}k\sum_{\begin{subarray}{c}|% s^{\prime}|\leqslant 4x/(p^{m}2^{k})\end{subarray}}r^{\prime}(P(s)),

where $P(s)=(4s+1)(p^{m}2^{k}(4s+1)+1)$ . Since $P(0)$ is odd, we can apply Nair’s result [30, Theorem, page 259] with $\delta=1$ . Indeed, we have $\|P\|\ll p^{m}2^{k}\ll(xp^{-m}2^{-k})$ because both $2^{k}$ and $p^{m}$ are at most $x^{1/4}$ . We obtain

\sum_{\begin{subarray}{c}|s^{\prime}|\leqslant 4x/(p^{m}2^{k})\end{subarray}}r% ^{\prime}(P(s))\ll\frac{x}{p^{m}2^{k}}\prod_{2<p\leqslant x}\bigg{(}1-\frac{2}% {p}\bigg{)}\exp\bigg{(}\sum_{p\leqslant x}\frac{2r^{\prime}(p)}{p}\bigg{)}\ll% \frac{x}{p^{m}2^{k}}

with an absolute implied constant. Thus, the overall contribution is

\ll x\frac{m}{p^{m}}\sum_{0\leqslant k\leqslant(\log x)/(4\log 2)}\frac{k}{2^{% k}},

which is sufficient. The contribution of the cases with $k=0$ can be dealt with in an analogous way by making substitutions in the term $r^{\prime}(n+1)$ . ∎

The proof of Theorem 1.3 is an application of Theorem 1.1 with

m_{p}(z)=1+\left[\frac{\log z}{\log p}\right],\quad k(m)=r(m)r(m+1)\mathds{1}(% p\leqslant z\Rightarrow v_{p}(m)<m_{p}(z)),

together with $s=1$ , $\omega(t)=\pi^{2}$ and

\rho(a,q)=\frac{\eta_{q}(a)\eta_{q}(a+1)}{q^{3}}\mathds{1}(a\equiv 0,1\left(% \textnormal{mod}\ 4\right))\prod_{p\nmid q}\left(1-\frac{1}{p^{2}}\right)

in (1.1). We have $W_{z}\leqslant\mathrm{e}^{3z}$ since

\log W_{z}=\sum_{p\leqslant z}\left(1+\left[\frac{\log z}{\log p}\right]\right% )\leqslant 3z

for all large enough $z$ by the Prime Number Theorem. Hence, by Theorem 3.1 we have

\frac{E(bP^{d},W_{z})W_{z}}{P^{d}}\ll\frac{(\log P)^{19}W_{z}^{5/2}\tau(W_{z})% ^{4}}{P^{d/6}}\ll\frac{(\log P)^{19}W_{z}^{3}}{P^{d/6}}\ll\frac{(\log P)^{19}}% {P^{d/12}}

as long as we assume $z=z(P):=\frac{d}{108}\log P$ and using the $q=1$ case. Thus, the error term in Theorem 1.1 is

\ll P^{-\delta}+\widetilde{\varepsilon}(z)+z^{-c}+\frac{(\log P)^{19}}{P^{d/12% }}\to 0\textrm{ as }P\to+\infty

since one can show that $\widetilde{\varepsilon}(z)\to 0$ in a similar manner as in §2.1. For the main term we note that

\lim_{P\to\infty}\int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}:sf(% \mathbf{t})>P^{-d}\end{subarray}}\omega(P^{d}sf(\mathbf{t}))\mathrm{d}\mathbf{% t}=\pi^{2}\mathrm{vol}\left(\{\mathbf{t}\in\mathscr{B}\,:\,f(\mathbf{t})>0\}% \right).

Since

\prod_{p\nmid T_{z}}\left(1-\frac{1}{p^{2}}\right)=\prod_{p>z}\left(1-\frac{1}% {p^{2}}\right)=1+O\left(\frac{1}{z}\right),

we see that

\sum_{\mathbf{t}\in(\mathbb{Z}/T_{z}\mathbb{Z})^{n}}\rho(f(\mathbf{t}),T_{z})=% \frac{1+O(1/z)}{{T_{z}}^{3}}\sum_{\begin{subarray}{c}\mathbf{t}\in(\mathbb{Z}/% T_{z}\mathbb{Z})^{n}\\ f(\mathbf{t})\equiv 0,1\left(\textnormal{mod}\ 4\right)\end{subarray}}\eta_{T_% {z}}(f(\mathbf{t}))\eta_{T_{z}}(f(\mathbf{t})+1).

The condition $f(\mathbf{t})\equiv 0,1\left(\textnormal{mod}\ 4\right)$ is implied by the fact that $4$ divides $T_{z}$ and the presence of the terms $\eta_{T_{z}}(f(\mathbf{t}))\eta_{T_{z}}(f(\mathbf{t})+1)$ . Hence,

		$\displaystyle\lim_{z\to\infty}\frac{\sharp\{(\mathbf{t},\mathbf{x},\mathbf{y})% \in(\mathbb{Z}/T_{z}\mathbb{Z})^{n+4}:x_{1}^{2}+x_{2}^{2}=f(\mathbf{t}),y_{1}^% {2}+y_{2}^{2}=1+f(\mathbf{t}),v_{p}(f(\mathbf{t}))<m_{p}(z)\,\,\forall p% \leqslant z\}}{T_{z}^{n+2}}$
	$\displaystyle=$	$\displaystyle\lim_{z\to\infty}T_{z}^{-n+1}\sum_{\begin{subarray}{c}\mathbf{t}% \in(\mathbb{Z}/T_{z}\mathbb{Z})^{n}\\ v_{p}(f(\mathbf{t}))<m_{p}(z)\,\,\forall p\leqslant z\end{subarray}}\rho(f(% \mathbf{t}),T_{z}),$

which can be shown to exist as in §2.1 by using the fact that $\rho(f(\mathbf{t}),T_{z})\geqslant 0$ . The condition $v_{p}(f(\mathbf{t}))<m_{p}(z)$ can be ignored since the number of $\mathbf{t}\in(\mathbb{Z}/p^{m_{p}(z)}\mathbb{Z})^{n}$ with $f(\mathbf{t})=0$ equals $p^{m_{p}(z)(n-1)}(1+O(p^{-1-c}))$ for some $c>0$ by (2.4) and (2.6) and the number of $\mathbf{x},\mathbf{y}\in\left(\mathbb{Z}/p^{m_{p}(z)\mathbb{Z}}\right)^{4}$ with $x_{1}^{2}+x_{2}^{2}=0,y_{1}^{2}+y_{2}^{2}=1$ is $O\left(m_{p}(z)p^{2m_{p}(z)}\right)$ by [25, pages 27–28]. Hence,

		$\displaystyle\prod_{p\leqslant z}\frac{\sharp\{(\mathbf{t},\mathbf{x},\mathbf{% y})\in(\mathbb{Z}/p^{m_{p}(z)}\mathbb{Z})^{n+4}:x_{1}^{2}+x_{2}^{2}=f(\mathbf{% t}),y_{1}^{2}+y_{2}^{2}=1+f(\mathbf{t}),v_{p}(f(\mathbf{t}))<m_{p}(z)\forall p% \leqslant z\}}{p^{(n+2)m_{p}(z)}}$
	$\displaystyle=$	$\displaystyle V(z)\prod_{p\leqslant z}\frac{\sharp\{(\mathbf{t},\mathbf{x},% \mathbf{y})\in(\mathbb{Z}/p^{m_{p}(z)}\mathbb{Z})^{n+4}:x_{1}^{2}+x_{2}^{2}=f(% \mathbf{t}),y_{1}^{2}+y_{2}^{2}=1+f(\mathbf{t})\}}{p^{(n+2)m_{p}(z)}},$

where $V(z)=\displaystyle\prod_{p\leqslant z}\left(1+O\left(m_{p}(z)p^{-m_{p}(z)}% \right)\right)$ . We have $\lim_{z\to\infty}V(z)=1$ since $\log V(z)$ is

\ll\sum_{p\leqslant z}\frac{m_{p}(z)}{p^{m_{p}(z)}}\leqslant(\log z)\sum_{p% \leqslant z}p^{-1-[\log z/\log p]},

which was shown to converge to $0$ in §2.1.

We have so far shown that

\lim_{P\to\infty}P^{-n}\sum_{\begin{subarray}{c}\mathbf{x}\in\mathbb{Z}^{n}% \cap P\mathscr{B}\\ v_{p}(f(\mathbf{x}))<m_{p}(z)\,\,\forall p\leqslant z\end{subarray}}r(f(% \mathbf{x}))r(1+f(\mathbf{x}))=\pi^{2}\mathrm{vol}\left(\{\mathbf{t}\in% \mathscr{B}\,:\,f(\mathbf{t})>0\}\right)\prod_{p}\sigma_{p}(f),

where $\sigma_{p}(f)$ are as in Theorem 1.3 and it remains to get rid of the condition $v_{p}(f(\mathbf{x}))<m_{p}(z)$ .

Letting $\mathscr{A}$ be the set of integers $\nu$ for which there exists $p\leqslant z$ with $p^{m_{p}(z)}\mid\nu$ we have

\sum_{\begin{subarray}{c}\mathbf{x}\in\mathbb{Z}^{n}\cap P\mathscr{B}\\ f(\mathbf{x})\in\mathscr{A}\end{subarray}}r(f(\mathbf{x}))r(1+f(\mathbf{x}))% \ll\sum_{\begin{subarray}{c}|\nu|\ll P^{d}\\ \nu\in\mathscr{A}\end{subarray}}r(\nu)r(1+\nu)\sum_{\begin{subarray}{c}\mathbf% {x}\in\mathbb{Z}^{n}\cap P\mathscr{B}\\ f(\mathbf{x})=\nu\end{subarray}}1.,

By [1, Lem.5.5] we can bound this by

\ll P^{n-d}\sum_{\begin{subarray}{c}|\nu|\ll P^{d}\\ \nu\in\mathscr{A}\end{subarray}}r(\nu)r(1+\nu)\mathfrak{S}(\nu)J(\nu P^{-d})% \ll P^{n-d}\sum_{\begin{subarray}{c}|\nu|\ll P^{d}\\ \nu\in\mathscr{A}\end{subarray}}r(\nu)r(1+\nu)

since the estimates $\mathfrak{S}(\nu),J(\nu P^{-d})=O(1)$ hold uniformly in $\nu$ by (2.3) and (2.8). Since $p^{m_{p}(z)}\leqslant z\leqslant P^{d/4}$ we can employ Lemma 3.5 to obtain

\sum_{\begin{subarray}{c}|\nu|\ll P^{d}\\ \nu\in\mathscr{A}\end{subarray}}r(\nu)r((1+\nu)\ll P^{d}\sum_{p\leqslant z}% \frac{m_{p}(z)}{p^{m_{p}(z)}}\leqslant P^{d}\frac{\log z}{\log 2}\widetilde{% \varepsilon}(z).

By the estimate $\widetilde{\varepsilon}(z)\ll(\log z)/\sqrt{z}$ proved in §2.1 we get $\ll P^{d}(\log z)^{2}/\sqrt{z}$ , which concludes the proof.

3.7. The proof of Theorem 1.4

We shall use Theorem 1.2 with $k$ being the indicator function of elements of $\mathscr{A}$ . By our assumption for $q=1$ we see that $\displaystyle\int_{1}^{T}\omega(t)\mathrm{d}t$ is asymptotic to $\sharp\{\mathscr{A}\cap[1,T]\}$ . By assumption, $\mathscr{A}$ is non-empty, hence $\rho(1,1)\displaystyle\int_{1}^{T}\omega(t)\mathrm{d}t$ is non-zero for all large $T$ . This verified the second assumption of Theorem 1.2. Note that $\rho(r,q)$ is non-negative, it being the limit of non-negative counting functions. The remaining assumptions are easy to verify.

3.8. The proof of Theorem 1.9

We apply Theorem 1.2 with $s=1$ and $\mathscr{A}_{k}$ being the set of $k$ -th powers of positive integers. For all $r,q$ and $x$ , one has

\sum_{\begin{subarray}{c}1\leqslant m\leqslant x\\ m\equiv r\left(\textnormal{mod}\ q\right)\end{subarray}}\mathds{1}_{\mathscr{A% }_{k}}(m)=\sum_{\begin{subarray}{c}1\leqslant t\leqslant x^{1/k}\\ t^{k}\equiv r\left(\textnormal{mod}\ q\right)\end{subarray}}1=\sum_{\begin{% subarray}{c}y\in\mathbb{Z}/q\mathbb{Z}\\ y^{k}\equiv r\left(\textnormal{mod}\ q\right)\end{subarray}}\left(\frac{x^{1/k% }}{q}+O(1)\right)=\rho(r,q)\int_{1}^{x}\frac{\mathrm{d}t}{kt^{1-1/k}}+O(1),

where $q\rho(r,q)=\sharp\{y\left(\textnormal{mod}\ q\right):y^{k}=r\}$ and the implied constant is independent of $x$ . Recall that we are applying Theorem 1.2 and hence only need to work with fixed modulus $q$ and don’t need to record the dependence on $q$ in the error terms. Thus, letting $\omega(t)=1/(kt^{1-1/k})$ we find that

\frac{\sharp\{a\in\mathscr{A}_{k}:a\leqslant x,a\equiv r\left(\textnormal{mod}% \ q\right)\}}{\displaystyle\int_{1}^{x}\omega(t)\mathrm{d}t}-\rho(r,q)\ll\frac% {1}{\displaystyle\int_{1}^{x}\omega(t)\mathrm{d}t}\to 0.

To verify the remaining assumption of Theorem 1.2, we use the assumption that $f$ assumes at least one strictly positive value in $\mathscr{B}$ so that for all large $P$ we have

{P^{d}\int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}\\ f(\mathbf{t})>P^{-d}\end{subarray}}\frac{\mathrm{d}\mathbf{t}}{k(P^{d}f(% \mathbf{t}))^{1-1/k}}=\frac{P^{d/k}}{k}\int\limits_{\begin{subarray}{c}\mathbf% {t}\in\mathscr{B}\\ f(\mathbf{t})>P^{-d}\end{subarray}}\frac{\mathrm{d}\mathbf{t}}{f(\mathbf{t})^{% 1-1/k}}\gg P^{d/k}\gg\int_{1}^{bP^{d}}\frac{\mathrm{d}t}{t^{1-1/k}}.}

(3.9)

Invoking Theorem 1.2 allows us to conclude the proof.

3.9. The proof of Theorem 1.10

A positive integer $k$ is $m$ -full equivalently when it has the shape $k=k_{1}^{m}k_{2}^{m+1}\cdots k_{m}^{2m-1}$ where $k_{2}\cdots k_{m}$ is square-free. Thus, for any $x,r$ and $q$ , we have

\sharp\{k\in\mathbb{N}\cap[1,x]:k\textrm{ is }m\textrm{-full },k\equiv r\left(% \textnormal{mod}\ q\right)\}=\sum_{\begin{subarray}{c}k_{1}^{m}k_{2}^{m+1}% \cdots k_{m}^{2m-1}\leqslant x\\ k_{1}^{m}k_{2}^{m+1}\cdots k_{m}^{2m-1}\equiv r\left(\textnormal{mod}\ q\right% )\end{subarray}}\mu(k_{2}\cdots k_{m})^{2}.

Fixing the values of $k_{2},\ldots,k_{m}$ , this becomes

\sum_{\begin{subarray}{c}k_{2},\ldots,k_{m}\in\mathbb{N}\\ k_{2}^{m+1}\cdots k_{m}^{2m-1}\leqslant x\end{subarray}}\mu(k_{2}\cdots k_{m})% ^{2}\operatorname*{\sum{}^{*}}_{\begin{subarray}{c}y\left(\textnormal{mod}\ q% \right)\end{subarray}}\sharp\left\{k_{1}\leqslant(xk_{2}^{-m-1}\cdots k_{m}^{-% 2m+1})^{1/m}:k_{1}\equiv y\left(\textnormal{mod}\ q\right)\right\},

where $\operatorname*{\sum{}^{*}}$ is subject to $y^{m}k_{2}^{m+1}\cdots k_{m}^{2m-1}\equiv r\left(\textnormal{mod}\ q\right)$ . We obtain

\sum_{\begin{subarray}{c}k_{2},\ldots,k_{m}\in\mathbb{N}\\ k_{2}^{m+1}\cdots k_{m}^{2m-1}\leqslant x\end{subarray}}\mu(k_{2}\cdots k_{m})% ^{2}\sum_{\begin{subarray}{c}y\left(\textnormal{mod}\ q\right)\end{subarray}}% \left(\frac{(xk_{2}^{-m-1}\cdots k_{m}^{-2m+1})^{1/m}}{q}+O(1)\right).

For any $\frac{1}{m(m+1)}>\varepsilon>0$ the error term contributes

\ll\sharp\left\{k_{2},\ldots,k_{m}\in\mathbb{N}:(k_{2}\cdots k_{m})^{m+1}% \leqslant x\right\}\leqslant\sum_{a\leqslant x^{1/(m+1)}}\tau_{m-1}(a)\ll_{% \varepsilon}x^{\frac{1}{m+1}+\varepsilon}=o(x^{1/m})

where $\tau_{m-1}(a)$ is the number of ways of writing $a$ as the product of $m-1$ positive integers and we used the standard bound $\tau_{m-1}(a)\ll_{\varepsilon}a^{\varepsilon}$ . Recall that we are applying Theorem 1.2 and hence only need to work with fixed modulus $q$ and don’t need to record the dependence on $q$ in the error terms. The main term equals

\frac{x^{1/m}}{q}\sum_{\begin{subarray}{c}k_{2}\ldots,k_{m}\in\mathbb{N}\\ k_{2}^{m+1}\cdots k_{m}^{2m-1}\leqslant x\end{subarray}}\mu(k_{2}\cdots k_{m})% ^{2}(k_{2}^{-m-1}\cdots k_{m}^{-2m+1})^{1/m}\operatorname*{\sum{}^{*}}_{\begin% {subarray}{c}y\left(\textnormal{mod}\ q\right)\end{subarray}}1.

By the trivial bound $\displaystyle\operatorname*{\sum{}^{*}}_{\begin{subarray}{c}y\left(\textnormal% {mod}\ q\right)\end{subarray}}1\leqslant q$ we see that the sum is convergent as $x\to\infty$ , hence we have obtained

\sharp\{k\in\mathbb{N}\cap[1,x]:k\textrm{ is }m\textrm{-full },k\equiv r\left(% \textnormal{mod}\ q\right)\}=\rho(r,q)x^{1/m}+o(x^{1/m}),

where

\rho(r,q)=\frac{1}{q}\sum_{\begin{subarray}{c}k_{2},\ldots,k_{m}\in\mathbb{N}% \end{subarray}}\frac{\mu(k_{2}\cdots k_{m})^{2}}{k_{2}^{1+1/m}k_{3}^{1+2/m}% \cdots k_{m}^{2-1/m}}\sharp\left\{y\in\mathbb{Z}/q\mathbb{Z}:y^{m}k_{2}^{m+1}% \cdots k_{m}^{2m-1}=r\right\}.

Letting $\omega(t)=\frac{1}{\displaystyle mt^{1/m-1}}$ we can then verify the remaining assumption of Theorem 1.2 as in (3.9).

3.10. The proof of Theorem 1.5

Let $k$ be the indicator function of the integers that are an integer power of $2$ . Then for all $q\in\mathbb{N}$ and $r\in[1,q]$ we have

\sum_{\begin{subarray}{c}m\leqslant x\\ m\equiv r\left(\textnormal{mod}\ q\right)\end{subarray}}k(m)=\sharp\left\{0% \leqslant t\leqslant\frac{\log x}{\log 2}:2^{t}\equiv r\left(\textnormal{mod}% \ q\right)\right\}.

Assume that $v_{2}(r)<v_{2}(q)$ . Then the congruence is equivalent to $2^{t-v_{2}(r)}\equiv r2^{-v_{2}(r)}\left(\textnormal{mod}\ q2^{-v_{2}(r)}\right)$ . If there exists $t>v_{2}(r)$ with such a property then since both $2^{t-v_{2}(r)},q2^{-v_{2}(r)}$ are even, one sees that $r2^{-v_{2}(r)}$ is even, which is a contradiction. Hence, the sum is

O(v_{2}(q))+\sharp\left\{v_{2}(q)\leqslant t\leqslant\frac{\log x}{\log 2}:2^{% t-v_{2}(q)}\equiv r2^{-v_{2}(q)}\left(\textnormal{mod}\ q2^{-v_{2}(q)}\right)% \right\}\mathds{1}(v_{2}(r)\geqslant v_{2}(q)).

We can now define

\langle 2\rangle=\left\{2^{j}\left(\textnormal{mod}\ q2^{-v_{2}(q)}\right):j% \in\mathbb{N}\right\}.

Then we must have $r2^{-v_{2}(q)}\in\langle 2\rangle$ , so we can write $r2^{-v_{2}(q)}\equiv 2^{\alpha}\left(\textnormal{mod}\ q2^{-v_{2}(q)}\right)$ for some integer $\alpha$ . Then the cardinality becomes

O(v_{2}(q))+\sharp\left\{v_{2}(q)\leqslant t\leqslant\frac{\log x}{\log 2}:2^{% t-v_{2}(q)-\alpha}\equiv 1\left(\textnormal{mod}\ q2^{-v_{2}(q)}\right)\right% \}\mathds{1}(v_{2}(r)\geqslant v_{2}(q),r2^{-v_{2}(q)}\in\langle 2\rangle).

Denoting the order of $2\left(\textnormal{mod}\ q2^{-v_{2}(q)}\right)$ by $g(q)$ this then becomes

\frac{\mathds{1}(v_{2}(r)\geqslant v_{2}(q),r2^{-v_{2}(q)}\in\langle 2\rangle)% }{g(q)}\frac{\log x}{\log 2}+O(v_{2}(q)).

Note that $v_{2}(q)\leqslant(\log q)/(\log 2)$ . Therefore, if we let

\omega(t)=\frac{1}{t\log 2},\quad\rho(r,q)=\frac{\mathds{1}(v_{2}(r)\geqslant v% _{2}(q),r2^{-v_{2}(q)}\in\langle 2\rangle)}{g(q)}

we have shown that

\sum_{\begin{subarray}{c}1\leqslant m\leqslant x\\ m\equiv r\left(\textnormal{mod}\ q\right)\end{subarray}}k(m)=\rho(r,q)\int_{1}% ^{x}\omega(t)\mathrm{d}t+O(\log q),

with an absolute implied constant. To use Theorem 1.2 we must verify the remaining assumption regarding $\int\omega$ . We will show something more, namely,

{\lim_{P\to\infty}\frac{\mathscr{J}(P)}{\log P}=d\sigma_{\infty},}

(3.10)

where

\mathscr{J}(P):=P^{d}\int\limits_{\begin{subarray}{c}\mathbf{t}\in\mathscr{B}% \\ F(\mathbf{t})>P^{-d}\end{subarray}}\omega(P^{d}f(\mathbf{t}))\mathrm{d}\mathbf% {t}.

By (1.3) we have

\mathscr{J}(P)=\frac{1}{\log 2}\int_{\gamma\in\mathbb{R}}I(\mathscr{B},\gamma)% \int_{P^{-d}}^{b}\frac{\mathrm{e}(-\gamma\mu)}{\mu}\mathrm{d}\mu\mathrm{d}\gamma.

By the trivial bound $|\mathrm{e}(-\gamma\mu)|\leqslant 1$ we deduce for $\gamma\in\mathbb{R}$ and all $\lambda\in(0,b)$ that

\int_{\lambda}^{b}\frac{\mathrm{e}(-\gamma\mu)}{\mu}\mathrm{d}\mu\ll\log\frac{% b}{\lambda}

with an absolute implied constant. Hence, by (2.8), the contribution of $|\gamma|>\sqrt{\log P}$ towards $\mathscr{J}(P)$ is

\ll(\log P)\int_{\gamma>\sqrt{\log P}}|I(\mathscr{B};\gamma)|\mathrm{d}\gamma=% o(\log P).

Similarly, the contribution of $\mu>1/\log P$ is

\ll\int_{|\gamma|\leqslant\sqrt{\log P}}|I(\mathscr{B};\gamma)|\int_{1/\log P}% ^{b}\frac{\mathrm{d}\mu}{\mu}\mathrm{d}\gamma\ll\log\log P.

Thus,

\mathscr{J}(P)=\frac{1}{\log 2}\int_{|\gamma|\leqslant\sqrt{\log P}}I(\mathscr% {B};\gamma)\int_{P^{-d}}^{1/\log P}\frac{\mathrm{e}(-\gamma\mu)}{\mu}\mathrm{d% }\mu\mathrm{d}\gamma+o(\log P).

In this range we have $|\gamma\mu|\leqslant 1/\sqrt{\log P}$ , hence,

\mathrm{e}(-\gamma\mu)=1+O\left(1/\sqrt{\log P}\right).

Substituting in the last expression for $\mathscr{J}$ leads us to

\mathscr{J}(P)=\frac{1}{\log 2}\int_{|\gamma|\leqslant\sqrt{\log P}}I(\mathscr% {B};\gamma)\int_{P^{-d}}^{1/\log P}\frac{\mathrm{d}\mu}{\mu}\mathrm{d}\gamma+O% (E),

where $E$ is given by

\frac{1}{\sqrt{\log P}}\int_{|\gamma|\leqslant\sqrt{\log P}}|I(\mathscr{B};% \gamma)|\int_{P^{-d}}^{1}\frac{\mathrm{d}\mu}{\mu}\mathrm{d}\gamma\ll\frac{1}{% \sqrt{\log P}}\int_{P^{-d}}^{1}\frac{\mathrm{d}\mu}{\mu}\ll\sqrt{\log P}=o(% \log P).

The main term is

		$\displaystyle\frac{1}{\log 2}\int_{\|\gamma\|\leqslant\sqrt{\log P}}I(\mathscr{B% };\gamma)(-\log\log P+d\log P)\mathrm{d}\gamma$
	$\displaystyle=$	$\displaystyle O(\log\log P)+\frac{d\log P}{\log 2}\int_{\|\gamma\|\leqslant\sqrt% {\log P}}I(\mathscr{B};\gamma)\mathrm{d}\gamma=\frac{d\sigma_{\infty}(f)}{\log 2% }\log P+o(\log P),$

which concludes the proof of (1.3).

We are now in position to apply Theorem 1.2. Before doing so, we simplify $\rho(r,q)$ by noting that

\frac{1}{\varphi\left(q2^{-v_{2}(q)}\right)}\sum_{1\leqslant j\leqslant\varphi% \left(q2^{-v_{2}(q)}\right)}\mathds{1}\left(r2^{-v_{2}(q)}\equiv 2^{j}\left(% \textnormal{mod}\ q2^{-v_{2}(q)}\right)\right)=\frac{\mathds{1}(r2^{-v_{2}(q)}% \in\langle 2\rangle)}{g(q)},

where $\varphi$ is Euler’s totient function. Hence, we may write

\rho(r,q)=\frac{\mathds{1}(v_{2}(q)\leqslant v_{2}(r))}{\varphi\left(q2^{-v_{2% }(q)}\right)}\sharp\left\{1\leqslant j\leqslant\varphi(q2^{-v_{2}(q)}):r2^{-v_% {2}(q)}\equiv 2^{j}\left(\textnormal{mod}\ q2^{-v_{2}(q)}\right)\right\}.

Let $t_{p}(z)$ and $T_{z}$ be as in Theorem 1.2 and let $w=T_{z}2^{-t_{2}(z)}$ . Each $\mathbf{t}\in(\mathbb{Z}/T_{z}\mathbb{Z})^{n}$ can be uniquely written as $\mathbf{t}=w\mathbf{x}+2^{t_{2}(z)}\mathbf{y}$ where $\mathbf{x}\in(\mathbb{Z}/2^{t_{2}(z)}\mathbb{Z})^{n}$ and $\mathbf{y}\in(\mathbb{Z}/w\mathbb{Z})^{n}$ . We find that

\sum_{\mathbf{t}\in(\mathbb{Z}/T_{z}\mathbb{Z})^{n}}\frac{\rho(f(\mathbf{t}),T% _{z})}{T_{z}^{n-1}}=\frac{\sharp\left\{\mathbf{x}\in(\mathbb{Z}/2^{t_{2}(z)}% \mathbb{Z})^{n}:f(\mathbf{x})=0\right\}}{2^{(n-1)t_{2}(z)}}\frac{\sharp\left\{% \mathbf{y}\in(\mathbb{Z}/w\mathbb{Z})^{n},h\in\mathbb{Z}/\varphi(w)\mathbb{Z}:% f(\mathbf{y})=2^{h}\right\}}{\varphi(w)w^{n-1}}.

The first fraction in the right hand side converges to $\sigma_{2}(f)$ as $z\to\infty$ . By the Chinese remainder theorem the second fraction is

\prod_{2<p\leqslant z}\frac{\sharp\left\{\mathbf{y}\in(\mathbb{Z}/p^{t_{p}(z)}% \mathbb{Z})^{n},h\in\mathbb{Z}/\varphi(p^{t_{p}(z)})\mathbb{Z}:f(\mathbf{y})=2% ^{h}\right\}}{(p-1)p^{nt_{p}(z)}}

and letting $z\to\infty$ completes the proof.

References

[1] B. J. Birch, Forms in many variables. Proc. R. Soc. Lond. Ser. A, Math. Phys. Sci. 265 (1962), 245–263.
[2] V. Blomer, Forms in many variables. Sums of Hecke eigenvalues over values of quadratic polynomials. Int. Math. Res. Not. IMRN 16 (2008), 245–263.
[3] R. de la Bretèche, T. D. Browning and E. Peyre, On Manin’s conjecture for a family of Châtelet surfaces. Ann. of Math. 175 (2012), 297–343.
[4] R. de la Bretèche and T. D. Browning, Binary forms as sums of two squares and Châtelet surfaces. Isr. J. Math. 191 (2012), 973–1012.
[5] R. de la Bretèche and G. Tenenbaum, Moyennes de fonctions arithmétiques de formes binaires. Mathematika 58(2) (2012), 290–304.
[6] R. de la Bretèche and G. Tenenbaum, Sur la conjecture de Manin pour certaines surfaces de Châtelet. Journal of the Institute of Mathematics of Jussieu 12(4) (2013), 759–819.
[7] M. Bright, S. Ramdorai, A. Skorobogatov, E. Sofos and P. Vishe, New Directions in Rational Points, BIRS Chennai, link to the report, (2024).
[8] T. D. Browning, Cubic forms and the circle method. Birkhäuser (2021), xiv–166.
[9] T. D. Browning, L. Matthiesen and A. N. Skorobogatov, Rational points on pencils of conics and quadrics with many degenerate fibers. Ann. of Math. 180 (2014), 381–402.
[10] T. D. Browning and S. M. Prendiville, Improvements in Birch’s theorem on forms in many variables. J. reine angew. Math. 731 (2017), 203–234.
[11] J. Brüdern and T. D. Wooley, The Hasse principle for pairs of diagonal cubic forms. Ann. of Math. 166 (2007), 865–895.
[12] Y. Bugeaud, M. Mignotte and S. Siksek, Classical and modular approaches to exponential Diophantine equations. I: Fibonacci and Lucas perfect powers. Ann. of Math. 163 (2006), 969–1018.
[13] F. Chamizo, Correlated sums of $r(n)$ . J. Math. Soc. Japan 51 (1999), 237–252.
[14] L. Chiriac and L. Yang, Summing Hecke eigenvalues over polynomials. Math. Z. 302 (2022), 643–662.
[15] S. Daniel, On the divisor-sum problem for binary forms. J. reine angew. Math. 507 (1999), 107–129.
[16] H. Davenport, On Some Infinite Series Involving Arithmetical Functions (II), The Quarterly Journal of Mathematics, os-8, (1937), 313–320.
[17] P. Deligne, La conjecture de Weil. I, Inst. Hautes Études Sci. Publ. Math. 43 (1974), 273–307.
[18] P. Deligne and J.-P. Serre, Formes modulaires de poids 1. Ann. Sci. École Norm. Sup. (4) 7 (1974), 507–530.
[19] K. Destagnol, La conjecture de Manin sur les surfaces de Châtelet. Acta Arithmetica. 174(1) (2016), 31–97.
[20] K. Destagnol and E. Sofos, Rational points and prime values of polynomials in moderately many variables. Bull. Sci. Math. 156 (2019), 102794, 33.
[21] T. Estermann, An asymptotic formula in the theory of numbers. Proc. London Math. Soc. 34 (1932), 280–292.
[22] G. Faltings, Endlichkeitssätze für abelsche Varietäten über Zahlkörpern. Invent. Math. 73 (1983), 349–366.
[23] J. Friedlander and H. Iwaniec, The polynomial $X^{4}+Y^{2}$ captures its primes. Ann. of Math. 148(3) (1998), 945–1040.
[24] R. Heath-Brown, Primes represented by $x^{3}+2y^{3}$ . Acta Math. 186(1) (2001), 1–84.
[25] R. Heath-Brown, Linear relations amongst sums of two squares. Number theory and algebraic geometry London Math. Soc. Lecture Note Ser., Cambridge Univ. Press, Cambridge 303 (2003), 133–176.
[26] C. Hooley, On the number of divisors of a quadratic polynomial. Acta Math. 110 (1963), 97–114.
[27] H. Iwaniec, Topics in Classical Automorphic Forms. Graduate Studies in Mathematics. American Mathematical Society Providence, RI 17 (1997).
[28] H. Iwaniec and E. Kowalski, Analytic number theory. Colloquium Publications. American Mathematical Society Providence, RI 53 (2004), xi–615.
[29] K. Matomäki, M. Radziwiłł and T. Tao, An averaged form of Chowla’s conjecture. Algebra Number Theory 9 (2015), 2167–2196.
[30] M. Nair, Multiplicative functions of polynomial values in short intervals. Acta Arith. 62 (1992), 257–269.
[31] M. Nair and G. Tenenbaum, Short sums of certain arithmetic functions. Acta Math. 180 (1998), 119–144.
[32] M. Northey and P. Vishe, On the Hasse principle for complete intersections. Compos. Math. 160 (2024), 771–835.
[33] H. Pasten, Powerful values of polynomials and a conjecture of Vojta. J. Number Theory 133 (2013), 2964–2998.
[34] M. Pieropan and D. Schindler, Hyperbola method on toric varieties. Journal de l’École Polytechnique 11 (2024), 107–157.
[35] P. Pongsriiam and R. C. Vaughan, The divisor function on residue classes I. Acta Arith. 168 (2015), 369–382.
[36] S.L. Rydin Myerson, Quadratic forms and systems of forms in many variables. Invent. math. 213 (2018), 205–235.
[37] W. Schmidt, Nothcott’s theorem on heights II. The quadratic case. Acta Arith. 70(4) (1995), 343–375.
[38] T. N. Shorey and R. Tijdeman, Exponential Diophantine equations. Cambridge Tracts in Mathematics Cambridge University Press 87 (1986).
[39] T. Tao and J. Teräväinen, The structure of logarithmically averaged correlations of multiplicative functions, with applications to the Chowla and Elliott conjectures. Duke Math. J. 168 (2019), 1977–2027.
[40] N. Templier, A nonsplit sum of coefficients of modular forms. Duke Math. J. 157 (2011), 109–165.
[41] N. Templier and J. Tsimerman, Non-split sums of coefficients of $GL(2)$ -automorphic forms. Israel J. Math. 195 (2013), 677–723.
[42] D. I. Tolev, On the remainder term in the circle problem in an arithmetic progression. Proc. Steklov Inst. Math. 276 (2012), 261–274.

Averages of arithmetic functions over polynomials in many variables

Abstract.

2020 Mathematics Subject Classification:

1. Introduction

1.1. The general result

Example 1.1.

Definition 1.1.

Definition 1.2.

Theorem 1.1.

Remark.

Remark.

1.2. A more accessible version

Theorem 1.2.

1.3. Analytic Hasse principle for intersections

Theorem 1.3.

1.4. Subsets of integers

Theorem 1.4.

1.5. Points on non-algebraic varieties

Theorem 1.5.

1.6. Chowla’s conjecture

Theorem 1.6.

1.7. Sums of Fourier coefficients of modular forms over values of multivariable polynomials

Theorem 1.7.

1.8. Divisor function over values of polynomials

Theorem 1.8.

1.9. Rational points on a hyperelliptic hypersurface

Theorem 1.9.

1.10. m𝑚mitalic_m-full numbers represented by polynomials

Theorem 1.10.

Acknowledgements.

2. Proof of Theorems 1.1 and 1.2

Lemma 2.1.

Proof.

Lemma 2.2.

Proof.

Lemma 2.3.

Proof.

Lemma 2.4.

Proof.

Lemma 2.5.

Proof.

Lemma 2.6.

Proof.

Lemma 2.7.

Proof.

Lemma 2.8.

Proof.

2.1. Proof of Theorem 1.2

3. Applications

3.1. The proof of Theorem 1.6

3.2. The proof of Theorem 1.7

3.3. The proof of Theorem 1.8

3.4. A shifted convolution problem in arithmetic progressions

Theorem 3.1.

Lemma 3.1.

Proof.

Lemma 3.2.

Proof.

Lemma 3.3.

Proof.

Lemma 3.4.

Proof.

3.5. The proof of Theorem 3.1

3.6. The proof of Theorem 1.3

Lemma 3.5.

Proof.

3.7. The proof of Theorem 1.4

3.8. The proof of Theorem 1.9

3.9. The proof of Theorem 1.10

3.10. The proof of Theorem 1.5

References

1.10. $m$ -full numbers represented by polynomials