Introduction

For a primitive Dirichlet character χ\chi of conductor qq, the associated Dirichlet LL-function is defined for Re⁡s>1\operatorname{Re} s > 1 by the formula

L(s,χ)=∑n≥1χ(n)ns,Re⁡s>1,L(s,\chi)=\sum_{n\ge1}\frac{\chi(n)}{n^s}, \qquad \operatorname{Re} s>1,

and admits a meromorphic continuation to s∈Cs \in\mathbb{C} [22] §4.6. When χ\chi is the trivial character, L(s,χ)L(s,\chi) recovers the Riemann zeta function ζ(s)\zeta(s) [41].

The zeros of Dirichlet LL-functions are of significant interest, such as for their role in governing the distribution of primes in arithmetic progressions. The Generalized Riemann Hypothesis predicts that all zeros of L(s,χ)L(s,\chi) in the critical strip 0<Re⁡s<10 < \operatorname{Re} s < 1 satisfy Re⁡s=1/2\operatorname{Re} s = 1/2. For ζ(s)\zeta(s), the weaker assertion that there exists ε>0\varepsilon> 0 such that ζ(s)\zeta(s) has no zeros in Re⁡s>1−ε\operatorname{Re} s > 1 - \varepsilon is called a quasi-Riemann Hypothesis (e.g., in [2] §1, [34], p. 274, and [4] §2). For Dirichlet LL-functions, we consider the analogous assertion with a single ε>0\varepsilon> 0 valid for every primitive character χ\chi, independently of both conductor and height.

Let K=Q(−3)K = \mathbb{Q}(\sqrt{-3}). The main result of this paper is the following.

Theorem 1.1. Every finite-order Hecke LL-function over KK has no zeros in the half-plane Re⁡s>11/12\operatorname{Re}s > 11/12. Thus every Dirichlet LL-function, including ζ(s)\zeta(s), has no zeros for Re⁡s>11/12\operatorname{Re}s > 11/12.

Section 3 proves the theorem assuming Proposition 3.1, whose proof is completed in Section 5.6.

Theorem 1.1 is a natural intermediate step toward the stronger zero-free region Re⁡s>7/8\operatorname{Re}s > 7/8 established in [36], Theorem 1.1. For simplicity, we have isolated Theorem 1.1 and its proof here. By a standard explicit-formula argument (see [6], Chapters 19–20), Theorem 1.1 gives the following quantitative form of the prime number theorem in arithmetic progressions.

Corollary 1.2. Let π(x;q,a)\pi(x;q,a) count the primes p≤xp \leq x with p≡a(modq)p \equiv a \pmod{q}. Write φ\varphi for Euler’s totient function. For x≥2x \geq2,

sup⁡1≤q≤x(a,q)=1∣π(x;q,a)−1φ(q)∫2x dtlog⁡t∣≪x11/12log⁡x,\sup_{\substack{1\le q\le x\\(a,q)=1}} \biggl|\pi(x;q,a)-\frac1{\varphi(q)} \int_2^x\frac{\,d t}{\log t}\biggr| \ll x^{11/12}\log x,

where the implied constant is absolute and effective.

Theorem 1.1 has a number of additional arithmetic consequences. By the zero-free-region method originating with Rodosskiĭ [42], in the form recorded in [33], Theorem 13.12, the least quadratic nonresidue modulo an odd prime pp is bounded by a fixed power of log⁡p\log p, which in particular proves Vinogradov’s conjecture [49] that this nonresidue is ≪εpε\ll_{\varepsilon}p^{\varepsilon} for every ε>0\varepsilon> 0. See also Bhargava, Ivanyos, Mittal, and Saxena [3], Theorem 6.7 for this consequence of a fixed zero-free half-plane. From a computational number theory perspective, this bound allows square roots modulo pp to be extracted deterministically in time polynomial in log⁡p\log p using the Tonelli–Shanks algorithm (see [11], §2.9). Theorem 1.1 also yields a deterministic polynomial-time implementation of Miller’s primality test [32]; note that polynomial-time primality testing was previously known unconditionally via the Agrawal–Kayal–Saxena algorithm [1]. For negative fundamental discriminants DD, Littlewood’s short Euler-product argument [29], applied to the fixed zero-free half-plane in Theorem 1.1, together with Dirichlet’s class-number formula (see [6], Chapter 6), yields the effective class-number bound h(D)≫∣D∣/log⁡log⁡∣D∣h(D)\gg\sqrt{|D|}/\log\log|D|, with an absolute computable implied constant. The class-group computations of Elsenhans, Klüners, and Nicolae [8], Theorem 2, building on Weinberger [50], give a complete list of imaginary quadratic fields with class-group exponent dividing two when there are no Landau–Siegel zeros. Together with Grube’s characterization of idoneal numbers and his reduction to fundamental discriminants [15] (see [23], Theorem 6 and §2.3), this shows that Theorem 1.1 confirms the completeness of Euler’s list of 65 idoneal numbers.

Prior work

There is a long line of work establishing zero-free regions for Dirichlet LL-functions. Hadamard and de la Vallée Poussin proved the Prime Number Theorem in 1896 by establishing that ζ(s)\zeta(s) has no zeros on the line Re⁡s=1\operatorname{Re}s=1 [16, 46]. De la Vallée Poussin subsequently obtained a quantitative zero-free region to the left of this line [47]. For nonprincipal primitive Dirichlet characters, the classical zero-free-region theorem of Grönwall and Titchmarsh [14, 45], in the modern form given in [22], Theorem 5.26 (see also [6], Chapter 14), gives an absolute, effective constant c0>0c_0>0 such that

L(σ+it,χ)≠0ifσ>1−c0log⁡(q(∣t∣+3)),L(\sigma+\mathrm i t,\chi)\ne0 \quad\text{if}\quad \sigma>1-\frac{c_0}{\log(q(|t|+3))},

with at most one exception for each primitive character. Any exception is a simple real zero, can occur only if χ\chi is quadratic, and is called a Landau–Siegel zero. Theorem 1.1 rules out the existence of such a zero.

There are two parameters in (1.1): the conductor qq and the height ∣t∣|t|. Even when qq is fixed, its width tends to zero with the height. Vinogradov and Korobov developed methods [48, 24] to prove that

ζ(σ+it)≠0ifσ>1−c1(log⁡∣t∣)2/3(log⁡log⁡∣t∣)1/3,∣t∣≥3,\zeta(\sigma+\mathrm i t)\ne0 \quad\text{if}\quad \sigma>1-\frac{c_1} { (\log|t|)^{2/3}(\log\log|t|)^{1/3}},\qquad |t|\ge3,

for an absolute c1>0c_1>0 (Ford [9] gives an explicit value of c1c_1). The assertion of Theorem 1.1 is a half-plane of fixed width, independent of both conductor and height.

As we will see in the outline, the proof of Theorem 1.1 draws on modern developments in character large sieves and metaplectic theta series. In particular, the connection between cubic Gauss sums and cubic theta coefficients [26, 38] links the argument to work on Patterson’s conjecture by Heath-Brown and Patterson [20], Heath-Brown [19], and Dunn and Radziwiłł [7]. The recursive large-sieve arguments also build on Heath-Brown’s proof of the quadratic large sieve inequality [18].

Organization

Section 2 gives a high-level outline of the argument. Section 3 deduces the zero-free region from a mean-square estimate for twisted Möbius sums. Section 4 reduces this estimate to a dual mean square with cubic Gauss-sum coefficients. Section 5 states the completed mean-square and transfer estimates, combines them into a recursive inequality, and proves the required mean-square bound. Sections 6 and 7 prove the completed estimate and the bounds for the remaining cube-divisor sums, respectively. Appendix A supplies the arithmetic identities and the detailed theta calculation. Appendix B records the smooth separation lemma used throughout.

Notation

Write O=OK=Z[ω]\mathcal{O}=\mathcal{O}_K=\mathbb{Z}[\omega], where ω=e2πi/3\omega=e^{2\pi\mathrm i/3}. For a∈Ka\in K, write NK/Q(a)=aa‾=∣a∣2\mathrm N_{K/\mathbb Q}(a)=a\overline a=|a|^2 for its norm. For a nonzero integral ideal a⊆O\mathfrak{a}\subseteq\mathcal{O}, write NK/Q(a)=#(O/a)\mathrm N_{K/\mathbb Q}(\mathfrak a)=\#(\mathcal O/\mathfrak a). We use 1C\mathbf1_{\mathcal C} for the indicator of a condition C\mathcal C. For a prime ideal represented by pp and a nonzero element or ideal aa, write vp(a)v_p(a) for its exponent in the prime factorization of aa.

We use standard asymptotic notation, writing f=O(g)f=O(g) or f≪gf\ll g when ∣f∣≤Cg|f|\le Cg, where g≥0g\ge0 and C>0C>0 is an absolute constant unless otherwise specified. Subscripts indicate possible dependence of the implied constant: for example, f≪ν,W,εgf\ll_{\nu,W,\varepsilon}g means ∣f∣≤Cν,W,εg|f|\le C_{\nu,W,\varepsilon}g, where the constant may depend on ν\nu, WW, ε\varepsilon but is uniform in all other varying parameters. For nonnegative f,gf,g, we write f≍gf\asymp g when f≪gf\ll g and g≪fg\ll f. A dyadic norm range is an interval R≤NK/Q(a)<2RR\le\mathrm N_{K/\mathbb Q}(a)<2R; a sum over dyadic scales uses R=2jR=2^j.

We use D≥2D\ge2 as an ambient size parameter; in the main argument, it is the original norm scale introduced in the outline below. Auxiliary scales may vary within ranges bounded by fixed powers of DD. We write

A≼BifA≪εDεBfor every ε>0,A\preccurlyeq B \quad\text{if} \quad A\ll_{\varepsilon}D^{\varepsilon}B \quad\text{for every } \varepsilon>0,

uniformly in the varying parameters over their stated ranges. Implied constants may also depend on the fixed data specified in each statement.

Outline of the argument

In this section, we sketch the proof of Theorem 1.1, suppressing various coprimality conditions, local factors, and details of smoothing. The precise statements appear in Sections 3–7.

Step 1: Reduction to power-saving estimates for twisted Möbius sums

Fix a finite-order Hecke character ν\nu of KK. Let LK(s,ν)L_K(s,\nu) be the corresponding Hecke LL-function. We extend ν\nu by zero to ideals not coprime to its conductor. We will deduce Theorem 1.1 from a power-saving estimate for ν\nu-twisted Möbius sums.

Let μ\mu denote the ideal Möbius function on KK. All original ideal sums run over nonzero integral ideals. As in Section 3, the ideals in these sums are understood to be prime to 2, 3, and the conductor of ν\nu. For a norm scale D>0D>0, consider

A1(D)=∑nμ(n)ν(n)W(NK/Q(n)/D),A_1(D)=\sum_{\mathfrak n} \mu(\mathfrak n)\nu(\mathfrak n) W(\mathrm N_{K/\mathbb Q}(\mathfrak n)/D),

where W∈Cc∞((0,∞);C)W\in C_c^\infty((0,\infty);\mathbb C) is a smooth cutoff function. Thus A1(D)A_1(D) is a smoothed ν\nu-twisted Möbius sum over ideals of norm comparable to DD, with roughly DD terms. We seek a power saving over this size: for a fixed δ>0\delta>0 and every ε>0\varepsilon>0,

A1(D)≪ν,W,εD1−δ+ε for every W∈Cc∞((0,∞);C).A_1(D)\ll_{\nu,W,\varepsilon}D^{1-\delta+\varepsilon} \text{ for every }W\in C_c^\infty((0,\infty);\mathbb C).

The implication from (2.1) to a zero-free half-plane is the smoothed Hecke version of the classical relation between Möbius sums and zero-free regions. In the zeta-function case, the equivalence between the Riemann Hypothesis and the bound ∑n≤xμ(n)≪εx1/2+ε\sum_{n\leq x}\mu(n)\ll_\varepsilon x^{1/2+\varepsilon} is due to Littlewood [28]; see also [31], §1. Section 3 gives the complete argument needed here.

Step 2: Embedding the sum in a family

In order to estimate A1(D)A_1(D), we embed it into a family of such sums, parametrized by u∈OKu\in\mathcal{O}_K, by introducing a sextic twist. We call an element of O\mathcal{O} primary if it is congruent to 1 modulo 3. The ring O=Z[ω]\mathcal{O}=\mathbb{Z}[\omega] is Euclidean, so it has unique factorization and every ideal is principal; for an ideal coprime to 3, multiplying any generator by a unique unit gives n≡1(mod3)n\equiv1\pmod3, since the six units represent the six invertible residue classes modulo 3. For a primary element nn, we write μ(n)=μ((n))\mu(n)=\mu((n)) and ν(n)=ν((n))\nu(n)=\nu((n)); divisor sums count each ideal divisor once, using its primary generator. For an ideal n\mathfrak n prime to 6, use its unique primary generator nn and write

χn(u)=χn(u)=(u/n)6.\chi_{\mathfrak n}(u)=\chi_n(u)=(u/n)_6.

Here (u/n)6(u/n)_6 is the sextic residue symbol, extended by zero when (u,n)≠1(u,n)\ne1.1 We also write (un)2\left(\frac{u}{n}\right)_2 and (un)3\left(\frac{u}{n}\right)_3 for the quadratic and cubic residue symbols, defined by the same convention with 6 replaced by 2 and 3, respectively. When (n,6)=1(n,6)=1, they equal χn(u)3\chi_n(u)^3 and χn(u)2\chi_n(u)^2. Restrictions that keep these symbols defined are understood when omitted in the outline. The family is

Au(D)=∑nμ(n)ν(n)χn(u)W(NK/Q(n)/D).A_u(D)=\sum_{\mathfrak n} \mu(\mathfrak n)\nu(\mathfrak n)\chi_{\mathfrak n}(u) W(\mathrm N_{K/\mathbb Q}(\mathfrak n)/D).

Fix 0<ϑ≤1/100<\vartheta\leq1/10. The crucial estimate, Proposition 3.1, is that

∑0<NK/Q(u)≤H∣Au(D)∣2≪ν,W,ϑ,εD1+εH,H=D1+ϑ.\sum_{0<\mathrm N_{K/\mathbb Q}(u)\le H}|A_u(D)|^2 \ll_{\nu,W,\vartheta,\varepsilon}D^{1+\varepsilon}H, \qquad H=D^{1+\vartheta}.

In a mean square, we call the variable in the outer sum the row and the variable in the inner sum the column. In (2.3) these are uu and n\mathfrak n, respectively. Up to the factor of DεD^\varepsilon, (2.3) can be thought of as saying that the family {Au(D)}\{A_u(D)\} exhibits “square-root cancellation” on average (in the L2L^2 sense) over uu.

The mean-square estimate (2.3) gives the desired power saving for A1(D)A_1(D) because Ap6(D)≈A1(D)A_{p^6}(D) \approx A_1(D) for many primes pp. For a primary prime pp with Y/2<NK/Q(p)≤YY/2<\mathrm N_{K/\mathbb Q}(p)\le Y, the identity χn(p6)=1p∤n\chi_n(p^6)=\mathbf1_{p\nmid n} leaves only OW(D/Y)O_W(D/Y) differing terms and hence gives

Ap6(D)=A1(D)+OW(D/Y).A_{p^6}(D) = A_1(D) + O_W(D/Y).

By Landau’s prime ideal theorem [27] (see [22], Theorem 5.33), there are ≍Y/log⁡Y\asymp Y/\log Y choices of such pp. Taking Y=H1/6Y = H^{1/6}, considering the contribution of such terms to (2.3) and using (2.4) gives

∣A1(D)∣2≪D1+εH5/6+D2H−1/3.|A_1(D)|^2 \ll D^{1+\varepsilon}H^{5/6} + D^2H^{-1/3}.

With H=D1+ϑH = D^{1+\vartheta}, this gives A1(D)≪D11/12+5ϑ/12+εA_1(D) \ll D^{11/12+5\vartheta/12+\varepsilon} for every fixed 0<ϑ≤1/100 < \vartheta\le1/10. Choosing ϑ\vartheta sufficiently small for each requested exponent loss yields A1(D)≪ν,W,εD11/12+εA_1(D) \ll_{\nu,W,\varepsilon} D^{11/12+\varepsilon}. Therefore it remains to establish (2.3).

Step 3: Poisson summation

We turn to the task of proving the mean-square estimate (2.3). Now that we have introduced the row variable uu, we can apply Poisson summation in this variable. This introduces sextic Gauss sums, arising from the Fourier transforms of the sextic characters. Upon application of the Gauss–Jacobi identities, these sextic Gauss sums absorb the factor of μ\mu and turn into cubic Gauss sums. We interpret the resulting cubic Gauss sums as coefficients of Kubota’s cubic theta function. In the next step, we use the automorphy of this theta function. For now, we explain the appearance of these Gauss sums in more detail.

Define the additive character

e(z)=exp⁡(4πi Im⁡z/3)(z∈C).e(z)=\exp(4\pi\mathrm i\,\operatorname{Im}z/\sqrt3)\quad(z\in\mathbb C).

For a squarefree Eisenstein integer n≡1(mod3)n \equiv1 \pmod{3} prime to 66 and an integer jj, define the normalized Gauss sums

γj(n)=1NK/Q(n)∑x mod nχn(x)je(x/n).\gamma_j(n)=\frac{1}{\sqrt{\mathrm N_{K/\mathbb Q}(n)}} \sum_{x\bmod n}\chi_n(x)^j e(x/n).

For j=−1j=-1, we interpret χn−1\chi_n^{-1} as the conjugate character χn‾\overline{\chi_n}.

Consider the classical identity

(−1m)=(1m∑x mod m(xm)e2πix/m)2,\left(\frac{-1}{m}\right) = \left(\frac{1}{\sqrt{m}}\sum_{x \bmod m}\left(\frac{x}{m}\right)e^{2\pi i x/m}\right)^2,

where mm is an odd squarefree positive integer and the symbols are Jacobi symbols. We would like a similar decomposition for μ\mu instead of (−1m)\left(\frac{-1}{m}\right). Following the work of Hasse [17], pp. 443–445 and Heath-Brown [19], (2), we derive the following identity in Appendix A.1, valid for squarefree primary nn away from a fixed set of excluded primes:

μ(n)γ−1(n)=χn(−1)G(n)−1α(n)‾γ2(n),α(n)=n∣n∣,G(n)=χn(4)‾γ3(n).\mu(n)\gamma_{-1}(n) = \chi_n(-1)G(n)^{-1}\overline{\alpha(n)}\gamma_2(n), \qquad\alpha(n) = \frac{n}{|n|}, \qquad G(n) = \overline{\chi_n(4)}\gamma_3(n).

Here G(n)G(n) is a fixed ray-class factor. The factor γ−1(n)\gamma_{-1}(n) comes from Poisson summation; as promised, it combines with μ(n)\mu(n) to produce the cubic coefficient α(n)‾γ2(n)\overline{\alpha(n)}\gamma_2(n) (up to ray-class factors). Further details are given in Section 4.

We call the sums obtained after applying Poisson summation dual sums.2 Rearranging these dual sums reduces the problem of estimating ∑u∣Au(D)∣2\sum_u |A_u(D)|^2 to estimating the following family of column sums (for clarity, we have suppressed auxiliary twists and coprimality conditions):

Bh(X)=∑n∈On≡1 (3)n squarefreeα(n)‾γ2(n)ξ(n)χn(h)U(NK/Q(n)/X).B_h(X)=\sum_{\substack{n\in\mathcal O\\n\equiv1\ (3)\\n\ {\rm squarefree}}} \overline{\alpha(n)}\gamma_2(n)\xi(n)\chi_n(h) U(\mathrm N_{K/\mathbb Q}(n)/X).

Here hh is the new row variable, Fourier dual to uu, while UU is a smooth weight restricting the column variable nn to norm comparable to XX. We choose ξ\xi to be either νη\nu\eta or ν‾η\overline{\nu}\eta, where η\eta ranges over a finite set of ray class characters of fixed modulus.3 For each choice of ξ\xi, this defines a family of sums Bh(X)B_h(X).

After treating the diagonal and separating the smooth weights, we roughly get that

∑0<NK/Q(u)≤H∣Au(D)∣2≼DH+HD∑0<NK/Q(h)≪H∣Bh(D)∣2.\sum_{0<\mathrm N_{K/\mathbb Q}(u)\le H}|A_u(D)|^2 \preccurlyeq DH+\frac HD \sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|B_h(D)|^2.

Here H≍D2/H\mathcal{H}\asymp D^2/H is the dual norm scale. This schematic comparison suppresses the common factors arising when the square is expanded. Thus the desired estimate is

∑0<NK/Q(h)≪H∣Bh(D)∣2≼D2.\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|B_h(D)|^2\preccurlyeq D^2.

Step 4: Relation to cubic Gauss sums and Kubota’s cubic theta function

To estimate (2.9), we first identify the cubic Gauss coefficients with Fourier coefficients of a theta function. Its transformation law applies to a sum with extra cube factors, which we will remove in Step 5.

Let θ(z,v)\theta(z,v) be Kubota’s cubic theta function on hyperbolic three-space (z∈Cz\in\mathbb{C}, v∈R>0v\in\mathbb{R}_{>0}), obtained as a residue of a cubic metaplectic Eisenstein series [26, 38]. We use the normalization of [7], §5.1, (5.6)–(5.8), recalled below in (6.1).

Put λ=1+2ω\lambda=1+2\omega. For primary elements n,b∈O=Z[ω]n,b\in\mathcal{O}=\mathbb{Z}[\omega], with nn squarefree and (nb,6)=1(nb,6)=1, let cθ(nb3)c_\theta(nb^3) denote the Fourier coefficient of θˉ\bar{\theta} indexed by λ−3nb3\lambda^{-3}nb^3 in its expansion in zz. Patterson’s formula [38], Theorem 8.1, recorded in [7], (5.7), gives

cθ(nb3)=35/2∣b∣ χn(λ)2‾γ2(n).c_\theta(nb^3) =3^{5/2}|b|\,\overline{\chi_n(\lambda)^2}\gamma_2(n).

Taking b=1b=1 in (2.10) and substituting into (2.7) yields

Bh(X)=3−5/2∑n∈On≡1 (3)n squarefreecθ(n)α(n)‾χn(λ)2ξ(n)χn(h)U(NK/Q(n)/X).B_h(X)=3^{-5/2}\sum_{\substack{n\in\mathcal O\\n\equiv1\ (3)\\n\ {\rm squarefree}}} c_\theta(n)\overline{\alpha(n)}\chi_n(\lambda)^2 \xi(n)\chi_n(h)U(\mathrm N_{K/\mathbb Q}(n)/X).

Thus (2.11) expresses Bh(X)B_h(X) as a smoothed, twisted sum of the squarefree-index coefficients of θˉ\bar{\theta}. This connection was used by Heath-Brown and Patterson to study Kummer sums [20].

To use the summation formula for theta coefficients, we embed the sum (2.11) in a completed sum, by including terms indexed by nb3nb^3. This use of cube completion follows Dunn and Radziwiłł [7], Lemma 5.4 and Proposition 5.3, with the underlying theta coefficients given by Patterson [38], Theorem 8.1. For fixed hh, ξ\xi, and UU, define

Th(X):=135/2X∑n,b∈On,b≡1 (3)n squarefreecθ(nb3)α(nb3)‾χnb3(λ)2ξ(nb3)χnb3(h)U(NK/Q(nb3)/X).T_h(X):=\frac{1}{3^{5/2}\sqrt X} \sum_{\substack{n,b\in\mathcal O\\n,b\equiv1\ (3)\\n\ {\rm squarefree}}} c_\theta(nb^3)\overline{\alpha(nb^3)}\chi_{nb^3}(\lambda)^2 \xi(nb^3)\chi_{nb^3}(h)U(\mathrm N_{K/\mathbb Q}(nb^3)/X).

The b=1b=1 terms are exactly X−1/2Bh(X)X^{-1/2}B_h(X); the other terms supply the cube indices in the theta expansion. The first goal is to bound the mean square of Th(X)T_h(X) over hh.

The automorphy of θ\theta gives a summation formula that transforms the completed column sum Th(X)T_h(X), with the row hh held fixed, into dual sums of theta coefficients with new smooth weights and character twists. A cusp is represented by a boundary point in K∪{∞}K \cup\{\infty\}. A cusp expansion is the Fourier expansion in the horizontal variable after a change of coordinates taking infinity to that point. We call its Fourier coefficients cusp coefficients. Our starting point is a variant (established in Appendix A.2) of the theta transformation of Dunn and Radziwiłł [7] §5, extending work of Patterson [38] and Yoshimoto [51]. The key point is how the character twists change. Suppose for illustration that hh is squarefree and primary, with (h,6)=1(h,6)=1. Then, by Proposition 6.2, the transformed expression is a finite linear combination of sums of the form

∑0≠m∈Odθ(m)α(m)NK/Q(m) χh(m)3V∗♯ ⁣(NK/Q(m)XNK/Q(h)2).\sum_{0\ne m\in\mathcal O} \frac{d_\theta(m)\alpha(m)}{\sqrt{\mathrm N_{K/\mathbb Q}(m)}}\, \chi_h(m)^3 V_*^\sharp\!\Bigl(\frac{\mathrm N_{K/\mathbb Q}(m)X}{\mathrm N_{K/\mathbb Q}(h)^2}\Bigr).

The dual column index is mm, while hh remains the row index. Here dθ(m)d_\theta(m) denotes the coefficient at the Fourier index λ−4m\lambda^{-4}m in a cusp expansion of θˉ\bar{\theta}, and V∗♯V_*^\sharp is the transform in (6.6), applied to V∗(y)=y1/2U(y)V_*(y)=y^{1/2}U(y). We have suppressed a finite sum over these cusp expansions, fixed periodic twists, bounded prefactors, and fixed scale constants. Within each fixed ray class of hh, the coefficient sequences and transformed weights are independent of hh by Lemma 6.3. Proposition 6.2 gives the precise formula, writing d(ℓ)d(\ell) for the cusp coefficient at ℓ=λ−4m\ell=\lambda^{-4}m.

Since V∗♯V_*^\sharp decays rapidly at infinity, the effective norm range in (2.13) is NK/Q(m)≪NK/Q(h)2/X\mathrm N_{K/\mathbb Q}(m)\ll\mathrm N_{K/\mathbb Q}(h)^2/X, in place of the original range NK/Q(nb3)≍X\mathrm N_{K/\mathbb Q}(nb^3)\asymp X. These dual sums arise from the theta transformation and involve theta coefficients on the transformed scale, now twisted by the quadratic character χh3\chi_h^3. This character arises because, at each prime p∣hp\mid h, the Fourier and theta factors combine as

χp−1χp−2=χp−3=χp3.\chi_p^{-1}\chi_p^{-2}=\chi_p^{-3}=\chi_p^3.

Now a key point is that since χp3\chi_p^3 is quadratic, Goldmakher and Louvel’s quadratic large sieve [13], Theorem 1.1 and Corollary 1.2 (a generalization of Heath-Brown’s quadratic large sieve [18] to number fields) bounds the mean square of the sums in (2.13) as hh varies. For squarefree quadratic families with row and column norm ranges MM, LL, the quadratic large sieve gives the factor M+LM+L, up to (ML)ε(ML)^\varepsilon. Crucially, this avoids the additional term (ML)2/3(ML)^{2/3} in Blomer, Goldmakher, and Louvel’s general higher-order large sieve [5], Theorem 1.3.

We then apply Cauchy–Schwarz in the cube variable and the quadratic large sieve in the squarefree column variable to estimate the mean square of Th(X)T_h(X). The details are given in the proof of Proposition 5.2, which gives, for H,X≥1\mathcal{H},X\ge1,

∑0<NK/Q(h)≪H∣Th(X)∣2≪ε,ξ,U(HX)ε(H+H2X).\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|T_h(X)|^2 \ll_{\varepsilon,\xi,U}(\mathcal H X)^\varepsilon \Bigl(\mathcal H+\frac{\mathcal H^2}{X}\Bigr).

At X=DX=D and H≤D\mathcal{H}\le D, the heuristic comparison Bh(D)≈D Th(D)B_h(D)\approx\sqrt D\,T_h(D) would therefore give the desired D2D^2 bound. The remaining task is to justify the corresponding mean-square bound by removing the cube factors.

Step 5: Removing the cube factors

We now pass from a mean-square bound for Th(X)T_h(X) to one for X−1/2Bh(X)X^{-1/2}B_h(X), with both averages taken over hh. Note that one cannot simply discard the terms with b≠1b\ne1, because the contributions from different bb can cancel. We begin by undoing the addition of cube factors using Möbius inversion. Related completions appear in Patterson [38], Theorem 6.1 and Heath-Brown [19], §3, with explicit removal of the cube factors by Möbius inversion in Dunn and Radziwiłł [7], Proposition 5.3 and (8.2).

Fix hh and ξ\xi, and write Ph(X)=X−1/2Bh(X)P_h(X)=X^{-1/2}B_h(X) for the b=1b=1 part of Th(X)T_h(X). Substituting the coefficient formula (2.10) into (2.12) gives

Th(X)=∑b∈Ob≡1 (3)wh(b)NK/Q(b) Ph ⁣(XNK/Q(b)3),T_h(X)=\sum_{\substack{b\in\mathcal O\\b\equiv1\ (3)}}\frac{w_h(b)}{\mathrm N_{K/\mathbb Q}(b)}\, P_h\!\Bigl(\frac{X}{\mathrm N_{K/\mathbb Q}(b)^3}\Bigr),

where

wh(b)=α(b)‾3ξ(b)3χb(h)3.w_h(b)=\overline{\alpha(b)}^3\xi(b)^3\chi_b(h)^3.

Each factor in (2.17) is completely multiplicative in bb. Thus wh(bc)=wh(b)wh(c)w_h(bc)=w_h(b)w_h(c) even when b,cb,c share prime factors, and ∣wh(b)∣≤1|w_h(b)|\le1. Möbius inversion (cf. [22] §1.3) gives

Ph(X)=∑d∈Od≡1 (3)μ(d)wh(d)NK/Q(d) Th ⁣(XNK/Q(d)3).P_h(X)=\sum_{\substack{d\in\mathcal O\\d\equiv1\ (3)}}\frac{\mu(d)w_h(d)}{\mathrm N_{K/\mathbb Q}(d)}\, T_h\!\Bigl(\frac{X}{\mathrm N_{K/\mathbb Q}(d)^3}\Bigr).

Indeed, substituting (2.16) into (2.18) and grouping by the total cube index bb gives the factor wh(b)∑d∣bμ(d)w_h(b)\sum_{d\mid b}\mu(d), which is 11 for b=1b=1 and 00 otherwise.

The objective is now an estimate for PhP_h, rather than for the completed sum ThT_h. For the simplified family and the parameter ranges arising from Step 3, the required bound is

∑0<NK/Q(h)≪H∣Ph(X)∣2≪ε,ξ,U(HX)εX.\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|P_h(X)|^2 \ll_{\varepsilon,\xi,U}(\mathcal H X)^\varepsilon X.

Since Bh(X)=XPh(X)B_h(X)=\sqrt{X}P_h(X), this is equivalent to a bound of size X2X^2, up to the same small power, for the mean square of Bh(X)B_h(X).

For 1≤Hc≤X1/31\le H_c\le X^{1/3}, let Ph,≤Hc(X)P_{h,\le H_c}(X) be the part of (2.18) with NK/Q(d)≤Hc\mathrm N_{K/\mathbb Q}(d)\le H_c. Applying weighted Cauchy–Schwarz for each fixed hh, then summing over hh and using (2.15), gives

∑0<NK/Q(h)≪H∣Ph,≤Hc(X)∣2≤(∑NK/Q(e)≤Hc1NK/Q(e))(∑NK/Q(d)≤Hc1NK/Q(d)∑0<NK/Q(h)≪H∣Th ⁣(X/NK/Q(d)3)∣2)≼∑NK/Q(d)≤Hc1NK/Q(d)(H+H2NK/Q(d)3X)≼H+H2Hc3X.\begin{aligned} &\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|P_{h,\le H_c}(X)|^2\\ &\quad\le\Bigl(\sum_{\mathrm N_{K/\mathbb Q}(e)\le H_c}\frac1{\mathrm N_{K/\mathbb Q}(e)}\Bigr) \Bigl(\sum_{\mathrm N_{K/\mathbb Q}(d)\le H_c}\frac1{\mathrm N_{K/\mathbb Q}(d)} \sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H} |T_h\!\bigl(X/\mathrm N_{K/\mathbb Q}(d)^3\bigr)|^2\Bigr)\\ &\quad\preccurlyeq\sum_{\mathrm N_{K/\mathbb Q}(d)\le H_c}\frac1{\mathrm N_{K/\mathbb Q}(d)} \Bigl(\mathcal H+\frac{\mathcal H^2\mathrm N_{K/\mathbb Q}(d)^3}{X}\Bigr) \preccurlyeq\mathcal H+\frac{\mathcal H^2H_c^3}{X}. \end{aligned}

Thus

∑0<NK/Q(h)≪H∣Ph,≤Hc(X)∣2≪ε,ξ,U(HX)ε(H+H2Hc3X).\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|P_{h,\le H_c}(X)|^2 \ll_{\varepsilon,\xi,U}(\mathcal H X)^\varepsilon \Bigl(\mathcal H+\frac{\mathcal H^2H_c^3}{X}\Bigr).

For H≤X\mathcal{H}\le X, the estimate (2.20) is within the target (2.19) provided H2Hc3/X≤X\mathcal{H}^2H_c^3/X\le X. Thus we apply the completed mean-square bound directly only up to the cutoff

Hc=min⁡ ⁣{X1/3,(XH)2/3}.H_c=\min\!\biggl\{X^{1/3}, \Bigl(\frac{X}{\mathcal H}\Bigr)^{2/3}\biggr\}.

At the basic scales from Step 3,

X≍D,H≍D1−ϑ,Hc≍D2ϑ/3.X\asymp D,\qquad\mathcal{H}\asymp D^{1-\vartheta},\qquad H_c\asymp D^{2\vartheta/3}.

The inverse sum can extend to NK/Q(d)≍D1/3\mathrm N_{K/\mathbb Q}(d)\asymp D^{1/3}, so we must still control the larger divisors.

Write τdiv(b)\tau_{\mathrm{div}}(b) for the number of nonzero integral ideal divisors of (b)(b). Let Ph,>Hc(X)P_{h,>H_c}(X) denote the terms with NK/Q(d)>Hc\mathrm N_{K/\mathbb Q}(d)>H_c in (2.18), and put Lb=X/NK/Q(b)3L_b=X/\mathrm N_{K/\mathbb Q}(b)^3. Substituting (2.16) into the truncated inversion formula for Ph,>Hc(X)P_{h,>H_c}(X) obtained from (2.18), and grouping by b=dcb=dc, gives

Ph,>Hc(X)=∑NK/Q(d)>Hc∑cμ(d)wh(d)wh(c)NK/Q(d)NK/Q(c)Ph ⁣(XNK/Q(dc)3)=∑NK/Q(b)>Hcβ0(b)χb(h)3NK/Q(b)Ph(Lb),\begin{aligned} P_{h,>H_c}(X) &=\sum_{\mathrm N_{K/\mathbb Q}(d)>H_c}\sum_c \frac{\mu(d)w_h(d)w_h(c)}{\mathrm N_{K/\mathbb Q}(d)\mathrm N_{K/\mathbb Q}(c)} P_h\!\Bigl(\frac{X}{\mathrm N_{K/\mathbb Q}(dc)^3}\Bigr)\\ &=\sum_{\mathrm N_{K/\mathbb Q}(b)>H_c} \frac{\beta_0(b)\chi_b(h)^3}{\mathrm N_{K/\mathbb Q}(b)}P_h(L_b), \end{aligned}

where all indices are primary and

β0(b)=α(b)‾ 3ξ(b)3∑d∣bNK/Q(d)>Hcμ(d),∣β0(b)∣≤τdiv(b).\beta_0(b)=\overline{\alpha(b)}^{\,3}\xi(b)^3 \sum_{\substack{d\mid b\\\mathrm N_{K/\mathbb Q}(d)>H_c}}\mu(d), \qquad |\beta_0(b)|\le\tau_{\mathrm{div}}(b).

The support of UU restricts NK/Q(b)≪UX1/3\mathrm N_{K/\mathbb Q}(b)\ll_U X^{1/3}. Weighted Cauchy–Schwarz for each hh, followed by summation, gives

∑0<NK/Q(h)≪H∣Ph,>Hc(X)∣2≤(∑bτdiv(b)NK/Q(b))∑bτdiv(b)NK/Q(b)∑0<NK/Q(h)≪H∣Ph(Lb)∣2≼sup⁡b∑0<NK/Q(h)≪H∣Ph(Lb)∣2.\begin{aligned} \sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|P_{h,>H_c}(X)|^2 &\le\Bigl(\sum_b\frac{\tau_{\mathrm{div}}(b)}{\mathrm N_{K/\mathbb Q}(b)}\Bigr) \sum_b\frac{\tau_{\mathrm{div}}(b)}{\mathrm N_{K/\mathbb Q}(b)} \sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|P_h(L_b)|^2\\ &\preccurlyeq\sup_b\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|P_h(L_b)|^2. \end{aligned}

Put aξ(n)=α(n)‾γ2(n)ξ(n)a_\xi(n)=\overline{\alpha(n)}\gamma_2(n)\xi(n). Write

E(H,X):=1X∑0<NK/Q(h)≪H∣∑n∗aξ(n)χn(h)U(NK/Q(n)/X)∣2.E(\mathcal H,X) :=\frac1X\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H} \Bigl|\sum_n^*a_\xi(n)\chi_n(h)U(\mathrm N_{K/\mathbb Q}(n)/X)\Bigr|^2.

A star restricts an index to squarefree primary elements. By (2.7) and the definition of PhP_h, we have

E(H,X)=1X∑0<NK/Q(h)≪H∣Bh(X)∣2=∑0<NK/Q(h)≪H∣Ph(X)∣2.E(\mathcal H,X) =\frac1X\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|B_h(X)|^2 =\sum_{0<\mathrm N_{K/\mathbb Q}(h)\ll\mathcal H}|P_h(X)|^2.

At the scales from Step 3, our goal is to prove

E(H,X)≼X.E(\mathcal H,X)\preccurlyeq X.

At X=DX=D, this is precisely the bound (2.9) for the mean square of Bh(D)B_h(D). Substitution into the Poisson comparison (2.8) then gives the required original mean-square estimate (2.3). Combining (2.20) for Ph,≤HcP_{h,\leq H_c} (the terms with NK/Q(d)≤Hc\mathrm N_{K/\mathbb Q}(d)\le H_c) and (2.22) for Ph,>HcP_{h,>H_c} (the terms with NK/Q(d)>Hc\mathrm N_{K/\mathbb Q}(d)>H_c), with the cutoff (2.21), gives

E(H,X)≼X+sup⁡b: 1≤Lb≤X/Hc3E(H,Lb).E(\mathcal H,X)\preccurlyeq X + \sup_{b:\,1\le L_b\le X/H_c^3}E(\mathcal H,L_b).

It therefore remains to prove E(H,Lb)≼XE(\mathcal H,L_b)\preccurlyeq X for 1<Lb≤X/Hc31<L_b\leq X/H_c^3. Here the row range H\mathcal{H} stays fixed, and the required bound is still of size XX even though the column scale has decreased to LbL_b.

We now expand the square and apply Poisson summation in hh as before; we now record the proof in terms of E(H,Lb)E(\mathcal{H},L_b). The Gauss-sum coefficients become Möbius coefficients, and the new row variable yy has norm at most about Lb2/HL_b^2/\mathcal{H}. After separating the weights, this gives schematically

E(H,Lb)≼X+HLb2∑0<NK/Q(y)≪Lb2/H∣M(y)∣2,E(\mathcal H,L_b) \preccurlyeq X+\frac{\mathcal H}{L_b^2} \sum_{0<\mathrm N_{K/\mathbb Q}(y)\ll L_b^2/\mathcal H}|M(y)|^2,
M(y)=∑NK/Q(n)≍Lb∗μ(n)ξ1(n)χn(y)‾U1(NK/Q(n)/Lb).M(y)=\sum_{\mathrm N_{K/\mathbb Q}(n)\asymp L_b}^* \mu(n)\xi_1(n)\overline{\chi_n(y)}U_1(\mathrm N_{K/\mathbb Q}(n)/L_b).

Here ξ1\xi_1 is another fixed ray class character and U1U_1 is a smooth compactly supported weight produced by separating the variables. The term XX includes the diagonal and zero-frequency contributions, using H≤X\mathcal{H}\le X.

Since Lb=X/NK/Q(b)3≤XL_b=X/\mathrm N_{K/\mathbb Q}(b)^3\le X, we may enlarge the nonnegative sum over yy to

NK/Q(y)≪Y,Y=XLbH≥Lb2H.\mathrm N_{K/\mathbb Q}(y)\ll Y,\qquad Y=\frac{XL_b}{\mathcal H}\ge\frac{L_b^2}{\mathcal H}.

The purpose of this enlargement is that a second application of Poisson summation gives a shorter row range:

H′=Lb2Y=HLbX.\mathcal{H}'=\frac{L_b^2}{Y}=\frac{\mathcal{H}L_b}{X}.

The coefficients return to cubic Gauss-sum coefficients, and

∑0<NK/Q(y)≪Y∣M(y)∣2≼YLb+YE′(H′,Lb).\sum_{0<\mathrm N_{K/\mathbb Q}(y)\ll Y}|M(y)|^2 \preccurlyeq YL_b+Y E'(\mathcal H',L_b).

Here E′E' has the same form as EE, with possibly different smooth weights and fixed characters; YLbYL_b accounts for the zero-frequency contribution. Consequently,

E(H,Lb)≼X+XLbE′ ⁣(HLbX,Lb),E(H,Lb)X≼1+E′(HLb/X,Lb)Lb.E(\mathcal H,L_b) \preccurlyeq X+\frac X{L_b} E'\!\Bigl(\frac{\mathcal HL_b}{X},L_b\Bigr), \qquad \frac{E(\mathcal H,L_b)}X \preccurlyeq 1+\frac{E'(\mathcal HL_b/X,L_b)}{L_b}.

Thus it suffices to prove

E′ ⁣(HLbX,Lb)≼Lb.E'\!\Bigl(\frac{\mathcal HL_b}{X},L_b\Bigr)\preccurlyeq L_b.

This is the original type of estimate at smaller parameters:

(H′,X′)=(HNK/Q(b)3,XNK/Q(b)3),H′X′=HX.(\mathcal H',X') =\Bigl(\frac{\mathcal H}{\mathrm N_{K/\mathbb Q}(b)^3}, \frac{X}{\mathrm N_{K/\mathbb Q}(b)^3}\Bigr), \qquad \frac{\mathcal H'}{X'}=\frac{\mathcal H}{X}.

Both scales decrease while their ratio stays fixed. The factor X/LbX/L_b in the preceding inequality is exactly what converts the new target bound LbL_b into the required bound XX. These two Poisson summations give the transfer estimate: a bound for the remaining mean square in terms of new mean squares of the same type. For related uses of an enlarged summation range, see Goldmakher–Louvel [13], Lemma 4.4 and the proof of Theorem 4.1, following Heath-Brown [18], Lemma 9.

Combining the bounds (2.20) and (2.22) from the first part of Step 5 with the transfer estimate above gives the recursive bound

E(H,X)X≼1+sup⁡NK/Q(b)>HcLb>1E′(H′,X′)X′,X′=Lb,H′=HLbX.\frac{E(\mathcal H,X)}{X} \preccurlyeq 1+\sup_{\substack{\mathrm N_{K/\mathbb Q}(b)>H_c\\L_b>1}}\frac{E'(\mathcal H',X')}{X'}, \qquad X'=L_b,\quad \mathcal H'=\frac{\mathcal H L_b}{X}.

For each bb in this supremum, the cutoff (2.21) gives H′<H(H/X)2≪D−2ϑH\mathcal{H}'<\mathcal{H}(\mathcal{H}/X)^2\ll D^{-2\vartheta}\mathcal{H}. Since H′/X′=H/X\mathcal{H}'/X'=\mathcal{H}/X, the same contraction applies at every step. After Oϑ(1)O_\vartheta(1) steps, every resulting mean square either has an empty remainder or has row parameter at most 1, where counting gives the desired bound at its reduced scales. Applying the recursive inequality back through these steps proves (2.23) for the original E(H,X)E(\mathcal H,X), completing the sketch of the proof.

The preceding sketch suppresses auxiliary twists and common factors for the purpose of illustration. To carry out this argument with the auxiliary twists included, we use the family

E(H,X,F)=1XF∑NK/Q(f)≍F∗∑0<NK/Q(k)≪H∣∑n∗aξ(n)χn(k)χn(f)4U(NK/Q(n)/X)∣2\mathcal E(\mathcal H,X,F) =\frac1{XF}\sum_{\mathrm N_{K/\mathbb Q}(f)\asymp F}^* \sum_{0<\mathrm N_{K/\mathbb Q}(k)\ll\mathcal H} \Bigl|\sum_n^*a_\xi(n)\chi_n(k)\chi_n(f)^4 U(\mathrm N_{K/\mathbb Q}(n)/X)\Bigr|^2

and prove E(H,X,F)≼XF\mathcal E(\mathcal H,X,F)\preccurlyeq XF in the parameter ranges of Proposition 5.1. Unlike the sketch above, the full argument must also handle the common factors and the resulting dyadic ranges. The precise family is defined in (4.9), and Proposition 5.4 states the transfer estimate. Combining it with the bounds for the two parts of the inverse sum in Step 5 reduces the row range at each step. Section 5 proves the desired bound by a finite iteration, following the admissible-exponent method of Heath-Brown [18], Lemma 8 and §8, [19], Lemma 9; see also [13], Theorem 4.1 and [5], §3.2.

From the mean-square estimate to the zero-free region

We first carry out Steps 1 and 2 of the outline: extract cancellation in A1(D)A_1(D) from a mean-square estimate for the family, then use a Mellin transform to deduce nonvanishing. Fix a finite-order Hecke character ν\nu of KK and a finite set SS of prime ideals that contains all prime ideals above 2 or 3, as well as all prime ideals dividing the conductor of ν\nu. For an ideal or element aa, write (a,S)=1(a,S)=1 if no prime ideal in SS divides aa. An ideal is supported on SS if all its prime factors belong to SS. For W∈Cc∞((0,∞);C)W\in C_c^\infty((0,\infty);\mathbb C), recall the family

Au(D):=∑(n,S)=1μ(n)ν(n)χn(u)W(NK/Q(n)/D),A_u(D):=\sum_{(n,S)=1}\mu(n)\nu(n)\chi_n(u)W(\mathrm N_{K/\mathbb Q}(n)/D),

introduced in (2.2). Here and throughout the original family, nn runs over ideals prime to SS, represented by their primary generators. The key estimate is the following.

Proposition 3.1. For every fixed 0<ϑ≤1/100<\vartheta\le1/10 and ε>0\varepsilon>0, there exists an integer k=k(ϑ,ε)≥1k=k(\vartheta,\varepsilon)\ge1 such that, for every compact interval I⊂(0,∞)I\subset(0,\infty), all smooth WW supported in II, and D≥2D\ge2,

∑0<NK/Q(u)≤D1+ϑ∣Au(D)∣2≪ν,S,I,ϑ,ε(max⁡0≤j≤k∥W(j)∥∞)2D2+ϑ+ε.\sum_{0<\mathrm N_{K/\mathbb Q}(u)\le D^{1+\vartheta}}|A_u(D)|^2 \ll_{\nu,S,I,\vartheta,\varepsilon} \Bigl(\max_{0\le j\le k}\|W^{(j)}\|_\infty\Bigr)^2 D^{2+\vartheta+\varepsilon}.

We first deduce Theorem 1.1 assuming Proposition 3.1. The proof of the proposition is completed in Section 5.6.

Proof of Theorem 1.1 from Proposition 3.1. Fix 0<ϑ≤1/100<\vartheta\le1/10, and put H=D1+ϑH=D^{1+\vartheta} and Y=H1/6=D(1+ϑ)/6Y=H^{1/6}=D^{(1+\vartheta)/6}. For prime ideals Y/2<NK/Q(p)≤YY/2<\mathrm N_{K/\mathbb Q}(\mathfrak p)\le Y, p∉S\mathfrak{p}\notin S, let pp be their primary generators, chosen with p≡1(mod3)p\equiv1\pmod3. Then χn(p6)=1p∤n\chi_n(p^6)=\mathbf1_{\mathfrak p\nmid n} and therefore

∣A1(D)−Ap6(D)∣≤∥W∥∞#{n:(n,S)=1, p∣n, NK/Q(n)≍D}≪WD/Y.|A_1(D)-A_{p^6}(D)| \le\|W\|_\infty \#\{n:(n,S)=1,\ \mathfrak p\mid n,\ \mathrm N_{K/\mathbb Q}(n)\asymp D\} \ll_W D/Y.

For this fixed field, Landau’s prime ideal theorem [27] (see [22], Theorem 5.33) gives J≍Y/log⁡YJ\asymp Y/\log Y such primes. Their sixth powers are distinct rows of norm at most HH. Apply (3.1) with loss ε/2\varepsilon/2 and use log⁡Y≪εDε/2\log Y\ll_\varepsilon D^{\varepsilon/2} to obtain

∣A1(D)∣2≤2J∑Y/2<NK/Q(p)≤Y(p,S)=1∣Ap6(D)∣2+OW(D2/Y2)\begin{aligned} |A_1(D)|^2 &\le \frac2J\sum_{\substack{Y/2<\mathrm N_{K/\mathbb Q}(p)\le Y\\(p,S)=1}}|A_{p^6}(D)|^2+O_W(D^2/Y^2)\end{aligned}
≪ν,S,W,ϑ,εD2+ϑ+ε/Y+D2/Y2=D11/6+5ϑ/6+ε+D5/3−ϑ/3.\ll_{\nu,S,W,\vartheta,\varepsilon}D^{2+\vartheta+\varepsilon}/Y+D^2/Y^2=D^{11/6+5\vartheta/6+\varepsilon}+D^{5/3-\vartheta/3}.

Thus A1(D)≪ν,S,W,ϑ,εD11/12+5ϑ/12+εA_1(D) \ll_{\nu,S,W,\vartheta,\varepsilon} D^{11/12+5\vartheta/12+\varepsilon}. Given any requested exponent loss, choose ϑ>0\vartheta>0 and then the loss in Proposition 3.1 sufficiently small. Renaming the resulting loss ε\varepsilon, we obtain

A1(D)≪ν,S,W,εD11/12+ε.A_1(D) \ll_{\nu,S,W,\varepsilon} D^{11/12+\varepsilon}.

We now use (3.3) to rule out a zero of LK(s,ν)L_K(s,\nu) in Re⁡s>11/12\operatorname{Re}s>11/12. Suppose such a zero ϱ\varrho exists. Choose 0≠ϕ∈Cc∞((1,2))0\ne\phi\in C_c^\infty((1,2)), ϕ≥0\phi\ge0, and W(y)=y−ϱϕ(y)W(y)=y^{-\varrho}\phi(y). With the Mellin convention W^(s)=∫0∞W(y)ys dy/y\widehat W(s)=\int_0^\infty W(y)y^s\,d y/y, we have W^(ϱ)=∫ϕ(y) dy/y>0\widehat W(\varrho)=\int\phi(y)\,d y/y>0. The function

MW(s)=∫0∞A1(D)D−s dDD\mathcal M_W(s)=\int_0^\infty A_1(D)D^{-s}\frac{\,d D}{D}

is holomorphic on Re⁡s>11/12\operatorname{Re}s>11/12: the estimate (3.3) controls the integral at infinity, and the compact support of WW makes A1(D)A_1(D) vanish for sufficiently small DD. Termwise integration for Re⁡s>1\operatorname{Re}s>1 gives

MW(s)=W^(s)∑(n,S)=1μ(n)ν(n)NK/Q(n)s=W^(s)LKS(s,ν),\mathcal M_W(s) =\widehat W(s)\sum_{(n,S)=1}\frac{\mu(n)\nu(n)}{\mathrm N_{K/\mathbb Q}(n)^s} =\frac{\widehat W(s)}{L_K^S(s,\nu)},

where LKSL_K^S is the Euler product outside SS:

LKS(s,ν):=LK(s,ν)∏p∈S(1−ν(p)NK/Q(p)−s).L_K^S(s,\nu):=L_K(s,\nu) \prod_{\mathfrak p\in S}(1-\nu(\mathfrak p)\mathrm N_{K/\mathbb Q}(\mathfrak p)^{-s}).

By Hecke’s meromorphic continuation theorem [21] (see [22], §5.10) and the identity theorem, the identity LKS(s,ν)MW(s)=W^(s)L_K^S(s,\nu)\mathcal M_W(s)=\widehat W(s) holds on Re⁡s>11/12\operatorname{Re}s>11/12. The omitted Euler factors are nonzero here, so evaluation at s=ϱs=\varrho gives 0=W^(ϱ)>00=\widehat{W}(\varrho)>0, a contradiction.

To deduce the assertion for Dirichlet LL-functions in Theorem 1.1, let χ−3\chi_{-3} be the nontrivial character modulo 3. For any Dirichlet character χ\chi, quadratic base change gives, up to Euler factors nonzero in Re⁡s>0\operatorname{Re}s>0,

LK(s,χ∘NK/Q)=L(s,χ)L(s,χχ−3).L_K(s,\chi\circ\mathrm N_{K/\mathbb Q}) =L(s,\chi)L(s,\chi\chi_{-3}).

The Hecke conclusion excludes zeros of either factor away from s=1s=1. At s=1s=1, the only possible pole–zero cancellation is ruled out by Dirichlet’s nonvanishing theorem, which gives L(1,χ−3)>0L(1,\chi_{-3})>0 (see [6], Chapters 4 and 6).

Poisson summation and the dual mean square

We now turn to the mean-square estimate in Proposition 3.1. Following Step 3 of the outline, we expand the square and apply Poisson summation in the row variable uu. The finite Fourier transforms of the characters supply Gauss sums. Arithmetic identities then combine these Gauss sums with the original Möbius coefficients to give the normalized cubic Gauss-sum coefficients that appear in the theta function.

The column indices produced by expanding the square need not be coprime. We first extract their common factor, then use Möbius inversion to separate the remaining coprimality condition. These operations introduce an auxiliary twisting index and an exclusion ideal. The comparison below removes the exclusion without changing the row range or the product of the column and auxiliary scales.

We record the arithmetic and Poisson identities first, then define this family. Proposition 4.5 states the estimate for it that suffices to prove (3.1); the rest of the section proves that implication.

The arithmetic identities

We first record the identities that convert between Möbius coefficients and normalized cubic Gauss sums. Recall that, for a squarefree primary element nn with (n,S)=1(n,S)=1,

α(n)=n∣n∣,γj(n)=1NK/Q(n)∑x mod nχn(x)je(x/n)(j∈Z),\alpha(n)=\frac{n}{|n|},\qquad \gamma_j(n)=\frac{1}{\sqrt{\mathrm N_{K/\mathbb Q}(n)}} \sum_{x\bmod n}\chi_n(x)^j e(x/n)\quad(j\in\mathbb Z),

where e(z)=exp⁡(4πi Im⁡z/3)e(z)=\exp(4\pi\mathrm i\,\operatorname{Im}z/\sqrt3), as in (2.6). Put

aξ(n)=α(n)‾γ2(n)ξ(n)(n squarefree),a_{\xi}(n)=\overline{\alpha(n)}\gamma_2(n)\xi(n)\qquad(n\text{ squarefree}),

where ξ\xi is a fixed ray class character. The character ξ\xi accounts for the original twist ν\nu and for the residue-class factors in the identities below.

By character orthogonality [22], we can absorb functions on a fixed ray class group into a finite sum of twists ξ\xi. Choose a fixed modulus, supported on SS, divisible by the conductor of ν\nu and sufficiently large that all reciprocity factors below depend only on the corresponding ray classes. We then expand these factors in characters of this finite group.

The following lemma collects standard consequences of the Gauss–Jacobi identities, quadratic Gauss-sum evaluations, and reciprocity.

Lemma 4.1. On a fixed ray class group, there are a function GG with values in {z∈C:∣z∣=1}\{z\in\mathbb{C}:|z|=1\} and a symmetric {±1}\{\pm1\}-valued bicharacter R\mathcal{R}. This means that R(a,b)=R(b,a)\mathcal{R}(a,b)=\mathcal{R}(b,a) and R\mathcal{R} is multiplicative in each argument separately:

R(aa′,b)=R(a,b)R(a′,b),\mathcal{R}(aa',b)=\mathcal{R}(a,b)\mathcal{R}(a',b),
R(a,bb′)=R(a,b)R(a,b′),\mathcal{R}(a,bb')=\mathcal{R}(a,b)\mathcal{R}(a,b'),

for all classes a,a′,b,b′a,a',b,b' in the ray class group. These functions have the following properties. Let a,b,na,b,n be primary elements prime to SS, with (a,b)=1(a,b)=1 and nn squarefree. Then

χb(a)=R(a,b)χa(b),G(ab)=G(a)G(b)R(a,b),\chi_b(a)=\mathcal R(a,b)\chi_a(b), G(ab)=G(a)G(b)\mathcal R(a,b),
γ2(n)3=μ(n)α(n),γ1(n)γ2(n)=μ(n)α(n)G(n),\gamma_2(n)^3=\mu(n)\alpha(n), \gamma_1(n)\gamma_2(n)=\mu(n)\alpha(n)G(n),
G(n)=χn(4)‾γ3(n),γ1(n)γ−1(n)=χn(−1).G(n)=\overline{\chi_n(4)}\gamma_3(n), \gamma_1(n)\gamma_{-1}(n)=\chi_n(-1).

Consequently

α(n)‾γ2(n)γ1(n)=μ(n)G(n),\overline{\alpha(n)}\gamma_2(n)\gamma_1(n)=\mu(n)G(n),
μ(n)γ−1(n)=χn(−1)G(n)−1α(n)‾γ2(n),\mu(n)\gamma_{-1}(n)=\chi_n(-1)G(n)^{-1}\overline{\alpha(n)}\gamma_2(n),
aξ(ab)=aξ(a)aξ(b)χb(a)4,a_{\xi}(ab)=a_{\xi}(a)a_{\xi}(b)\chi_b(a)^4,
χa(−1)G(a)‾G(b)R(a,b)=G(ba−1).\chi_a(-1)\overline{G(a)}G(b)\mathcal{R}(a,b)=G(ba^{-1}).

In (4.6), a,ba,b are also squarefree. The identity for G(ab)G(ab) in (4.1) extends to all classes of the fixed ray class group; the quotient in (4.7) is taken in that group.

The proof, including the dependence of GG and R\mathcal{R} on fixed ray classes, is given in Appendix A.1. Under Poisson summation, (4.5) converts the Möbius coefficients to aξ(n)a_{\xi}(n), up to fixed ray class factors, while (4.4) converts them back.

Poisson summation with excluded primes

The row sum to which we apply Poisson summation will have an additional coprimality restriction. We record the formula with that restriction included, so that its effect on the Fourier frequencies and the normalization is explicit.

For a nonzero ideal r\mathfrak{r}, write rad⁡r\operatorname{rad}\mathfrak{r} for the product of its distinct prime divisors. We use the additive character ee introduced in Step 3. In quotients and Gauss sums, use a fixed generator for each ideal, chosen primary when the ideal is prime to 33.

Lemma 4.2. Let H>0\mathcal{H} > 0, let r\mathfrak{r} be a nonzero ideal, and let Φ(NK/Q(k)/H)\Phi(\mathrm N_{K/\mathbb Q}(k)/\mathcal H) be a smooth radial Schwartz weight on the row lattice. Let χ\chi be a primitive multiplicative character of (O/m)×(\mathcal{O}/\mathfrak{m})^\times, viewed as a function on O\mathcal{O} by reduction modulo m\mathfrak{m} and extension by zero on nonunits. Write

γ(χ)=NK/Q(m)−1/2∑x mod mχ(x)e(x/m)\gamma(\chi)=\mathrm N_{K/\mathbb Q}(\mathfrak m)^{-1/2} \sum_{x\bmod\mathfrak m}\chi(x)e(x/\mathfrak m)

for its normalized Gauss sum. Since χ\chi is primitive, its modulus m\mathfrak{m} is determined by χ\chi and is suppressed in the notation. Then

∑kχ(k)1(k,r)=1Φ(NK/Q(k)/H)=Hγ(χ)NK/Q(m)∑d∣rad⁡rμ(d)χ(d)NK/Q(d)∑hχ(h)‾Φ^ ⁣(HNK/Q(h)NK/Q(d)NK/Q(m)).\begin{aligned} &\sum_k\chi(k)\mathbf1_{(k,\mathfrak r)=1}\Phi(\mathrm N_{K/\mathbb Q}(k)/\mathcal{H})\\ &\quad=\frac{\mathcal{H}\gamma(\chi)}{\sqrt{\mathrm N_{K/\mathbb Q}(\mathfrak m)}} \sum_{d\mid\operatorname{rad} \mathfrak r}\frac{\mu(d)\chi(d)}{\mathrm N_{K/\mathbb Q}(d)} \sum_h\overline{\chi(h)} \widehat\Phi\!\Bigl(\frac{\mathcal{H}\mathrm N_{K/\mathbb Q}(h)}{\mathrm N_{K/\mathbb Q}(d)\mathrm N_{K/\mathbb Q}(\mathfrak m)}\Bigr). \end{aligned}

Here Φ^\widehat{\Phi} is defined by

Φ^(∣w∣2)=23∫CΦ(∣z∣2)e(−zw) dx dy,z=x+iy.\widehat\Phi(|w|^2)=\frac2{\sqrt3}\int_{\mathbb C} \Phi(|z|^2)e(-zw)\,dx\,dy,\qquad z=x+\mathrm i y.

The measure is normalized so that O\mathcal{O} has covolume one. Here k,h∈Ok,h\in\mathcal{O}, and hh is the Fourier frequency. We have

χ nonprincipal  ⟹  χ(0)‾=0.\chi\ \text{nonprincipal} \implies\overline{\chi(0)}=0.

For the principal primitive character (χ=1,m=1)(\chi=\mathbf1,\mathfrak m=1), the zero-frequency contribution is

HΦ^(0)∏p∣r(1−1NK/Q(p)).\mathcal{H}\widehat{\Phi}(0)\prod_{\mathfrak{p}\mid\mathfrak{r}}\left(1-\frac{1}{N_{K/\mathbb{Q}}(\mathfrak{p})}\right).

Proof. Inclusion–exclusion followed by k=dℓk=d\ell gives

∑kχ(k)1(k,r)=1Φ(NK/Q(k)/H)=∑d∣rad⁡rμ(d)χ(d)∑ℓχ(ℓ)Φ ⁣(NK/Q(ℓ)H/NK/Q(d)).\sum_k\chi(k)\mathbf1_{(k,\mathfrak r)=1}\Phi(\mathrm N_{K/\mathbb Q}(k)/\mathcal{H}) =\sum_{d\mid\operatorname{rad} \mathfrak r}\mu(d)\chi(d) \sum_\ell\chi(\ell)\Phi\!\Bigl(\frac{\mathrm N_{K/\mathbb Q}(\ell)}{\mathcal{H}/\mathrm N_{K/\mathbb Q}(d)}\Bigr).

Apply lattice Poisson summation [22] in ℓ\ell on residue classes modulo m\mathfrak{m}, at scale H/NK/Q(d)\mathcal H/\mathrm N_{K/\mathbb Q}(d). The primitive Gauss-sum identity [22] evaluates the finite Fourier transform as NK/Q(m)γ(χ)χ(h)‾\sqrt{\mathrm N_{K/\mathbb Q}(\mathfrak m)}\gamma(\chi)\overline{\chi(h)} for every hh, giving (4.8). For a nonprincipal primitive character, χ(0)‾=0\overline{\chi(0)}=0, so the zero-frequency term vanishes. For the principal character, sum μ(d)/NK/Q(d)\mu(d)/\mathrm N_{K/\mathbb Q}(d) over d∣rad⁡rd\mid\operatorname{rad}\mathfrak{r} to obtain the displayed zero-frequency contribution.

The dual mean squares

We use the following mean square of the column sums with coefficients aξ(n)a_\xi(n).

Definition 4.3. The parameters H\mathcal{H}, XX, and FF are the norm scales of kk, nn, and ff, respectively. Here kk ranges over elements of O\mathcal{O}, and a star restricts a sum to squarefree ideals prime to SS, represented by their primary generators. For a smooth compactly supported weight WW on (0,∞)(0,\infty), define the dual mean square by

E(H,X,F;ξ,W)=1XF∑F≤NK/Q(f)<2F(f,S)=1∗∑0<NK/Q(k)≤H∣∑(n,S)=1∗aξ(n)χn(k)χn(f)4W(NK/Q(n)/X)∣2.\begin{aligned} &\mathcal E(\mathcal{H},X,F;\xi,W)\\ &\qquad=\frac1{XF} \sum_{\substack{F\le\mathrm N_{K/\mathbb Q}(f)<2F\\(f,S)=1}}^*\sum_{0<\mathrm N_{K/\mathbb Q}(k)\le \mathcal{H}} \Bigl|\sum_{(n,S)=1}^*a_\xi(n)\chi_n(k)\chi_n(f)^4 W(\mathrm N_{K/\mathbb Q}(n)/X)\Bigr|^2. \end{aligned}

Below, ξ\xi ranges over all characters of the fixed ray class group chosen above, through which ν\nu, GG, and R\mathcal{R} factor.

Write Er\mathcal{E}_{\mathfrak r} for the same expression with the additional restriction (n,r)=1(n,\mathfrak r)=1. This notation is only needed in the Poisson reductions; the following comparison returns to E\mathcal{E}.

Lemma 4.4. Let r0r_0 be the product of the primes dividing r\mathfrak r outside SS. For H,X>0\mathcal{H},X>0, F≥1F\geq1, and smooth compactly supported WW,

Er(H,X,F;ξ,W)≤τdiv(r0)∑d∣r0E ⁣(H,XNK/Q(d),FNK/Q(d);ξ,W).\mathcal E_{\mathfrak r}(\mathcal H,X,F;\xi,W) \le\tau_{\mathrm{div}}(r_0) \sum_{d\mid r_0}\mathcal E\!\Bigl(\mathcal H,\frac X{\mathrm N_{K/\mathbb Q}(d)},F\mathrm N_{K/\mathbb Q}(d);\xi,W\Bigr).

Proof. Let Cr(X;k,f)C_{\mathfrak r}(X;k,f) denote the inner column sum defining Er\mathcal{E}_{\mathfrak r}, with ξ,W\xi,W fixed. Inclusion–exclusion, n=dmn=dm, and (4.6) give

Cr(X;k,f)=∑d∣r0μ(d)∑(n,S)=1d∣n∗aξ(n)χn(k)χn(f)4W(NK/Q(n)/X)=∑d∣r0(d,f)=1μ(d)aξ(d)χd(k)χd(f)4C1(X/NK/Q(d);k,df).\begin{aligned} C_{\mathfrak r}(X;k,f) &=\sum_{d\mid r_0}\mu(d) \sum_{\substack{(n,S)=1\\d\mid n}}^* a_\xi(n)\chi_n(k)\chi_n(f)^4W(\mathrm N_{K/\mathbb Q}(n)/X)\\ &=\sum_{\substack{d\mid r_0\\(d,f)=1}} \mu(d)a_\xi(d)\chi_d(k)\chi_d(f)^4 C_1(X/\mathrm N_{K/\mathbb Q}(d);k,df). \end{aligned}

Indeed, χm(d)4\chi_m(d)^4 enforces (m,d)=1(m,d)=1 and combines with χm(f)4\chi_m(f)^4; terms with (d,f)≠1(d,f)\ne1 vanish. Each exterior coefficient has modulus at most one. Apply Cauchy–Schwarz in dd, then enlarge the injective image f↦dff\mapsto df to the squarefree range FNK/Q(d)≤NK/Q(df)<2FNK/Q(d)F\mathrm N_{K/\mathbb Q}(d)\le\mathrm N_{K/\mathbb Q}(df)<2F\mathrm N_{K/\mathbb Q}(d). The normalizing product is unchanged: (X/NK/Q(d))(FNK/Q(d))=XF(X/\mathrm N_{K/\mathbb Q}(d))(F\mathrm N_{K/\mathbb Q}(d))=XF. □

Reduction to the dual mean square

The following proposition gives the dual estimate sufficient for (3.1). Its proof applies Poisson summation directly to the original Möbius sums.

Proposition 4.5. Fix 0<ϑ≤1/100<\vartheta\leq1/10 and put H=D1+ϑH=D^{1+\vartheta}. For each fixed C≥1C\geq1 and all real BB, F≥1F\geq1, consider the ranges

X=DBF,0<H≤CD2HB2,XF=DB.X=\frac D{BF},\qquad 0<\mathcal{H}\le\frac{CD^2}{HB^2},\qquad XF=\frac DB.

Suppose that for every ε>0\varepsilon>0 there exists an integer J=J(ϑ,ε)≥1J=J(\vartheta,\varepsilon)\geq1 such that, for every compact interval I⊂(0,∞)I\subset(0,\infty) and every smooth WW supported in II,

E(H,X,F;ξ,W)≪ν,S,I,C,ϑ,ε∥W∥CJ(I)2DεXF,\mathcal E(\mathcal H,X,F;\xi,W) \ll_{\nu,S,I,C,\vartheta,\varepsilon}\|W\|_{C^J(I)}^2D^\varepsilon XF,

where

∥W∥CJ(I)=max⁡0≤j≤Jsup⁡x∈I∣W(j)(x)∣.\|W\|_{C^J(I)}=\max_{0\le j\le J}\sup_{x\in I}|W^{(j)}(x)|.

For functions of several variables, the norm uses all partial derivatives of total order at most JJ. Then the original mean-square estimate in Proposition 3.1 holds.

Proof. Write I=[a0,b0]I=[a_0,b_0]. Choose a nonnegative radial Schwartz weight Φ\Phi such that Φ(x)≥1\Phi(x)\geq1 on [0,1][0,1] and supp⁡Φ^⊂[0,CΦ]\operatorname{supp}\widehat{\Phi}\subset[0,C_{\Phi}] for some fixed CΦ>0C_{\Phi}>0. It suffices to prove MD≪ν,S,I,ϑ,ε∥W∥CJ′(I)2HDε\mathcal M_D\ll_{\nu,S,I,\vartheta,\varepsilon} \|W\|_{C^{J'}(I)}^2HD^\varepsilon for some J′=J′(ϑ,ε)J'=J'(\vartheta,\varepsilon), where

MD:=1D∑u∈OΦ(NK/Q(u)/H)∣∑(n,S)=1μ(n)ν(n)χn(u)W(NK/Q(n)/D)∣2.\mathcal M_D:=\frac1D\sum_{u\in\mathcal O}\Phi(\mathrm N_{K/\mathbb Q}(u)/H) \Bigl|\sum_{(n,S)=1}\mu(n)\nu(n)\chi_n(u)W(\mathrm N_{K/\mathbb Q}(n)/D)\Bigr|^2.

Conjugate the inner sum, expand the square, and put g=(n1,n2)g=(n_1,n_2) and nj=gzjn_j=gz_j. The squarefree indices satisfy (z1,z2)=(z1z2,g)=1(z_1,z_2)=(z_1z_2,g)=1, and

χn1(u)‾χn2(u)=1(u,g)=1χz1(u)‾χz2(u).\overline{\chi_{n_1}(u)}\chi_{n_2}(u) =\mathbf1_{(u,g)=1}\overline{\chi_{z_1}(u)}\chi_{z_2}(u).

Apply Lemma 4.2 with χ=χ‾z1χz2\chi=\overline{\chi}_{z_1}\chi_{z_2} and exclusion gg:

∑uΦ(NK/Q(u)/H)1(u,g)=1χ(u)=∑e∣gHμ(e)χ(e)γ(χ)NK/Q(e)NK/Q(z1)NK/Q(z2)∑hχ(h)‾Φ^ ⁣(HNK/Q(h)NK/Q(e)NK/Q(z1)NK/Q(z2)).\begin{aligned} &\sum_u\Phi(\mathrm N_{K/\mathbb Q}(u)/H)\mathbf1_{(u,g)=1}\chi(u)\\ &\quad=\sum_{e\mid g} \frac{H\mu(e)\chi(e)\gamma(\chi)} {\mathrm N_{K/\mathbb Q}(e)\sqrt{\mathrm N_{K/\mathbb Q}(z_1)\mathrm N_{K/\mathbb Q}(z_2)}} \sum_h\overline{\chi(h)} \widehat\Phi\!\Bigl(\frac{H\mathrm N_{K/\mathbb Q}(h)}{\mathrm N_{K/\mathbb Q}(e)\mathrm N_{K/\mathbb Q}(z_1)\mathrm N_{K/\mathbb Q}(z_2)}\Bigr). \end{aligned}

The zero frequency occurs only when z1=z2=1z_1=z_2=1 and contributes

Z=HDΦ^(0)∑(g,S)=1∗∣W(NK/Q(g)/D)∣2∏p∣g(1−1NK/Q(p))≪Φ,IH∥W∥∞2.Z=\frac HD\widehat\Phi(0) \sum_{(g,S)=1}^*|W(\mathrm N_{K/\mathbb Q}(g)/D)|^2 \prod_{p\mid g}\Bigl(1-\frac1{\mathrm N_{K/\mathbb Q}(p)}\Bigr) \ll_{\Phi,I}H\|W\|_\infty^2.

Retain all h≠0h\ne0, including those with z1=z2=1z_1=z_2=1.

Expand the fixed ray class function ν(t)‾G(t−1)=∑ξcξξ(t)\overline{\nu(t)}G(t^{-1})=\sum_\xi c_\xi\xi(t). The Chinese remainder theorem and (4.5)–(4.7) give

μ(z1)μ(z2)ν‾(z1)ν(z2)γ(χ‾z1χz2)=∑ξcξaξ(z1)aξ‾(z2).\mu(z_1)\mu(z_2)\overline{\nu}(z_1)\nu(z_2)\gamma(\overline{\chi}_{z_1}\chi_{z_2})=\sum_\xi c_\xi a_\xi(z_1)\overline{a_\xi}(z_2).

Write N=NK/QN=\mathrm N_{K/\mathbb Q} for the remainder of this proof, and put W0(x)=x−1/2W‾(x)W_0(x)=x^{-1/2}\overline{W}(x). Substituting (4.15) into (4.14), including the normalization 1/D1/D in (4.13), gives

MD−Z=∑ξcξSξ,∣MD−Z∣≪max⁡ξ∣Sξ∣,\mathcal M_D-Z=\sum_\xi c_\xi\mathcal S_\xi, \qquad |\mathcal M_D-Z|\ll\max_\xi|\mathcal S_\xi|,

since the character sum is fixed and finite. Fix ξ\xi. Using χz(e)‾=χz(e5)\overline{\chi_z(e)}=\chi_z(e^5), its contribution is

Sξ=HD2∑g∗∑e∣gμ(e)N(g)N(e)∑h≠0∑z1,z2(z1,z2)=1(z1z2,g)=1∗aξ(z1)aξ(z2)‾χz1(he5)χz2(he5)‾×W0 ⁣(N(gz1)D)W0 ⁣(N(gz2)D)‾Φ^ ⁣(HN(h)N(e)N(z1)N(z2)).\begin{aligned} \mathcal S_\xi ={}&\frac H{D^2} \sum_g^*\sum_{e\mid g}\frac{\mu(e)N(g)}{N(e)} \sum_{h\ne0} \sum_{\substack{z_1,z_2\\(z_1,z_2)=1\\(z_1z_2,g)=1}}^* a_\xi(z_1)\overline{a_\xi(z_2)} \chi_{z_1}(he^5)\overline{\chi_{z_2}(he^5)} \\ &\quad{}\times W_0\!\Bigl(\frac{N(gz_1)}D\Bigr) \overline{W_0\!\Bigl(\frac{N(gz_2)}D\Bigr)} \widehat\Phi\!\Bigl( \frac{HN(h)}{N(e)N(z_1)N(z_2)} \Bigr). \end{aligned}

Here and below every starred variable is squarefree, primary, and prime to SS, while h,kh,k range over nonzero elements of O\mathcal O. We must show ∣Sξ∣≪ν,S,I,ϑ,ε∥W∥CJ′(I)2HDε|\mathcal S_\xi|\ll_{\nu,S,I,\vartheta,\varepsilon} \|W\|_{C^{J'}(I)}^2HD^\varepsilon.

Insert the coprimality identity and change variables:

1(z1,z2)=1=∑v∣z1, v∣z2μ(v),zj=vmj.\mathbf1_{(z_1,z_2)=1} =\sum_{v\mid z_1,\ v\mid z_2}\mu(v), \qquad z_j=vm_j.

By (4.6), the common factor from the two columns satisfies

aξ(vm)=aξ(v)aξ(m)χm(v)4,∣aξ(v)χv(he5)∣2=1(v,h)=1,a_\xi(vm)=a_\xi(v)a_\xi(m)\chi_m(v)^4, \qquad |a_\xi(v)\chi_v(he^5)|^2=\mathbf1_{(v,h)=1},

where (v,e)=1(v,e)=1. Thus the same sum becomes

Sξ=HD2∑g,v(g,v)=1∗∑e∣gμ(e)μ(v)N(g)N(e)∑h≠0(h,v)=1∑m1,m2(m1m2,gv)=1∗aξ(m1)aξ(m2)‾×χm1(he5v4)χm2(he5v4)‾W0 ⁣(N(gvm1)D)W0 ⁣(N(gvm2)D)‾×Φ^ ⁣(HN(h)N(e)N(v)2N(m1)N(m2)).\begin{aligned} \mathcal S_\xi ={}&\frac H{D^2} \sum_{\substack{g,v\\(g,v)=1}}^* \sum_{e\mid g}\frac{\mu(e)\mu(v)N(g)}{N(e)} \sum_{\substack{h\ne0\\(h,v)=1}} \sum_{\substack{m_1,m_2\\(m_1m_2,gv)=1}}^* a_\xi(m_1)\overline{a_\xi(m_2)} \\ &\quad{}\times \chi_{m_1}(he^5v^4)\overline{\chi_{m_2}(he^5v^4)} W_0\!\Bigl(\frac{N(gvm_1)}D\Bigr) \overline{W_0\!\Bigl(\frac{N(gvm_2)}D\Bigr)} \\ &\quad{}\times \widehat\Phi\!\Bigl( \frac{HN(h)}{N(e)N(v)^2N(m_1)N(m_2)} \Bigr). \end{aligned}

There is now no restriction (m1,m2)=1(m_1,m_2)=1. Make the bijective change of variables

b=g/e,f=ev,k=eh,b=g/e,\qquad f=ev,\qquad k=eh,
e=(f,k),v=f/e,g=be,h=k/e.e=(f,k),\qquad v=f/e,\qquad g=be,\qquad h=k/e.

The resulting b,fb,f are coprime and squarefree, k≠0k\ne0 is arbitrary, and

gv=bf,he5v4=kf4,μ(e)μ(v)=μ(f),N(g)N(e)=N(b).gv=bf,\qquad he^{5}v^{4}=kf^{4},\qquad\mu(e)\mu(v)=\mu(f),\qquad\frac{N(g)}{N(e)}=N(b).

Consequently,

Sξ=HD2∑b,f(b,f)=1∗μ(f)N(b)∑k≠0∑m1,m2(m1m2,b)=1∗aξ(m1)aξ(m2)‾χm1(kf4)χm2(kf4)‾\begin{aligned} \mathcal S_\xi ={}&\frac H{D^2} \sum_{\substack{b,f\\(b,f)=1}}^*\mu(f)N(b) \sum_{k\ne0} \sum_{\substack{m_1,m_2\\(m_1m_2,b)=1}}^* a_\xi(m_1)\overline{a_\xi(m_2)} \chi_{m_1}(kf^4)\overline{\chi_{m_2}(kf^4)} \end{aligned}
×W0 ⁣(N(bfm1)D)W0 ⁣(N(bfm2)D)‾Φ^ ⁣(HN(k)N(f)2N(m1)N(m2)).\quad{}\times W_0\!\Bigl(\frac{N(bfm_1)}D\Bigr) \overline{W_0\!\Bigl(\frac{N(bfm_2)}D\Bigr)} \widehat\Phi\!\Bigl( \frac{HN(k)}{N(f)^2N(m_1)N(m_2)} \Bigr).

The characters enforce (mj,f)=1(m_j,f)=1, leaving only the displayed exclusion (mj,b)=1(m_j,b)=1. The supports of W0W_0 and Φ^\widehat{\Phi} give

N(b)N(f)≤b0D,0<N(k)≤CΦb02D2HN(b)2.N(b)N(f)\le b_0D,\qquad0<N(k)\le\frac{C_{\Phi}b_0^2D^2}{HN(b)^2}.

Partition B≤N(b)<2BB\le N(b)<2B and F≤N(f)<2FF\le N(f)<2F into dyadic ranges. Denote the corresponding contribution to (4.16) by Sξ;B,F\mathcal S_{\xi;B,F}, and put

X=DBF,H=CI,ΦD2HB2,CI,Φ=max⁡(2,CΦb02).X=\frac{D}{BF},\qquad\mathcal{H}=\frac{C_{I,\Phi}D^2}{HB^2},\qquad C_{I,\Phi}=\max(2,C_{\Phi}b_0^2).

These satisfy (4.11). In the variables xj=N(mj)/Xx_j=N(m_j)/X, the coupled smooth weight is

Kb,f,k(x1,x2)=W0(rx1)W0(rx2)‾Φ^ ⁣(ax1x2),\mathcal K_{b,f,k}(x_1,x_2) =W_0(r x_1)\overline{W_0(r x_2)} \widehat\Phi\!\Bigl(\frac{a}{x_1x_2}\Bigr),
r=N(b)N(f)BF∈[1,4),a=HN(k)N(f)2X2≤CI,Φ.r=\frac{N(b)N(f)}{BF}\in[1,4),\qquad a=\frac{HN(k)}{N(f)^2X^2}\le C_{I,\Phi}.

For every integer q≥0q\ge0, these kernels satisfy

sup⁡b,f,k∥Kb,f,k∥Cq([a0/4,b0]2)≪I,Φ,q∥W∥Cq(I)2.\sup_{b,f,k}\|\mathcal K_{b,f,k}\|_{C^q([a_0/4,b_0]^2)} \ll_{I,\Phi,q}\|W\|_{C^q(I)}^2.

For U∈Cc∞(I∗)U\in C_c^\infty(I_*), where I∗=[a0/8,2b0]I_*=[a_0/8,2b_0], put

SU(b,f,k)=∑(n,S)=1(n,b)=1∗aξ(n)χn(k)χn(f)4U(N(n)/X).S_U(b,f,k)=\sum_{\substack{(n,S)=1\\(n,b)=1}}^{*}a_{\xi}(n)\chi_n(k)\chi_n(f)^4U(N(n)/X).

The shifts (X,F)↦(X/N(d),FN(d))(X,F)\mapsto(X/N(d),FN(d)) preserve (4.11). Thus Lemma 4.4 and (4.12), with J=J(ϑ,ε/4)J=J(\vartheta,\varepsilon/4), give

∑f,k∣SU(b,f,k)∣2=XF E(b)(H,X,F;ξ,U)≤XF τdiv(b)∑d∣bE(H,X/N(d),FN(d);ξ,U)\begin{aligned} \sum_{f,k}|S_U(b,f,k)|^2 &=XF\,\mathcal E_{(b)}(\mathcal H,X,F;\xi,U)\\ &\le XF\,\tau_{\mathrm{div}}(b) \sum_{d\mid b}\mathcal E(\mathcal H,X/N(d),FN(d);\xi,U)\end{aligned}
≪ν,S,I,ϑ,εDε/2(XF)2∥U∥CJ(I∗)2.\ll_{\nu,S,I,\vartheta,\varepsilon} D^{\varepsilon/2}(XF)^2\|U\|_{C^J(I_*)}^2.

Here N(b)≪IDN(b) \ll_I D, so the divisor factors are absorbed in Dε/4D^{\varepsilon/4}. We use the following smooth-weight principle, stated and proved in Lemma B.2: a mean-square bound valid for every common test function, with a CJC^J norm, also bounds the corresponding quadratic sum with a kernel depending on the row, at a cost given by the kernel’s C2J+4C^{2J+4} norm. For each fixed bb, apply Lemma B.2 with row index (f,k)(f,k), coefficient μ(f)1(f,b)=1\mu(f)\mathbf1_{(f,b)=1}, and kernel Kb,f,k\mathcal{K}_{b,f,k}. The larger row range 0<N(k)≤H0<N(k)\leq\mathcal{H} adds only terms whose original kernel vanishes. With q=2J+4q=2J+4, this gives

∣Sξ;B,F∣≪ν,S,I,ϑ,εDε/2∥W∥Cq(I)2∑B≤N(b)<2B∗HN(b)D2(XF)2≪Dε/2∥W∥Cq(I)2HBD2 B (XF)2=Dε/2∥W∥Cq(I)2H,\begin{aligned} |\mathcal S_{\xi;B,F}| &\ll_{\nu,S,I,\vartheta,\varepsilon} D^{\varepsilon/2}\|W\|_{C^q(I)}^2 \sum_{B\le N(b)<2B}^*\frac{HN(b)}{D^2}(XF)^2\\ &\ll D^{\varepsilon/2}\|W\|_{C^q(I)}^2 \frac{HB}{D^2}\,B\,(XF)^2 =D^{\varepsilon/2}\|W\|_{C^q(I)}^2H, \end{aligned}

where XF=D/BXF=D/B and ideal counting gives O(B)O(B) choices of bb. Summing the O((log⁡D)2)O((\log D)^2) nonempty dyadic ranges bounds Sξ\mathcal S_\xi as required. Sum over the fixed finite set of ξ\xi, include ZZ, and restore the factor DD to obtain (3.1). The required derivative order depends only on ϑ,ε\vartheta,\varepsilon.

Iteration of the dual mean-square estimate

Put Σ=XF\Sigma=XF. We prove the following estimate for the family (4.9), with exclusions removed by Lemma 4.4. The two inputs are proved in Sections 6 and 7.

Proposition 5.1. Fix κ>0\kappa>0 and C0≥1C_0\geq1. Suppose that

H,X,F≥1,Σ=XF≤DC0,H≤ΣD−κ.\mathcal{H},X,F\geq1,\qquad\Sigma=XF\leq D^{C_0},\qquad\mathcal{H}\leq\Sigma D^{-\kappa}.

For every ε>0\varepsilon>0 there is an integer J=J(κ,C0,ε)≥1J=J(\kappa,C_0,\varepsilon)\geq1 such that

E(H,X,F;ξ,W)≪I,ν,S,κ,C0,ε∥W∥CJ(I)2DεΣ\mathcal E(\mathcal H,X,F;\xi,W) \ll_{I,\nu,S,\kappa,C_0,\varepsilon} \|W\|_{C^J(I)}^2D^\varepsilon\Sigma

for every compact interval I⊂(0,∞)I\subset(0,\infty) and smooth WW supported in II.

The proof is given in Section 5.5, using Lemma 5.3 and Proposition 5.4.

The completed sums

Let Ψ\Psi be a completely multiplicative C\mathbb{C}-valued function on the nonzero integral ideals of OK=Z[ω]\mathcal{O}_K=\mathbb{Z}[\omega] coprime to 3, vanishing on ideals divisible by a prime in SS. For a primary element nn, write Ψ(n)=Ψ((n))\Psi(n)=\Psi((n)). Define

T(X;Ψ)=∑(n,S)=1∗∑b∈Ob≡1 (3)(b,S)=1α(n)‾γ2(n)Ψ(n)α(b)‾ 3Ψ(b)3NK/Q(n) NK/Q(b) V∗(NK/Q(n)NK/Q(b)3/X),T(X;\Psi)= \sum_{(n,S)=1}^*\sum_{\substack{b\in\mathcal O\\b\equiv1\ (3)\\(b,S)=1}} \frac{\overline{\alpha(n)}\gamma_2(n)\Psi(n) \overline{\alpha(b)}^{\,3}\Psi(b)^3} {\sqrt{\mathrm N_{K/\mathbb Q}(n)}\,\mathrm N_{K/\mathbb Q}(b)}\,V_*(\mathrm N_{K/\mathbb Q}(n)\mathrm N_{K/\mathbb Q}(b)^3/X),

where V∗(y)=yW(y)V_*(y)=\sqrt{y}W(y). The extra index bb supplies the cubes in the Fourier expansion of the Kubota theta function. Here n,bn,b are primary elements of O\mathcal{O}; only nn is required to be squarefree. The b=1b=1 part is exactly X−1/2∑(n,S)=1∗α(n)‾γ2(n)Ψ(n)W(NK/Q(n)/X)X^{-1/2}\sum_{(n,S)=1}^*\overline{\alpha(n)}\gamma_2(n)\Psi(n)W(\mathrm N_{K/\mathbb Q}(n)/X). To identify this with the normalized column sum in (4.9), fix a ray class character ξ\xi. For a nonzero row k∈Ok\in\mathcal{O} and a squarefree primary ff prime to SS, set

Ψk(n):=ξ(n)χn(k)χn(f)4,T(X;k,f):=T(X;Ψk).\Psi_k(n):=\xi(n)\chi_n(k)\chi_n(f)^4,\qquad T(X;k,f):=T(X;\Psi_k).

The displayed product defines Ψk(n)\Psi_k(n) for (n,S)=1(n,S)=1; set Ψk(n)=0\Psi_k(n)=0 otherwise. Thus the zero extension required in (5.3) is part of this definition.

The completed mean-square estimate

For the twist (5.4), the theta transformation converts the relevant sextic twists into quadratic characters. Combining it with the quadratic large sieve gives the following estimate, proved in Section 6.

Proposition 5.2. Fix ε>0\varepsilon> 0, C0≥1C_0 \ge1, and a character ξ\xi of the fixed ray class group. There is J=J(ε,C0)≥1J = J(\varepsilon,C_0) \ge1 such that

∑0<NK/Q(k)≪H∣T(X;k,f)∣2≪I,ξ,S,ε,C0Dε∥W∥CJ(I)2(H+H2NK/Q(f)X)\sum_{0<\mathrm N_{K/\mathbb Q}(k)\ll\mathcal H}|T(X;k,f)|^2 \ll_{I,\xi,S,\varepsilon,C_0} D^\varepsilon\|W\|_{C^J(I)}^2 \Bigl(\mathcal H+\frac{\mathcal H^2\mathrm N_{K/\mathbb Q}(f)}{X}\Bigr)

whenever

1≤H,X,NK/Q(f)≤DC0.1\le\mathcal H,X,\mathrm N_{K/\mathbb Q}(f)\le D^{C_0}.

Here ff is squarefree and primary with (f,S)=1(f,S)=1, WW is smooth and supported in a compact interval I⊂(0,∞)I \subset(0,\infty), and TT is defined by (5.4).

Cube inversion and the remaining sums

Set

Hc3=min⁡ ⁣(X,X2H2),Lb=XNK/Q(b)3.H_c^3=\min\!\Bigl(X,\frac{X^2}{\mathcal H^2}\Bigr), \qquad L_b=\frac{X}{\mathrm N_{K/\mathbb Q}(b)^3}.

Lemma 5.3. Fix C0≥1C_0 \ge1 and ε>0\varepsilon> 0. Suppose 1≤H,X,F≤DC01 \le\mathcal{H}, X, F \le D^{C_0} and H≤Σ=XF\mathcal{H} \le\Sigma= XF. There is an integer J=J(ε,C0)≥1J = J(\varepsilon,C_0) \ge1 such that

E(H,X,F;ξ,W)≪I,ν,S,C0,εDε(Σ∥W∥CJ(I)2+sup⁡b≡1 (3), (b,S)=1NK/Q(b)>Hc, Lb>1E(H,Lb,F;ξ,W))\mathcal E(\mathcal H,X,F;\xi,W) \ll_{I,\nu,S,C_0,\varepsilon}D^\varepsilon \biggl(\Sigma\|W\|_{C^J(I)}^2+ \sup_{\substack{b\equiv1\ (3),\ (b,S)=1\\ \mathrm N_{K/\mathbb Q}(b)>H_c,\ L_b>1}} \mathcal E(\mathcal H,L_b,F;\xi,W)\biggr)

for every compact interval I⊂(0,∞)I \subset(0,\infty) and W∈Cc∞(I)W\in C_c^\infty(I). An empty supremum is zero; in particular it is empty when H2≤X\mathcal{H}^2 \le X.

Proof. Write N=NK/QN=\mathrm N_{K/\mathbb Q} and I=[u,v]I=[u,v]. Complete multiplicativity in (5.3), including at the zeros of Ψk\Psi_k, gives

X−1/2∑(n,S)=1∗aξ(n)χn(k)χn(f)4W(N(n)/X)=∑(h,S)=1μ(h)α(h)‾ 3Ψk(h)3N(h)T(X/N(h)3;k,f).\begin{aligned} &X^{-1/2}\sum_{(n,S)=1}^*a_\xi(n)\chi_n(k)\chi_n(f)^4 W(N(n)/X)\\ &\qquad=\sum_{(h,S)=1} \frac{\mu(h)\overline{\alpha(h)}^{\,3}\Psi_k(h)^3}{N(h)} T(X/N(h)^3;k,f). \end{aligned}

Indeed, after substitution of (5.3), the coefficient of the total cube index bb contains ∑h∣bμ(h)=1b=1\sum_{h\mid b}\mu(h)=\mathbf{1}_{b=1}. This is Möbius inversion; compare [22] (1.18) and [7] (8.2).

Split the right-hand side of (5.8) into Pshort(k,f)+Plong(k,f)P_{\rm short}(k,f)+P_{\rm long}(k,f), where

Pshort(k,f)=∑(h,S)=1N(h)≤Hcμ(h)α(h)‾ 3Ψk(h)3N(h)T(X/N(h)3;k,f),\begin{aligned} P_{\rm short}(k,f)&= \sum_{\substack{(h,S)=1\\N(h)\le H_c}} \frac{\mu(h)\overline{\alpha(h)}^{\,3}\Psi_k(h)^3}{N(h)} T(X/N(h)^3;k,f),\end{aligned}
Eshort=1F∑F≤N(f)<2F(f,S)=1∗∑0<N(k)≤H∣Pshort(k,f)∣2.\begin{aligned} \mathcal E_{\rm short}&=\frac1F\sum_{\substack{F\le N(f)<2F\\(f,S)=1}}^* \sum_{0<N(k)\le\mathcal H}|P_{\rm short}(k,f)|^2. \end{aligned}

Define Elong\mathcal E_{\rm long} in the same way using PlongP_{\rm long}. Weighted Cauchy–Schwarz for each k,fk,f, followed by Proposition 5.2, gives

Eshort≤(∑N(h)≤Hc1N(h))∑N(h)≤Hc1N(h)1F∑F≤N(f)<2F(f,S)=1∗∑0<N(k)≤H∣T(X/N(h)3;k,f)∣2\begin{aligned} \mathcal E_{\rm short} &\le\Bigl(\sum_{N(h)\le H_c}\frac1{N(h)}\Bigr) \sum_{N(h)\le H_c}\frac1{N(h)} \frac1F\sum_{\substack{F\le N(f)<2F\\(f,S)=1}}^* \sum_{0<N(k)\le\mathcal H}|T(X/N(h)^3;k,f)|^2\end{aligned}
≪I,ν,S,C0,εDε/2∥W∥CJ(I)2(H+H2FHc3X)≪Dε/2Σ∥W∥CJ(I)2.\ll_{I,\nu,S,C_0,\varepsilon} D^{\varepsilon/2}\|W\|_{C^J(I)}^2 \Bigl(\mathcal H+\frac{\mathcal H^2FH_c^3}{X}\Bigr) \ll D^{\varepsilon/2}\Sigma\|W\|_{C^J(I)}^2.

Here X/N(h)3≥1X/N(h)^3\ge1, and the reciprocal-norm sum contributes only a logarithm. We use Proposition 5.2 with exponent ε/4\varepsilon/4 and range C0+1C_0+1, since N(f)<2DC0N(f)<2D^{C_0}. If Hc<1H_c<1, the sum is empty.

Expand TT using (5.3) in the remaining terms, with cube index cc:

Plong(k,f)=1X∑(h,S)=1N(h)>Hc∑(c,S)=1μ(h)N(hc) α(hc)‾ 3Ψk(hc)3×∑(n,S)=1∗aξ(n)χn(k)χn(f)4W ⁣(N(n)N(hc)3X).\begin{aligned} P_{\rm long}(k,f) ={}&\frac1{\sqrt X} \sum_{\substack{(h,S)=1\\N(h)>H_c}}\sum_{(c,S)=1} \mu(h)\sqrt{N(hc)}\,\overline{\alpha(hc)}^{\,3}\Psi_k(hc)^3\\ &\quad\times\sum_{(n,S)=1}^*a_\xi(n)\chi_n(k)\chi_n(f)^4 W\!\Bigl(\frac{N(n)N(hc)^3}{X}\Bigr). \end{aligned}

Put b=hcb=hc; the pairs giving bb are exactly (h,b/h)(h,b/h) with h∣bh\mid b and N(h)>HcN(h)>H_c. Since Ψk(b)3=ξ(b)3χb(k)31(b,f)=1\Psi_k(b)^3=\xi(b)^3\chi_b(k)^3\mathbf{1}_{(b,f)=1}, this gives

Plong(k,f)=∑(b,S)=1N(b)>Hcβ0(b,f)χb(k)3N(b)1Lb∑(n,S)=1∗aξ(n)χn(k)χn(f)4W(N(n)/Lb),\begin{aligned} P_{\rm long}(k,f) ={}&\sum_{\substack{(b,S)=1\\N(b)>H_c}} \frac{\beta_0(b,f)\chi_b(k)^3}{N(b)} \frac1{\sqrt{L_b}} \sum_{(n,S)=1}^*a_\xi(n)\chi_n(k)\chi_n(f)^4 W(N(n)/L_b),\end{aligned}
β0(b,f):=α(b)‾ 3ξ(b)31(b,f)=1∑h∣bN(h)>Hcμ(h),∣β0(b,f)∣≤τdiv(b).\begin{aligned} \beta_0(b,f):={}&\overline{\alpha(b)}^{\,3}\xi(b)^3 \mathbf1_{(b,f)=1}\sum_{\substack{h\mid b\\N(h)>H_c}}\mu(h), \qquad |\beta_0(b,f)|\le\tau_{\mathrm{div}}(b). \end{aligned}

The divisor function τdiv(b)\tau_{\mathrm{div}}(b) counts ideal divisors. All these indices are primary and prime to SS; bb need not be squarefree. The support of WW restricts the sum to N(b)3≤vXN(b)^3\le vX. Weighted Cauchy–Schwarz now gives

Elong≤(∑Hc<N(b)≤(vX)1/3(b,S)=1τdiv(b)N(b))∑Hc<N(b)≤(vX)1/3(b,S)=1τdiv(b)N(b)E(H,Lb,F;ξ,W).\begin{aligned} \mathcal E_{\rm long} &\le\Bigl(\sum_{\substack{H_c<N(b)\le(vX)^{1/3}\\(b,S)=1}} \frac{\tau_{\mathrm{div}}(b)}{N(b)}\Bigr) \sum_{\substack{H_c<N(b)\le(vX)^{1/3}\\(b,S)=1}} \frac{\tau_{\mathrm{div}}(b)}{N(b)} \mathcal E(\mathcal H,L_b,F;\xi,W). \end{aligned}

The reciprocal divisor sum is ≪Ilog⁡2(2+X)\ll_I\log^2(2+X), since

∑N(b)≤Rτdiv(b)N(b)=∑N(cd)≤R1N(c)N(d)≤(∑N(c)≤R1N(c))2≪log⁡2(2R)(R≥1).\sum_{N(b)\le R}\frac{\tau_{\mathrm{div}}(b)}{N(b)}=\sum_{N(cd)\le R}\frac{1}{N(c)N(d)}\le\left(\sum_{N(c)\le R}\frac{1}{N(c)}\right)^2\ll\log^2(2R)\qquad(R\ge1).

If Lb≤1L_b\le1, the column sum is empty unless Lb≥1/vL_b\ge1/v; otherwise it has OI(1)O_I(1) terms, so

E(H,Lb,F;ξ,W)≪IHLb∥W∥∞2≪IΣ∥W∥∞2.\mathcal E(\mathcal H,L_b,F;\xi,W) \ll_I\frac{\mathcal H}{L_b}\|W\|_\infty^2 \ll_I\Sigma\|W\|_\infty^2.

Absorb the logarithms in DεD^\varepsilon and combine the two parts using E≤2Eshort+2Elong\mathcal E\le2\mathcal E_{\rm short}+2\mathcal E_{\rm long}. This proves (5.7). If H2≤X\mathcal{H}^2\le X, then Hc=X1/3H_c=X^{1/3}, so N(b)>HcN(b)>H_c implies Lb<1L_b<1.

The transfer estimate

Fix a nonnegative radial Schwartz weight Φ\Phi with Φ(t)≥1\Phi(t) \ge1 for 0≤t≤10 \le t \le1 and supp⁡Φ^⊂[0,CΦ]\operatorname{supp}\widehat{\Phi} \subset[0,C_{\Phi}]. Such a weight is obtained by squaring and rescaling a real radial Schwartz function with compactly supported Fourier transform. Write N=NK/QN=\mathrm N_{K/\mathbb Q}. For fixed H,L,F,ξ\mathcal{H},L,F,\xi, put

A(W):=1LF∑F≤N(f)<2F(f,S)=1∗∑k∈OΦ(N(k)/H)∣∑(n,S)=1∗aξ(n)χn(k)χn(f)4W(N(n)/L)∣2.\mathcal A(W):=\frac1{LF} \sum_{\substack{F\le N(f)<2F\\(f,S)=1}}^* \sum_{k\in\mathcal O}\Phi(N(k)/\mathcal H) \Bigl|\sum_{(n,S)=1}^* a_\xi(n)\chi_n(k)\chi_n(f)^4W(N(n)/L)\Bigr|^2.

Thus E(H,L,F;ξ,W)≤A(W)\mathcal{E}(\mathcal{H},L,F;\xi,W)\le\mathcal{A}(W). In the following proposition Σ\Sigma is an independent target scale; we will apply it at column scale LbL_b with Σ=XF\Sigma=XF.

Proposition 5.4. Assume 1≤H,L,F,Σ≤DC01\le\mathcal{H},L,F,\Sigma\le D^{C_0} and max⁡{H,LF}≤Σ\max\{\mathcal{H},LF\}\le\Sigma, for fixed C0≥1C_0\ge1. For every integer m≥0m\ge0 and ε>0\varepsilon>0,

A(W)≪DεΣ∥W∥C4m+12(I)2(1+sup⁡E(H′,X′,F′;ξ′,U)Σ′)\mathcal A(W)\ll D^\varepsilon\Sigma\|W\|_{C^{4m+12}(I)}^2 \Bigl(1+\sup\frac{\mathcal E(\mathcal H',X',F';\xi',U)}{\Sigma'}\Bigr)

where the supremum is over the family (4.9), with Σ′=X′F′\Sigma'=X'F', H′,X′,F′≥1\mathcal{H}',X',F'\ge1, and

H′≤HLΣF,H′Σ′≤HΣ,Σ′≤L.\mathcal H'\le\frac{\mathcal HL}{\Sigma F},\qquad \frac{\mathcal H'}{\Sigma'}\le \frac{\mathcal H}{\Sigma},\qquad \Sigma'\le L.

For I=[u,v]I=[u,v], the tests satisfy U∈Cc∞(I′)U\in C_c^\infty(I'), I′=[u/16,4v]I'=[u/16,4v], and ∥U∥Cm(I′)≤1\|U\|_{C^m(I')}\le1; an empty supremum is zero. The estimate holds for every compact I⊂(0,∞)I\subset(0,\infty) and W∈Cc∞(I)W\in C_c^\infty(I), with implied constant depending on m,C0,ε,I,ν,Sm,C_0,\varepsilon,I,\nu,S and the fixed ray class group.

The proof is given in Section 7, by combining Lemmas 7.1 and 7.3.

Proof of Proposition 5.1

Proof of Proposition 5.1. Fix κ>0\kappa>0 and C0≥1C_0\ge1. We prove by induction on j≥0j\ge0 that the proposition holds under the additional restriction H≤Djκ\mathcal{H}\le D^{j\kappa}, with the derivative order and implied constant allowed to depend on jj. More precisely, for every ε>0\varepsilon>0 there is an integer J=J(j,κ,C0,ε)≥1J=J(j,\kappa,C_0,\varepsilon)\ge1 such that

E(H,X,F;ξ,W)≪I,ν,S,j,κ,C0,εDεΣ∥W∥CJ(I)2\mathcal E(\mathcal H,X,F;\xi,W) \ll_{I,\nu,S,j,\kappa,C_0,\varepsilon} D^\varepsilon\Sigma\|W\|_{C^J(I)}^2

for every compact interval I=[u,v]⊂(0,∞)I=[u,v]\subset(0,\infty) and W∈Cc∞(I)W\in C_c^\infty(I), under the hypotheses of the proposition. The derivative order is independent of II; this allows us to apply the induction hypothesis on the enlarged interval in the transfer estimate. For the base case j=0j=0, we have H≤1\mathcal{H}\le1. In (4.9), the support of WW restricts nn to the OI(X)O_I(X) ideals with NK/Q(n)∈[uX,vX]\mathrm N_{K/\mathbb Q}(n)\in[uX,vX]. The factors aξ(n)a_{\xi}(n), χn(k)\chi_n(k), and χn(f)4\chi_n(f)^4 have absolute value at most one, so the triangle inequality bounds each inner sum by OI(X∥W∥∞)O_I(X\|W\|_\infty). There are O(F)O(F) choices of ff and O(H)O(\mathcal{H}) choices of kk. Consequently,

E(H,X,F;ξ,W)≪IFHXF(X∥W∥∞)2=HX∥W∥∞2.\mathcal E(\mathcal H,X,F;\xi,W) \ll_I\frac{F\mathcal H}{XF} \bigl(X\|W\|_\infty\bigr)^2 =\mathcal H X\|W\|_\infty^2.

Dividing by Σ=XF\Sigma=XF and using H≤1\mathcal{H}\le1 and F≥1F\ge1 proves the assertion for j=0j=0, with J=1J=1.

Suppose the assertion holds for jj, and fix ε>0\varepsilon>0. Let mm be the derivative order supplied by the induction hypothesis with exponent ε/3\varepsilon/3. For H≤D(j+1)κ\mathcal{H}\le D^{(j+1)\kappa}, apply Lemma 5.3 with exponent ε/3\varepsilon/3. We must bound the mean squares in its supremum uniformly in bb. If this supremum is nonempty, then H2>X\mathcal{H}^2>X and Hc3=X2/H2H_c^3=X^2/\mathcal H^2. For each such bb, recall that Lb=X/NK/Q(b)3>1L_b=X/\mathrm N_{K/\mathbb Q}(b)^3>1. Since Lb≤XL_b\le X and Σ≥max⁡{H,LbF}\Sigma\ge\max\{\mathcal{H},L_bF\}, Proposition 5.4 applies at column scale LbL_b, with derivative order mm, and exponent ε/3\varepsilon/3. Its estimate also bounds E(H,Lb,F;ξ,W)\mathcal{E}(\mathcal{H},L_b,F;\xi,W) because E(H,Lb,F;ξ,W)≤A(W)\mathcal{E}(\mathcal{H},L_b,F;\xi,W)\le\mathcal{A}(W). Let H′\mathcal{H}', X′X', F′F' be any parameters in the supremum in (5.11), with Σ′=X′F′\Sigma'=X'F' as in that proposition. By (5.12) and NK/Q(b)>Hc\mathrm N_{K/\mathbb Q}(b)>H_c,

H′≤HLbΣF=HNK/Q(b)3F2<H(HΣ)2≤D−2κH≤Djκ,H′Σ′≤HΣ≤D−κ,Σ′≤Lb≤Σ≤DC0.\begin{aligned} \mathcal H' &\le\frac{\mathcal H L_b}{\Sigma F} =\frac{\mathcal H}{\mathrm N_{K/\mathbb Q}(b)^3F^2} <\mathcal H\Bigl(\frac{\mathcal H}{\Sigma}\Bigr)^2 \le D^{-2\kappa}\mathcal H\le D^{j\kappa},\\ \frac{\mathcal H'}{\Sigma'} &\le\frac{\mathcal H}{\Sigma}\le D^{-\kappa}, \qquad \Sigma'\le L_b\le\Sigma\le D^{C_0}. \end{aligned}

Thus each of these mean squares is covered by the induction hypothesis. Its weight UU is supported in [u/16,4v][u/16,4v] and has CmC^m norm at most one. Applying the induction hypothesis on this interval gives

E(H′,X′,F′;ξ′,U)Σ′≪I,ν,S,j,κ,C0,εDε/3.\frac{\mathcal{E}(\mathcal{H}',X',F';\xi',U)}{\Sigma'}\ll_{I,\nu,S,j,\kappa,C_0,\varepsilon}D^{\varepsilon/3}.

The enlarged interval is determined by II, so the implied constant has only the permitted dependence on the support. Substituting into (5.11) yields, uniformly in bb,

E(H,Lb,F;ξ,W)≪I,ν,S,j,κ,C0,εD2ε/3Σ∥W∥C4m+12(I)2.\mathcal E(\mathcal H,L_b,F;\xi,W) \ll_{I,\nu,S,j,\kappa,C_0,\varepsilon} D^{2\varepsilon/3}\Sigma\|W\|_{C^{4m+12}(I)}^2.

Substitution into (5.7) contributes the remaining factor Dε/3D^{\varepsilon/3}. Choose JJ at least 4m+124m+12 and at least the derivative order required by Lemma 5.3. We obtain

E(H,X,F;ξ,W)≪I,ν,S,j,κ,C0,εDεΣ∥W∥CJ(I)2.\mathcal E(\mathcal H,X,F;\xi,W) \ll_{I,\nu,S,j,\kappa,C_0,\varepsilon} D^\varepsilon\Sigma\|W\|_{C^J(I)}^2.

The choice of JJ depends only on jj, κ\kappa, C0C_0, ε\varepsilon. If the supremum in (5.7) is empty, the same bound follows directly from that lemma. This completes the induction.

Finally, take j=⌈C0/κ⌉j=\lceil C_0/\kappa\rceil. The hypothesis H≤Σ≤DC0\mathcal{H}\leq\Sigma\leq D^{C_0} ensures H≤Djκ\mathcal{H}\leq D^{j\kappa}, so the induction gives the proposition.

The original mean square and the exponent 11/1211/12

Proof of Proposition 3.1. Fix 0<ϑ≤1/100<\vartheta\leq1/10 and put H=D1+ϑH=D^{1+\vartheta}. The ranges in (4.11) give

H≤CD1−ϑB2,Σ=XF=DB,HΣ≪D−ϑ.\mathcal{H}\leq\frac{CD^{1-\vartheta}}{B^2}, \qquad \Sigma=XF=\frac{D}{B}, \qquad \frac{\mathcal{H}}{\Sigma}\ll D^{-\vartheta}.

For large DD, Proposition 5.1 applies with κ=ϑ/2\kappa=\vartheta/2 and C0=2C_0=2. If X<1X<1, nonempty support forces X≫I1X\gg_I 1, and counting gives E≪IH∥W∥∞2≪Σ∥W∥∞2\mathcal E\ll_I\mathcal H\|W\|_\infty^2 \ll\Sigma\|W\|_\infty^2; if H<1\mathcal{H}<1, the sum is empty. Thus (4.12) holds, with derivative order depending only on ϑ\vartheta, ε\varepsilon, and Proposition 4.5 proves (3.1). Bounded DD is again covered by counting.

Proof of the completed mean-square estimate

We prove Proposition 5.2. We first express the completed sum (5.3) using cubic theta coefficients and state the transformation formula. We then apply the quadratic large sieve and account for repeated prime factors in kk.

Realization by the cubic theta function

With the row kk and auxiliary index ff fixed, we express T(X;k,f)T(X;k,f) as a weighted sum of Fourier coefficients of the cubic theta function. Its automorphy then expresses this sum in terms of coefficients at other cusps, as in [7], §5 and Appendix A. We use Ψk\Psi_k from (5.4), suppressing its dependence on the fixed ff and ξ\xi.

Write (z,v)∈C×R>0(z,v) \in\mathbb{C} \times\mathbb{R}_{>0} for upper half-space coordinates, with horizontal coordinate zz and height vv. We use Kubota’s cubic theta function in the normalization of [7]:

θ(z,v)=35/22v2/3+∑ℓ∈λ−3Oℓ≠0τ(ℓ)vK1/3(4π∣ℓ∣v)exp⁡ ⁣(2πi(ℓz+ℓz‾)).\theta(z,v)=\frac{3^{5/2}}2v^{2/3} +\sum_{\substack{\ell\in\lambda^{-3}\mathcal O\\\ell\ne0}} \tau(\ell)vK_{1/3}(4\pi|\ell|v) \exp\!\bigl(2\pi\mathrm i(\ell z+\overline{\ell z})\bigr).

Here K1/3K_{1/3} is the modified Bessel function of the second kind, and τ(ℓ)\tau(\ell) is the coefficient sequence given explicitly in [7], following Patterson’s calculation [38]. We define Θk(z,v)\Theta_k(z,v) by twisting the Fourier coefficients of θˉ\bar{\theta}. First define ϕk:O→C\phi_k:\mathcal{O}\to\mathbb{C} by

ϕk(n)={χn(λ)2Ψk(n),n≡1(mod3), (n,S)=1,0,otherwise.\begin{aligned}\phi_k(n)= \begin{cases} \chi_n(\lambda)^2\Psi_k(n), & n\equiv1\pmod{3},\ (n,S)=1,\\ 0, & \text{otherwise}. \end{cases}\end{aligned}

Then we define Θk(z,v)\Theta_k(z,v) by multiplying the Fourier coefficient of θˉ\bar{\theta} at ℓ\ell, which is τ(−ℓ)‾\overline{\tau(-\ell)}, by ϕk(λ3ℓ)\phi_k(\lambda^3\ell):

Θk(z,v)=∑ℓ∈λ−3Oℓ≠0τ(−ℓ)‾ϕk(λ3ℓ)vK1/3(4π∣ℓ∣v)exp⁡ ⁣(2πi(ℓz+ℓz‾)).\Theta_k(z,v)= \sum_{\substack{\ell\in\lambda^{-3}\mathcal O\\\ell\ne0}} \overline{\tau(-\ell)}\phi_k(\lambda^3\ell) vK_{1/3}(4\pi|\ell|v) \exp\!\bigl(2\pi\mathrm i(\ell z+\overline{\ell z})\bigr).

Recall that cθ(nb3)c_\theta(nb^3) is defined in (2.10). Choose a nonzero q∈Oq\in\mathcal{O} such that ϕk\phi_k is periodic modulo (q)(q), and put

ϕ^k(h)=1NK/Q(q)∑n mod qϕk(n)e(−hn/q).\widehat\phi_k(h)=\frac1{\mathrm N_{K/\mathbb Q}(q)}\sum_{n\bmod q}\phi_k(n)e(-hn/q).

Lemma 6.1. With these definitions,

T(X;k,f)=135/2X∑n,b∈On,b≡1 (3)n squarefree(nb,S)=1cθ(nb3)ϕk(nb3)α(nb3)‾W ⁣(NK/Q(nb3)X).T(X;k,f)=\frac1{3^{5/2}\sqrt X} \sum_{\substack{n,b\in\mathcal O\\n,b\equiv1\ (3)\\n\ {\rm squarefree}\\(nb,S)=1}} c_\theta(nb^3)\phi_k(nb^3)\overline{\alpha(nb^3)} W\!\Bigl(\frac{\mathrm N_{K/\mathbb Q}(nb^3)}X\Bigr).

Moreover,

Θk(z,v)=∑h∈O/(q)ϕ^k(h)θ(z+λ2h/q,v)‾.\Theta_k(z,v)=\sum_{h\in\mathcal O/(q)}\widehat\phi_k(h) \overline{\theta(z+\lambda^2h/q,v)}.

Proof. The coefficient formula (2.10) gives

cθ(nb3)ϕk(nb3)=35/2∣b∣γ2(n)Ψk(n)Ψk(b)3.c_\theta(nb^3)\phi_k(nb^3)=3^{5/2}|b|\gamma_2(n)\Psi_k(n)\Psi_k(b)^3.

Substitution into (5.3), using V∗(y)=yW(y)V_*(y)=\sqrt{y}W(y), proves (6.3). Finite Fourier inversion gives

ϕk(n)=∑h∈O/(q)ϕk^(h)e(nh/q),∑h∈O/(q)ϕk^(h)=ϕk(0)=0.\phi_k(n)=\sum_{h\in\mathcal{O}/(q)}\widehat{\phi_k}(h)e(nh/q),\qquad\sum_{h\in\mathcal{O}/(q)}\widehat{\phi_k}(h)=\phi_k(0)=0.

Multiplication of the ℓ\ellth Fourier mode by e(λ3ℓh/q)e(\lambda^3\ell h/q) translates zz to z+λ2h/qz+\lambda^2h/q. The second identity cancels the constant term, proving (6.4). □

The resulting transformation formula is stated in the next subsection; its automorphy calculation is given in Appendix A.2.

The theta transformation

Recall that, for a primary prime p∉Sp \notin S, χp(x)=(x/p)6\chi_p(x) = (x/p)_6 is the sextic residue character on (O/(p))×(\mathcal{O}/(p))^\times. For every integer jj, write χpj\chi_p^j for its jjth power on this group, extended by zero on multiples of pp; in particular, χp0(x)=1p∤x\chi_p^0(x)=\mathbf1_{p\nmid x}. Sextic reciprocity (4.1) expresses the factors of Ψk\Psi_k at primes outside SS as such powers, with 0≤j≤50 \leq j \leq5.

For a primary prime p∉Sp \notin S, j∈{0,…,5}j \in\{0,\ldots,5\}, and x∈Ox \in\mathcal{O}, define the local factor

Bp,j(x)={χp(x)−j−2,j≠0,4,NK/Q(p)−1/2(−1+NK/Q(p)1p∣x),j=4,NK/Q(p)−1/2χp(x)−2,j=0.\begin{aligned} B_{p,j}(x)= \begin{cases} \chi_p(x)^{-j-2},&j\ne0,4,\\ \mathrm N_{K/\mathbb Q}(p)^{-1/2}(-1+\mathrm N_{K/\mathbb Q}(p)\mathbf1_{p\mid x}),&j=4,\\ \mathrm N_{K/\mathbb Q}(p)^{-1/2}\chi_p(x)^{-2},&j=0. \end{cases} \end{aligned}

The key case is the quadratic character

Bp,1(x)=χp(x)3.B_{p,1}(x)=\chi_p(x)^3.

For the smooth compactly supported weight V∗V_* in (5.3), put V^∗(s)=∫0∞V∗(x)xs dx/x\widehat{V}_*(s)=\int_0^\infty V_*(x)x^s\,dx/x. We write ∫(σ)\int_{(\sigma)} for integration upwards along the vertical line with real part σ\sigma. With Γ\Gamma denoting Euler’s gamma function, the accompanying transform of the weight is

V∗♯(x)=12πi∫(0)V^∗(−t)Γ(7/6+t)Γ(5/6+t)Γ(7/6−t)Γ(5/6−t)((2π)4x27)−t dt,x>0.V_*^\sharp(x)=\frac1{2\pi\mathrm i}\int_{(0)} \widehat V_*(-t) \frac{\Gamma(7/6+t)\Gamma(5/6+t)} {\Gamma(7/6-t)\Gamma(5/6-t)} \Bigl(\frac{(2\pi)^4x}{27}\Bigr)^{-t}\,d t,\qquad x>0.

Here θ\theta is Kubota’s cubic theta function, normalized in (6.1). The dual sums use the Fourier coefficients of θ‾\overline{\theta} at three cusps. Take the representatives

γ0=I,γ+=(10ω1),γ−=(10ω21),\begin{aligned} \gamma_0=I,\qquad \gamma_+=\Bigl(\begin{matrix}1&0\\\omega&1\end{matrix}\Bigr),\qquad \gamma_-=\Bigl(\begin{matrix}1&0\\\omega^2&1\end{matrix}\Bigr), \end{aligned}

acting on upper half-space. For σ∈{0,+,−}\sigma\in\{0,+,-\}, define dσ(ℓ)d_\sigma(\ell) by the Fourier expansion [7], (5.9), (5.15)

θ(γσ(z,v))‾=35/221σ=0v2/3+∑0≠ℓ∈λ−4Odσ(ℓ)vK1/3(4π∣ℓ∣v)exp⁡ ⁣(2πi(ℓz+ℓz‾)).\overline{\theta(\gamma_\sigma(z,v))} =\frac{3^{5/2}}2\mathbf1_{\sigma=0}v^{2/3} +\sum_{0\ne\ell\in\lambda^{-4}\mathcal O} d_\sigma(\ell)vK_{1/3}(4\pi|\ell|v) \exp\!\bigl(2\pi\mathrm i(\ell z+\overline{\ell z})\bigr).

Thus d0(ℓ)=τ(−ℓ)‾d_0(\ell)=\overline{\tau(-\ell)}, with τ\tau extended by zero outside λ−3O\lambda^{-3}\mathcal{O}. In particular, cθ(nb3)=d0(λ−3nb3)c_\theta(nb^3)=d_0(\lambda^{-3}nb^3) for the indices in (2.10); the outline uses dθ(m)d_\theta(m) for dσ(λ−4m)d_\sigma(\lambda^{-4}m) with one of these three choices of σ\sigma The arithmetic formulas for all three sequences are given by (A.6) and (A.7) in Appendix A.2.

We state the formula for a general product of local twists. Fix a ray class character whose conductor is supported on SS, and let Ψ0\Psi_0 be its extension by zero at every prime in SS. Let P\mathcal{P} be a finite set of primes outside SS, each represented by its primary generator in O\mathcal{O}, and choose integers 0≤jp≤50\leq j_p\leq5 for each p∈Pp\in\mathcal{P}. For primary n∈On\in\mathcal{O}, set

Ψ(n)=Ψ0(n)∏p∈Pχpjp(n).\Psi(n)=\Psi_0(n)\prod_{p\in\mathcal{P}}\chi_p^{j_p}(n).

In the transformed sums, A\mathcal{A} will range over subsets satisfying

{p∈P:jp≠0}⊆A⊆P.\{p\in\mathcal{P}:j_p\ne0\}\subseteq\mathcal{A}\subseteq\mathcal{P}.

Thus only primes with jp=0j_p=0 may be omitted from A\mathcal{A}. We call the primes in A\mathcal{A} active and those in P∖A\mathcal{P}\setminus\mathcal{A} inactive. For such a subset and c0∈O∖{0}c_0\in\mathcal{O}\setminus\{0\}, write

c=c0∏p∈Ap.c=c_0\prod_{p\in\mathcal{A}}p.

Proposition 6.2. The quantity T(X;Ψ)T(X;\Psi) defined in (5.3) is a sum of OΨ0,S(2∣P∣)O_{\Psi_0,S}(2^{|\mathcal{P}|}) terms of the form

C∑0≠ℓ∈λ−4Od(ℓ)α(ℓ)NK/Q(ℓ) ψ(λ4ℓ)∏p∈ABp,jp(λ4ℓ)V∗♯ ⁣(NK/Q(ℓ)XNK/Q(c)2).C\sum_{0\ne\ell\in\lambda^{-4}\mathcal O} \frac{d(\ell)\alpha(\ell)}{\sqrt{\mathrm N_{K/\mathbb Q}(\ell)}}\, \psi(\lambda^4\ell) \prod_{p\in\mathcal A}B_{p,j_p}(\lambda^4\ell) V_*^\sharp\!\Bigl(\frac{\mathrm N_{K/\mathbb Q}(\ell) X}{\mathrm N_{K/\mathbb Q}(c)^2}\Bigr).

Here A\mathcal{A} and cc are as above, and ∣C∣≪Ψ0,S1|C|\ll_{\Psi_0,S}1. The triples (d,ψ,c0)(d,\psi,c_0) belong to a fixed finite family depending only on Ψ0,S\Psi_0,S, where d∈{d0,d+,d−}d\in\{d_0,d_+,d_-\}, c0∈O∖{0}c_0\in\mathcal{O}\setminus\{0\}, and ψ\psi is a unit-modulus additive character on O\mathcal{O}.

To average the transformed sums over k0k_0, we need to choose their cusp coefficients and additive characters consistently as k0k_0 varies. Fix a finite set Pfix\mathcal P_{\rm fix} of primary primes outside SS and exponents jp∈{0,…,5}j_p\in\{0,\ldots,5\}, and put

Ψk0(n)=Ψ0(n)∏p∈Pfixχpjp(n)∏p∣k0χp(n),\Psi_{k_0}(n)=\Psi_0(n) \prod_{p\in\mathcal P_{\rm fix}}\chi_p^{j_p}(n) \prod_{p\mid k_0}\chi_p(n),

where k0k_0 is primary and squarefree, with (k0,S∏p∈Pfixp)=1(k_0,S\prod_{p\in\mathcal P_{\rm fix}}p)=1. Let Afix\mathcal A_{\rm fix} range over the subsets satisfying

{p∈Pfix:jp≠0}⊆Afix⊆Pfix,A=Afix∪{p:p∣k0}.\{p\in\mathcal P_{\rm fix}:j_p\ne0\} \subseteq\mathcal A_{\rm fix}\subseteq\mathcal P_{\rm fix}, \qquad \mathcal A=\mathcal A_{\rm fix}\cup\{p:p\mid k_0\}.

Every prime dividing k0k_0 belongs to A\mathcal{A}, since its exponent is 1.

Lemma 6.3 (Uniformity in the twist). For the family Ψk0\Psi_{k_0} above, the summands in Proposition 6.2 may be indexed by (h,Afix)(h,\mathcal A_{\rm fix}), with hh in a fixed finite set depending only on Ψ0,S\Psi_0,S. Zero scalar coefficients are permitted. For each fixed index, the triple (d,ψ,c0)(d,\psi,c_0) depends on k0k_0 only through its ray class modulo a fixed ideal supported on SS and depending only on Ψ0,S\Psi_0,S.

To estimate the dual sums, we need bounds for the cusp coefficients and the transformed weight.

Lemma 6.4. For each d∈{d0,d+,d−}d\in\{d_0,d_+,d_-\}, the coefficient d(ℓ)d(\ell) vanishes unless ℓ\ell can be written as

ℓ=uλmnb3,\ell=u\lambda^mnb^3,

where u∈O×u\in\mathcal{O}^{\times}, m∈Zm\in\mathbb{Z} with m≥−4m\ge-4, and n,b∈On,b\in\mathcal{O} are primary with nn squarefree.

∣d(ℓ)∣=∣d(uλmnb3)∣≤27⋅3m/6∣b∣.|d(\ell)|=|d(u\lambda^mnb^3)|\le27\cdot3^{m/6}|b|.

For A>0A>0, integers j≥0j\ge0, and V∗V_* supported in a fixed compact interval I⊂(0,∞)I\subset(0,\infty), there is J=J(A,j)J=J(A,j) such that

∣(x∂x)jV∗♯(x)∣≪A,j,I∥V∗∥CJ(I)min⁡(x1/4,x−A),x>0.|(x\partial_x)^jV_*^\sharp(x)|\ll_{A,j,I}\|V_*\|_{C^J(I)}\min(x^{1/4},x^{-A}),\qquad x>0.

Proposition 6.2 and Lemmas 6.3 and 6.4 are proved in Appendix A.2.

The quadratic mean square on the dual side

We record the quadratic large-sieve estimate needed below.

Lemma 6.5. Let H,U≥1\mathcal{H},U\ge1, and let β(n)\beta(n) be complex coefficients on squarefree primary nn with NK/Q(n)≍U\mathrm N_{K/\mathbb Q}(n)\asymp U. For every ε>0\varepsilon>0,

∑k≡1 (3), (k,S)=1NK/Q(k)≤H∗∣∑n≡1 (3), n squarefreeNK/Q(n)≍Uβ(n)χk(n)3∣2≪S,ε(HU)ε(H+U)∑n∣β(n)∣2.\sum_{\substack{k\equiv1\ (3),\ (k,S)=1\\\mathrm N_{K/\mathbb Q}(k)\le\mathcal H}}^* \Bigl|\sum_{\substack{n\equiv1\ (3),\ n\ {\rm squarefree}\\ \mathrm N_{K/\mathbb Q}(n)\asymp U}} \beta(n)\chi_k(n)^3\Bigr|^2 \ll_{S,\varepsilon}(\mathcal HU)^\varepsilon(\mathcal H+U) \sum_n|\beta(n)|^2.

The star restricts kk to squarefree primary elements.

Proof. This is Goldmakher and Louvel’s quadratic large sieve [13], Theorem 1.1, after fixing the product of the prime factors of nn lying in SS, and finitely many ray classes. For completeness, let ek∈{0,1}e_k \in\{0,1\} according as NK/Q(k)≡1,3(mod4)\mathrm N_{K/\mathbb Q}(k)\equiv1,3\pmod4, and let κλ\kappa_\lambda be the nontrivial character modulo λ\lambda. The character x↦(x/k)2κλ(x)ekx \mapsto(x/k)_2\kappa_\lambda(x)^{e_k} is trivial on units and has primitive conductor kλekk\lambda^{e_k}. When ek=0e_k=0, the factor κλ(x)ek\kappa_\lambda(x)^{e_k} is omitted. Classes modulo 24O24\mathcal{O} fix the supplementary characters and reciprocity factors; for coprime k1,k2k_1,k_2 in the same class, the product character has conductor k1k2k_1k_2. These are the hypotheses in [13], Definition 1 and §2.

Squarefree rows and the completed bound

Write N(a)=NK/Q(a)N(a)=\mathrm N_{K/\mathbb Q}(a). In (5.4), allow any g∈O∖{0}g\in\mathcal{O}\setminus\{0\} in place of ff, so Ψk(n)=ξ(n)χn(k)χn(g)4\Psi_k(n)=\xi(n)\chi_n(k)\chi_n(g)^4. The zero extension at SS is retained. Choose prime-ideal generators, primary away from 33, and extend multiplicatively to all ideals. We first bound squarefree rows.

Lemma 6.6. For every ε>0\varepsilon>0 and C0≥1C_0\geq1 there is J=J(ε,C0)≥1J=J(\varepsilon,C_0)\geq1 such that

∑0<N(s)≤Hs squarefree∣T(X;u0s,g)∣2≪I,ξ,S,ε,C0Dε∥W∥CJ(I)2(H+H2N(g)X)\sum_{\substack{0<N(s)\le\mathcal H\\s\ {\rm squarefree}}} |T(X;u_0s,g)|^2 \ll_{I,\xi,S,\varepsilon,C_0}D^\varepsilon\|W\|_{C^J(I)}^2 \Bigl(\mathcal H+\frac{\mathcal H^2N(g)}X\Bigr)

for 1≤H,X,N(g)≤DC01\leq\mathcal{H},X,N(g)\leq D^{C_0}, u0∈O×u_0\in\mathcal{O}^{\times}, every compact interval I⊂(0,∞)I\subset(0,\infty), and W∈Cc∞(I)W\in C_c^\infty(I). The sum uses the chosen generators of squarefree ideals, including those meeting SS.

Proof. By homogeneity assume ∥W∥CJ(I)≤1\|W\|_{C^J(I)}\le1, with JJ chosen below. Write

s=tk0,t=∏p∣sp∣g or p∈Sp,N(k0)≤H0:=HN(t).s=t k_0,\qquad t=\prod_{\substack{p\mid s\\p\mid g\ \text{or}\ p\in S}}p,\qquad N(k_0)\le\mathcal H_0:=\frac{\mathcal H}{N(t)}.

Fix tt; then k0k_0 is squarefree and primary, with (k0,g)=1(k_0,g)=1 and (k0,S)=1(k_0,S)=1. Discard empty ranges, so H0≥1\mathcal{H}_0\geq1, and retain these restrictions below. For a coefficient function A(n,b)A(n,b), a positive scale YY, and an integer m≥−4m\geq-4, put

Sm[A,Y](k0):=∑n,b≡1 (3)n squarefreeA(n,b)χk0(nb)33m/3N(n) N(b)V∗♯ ⁣(3mN(n)N(b)3H02YN(k0)2).\mathscr S_m[A,Y](k_0):= \sum_{\substack{n,b\equiv1\ (3)\\n\ {\rm squarefree}}} \frac{A(n,b)\chi_{k_0}(nb)^3}{3^{m/3}\sqrt{N(n)}\,N(b)} V_*^\sharp\!\Bigl(\frac{3^mN(n)N(b)^3\mathcal H_0^2} {YN(k_0)^2}\Bigr).

We shall obtain, after partitioning k0k_0 into fixed ray classes,

T(X;u0tk0,g)=∑ι∈I∑m≥−4cι,m(k0)Sm[Aι,m,Yι](k0),T(X;u_0tk_0,g) =\sum_{\iota\in\mathcal I}\sum_{m\ge-4} c_{\iota,m}(k_0)\mathscr S_m[A_{\iota,m},Y_\iota](k_0),

where I\mathcal{I} is independent of k0k_0, ∣I∣≼1|\mathcal I|\preccurlyeq1, ∣cι,m(k0)∣≪1|c_{\iota,m}(k_0)|\ll1, and Aι,mA_{\iota,m} is independent of k0k_0. For some aι,Yι>0a_\iota,Y_\iota>0, the coefficients and lengths satisfy

∣Aι,m(n,b)∣≪aι,aι2≤1,aι2Yι≪H02XN(t)2N(g).|A_{\iota,m}(n,b)|\ll a_\iota,\qquad a_\iota^2\leq1,\qquad a_\iota^2Y_\iota\ll\frac{\mathcal{H}_0^2}{X}N(t)^2N(g).

All column restrictions are included by extending Aι,mA_{\iota,m} by zero. Put k=u0tk0k=u_0tk_0. For P={p∉S:p∣kg}\mathcal{P}=\{p\notin S:p\mid kg\}, sextic reciprocity (4.1) gives

Ψk(n)=Ψ0(n)∏p∈Pχpjp(n),jp≡vp(k)+4vp(g)(mod6),0≤jp≤5.\Psi_k(n)=\Psi_0(n)\prod_{p\in\mathcal{P}}\chi_p^{j_p}(n),\qquad j_p\equiv v_p(k)+4v_p(g)\pmod{6},\qquad 0\leq j_p\leq5.

Even when jp=0j_p=0, the factor χp0(n)=1p∤n\chi_p^0(n)=\mathbf{1}_{p\nmid n} retains the zero extension at p∣kgp\mid kg. The fixed factor Ψ0\Psi_0 contains ξ\xi, the factors at SS, the unit factors, and n↦R(n,∏p∉Spvp(k))n\mapsto\mathcal{R}\left(n,\prod_{p\notin S}p^{v_p(k)}\right). The fourth power at gg contributes no reciprocity sign, and the factors at SS depend only on exponents modulo six. Thus Ψ0\Psi_0 ranges over a fixed finite family. In the notation of Lemma 6.3, take Pfix={p∉S:p∣tg}\mathcal P_{\rm fix}=\{p\notin S:p\mid tg\}, so that P=Pfix⊔{p:p∣k0}\mathcal P=\mathcal P_{\rm fix}\sqcup\{p:p\mid k_0\}; the exponents on Pfix\mathcal P_{\rm fix} are fixed with t,gt,g. Partitioning k0k_0 into fixed ray classes fixes Ψ0\Psi_0 and, by Lemma 6.3, the data d,ψ,c0d,\psi,c_0 in each transformed term.

Fix one transformed term, with active primes Afix\mathcal A_{\rm fix} away from k0k_0, and put c∗=c0∏p∈Afixpc_*=c_0\prod_{p\in\mathcal A_{\rm fix}}p. Thus c=c∗k0c=c_*k_0 in (6.9). By (6.10), write ℓ=uλmnb3\ell=u\lambda^mnb^3, with u,m,n,bu,m,n,b as there. Every prime dividing k0k_0 has jp=1j_p=1, and hence

Bp,1(uλm+4nb3)=χp(uλm+4)3χp(nb)3.B_{p,1}(u\lambda^{m+4}nb^3) =\chi_p(u\lambda^{m+4})^3\chi_p(nb)^3.

Indeed, the transformed exponent is −1−2≡3(mod6)-1-2\equiv3\pmod{6} and χp(b)9=χp(b)3\chi_p(b)^9=\chi_p(b)^3, also when p∣bp\mid b by zero extension. The coefficient bound and rapid decay in Lemma 6.4 give absolute convergence of the transformed series, so we may regroup its terms below. Set

Au,m(n,b)=d(uλmnb3)α(uλmnb3)ψ(uλm+4nb3)3m/6N(b),∣Au,m(n,b)∣≤27.A_{u,m}(n,b)=\frac{d(u\lambda^mnb^3)\alpha(u\lambda^mnb^3)\psi(u\lambda^{m+4}nb^3)}{3^{m/6}\sqrt{N(b)}},\qquad|A_{u,m}(n,b)|\le27.

The transformed term is therefore

C(k0)∑u∈O×∑m≥−4χk0(uλm+4)3∑n,b≡1 (3)n squarefreeAu,m(n,b)χk0(nb)33m/3N(n) N(b)×∏p∈AfixBp,jp(uλm+4nb3)V∗♯ ⁣(3mN(n)N(b)3XN(c∗)2N(k0)2),∣C(k0)∣≪1.\begin{aligned} &C(k_0)\sum_{u\in\mathcal O^\times}\sum_{m\ge-4} \chi_{k_0}(u\lambda^{m+4})^3 \sum_{\substack{n,b\equiv1\ (3)\\n\ {\rm squarefree}}} \frac{A_{u,m}(n,b)\chi_{k_0}(nb)^3} {3^{m/3}\sqrt{N(n)}\,N(b)}\\ &\qquad\times\prod_{p\in\mathcal A_{\rm fix}} B_{p,j_p}(u\lambda^{m+4}nb^3) V_*^\sharp\!\Bigl(\frac{3^mN(n)N(b)^3X} {N(c_*)^2N(k_0)^2}\Bigr),\qquad |C(k_0)|\ll1. \end{aligned}

The representation is unique: away from λ\lambda, the prime exponents of ℓ\ell are 3vp(b)3v_p(b) or 1+3vp(b)1+3v_p(b). Apart from the quadratic character and weight, the k0k_0-dependence is a bounded scalar.

Fix u,mu,m and write q=N(p)q=N(p) for an active prime p∤k0p\nmid k_0. For jp=4j_p=4, the local identity is

Bp,4(uλm+4nb3)=−q−1/2+q1/21p∣n+q1/21p∤n, p∣b.B_{p,4}(u\lambda^{m+4}nb^3)=-q^{-1/2}+q^{1/2}\mathbf{1}_{p\mid n}+q^{1/2}\mathbf{1}_{p\nmid n,\ p\mid b}.

To display the reindexing in (6.15), let A(p)(n,b)A^{(p)}(n,b) include all the other local factors and set

A(p,n)(n,b)=1p∤nA(p)(pn,b),A(p,b)(n,b)=1p∤nA(p)(n,pb).A^{(p,n)}(n,b)=\mathbf{1}_{p\nmid n}A^{(p)}(pn,b),\qquad A^{(p,b)}(n,b)=\mathbf{1}_{p\nmid n}A^{(p)}(n,pb).

Then the changes n=pn′n=pn' and b=pb′b=pb' give the exact identity

Sm[A(p)(n,b)Bp,4(uλm+4nb3),q2Y](k0)=−q−1/2Sm[A(p),q2Y](k0)+χk0(p)3Sm[A(p,n),qY](k0)+q−1/2χk0(p)3Sm[A(p,b),Y/q](k0).\begin{aligned} &\mathscr S_m[A^{(p)}(n,b)B_{p,4}(u\lambda^{m+4}nb^3),q^2Y](k_0)\\ &\quad=-q^{-1/2}\mathscr S_m[A^{(p)},q^2Y](k_0) +\chi_{k_0}(p)^3\mathscr S_m[A^{(p,n)},qY](k_0)\\ &\qquad\quad+q^{-1/2}\chi_{k_0}(p)^3 \mathscr S_m[A^{(p,b)},Y/q](k_0). \end{aligned}

In the second term p∤n′p\nmid n' preserves squarefreeness; in the third, p∤np\nmid n is retained and b′b' is unrestricted at pp. Both extracted phases have modulus one. The factors q1/2q^{1/2} in the local identity cancel against N(pn′)\sqrt{N(pn')} or leave q−1/2q^{-1/2} after division by N(pb′)N(pb').

Let aa track a bound ∣A(n,b)∣≤27a|A(n,b)|\le27a for the coefficient in (6.15), and let YY be its scale. Before inserting the active primes, these are a=1a=1 and Y=N(c0)2H02/XY=N(c_0)^2\mathcal H_0^2/X. Each active prime contributes q2q^2 to the squared conductor norm. Including this factor, the updates are

(a2,Y)⟼{(a2,q2Y),jp∉{0,4},(a2/q,q2Y),jp=0,(a2/q,q2Y), (a2,qY), (a2/q,Y/q),jp=4.\begin{aligned}(a^2,Y)\longmapsto \begin{cases} (a^2,q^2Y), & j_p\notin\{0,4\},\\ (a^2/q,q^2Y), & j_p=0,\\ (a^2/q,q^2Y),\ (a^2,qY),\ (a^2/q,Y/q), & j_p=4. \end{cases}\end{aligned}

Thus a2a^2 never increases. If p∣tp \mid t and p∉Sp \notin S, then jpj_p is odd, so a2Ya^2Y costs at most N(p)2N(p)^2. If p∣gp \mid g and p∤tp \nmid t, the cost is N(p)N(p) for jp=0,4j_p=0,4 and N(p)2N(p)^2 for jp=2j_p=2; the latter case requires vp(g)≥2v_p(g) \ge2. Hence

a2Y≪H02X∏p∣tp∉SN(p)2∏p∣gp∤t, p∉SN(p)vp(g)≤H02N(t)2N(g)X.a^2Y \ll\frac{\mathcal{H}_0^2}{X}\prod_{\substack{p\mid t\\p\notin S}}N(p)^2\prod_{\substack{p\mid g\\p\nmid t,\ p\notin S}}N(p)^{v_p(g)} \le\frac{\mathcal{H}_0^2N(t)^2N(g)}{X}.

Reindexing at distinct primes preserves the form and zero extensions. By Lemma 6.3, the transformed terms, units, and at most three choices per prime p∈Afixp\in\mathcal A_{\rm fix} with jp=4j_p=4 form an index set of divisor-bounded size in tgtg. Here and below, divisor-bounded in aa means bounded by Cτdiv(a)AC\tau_{\mathrm{div}}(a)^A for fixed constants C,AC,A; in particular, this is ≪εNK/Q(a)ε\ll_\varepsilon\mathrm N_{K/\mathbb Q}(a)^\varepsilon for every ε>0\varepsilon>0. For one index ι\iota, abbreviate a=aιa=a_\iota, Y=YιY=Y_\iota, and Am=Aι,mA_m=A_{\iota,m}. The local updates give Y≪DC1Y\ll D^{C_1} for some fixed C1=C1(C0)C_1=C_1(C_0). Split bb into B≤N(b)<2BB\le N(b)<2B, and nn by a smooth dyadic partition V(N(n)/U)V(N(n)/U), with U,B≥1U,B\ge1 and uniformly bounded cutoffs. For this part of (6.15), write

Sm;U,B(k0)=3−m/3∑B≤N(b)<2Bχk0(b)3N(b)Fb(k0),\mathscr S_{m;U,B}(k_0) =3^{-m/3}\sum_{B\le N(b)<2B}\frac{\chi_{k_0}(b)^3}{N(b)}F_b(k_0),
Fb(k0)=∑n squarefreeAm(n,b)χk0(n)3N(n)V(N(n)/U)V∗♯ ⁣(3mN(n)N(b)3H02YN(k0)2).F_b(k_0) =\sum_{n\ {\rm squarefree}} \frac{A_m(n,b)\chi_{k_0}(n)^3}{\sqrt{N(n)}}V(N(n)/U) V_*^\sharp\!\Bigl(\frac{3^mN(n)N(b)^3\mathcal H_0^2} {YN(k_0)^2}\Bigr).

Here n,bn,b remain primary. Put z∗=3mUB3/Yz_* = 3^mUB^3/Y and fix a decay exponent A>0A>0. With m,U,B,Ym,U,B,Y fixed, Lemma B.1, applied only in N(n)/UN(n)/U, gives the following representation in a real Mellin variable ss:

Fb(k0)=∫Rcb,k0(s)Gb,s(k0) ds,F_b(k_0)=\int_{\mathbb R}c_{b,k_0}(s)G_{b,s}(k_0)\,ds,
Gb,s(k0)=∑n squarefreeN(n)≍UAm(n,b)χk0(n)3N(n)(N(n)/U)is,\begin{aligned} G_{b,s}(k_0)&=\sum_{\substack{n\ {\rm squarefree}\\N(n)\asymp U}} \frac{A_m(n,b)\chi_{k_0}(n)^3}{\sqrt{N(n)}} (N(n)/U)^{\mathrm i s},\end{aligned}
∣cb,k0(s)∣≪I,A(1+z∗)−A(1+∣s∣)−2.|c_{b,k_0}(s)|\ll_{I,A}(1+z_*)^{-A}(1+|s|)^{-2}.

Indeed, the scale RR in (B.2) is 3mUN(b)3H02/(YN(k0)2)≥z∗3^mUN(b)^3\mathcal{H}_0^2/(YN(k_0)^2)\ge z_*. Lemma 6.4 supplies the required derivative bounds, independently of b,k0,U,B,mb,k_0,U,B,m. Since ∑N(n)≍UN(n)−1≪1\sum_{N(n)\asymp U}N(n)^{-1}\ll1, Lemma 6.5 gives, for any σ>0\sigma>0,

∑N(k0)≤H0∗∣Gb,s(k0)∣2≪S,σa2(H0U)σ(H0+U).\sum_{N(k_0)\le\mathcal H_0}^*|G_{b,s}(k_0)|^2 \ll_{S,\sigma}a^2(\mathcal H_0U)^\sigma(\mathcal H_0+U).

Apply weighted Cauchy–Schwarz to the Mellin integral, as recorded in (B.6), with the common majorant in (6.19), then weighted Cauchy–Schwarz in bb, using ∑B≤N(b)<2BN(b)−1≪1\sum_{B\le N(b)<2B}N(b)^{-1}\ll1:

∑N(k0)≤H0∗∣Sm;U,B(k0)∣2≤3−2m/3(∑B≤N(b)<2B1N(b))∑B≤N(b)<2B1N(b)∑N(k0)≤H0∗∣Fb(k0)∣2\sum_{N(k_0)\le\mathcal H_0}^*|\mathscr S_{m;U,B}(k_0)|^2 \le3^{-2m/3}\Bigl(\sum_{B\le N(b)<2B}\frac1{N(b)}\Bigr) \sum_{B\le N(b)<2B}\frac1{N(b)} \sum_{N(k_0)\le\mathcal H_0}^*|F_b(k_0)|^2
≪I,S,A,σa23−2m/3(H0U)σ(H0+U)(1+z∗)−2A.\ll_{I,S,A,\sigma}a^2 3^{-2m/3} (\mathcal H_0U)^\sigma(\mathcal H_0+U)(1+z_*)^{-2A}.

Since m≥−4m\ge-4, we have U≤81Yz∗U\le81Yz_* and hence H0U≪DC0+C1(1+z∗)\mathcal{H}_0U\ll D^{C_0+C_1}(1+z_*). Choose 0<σ≤10<\sigma\le1 small in terms of ε1>0\varepsilon_1>0, C0C_0 and take A=3A=3. Taking square roots in (6.20) gives

(∑N(k0)≤H0∗∣Sm;U,B(k0)∣2)1/2≪aDε13−m/3(H0+U)(1+3mUB3/Y)−2.\Bigl(\sum_{N(k_0)\le\mathcal H_0}^*|\mathscr S_{m;U,B}(k_0)|^2\Bigr)^{1/2} \ll aD^{\varepsilon_1}3^{-m/3} (\sqrt{\mathcal H_0}+\sqrt U)(1+3^mUB^3/Y)^{-2}.

The infinite dyadic sums converge: for dyadic U,B≥1U,B\ge1 and every Z>0Z>0,

∑U,B dyadic(1+UB3/Z)−2≪log⁡2(2+Z),∑U,B dyadicU(1+UB3/Z)−2≪Z.\sum_{U,B\ {\rm dyadic}}(1+UB^3/Z)^{-2}\ll\log^2(2+Z), \qquad \sum_{U,B\ {\rm dyadic}}\sqrt U(1+UB^3/Z)^{-2}\ll\sqrt Z.

These estimates make the sum of ℓ2(k0)\ell^2(k_0) norms finite. The triangle inequality and the geometric sum over m≥−4m \ge-4 therefore give

(∑N(k0)≤H0∗∣∑m≥−4cι,m(k0)Sm[Am,Y](k0)∣2)1/2≪aDε1(H0log⁡2(2+Y)+Y).\Bigl(\sum_{N(k_0)\le\mathcal H_0}^* \Bigl|\sum_{m\ge-4}c_{\iota,m}(k_0) \mathscr S_m[A_m,Y](k_0)\Bigr|^2\Bigr)^{1/2} \ll aD^{\varepsilon_1} \bigl(\sqrt{\mathcal H_0}\log^2(2+Y)+\sqrt Y\bigr).

Combine the divisor-bounded choices of ι\iota and use (6.16) to obtain

∑N(k0)≤H0∗∣T(X;u0tk0,g)∣2≼max⁡ιaι2(H0+Yι)≪H0+H02N(t)2N(g)X≤H+H2N(g)X.\sum_{N(k_0)\le\mathcal H_0}^*|T(X;u_0tk_0,g)|^2 \preccurlyeq\max_\iota a_\iota^2(\mathcal H_0+Y_\iota) \ll\mathcal H_0+\frac{\mathcal H_0^2N(t)^2N(g)}X \le\mathcal H+\frac{\mathcal H^2N(g)}X.

Sum the fixed ray classes and the divisor-bounded choices t∣rad⁡(g∏p∈Sp)t\mid\operatorname{rad}(g\prod_{\mathfrak p\in S}\mathfrak p), choosing ε1\varepsilon_1 and the other small-power losses in terms of ε\varepsilon. The weight estimates above require a fixed number J=J(ε,C0)J=J(\varepsilon,C_0) of derivatives. Homogeneity restores ∥W∥CJ(I)2\|W\|_{C^J(I)}^2 and proves (6.14). □

Proof of Proposition 5.2. Write uniquely k=u0sv2k=u_0sv^2, with ss squarefree and s,vs,v among the chosen ideal generators. No condition (s,v)=1(s,v)=1 is imposed. Since 8≡2(mod6)8\equiv2\pmod6, including the zero extensions,

χn(u0sv2)χn(f)4=χn(u0s)χn(fv2)4,T(X;u0sv2,f)=T(X;u0s,fv2).\chi_n(u_0sv^2)\chi_n(f)^4=\chi_n(u_0s)\chi_n(fv^2)^4,\qquad T(X;u_0sv^2,f)=T(X;u_0s,fv^2).

Apply Lemma 6.6 with row bound H/N(v)2\mathcal{H}/N(v)^2 and auxiliary twist g=fv2g=fv^2. Here N(g)≤HN(f)≤D2C0N(g)\le\mathcal{H}N(f)\le D^{2C_0}, so use that lemma with 2C02C_0. Summing its bounds gives

∑0<N(k)≤H∣T(X;k,f)∣2=∑u0∈O×∑N(v)2≤H∑N(s)≤H/N(v)2s squarefree∣T(X;u0s,fv2)∣2≪Dε∥W∥CJ(I)2(H+H2N(f)X)∑v1N(v)2.\begin{aligned} \sum_{0<N(k)\le\mathcal H}|T(X;k,f)|^2 &=\sum_{u_0\in\mathcal O^\times}\sum_{N(v)^2\le\mathcal H} \sum_{\substack{N(s)\le\mathcal H/N(v)^2\\s\ {\rm squarefree}}} |T(X;u_0s,fv^2)|^2\\ &\ll D^\varepsilon\|W\|_{C^J(I)}^2 \Bigl(\mathcal H+\frac{\mathcal H^2N(f)}X\Bigr) \sum_v\frac1{N(v)^2}. \end{aligned}

The last sum is ζK(2)<∞\zeta_K(2)<\infty, where ζK(s)=LK(s,1)\zeta_K(s)=L_K(s,\mathbf1) is the Dedekind zeta function of KK. This proves (5.5). □

Proof of the transfer proposition

We prove Proposition 5.4 for A(W)\mathcal{A}(W) defined in (5.10). Throughout this section, H,L,F,Σ,ξ\mathcal{H},L,F,\Sigma,\xi satisfy the hypotheses of that proposition. We retain the fixed weight Φ\Phi and the support bound CΦC_\Phi chosen before (5.10). The enlargement by Σ/(LF)≥1\Sigma/(LF)\ge1 in the intermediate mean square makes the second Poisson summation return to (4.9) with row range at most HL/(ΣF)\mathcal{H}L/(\Sigma F). Lemma B.2 separates the weights, and Lemma 4.4 removes the remaining exclusion.

First application of Poisson summation

Write N=NK/QN=\mathrm N_{K/\mathbb Q}, I=[u,v]I=[u,v], and I∗=[u/2,2v]I_*=[u/2,2v]. All ideal indices below are prime to SS and represented by their primary generators; a star additionally requires squarefreeness. The variables k,h,yk,h,y range over O\mathcal{O}. For squarefree CC, tt with (C,t)=1(C,t)=1, and for d∣Cd\mid C, put

ℓ=LN(C)N(t),wC,d=HN(C)N(d)L2F,YC,d=cIΣLFN(d)HN(C)2,\ell=\frac{L}{N(C)N(t)},\qquad w_{C,d}=\frac{\mathcal{H}N(C)}{N(d)L^2F},\qquad Y_{C,d}=c_I\frac{\Sigma LF N(d)}{\mathcal{H}N(C)^2},

where cI≥max⁡(1,4CΦv2)c_I\ge\max(1,4C_\Phi v^2) is fixed. Retain only ℓ≥1/(2v)\ell\ge1/(2v). For U∈Cc∞(I∗)U\in C_c^\infty(I_*) and a character ξ1\xi_1 of the fixed ray class group, define

py(n)=μ(n)ξ1(n)1(n,t)=1χn(y)‾χn(C)4χn(d),p_y(n)=\mu(n)\xi_1(n)\mathbf1_{(n,t)=1} \overline{\chi_n(y)}\chi_n(C)^4\chi_n(d),
PC,d,t(y;U)=∑n∗py(n)U(N(n)/ℓ).P_{C,d,t}(y;U)=\sum_n^*p_y(n)U(N(n)/\ell).

The notation pyp_y suppresses its dependence on C,d,t,ξ1C,d,t,\xi_1. The nonnegative form required in the second Poisson calculation is

Qξ1(U)=∑C,t squarefree(Ct,S)=1,(C,t)=1N(C)N(t)≤2vL∑d∣CwC,d∑yΦ(N(y)/YC,d)∣PC,d,t(y;U)∣2.\mathcal Q_{\xi_1}(U)= \sum_{\substack{C,t\text{ squarefree}\\(Ct,S)=1, (C,t)=1\\ N(C)N(t)\le2vL}} \sum_{d\mid C}w_{C,d} \sum_y\Phi(N(y)/Y_{C,d})|P_{C,d,t}(y;U)|^2.

The outer domain has no additional restriction coming from the support of Φ^\widehat{\Phi}. Here YC,d>0Y_{C,d}>0 may be smaller than 1; the outer sums are finite, and the yy-sum converges absolutely.

Lemma 7.1. Under the hypotheses of Proposition 5.4, fix an integer j≥0j\ge0 and ε0>0\varepsilon_0>0. If M≥0M\ge0 satisfies Qξ1(U)≤M∥U∥Cj(I∗)2\mathcal Q_{\xi_1}(U)\le M\|U\|_{C^j(I_*)}^2 for every U∈Cc∞(I∗)U\in C_c^\infty(I_*) and every ξ1\xi_1, then

A(W)≪Dε0(Σ+M)∥W∥C2j+4(I)2.\mathcal{A}(W)\ll D^{\varepsilon_0}(\Sigma+M)\|W\|_{C^{2j+4}(I)}^2.

Proof. Expand the square defining A(W)\mathcal{A}(W), and write ni=Cuin_i=Cu_i, where C=(n1,n2)C=(n_1,n_2) and (u1,u2)=(u1u2,C)=1(u_1,u_2)=(u_1u_2,C)=1. The row character is χu1χu2‾\chi_{u_1}\overline{\chi_{u_2}}, primitive modulo u1u2u_1u_2, with the extra restriction (k,C)=1(k,C)=1. Lemma 4.2 introduces d∣Cd\mid C and a frequency hh. Let ZZ denote the zero-frequency contribution. It requires u1=u2=1u_1=u_2=1, so

∣Z∣≪HLF∑F≤N(f)<2F∗∑C∗∣W(N(C)/L)∣2≪IH∥W∥∞2.|Z|\ll\frac{\mathcal{H}}{LF}\sum_{F\le N(f)<2F}^{*}\sum_C^{*}|W(N(C)/L)|^2\ll_I\mathcal{H}\|W\|_\infty^2.

For the other frequencies, the Chinese remainder theorem, (4.4), and (4.7) give the paired identity

aξ(u1)aξ(u2)‾γ(χu1χu2‾)=μ(u1)μ(u2)(ξG)(u1u2−1).a_{\xi}(u_1)\overline{a_{\xi}(u_2)}\gamma(\chi_{u_1}\overline{\chi_{u_2}})=\mu(u_1)\mu(u_2)(\xi G)(u_1u_2^{-1}).

Expand (ξG)(z)=∑ξ1ℓξ1ξ1(z)(\xi G)(z)=\sum_{\xi_1}\ell_{\xi_1}\xi_1(z) on the fixed ray class group. Insert 1(u1,u2)=1=∑t∣u1, t∣u2μ(t)\mathbf{1}_{(u_1,u_2)=1}=\sum_{t\mid u_1,\ t\mid u_2}\mu(t) and put ui=txiu_i=tx_i. The inverse is xi=ui/tx_i=u_i/t; squarefreeness imposes (xi,t)=1(x_i,t)=1, but no condition (x1,x2)=1(x_1,x_2)=1 remains. The factors at C,tC,t cancel against their conjugates, leaving μ(d)μ(t)\mu(d)\mu(t) and the restrictions (f,Ct)=(h,t)=1(f,Ct)=(h,t)=1. Thus

A(W)−Z=∑ξ1ℓξ1∑C,t squarefree(Ct,S)=1,(C,t)=1N(C)N(t)≤2vL∑d∣Cμ(d)μ(t)wC,d×∑F≤N(f)<2F(f,Ct)=1∗∑h≠0,(h,t)=1N(h)≤CΦv2L2N(d)/(HN(C)2)∑x1,x2∗phf2(x1)phf2(x2)‾Kh(N(x1)/ℓ,N(x2)/ℓ),\begin{aligned} &\mathcal A(W)-Z =\sum_{\xi_1}\ell_{\xi_1} \sum_{\substack{C,t\text{ squarefree}\\(Ct,S)=1, (C,t)=1\\ N(C)N(t)\le2vL}} \sum_{d\mid C}\mu(d)\mu(t)w_{C,d}\\ &\quad\times\sum_{\substack{F\le N(f)<2F\\(f,Ct)=1}}^* \sum_{\substack{h\ne0, (h,t)=1\\ N(h)\le C_\Phi v^2L^2N(d)/(\mathcal H N(C)^2)}} \sum_{x_1,x_2}^* p_{hf^2}(x_1)\overline{p_{hf^2}(x_2)} \mathscr K_h(N(x_1)/\ell,N(x_2)/\ell), \end{aligned}

where, for each integer j0≥0j_0\ge0,

Kh(z1,z2)=(z1z2)−1/2W(z1)W(z2)‾Φ^ ⁣(HN(h)N(C)2L2N(d)z1z2),\mathscr K_h(z_1,z_2) =(z_1z_2)^{-1/2}W(z_1)\overline{W(z_2)} \widehat\Phi\!\Bigl( \frac{\mathcal H N(h)N(C)^2}{L^2N(d)z_1z_2}\Bigr),
sup⁡C,d,h∥Kh∥Cj0(I2)≪I,Φ,j0∥W∥Cj0(I)2.\sup_{C,d,h}\|\mathscr K_h\|_{C^{j_0}(I^2)} \ll_{I,\Phi,j_0}\|W\|_{C^{j_0}(I)}^2.

Here χxi(C)4\chi_{x_i}(C)^4 supplies (xi,C)=1(x_i,C)=1; the other original zeros are supplied by phf2(xi)p_{hf^2}(x_i). The displayed frequency range follows from the support of Φ^\widehat{\Phi}.

For fixed C,d,tC,d,t, the map (f,h)↦y=hf2(f,h)\mapsto y=hf^2 has divisor-bounded multiplicity in yy. Since Σ/(LF)≥1\Sigma/(LF)\ge1, its nonzero image satisfies N(y)≤4CΦv2L2F2N(d)/(HN(C)2)≤YC,dN(y)\le4C_{\Phi}v^2L^2F^2N(d)/(\mathcal{H}N(C)^2)\le Y_{C,d}. Since Φ≥1\Phi\ge1 on [0,1][0,1] and is nonnegative,

∑C,t∑d∣CwC,d∑f,h∣PC,d,t(hf2;U)∣2≼Qξ1(U)≤M∥U∥Cj(I∗)2,\sum_{C,t}\sum_{d\mid C}w_{C,d}\sum_{f,h}|P_{C,d,t}(hf^2;U)|^2 \preccurlyeq\mathcal Q_{\xi_1}(U) \le M\|U\|_{C^j(I_*)}^2,

with the preceding ranges on the left. Apply Lemma B.2 to the full displayed expression for A(W)−Z\mathcal{A}(W)-Z. The kernel bound gives

∣A(W)−Z∣≪Dε0/2M∥W∥C2j+4(I)2.|\mathcal A(W)-Z| \ll D^{\varepsilon_0/2}M\|W\|_{C^{2j+4}(I)}^2.

Together with H≤Σ\mathcal{H}\leq\Sigma, this proves the result after allocating the divisor losses within ε0\varepsilon_0.

Second application of Poisson summation

Let C,tC,t be coprime squarefree primary elements prime to SS, and let d∣Cd\mid C. Fix a character ξ1\xi_1 of the ray class group. For arbitrary scales ℓ,Y>0\ell,Y>0 and U∈Cc∞(I∗)U\in C_c^\infty(I_*), with I∗=[a,b]⊂(0,∞)I_*=[a,b]\subset(0,\infty), put

P(y)=∑(n,S)=1∗μ(n)ξ1(n)1(n,t)=1χn(y)‾χn(C)4χn(d)U(N(n)/ℓ),M=∑y∈OΦ(N(y)/Y)∣P(y)∣2.\begin{aligned} P(y)&=\sum_{(n,S)=1}^*\mu(n)\xi_1(n)\mathbf1_{(n,t)=1} \overline{\chi_n(y)}\chi_n(C)^4\chi_n(d)U(N(n)/\ell),\\ \mathcal M&=\sum_{y\in\mathcal O}\Phi(N(y)/Y)|P(y)|^2. \end{aligned}

Here P(y)P(y) abbreviates PC,d,t(y;U)P_{C,d,t}(y;U) from (7.1), with the scale ℓ\ell now arbitrary.

The second Poisson summation turns the Möbius coefficients back into cubic Gauss-sum coefficients. To state the identity, put U0(x)=x−1/2U(x)U_0(x)=x^{-1/2}U(x) and expand

ξ1(z)G(z−1)=∑ξ′cξ′ξ′(z)\xi_1(z)G(z^{-1})=\sum_{\xi'}c_{\xi'}\xi'(z)

on the fixed ray class group.

On the transformed side, gg is the common divisor of the original columns, e∣ge\mid g comes from the row exclusion, and ww removes the remaining coprimality condition. The index hh is the nonzero Fourier frequency. Let g,wg,w range over squarefree primary elements prime to SS, ee over divisors of gg, and hh over O∖{0}\mathcal{O}\setminus\{0\}, subject to

(g,Ct)=(w,gCth)=1,N(gw)≤bℓ,N(h)≤CΦb2ℓ2N(e)YN(g)2.(g,Ct)=(w,gCth)=1,\qquad N(gw)\le b\ell,\qquad N(h)\le\frac{C_\Phi b^2\ell^2N(e)}{YN(g)^2}.

These conditions ensure that tg/etg/e and CewCew are coprime and squarefree. For each such choice and each character ξ′\xi', the new column scale, coefficient, and coupled smooth weight are

X′=ℓN(gw),X'=\frac{\ell}{N(gw)},
a♯(n)=aξ′(n)1(n,tg/e)=1χn(deh)χn(Cew)4,K(x1,x2)=U0(x1)U0(x2)‾Φ^(YN(h)N(g)2N(e)ℓ2x1x2).a^\sharp(n)=a_{\xi'}(n)\mathbf1_{(n,tg/e)=1}\chi_n(deh)\chi_n(Cew)^4,\qquad\mathscr K(x_1,x_2)=U_0(x_1)\overline{U_0(x_2)}\widehat\Phi\left(\frac{YN(h)N(g)^2}{N(e)\ell^2x_1x_2}\right).

Here K\mathscr K depends on U,Y,ℓ,g,e,hU,Y,\ell,g,e,h. A star on a sum over n1,n2n_1,n_2 restricts each index separately to squarefree primary elements. Define also

Z=YΦ^(0)∑(g,Ct)=1(g,S)=1∗∣U(N(g)/ℓ)∣2∏p∣g(1−1N(p)).Z=Y\widehat{\Phi}(0)\sum_{\substack{(g,Ct)=1\\(g,S)=1}}^{*}|U(N(g)/\ell)|^2\prod_{p\mid g}\left(1-\frac{1}{N(p)}\right).

Lemma 7.2 (Second Poisson identity). With the notation and index ranges above,

M=Z+∑g,e,w,h,ξ′Yμ(e)μ(w)cξ′N(g)N(e)ℓ∑(n1n2,S)=1∗a♯(n1)a♯(n2)‾K ⁣(N(n1)X′,N(n2)X′).\mathcal M=Z+ \sum_{g,e,w,h,\xi'}\frac{Y\mu(e)\mu(w)c_{\xi'}N(g)}{N(e)\ell} \sum_{(n_1n_2,S)=1}^*a^\sharp(n_1)\overline{a^\sharp(n_2)} \mathscr K\!\Bigl(\frac{N(n_1)}{X'},\frac{N(n_2)}{X'}\Bigr).

The zero-frequency term satisfies

∣Z∣≪I∗,ΦYℓ∥U∥∞2.|Z|\ll_{I_*,\Phi}Y\ell\|U\|_\infty^2.

Proof. Expand the square, put g=(n1,n2)g=(n_1,n_2) and ni=gzin_i=gz_i. The nonzero terms satisfy

(z1,z2)=(z1z2,gCt)=(g,Ct)=1,χn1(y)‾χn2(y)=1(y,g)=1χz1(y)‾χz2(y).(z_1,z_2)=(z_1z_2,gC t)=(g,C t)=1,\qquad\overline{\chi_{n_1}(y)}\chi_{n_2}(y)=\mathbf{1}_{(y,g)=1}\overline{\chi_{z_1}(y)}\chi_{z_2}(y).

Apply Lemma 4.2 to this last character, with exclusion gg and divisor e∣ge\mid g. Its zero frequency occurs exactly when z1=z2=1z_1=z_2=1 and gives the displayed ZZ; counting gg gives its bound. For the nonzero frequencies, (4.5)–(4.7) and the Chinese remainder theorem give

μ(z1)μ(z2)ξ1(z1)ξ1(z2)‾γ(χz1χz2‾)=∑ξ′cξ′aξ′(z1)aξ′(z2)‾.\mu(z_1)\mu(z_2)\xi_1(z_1)\overline{\xi_1(z_2)}\gamma(\overline{\chi_{z_1}\chi_{z_2}})=\sum_{\xi'}c_{\xi'}a_{\xi'}(z_1)\overline{a_{\xi'}(z_2)}.

Use also

χz(dh)χz(e)‾χz(C)4=χz(deh)χz(Ce)4,\chi_z(dh)\overline{\chi_z(e)}\chi_z(C)^4 =\chi_z(deh)\chi_z(Ce)^4,
U(N(gz1)/ℓ)U(N(gz2)/ℓ)‾N(z1)N(z2)=N(g)ℓU0(N(gz1)/ℓ)U0(N(gz2)/ℓ)‾.\frac{U(N(gz_1)/\ell)\overline{U(N(gz_2)/\ell)}} {\sqrt{N(z_1)N(z_2)}} =\frac{N(g)}\ell U_0(N(gz_1)/\ell)\overline{U_0(N(gz_2)/\ell)}.

Thus the full nonzero contribution is

M−Z=∑ξ′cξ′∑g(g,Ct)=1∗∑e∣gYμ(e)N(g)N(e)ℓ∑h≠0\begin{aligned} \mathcal M-Z={}&\sum_{\xi'}c_{\xi'} \sum_{\substack{g\\(g,Ct)=1}}^*\sum_{e\mid g} \frac{Y\mu(e)N(g)}{N(e)\ell}\sum_{h\ne0}\end{aligned}
×∑z1,z2(z1,z2)=1(z1z2,gCt)=1∗aξ′(z1)aξ′(z2)‾χz1(deh)χz2(deh)‾χz1(Ce)4χz2(Ce)4‾\begin{aligned} &\quad\times\sum_{\substack{z_1,z_2\\(z_1,z_2)=1\\(z_1z_2,gCt)=1}}^* a_{\xi'}(z_1)\overline{a_{\xi'}(z_2)} \chi_{z_1}(deh)\overline{\chi_{z_2}(deh)} \chi_{z_1}(Ce)^4\overline{\chi_{z_2}(Ce)^4}\end{aligned}
×U0 ⁣(N(gz1)ℓ)U0 ⁣(N(gz2)ℓ)‾Φ^ ⁣(YN(h)N(e)N(z1)N(z2)).\quad\times U_0\!\Bigl(\frac{N(gz_1)}\ell\Bigr) \overline{U_0\!\Bigl(\frac{N(gz_2)}\ell\Bigr)} \widehat\Phi\!\Bigl(\frac{YN(h)}{N(e)N(z_1)N(z_2)}\Bigr).

All starred ideal indices are prime to SS. Insert

1(z1,z2)=1=∑w∣z1, w∣z2μ(w),zi=wni.\mathbf1_{(z_1,z_2)=1}=\sum_{w\mid z_1,\ w\mid z_2}\mu(w), \qquad z_i=wn_i.

Here (w,gCt)=1(w,gCt)=1. Equation (4.6) and the zero-extended characters give

aξ′(wn)=aξ′(w)aξ′(n)χn(w)4((w,n)=1),a_{\xi'}(wn)=a_{\xi'}(w)a_{\xi'}(n)\chi_n(w)^4 \quad((w,n)=1),
∣aξ′(w)χw(deh)χw(Ce)4∣2=1(w,h)=1,|a_{\xi'}(w)\chi_w(deh)\chi_w(Ce)^4|^2=\mathbf1_{(w,h)=1},
1(n,gCtw)=1χn(deh)χn(Cew)4=1(n,tg/e)=1χn(deh)χn(Cew)4.\mathbf1_{(n,gCtw)=1}\chi_n(deh)\chi_n(Cew)^4 =\mathbf1_{(n,tg/e)=1}\chi_n(deh)\chi_n(Cew)^4.

Consequently the preceding full sum becomes

M−Z=∑ξ′cξ′∑g,w(g,Ct)=(w,gCt)=1∗∑e∣gYμ(e)μ(w)N(g)N(e)ℓ∑h≠0(h,w)=1\begin{aligned} &\mathcal M-Z=\sum_{\xi'}c_{\xi'} \sum_{\substack{g,w\\(g,Ct)=(w,gCt)=1}}^* \sum_{e\mid g}\frac{Y\mu(e)\mu(w)N(g)}{N(e)\ell} \sum_{\substack{h\ne0\\(h,w)=1}}\end{aligned}
×∑(n1n2,S)=1∗a♯(n1)a♯(n2)‾U0 ⁣(N(gwn1)ℓ)U0 ⁣(N(gwn2)ℓ)‾Φ^ ⁣(YN(h)N(e)N(w)2N(n1)N(n2)).\quad\times\sum_{(n_1n_2,S)=1}^*a^\sharp(n_1)\overline{a^\sharp(n_2)} U_0\!\Bigl(\frac{N(gwn_1)}\ell\Bigr) \overline{U_0\!\Bigl(\frac{N(gwn_2)}\ell\Bigr)} \widehat\Phi\!\Bigl( \frac{YN(h)}{N(e)N(w)^2N(n_1)N(n_2)}\Bigr).

There is no remaining condition (n1,n2)=1(n_1,n_2)=1. The column support gives N(gw)≤bℓN(gw)\le b\ell and N(z1z2)≤b2ℓ2/N(g)2N(z_1z_2)\le b^2\ell^2/N(g)^2; the support of Φ^\widehat{\Phi} therefore imposes the stated bound on N(h)N(h). Substituting (7.5) gives (7.6), with only finitely many terms.

Lemma 7.3. Return to I=[u,v]I=[u,v] and I∗=[u/2,2v]I_*=[u/2,2v]. For an integer m≥0m\ge0, let SmS_m be the supremum in (5.11). Then, for every ε0>0\varepsilon_0>0,

Qξ1(U)≪Dε0Σ(1+Sm)∥U∥C2m+4(I∗)2(U∈Cc∞(I∗)).\mathcal Q_{\xi_1}(U)\ll D^{\varepsilon_0}\Sigma(1+S_m) \|U\|_{C^{2m+4}(I_*)}^2 \qquad(U\in C_c^\infty(I_*)).

The test functions in SmS_m are supported in I′=[u/16,4v]I'=[u/16,4v] and have Cm(I′)C^m(I') norm at most one; an empty supremum is zero.

Proof. Write N=NK/QN=\mathrm N_{K/\mathbb Q}. Apply Lemma 7.2 with input length ℓ\ell, row scale YC,dY_{C,d}, and the same C,d,tC,d,t, using the character ξ1\xi_1. Here hh is the new Fourier variable for the sum over yy. Its nonzero output has squarefree g,w,e∣gg,w,e\mid g, and

(g,Ct)=(w,gCth)=1,r=tg/e,f′=Cew,k′=deh.(g,Ct)=(w,gCth)=1,\qquad r=tg/e,\qquad f'=Cew,\qquad k'=deh.

Thus r,f′r,f' are coprime and squarefree, and the coefficient a♯a^\sharp from (7.5) becomes

a♯(n)=aξ′(n)1(n,r)=1χn(k′)χn(f′)4.a^\sharp(n)=a_{\xi'}(n)\mathbf{1}_{(n,r)=1}\chi_n(k')\chi_n(f')^4.

Write X0′X'_0 for the individual column scale X′X' in (7.5); below, X′X' will denote a common dyadic scale. The length, coefficient, and Fourier parameter simplify to

X0′=ℓN(g)N(w)=LN(r)N(f′),X'_0=\frac{\ell}{N(g)N(w)}=\frac L{N(r)N(f')},
wC,dYC,dN(g)N(e)ℓ=cIΣN(r)L2,w_{C,d}\frac{Y_{C,d}N(g)}{N(e)\ell} =\frac{c_I\Sigma N(r)}{L^2},
YC,dN(h)N(g)2N(e)ℓ2=cIΣFN(k′)N(r)2HL.\frac{Y_{C,d}N(h)N(g)^2}{N(e)\ell^2} =\frac{c_I\Sigma F N(k')N(r)^2}{\mathcal H L}.

Put U0(x)=x−1/2U(x)U_0(x)=x^{-1/2}U(x). The coupled kernel is therefore

Kr,f′,k′(x1,x2)=U0(x1)U0(x2)‾Φ^ ⁣(cIΣFN(k′)N(r)2HLx1x2),\mathscr K_{r,f',k'}(x_1,x_2) =U_0(x_1)\overline{U_0(x_2)} \widehat\Phi\!\Bigl(\frac{c_I\Sigma F N(k')N(r)^2}{\mathcal H Lx_1x_2}\Bigr),

with both columns evaluated at xj=N(nj)/X0′x_j=N(n_j)/X'_0. In particular its nonzero support requires

N(r)N(f′)≤2vL,0<N(k′)≤4CΦv2cIHLΣFN(r)2≤HLΣFN(r)2,N(r)N(f')\le2vL,\qquad0<N(k')\le\frac{4C_{\Phi}v^2}{c_I}\frac{\mathcal{H}L}{\Sigma F N(r)^2}\le\frac{\mathcal{H}L}{\Sigma F N(r)^2},

since CΦ(2v)2/cI≤1C_{\Phi}(2v)^2/c_I\le1.

For fixed r,f′,k′r,f',k', all preimages are obtained by choosing

t∣r,f′=Cew,e∣k′,(w,k′)=1,d∣(C,k′),t\mid r,\qquad f'=Cew,\qquad e\mid k',\qquad(w,k')=1,\qquad d\mid(C,k'),

and setting g=e(r/t)g=e(r/t), h=k′/(de)h=k'/(de). The factorization f′=Cewf'=Cew is disjoint and squarefree. All these preimages have the same kernel: Ctgw=rf′C t g w=rf', and the Fourier parameter displayed above depends only on r,k′r,k'. Thus the support restrictions discard only zero kernels. The choices t∣rt\mid r contribute τdiv(r)\tau_{\mathrm{div}}(r). At a prime p∣f′p\mid f', the choices p∣Cp\mid C, p∣ep\mid e, and p∣wp\mid w, respectively, contribute

(1+1p∣k′)−1p∣k′−1p∤k′=1p∣k′(1+\mathbf{1}_{p\mid k'})-\mathbf{1}_{p\mid k'}-\mathbf{1}_{p\nmid k'}=\mathbf{1}_{p\mid k'}

to the sum of μ(e)μ(w)\mu(e)\mu(w) over the preimages; the two choices inside the first term are p∤dp\nmid d and p∣dp\mid d. Let Z\mathcal Z be the total zero-frequency contribution: the sum of the terms ZZ from Lemma 7.2, with the outer weights wC,dw_{C,d} and ranges in (7.2). Consequently, regrouping before taking absolute values gives

Qξ1(U)−Z=cIΣL2∑ξ′cξ′∑r,f′ squarefree(r,f′)=1N(r)τdiv(r)∑k′≠0f′∣k′×∑n1,n2∗a♯(n1)a♯(n2)‾Kr,f′,k′ ⁣(N(n1)X0′,N(n2)X0′).\begin{aligned} \mathcal Q_{\xi_1}(U)-\mathcal Z ={}&\frac{c_I\Sigma}{L^2} \sum_{\xi'}c_{\xi'} \sum_{\substack{r,f'\ \mathrm{squarefree}\\(r,f')=1}} N(r)\tau_{\mathrm{div}}(r)\sum_{\substack{k'\ne0\\f'\mid k'}}\\ &\quad\times \sum_{n_1,n_2}^*a^\sharp(n_1)\overline{a^\sharp(n_2)} \mathscr K_{r,f',k'}\!\Bigl(\frac{N(n_1)}{X'_0}, \frac{N(n_2)}{X'_0}\Bigr). \end{aligned}

Fix R≤N(r)<2RR\le N(r)<2R, F′≤N(f′)<2F′F'\le N(f')<2F' dyadically, and put

X′=LRF′,Σ′=X′F′=LR,H′=HLΣFR2.X'=\frac{L}{RF'},\qquad\Sigma'=X'F'=\frac{L}{R},\qquad\mathcal{H}'=\frac{\mathcal{H}L}{\Sigma FR^2}.

Retain only those rr for which a nonzero term occurs in these ranges. Since f′∣k′f'\mid k', each such rr satisfies

F′≤N(f′)≤N(k′)≤HLΣFN(r)2.F' \le N(f') \le N(k') \le\frac{\mathcal{H}L}{\Sigma F N(r)^2}.

In particular H′≥1\mathcal{H}' \ge1. For every j∣rj\mid r, using H≤Σ\mathcal{H}\le\Sigma, we also have

X′N(j)≥ΣFN(r)2HRN(j)≥FΣH≥1.\frac{X'}{N(j)} \ge\frac{\Sigma F N(r)^2}{\mathcal{H}R N(j)} \ge\frac{F\Sigma}{\mathcal{H}} \ge1.

Moreover 1≤X′/X0′<41\le X'/X'_0<4. The rescaled kernels

Kr,f′,k′(x1,x2)=Kr,f′,k′ ⁣(X′X0′x1,X′X0′x2)\mathcal K_{r,f',k'}(x_1,x_2) =\mathscr K_{r,f',k'}\!\Bigl(\frac{X'}{X'_0}x_1, \frac{X'}{X'_0}x_2\Bigr)

have support in [u/8,2v]2⊂int⁡(I′2)[u/8,2v]^2\subset\operatorname{int}(I'^2) and satisfy

sup⁡r,f′,k′∥Kr,f′,k′∥C2m+4(I′2)≪I∗,m,Φ∥U∥C2m+4(I∗)2,\sup_{r,f',k'}\|\mathcal K_{r,f',k'}\|_{C^{2m+4}(I'^2)} \ll_{I_*,m,\Phi}\|U\|_{C^{2m+4}(I_*)}^2,

where the supremum is over the retained indices in the dyadic block. This follows from the Schwartz bounds for Φ^\widehat{\Phi}. The common parameters satisfy exactly

H′Σ′=HΣFR≤HΣ,Σ′≤L,H′≤HLΣF.\frac{\mathcal H'}{\Sigma'}= \frac{\mathcal H}{\Sigma FR}\le\frac{\mathcal H}{\Sigma},\qquad \Sigma'\le L,\qquad \mathcal H'\le\frac{\mathcal H L}{\Sigma F}.

Fix δ>0\delta>0, to be chosen in terms of ε0\varepsilon_0 at the end. For each retained rr, enlarge the f′,k′f',k' ranges by positivity, keeping this set of rr fixed. Lemma 4.4 then gives, for V∈Cc∞(I′)V\in C_c^\infty(I'),

∑F′≤N(f′)<2F′∗∑0<N(k′)≤H′∣∑n∗a♯(n)V(N(n)/X′)∣2=Σ′E(r)(H′,X′,F′;ξ′,V)≪δDδ(Σ′)2Sm∥V∥Cm(I′)2.\begin{aligned} \sum_{F'\le N(f')<2F'}^*\sum_{0<N(k')\le\mathcal H'} \Bigl|\sum_n^*a^\sharp(n)V(N(n)/X')\Bigr|^2 &=\Sigma'\mathcal E_{(r)}(\mathcal H',X',F';\xi',V)\\ &\ll_\delta D^\delta(\Sigma')^2S_m \|V\|_{C^m(I')}^2. \end{aligned}

Indeed exclusion removal replaces (X′,F′)(X',F') by (X′/N(j),F′N(j))(X'/N(j),F'N(j)), j∣rj\mid r, preserving H′\mathcal{H}' and Σ′\Sigma'. Equation (7.7) gives X′/N(j)≥1X'/N(j)\ge1, and F′N(j)≥1F'N(j)\ge1 as well. Thus every resulting mean square lies in the defining supremum.

There are O(R)O(R) ideals rr in the dyadic range. Apply Lemma B.2 with row index (r,f′,k′)(r,f',k'), the preceding mean-square bound, and the uniform kernel bound. The absolute contribution of this dyadic block is at most

D2δΣRL2R(Σ′)2Sm∥U∥C2m+4(I∗)2=D2δΣSm∥U∥C2m+4(I∗)2,D^{2\delta}\frac{\Sigma R}{L^2}R(\Sigma')^2S_m\|U\|_{C^{2m+4}(I_*)}^2 =D^{2\delta}\Sigma S_m\|U\|_{C^{2m+4}(I_*)}^2,

because (ΣR/L2)R(L/R)2=Σ(\Sigma R/L^2)R(L/R)^2=\Sigma.

Finally, the zero-frequency bound in Lemma 7.2 gives

∣Z∣≪I∗,Φ∥U∥∞2∑N(Ct)≤2vL∑d∣CwC,dYC,dℓ|\mathcal Z| \ll_{I_*,\Phi}\|U\|_\infty^2 \sum_{N(Ct)\le 2vL}\sum_{d\mid C}w_{C,d}Y_{C,d}\ell
≪I∗Σ∥U∥∞2∑Cτdiv(C)N(C)2∑N(t)≤2vL1N(t)≪I∗,δDδΣ∥U∥∞2.\ll_{I_*}\Sigma\|U\|_\infty^2 \sum_C\frac{\tau_{\mathrm{div}}(C)}{N(C)^2} \sum_{N(t)\le 2vL}\frac1{N(t)} \ll_{I_*,\delta}D^\delta\Sigma\|U\|_\infty^2.

This uses only YC,d>0Y_{C,d}>0, not YC,d≥1Y_{C,d}\ge1. Sum the OI∗,C0((log⁡(2D))2)O_{I_*,C_0}((\log(2D))^2) dyadic blocks and the fixed finite character set, and take δ=ε0/4\delta=\varepsilon_0/4. The stated bound follows, with implied constant depending only on I∗,m,C0,ε0I_*,m,C_0,\varepsilon_0, the fixed cutoffs and arithmetic data.

Proof of Proposition 5.4. Use SmS_m as in Lemma 7.3, with mm as in the proposition. Lemma 7.3, with loss Dε/4D^{\varepsilon/4}, supplies the hypothesis of Lemma 7.1 with j=2m+4j=2m+4 and M≪Dε/4Σ(1+Sm)M\ll D^{\varepsilon/4}\Sigma(1+S_m). Applying that lemma with the same loss gives

A(W)≪Dε/2Σ(1+Sm)∥W∥C4m+12(I)2,\mathcal A(W)\ll D^{\varepsilon/2}\Sigma(1+S_m) \|W\|_{C^{4m+12}(I)}^2,

which proves (5.11).

Arithmetic identities and theta calculations

Reciprocity and Gauss sums

Proof of Lemma 4.1. Proof of (4.2). For squarefree nn, set G(n)=χn(4)‾γ3(n)G(n)=\overline{\chi_n(4)}\gamma_3(n). Its dependence on a fixed ray class will be verified below. Fix a primary prime p≡1(mod3)p \equiv1 \pmod{3} outside SS. For multiplicative characters AA, BB of (O/(p))×(\mathcal{O}/(p))^\times, extended by zero at 00, write

J(A,B)=∑x mod pA(x)B(1−x).J(A,B) = \sum_{x \bmod p} A(x)B(1-x).

For each y∈O/(p)y \in\mathcal{O}/(p), the substitution u=2x−1u=2x-1 gives

#{x mod p:4x(1−x)=y}=#{u mod p:u2=1−y}=1+χp3(1−y),\#\{x \bmod p : 4x(1-x)=y\} = \#\{u \bmod p : u^2=1-y\} = 1+\chi_p^3(1-y),

since p∤2p\nmid2 and χp3\chi_p^3 is the quadratic character, extended by zero at 00. Hence

χp(4)J(χp,χp)=∑x mod pχp(4x(1−x))=∑y mod pχp(y)(1+χp3(1−y))=J(χp,χp3).\chi_p(4)J(\chi_p,\chi_p)=\sum_{x \bmod p}\chi_p(4x(1-x))=\sum_{y \bmod p}\chi_p(y)(1+\chi_p^3(1-y))=J(\chi_p,\chi_p^3).

The Gauss–Jacobi identity (cf. [22], (3.18)) gives

γ1(p)γ3(p)γ4(p)=J(χp,χp3)NK/Q(p)=χp(4)J(χp,χp)NK/Q(p)=χp(4)γ1(p)2γ2(p).\frac{\gamma_1(p)\gamma_3(p)}{\gamma_4(p)} =\frac{J(\chi_p,\chi_p^3)}{\sqrt{\mathrm N_{K/\mathbb Q}(p)}} =\chi_p(4)\frac{J(\chi_p,\chi_p)}{\sqrt{\mathrm N_{K/\mathbb Q}(p)}} =\chi_p(4)\frac{\gamma_1(p)^2}{\gamma_2(p)}.

For j≢0(mod6)j\not\equiv0 \pmod{6}, the normalized Gauss sums satisfy

γj(p)γ−j(p)=χp(−1)j.\gamma_j(p)\gamma_{-j}(p)=\chi_p(-1)^j.

Taking j=2j=2 gives γ2(p)γ4(p)=χp2(−1)=1\gamma_2(p)\gamma_4(p)=\chi_p^2(-1)=1, since χp2\chi_p^2 is cubic. Combining this with the Gauss–Jacobi relation above gives, after cancelling γ1(p)\gamma_1(p) and rearranging,

γ1(p)γ2(p)=χp(4)‾γ3(p)γ2(p)3.\gamma_1(p)\gamma_2(p)=\overline{\chi_p(4)}\gamma_3(p)\gamma_2(p)^3.

To finish the prime case of (4.2), it remains to evaluate γ2(p)3\gamma_2(p)^3. We do so by computing J(χp2,χp2)J(\chi_p^2,\chi_p^2). Put q=NK/Q(p)q=\mathrm N_{K/\mathbb Q}(p). We have J(χp2,χp2)∈OJ(\chi_p^2,\chi_p^2)\in\mathcal{O} and NK/Q(J(χp2,χp2))=q\mathrm N_{K/\mathbb Q}\bigl(J(\chi_p^2,\chi_p^2)\bigr)=q. For 0≤ℓ≤(q−1)/30\leq\ell\leq(q-1)/3, the exponent (q−1)/3+ℓ(q-1)/3+\ell lies strictly between 00 and q−1q-1, so ∑x mod px(q−1)/3+ℓ=0\sum_{x\bmod p}x^{(q-1)/3+\ell}=0 in O/(p)\mathcal{O}/(p). Reducing modulo pp and expanding the binomial therefore gives

J(χp2,χp2)≡∑x mod px(q−1)/3(1−x)(q−1)/3=∑ℓ=0(q−1)/3(−1)ℓ((q−1)/3ℓ)∑x mod px(q−1)/3+ℓ=0(modp).\begin{aligned} J(\chi_p^2,\chi_p^2) &\equiv\sum_{x\bmod p}x^{(q-1)/3}(1-x)^{(q-1)/3}\\ &=\sum_{\ell=0}^{(q-1)/3}(-1)^\ell\binom{(q-1)/3}{\ell} \sum_{x\bmod p}x^{(q-1)/3+\ell} =0\pmod p. \end{aligned}

Thus J(χp2,χp2)/p∈OJ(\chi_p^2,\chi_p^2)/p\in\mathcal{O} has norm one, so J(χp2,χp2)J(\chi_p^2,\chi_p^2) is a unit multiple of pp. To determine the unit, write χp2(x(1−x))=ωjx\chi_p^2(x(1-x))=\omega^{j_x} with jx∈{0,1,2}j_x\in\{0,1,2\} for x∈O/(p)∖{0,1}x\in\mathcal{O}/(p)\setminus\{0,1\}. Then

∏x≠0,1x(1−x)=1⟹∑x≠0,1jx≡0(mod3).\prod_{x\ne0,1}x(1-x)=1\quad\Longrightarrow\quad\sum_{x\ne0,1}j_x\equiv0\pmod3.

Since (ω−1)2(\omega-1)^2 generates (3)(3), expansion modulo this ideal gives

J(χp2,χp2)=∑x≠0,1ωjx≡(q−2)+(ω−1)∑x≠0,1jx≡−1(mod3).J(\chi_p^2,\chi_p^2)=\sum_{x\ne0,1}\omega^{j_x}\equiv(q-2)+(\omega-1)\sum_{x\ne0,1}j_x\equiv-1\pmod3.

Together with p≡1(mod3)p \equiv1 \pmod{3}, this fixes the unit:

J(χp2,χp2)=−p.J(\chi_p^2,\chi_p^2)=-p.

Consequently

γ2(p)3=γ2(p)2γ4(p)=J(χp2,χp2)NK/Q(p)=−α(p).\gamma_2(p)^3 =\frac{\gamma_2(p)^2}{\gamma_4(p)} =\frac{J(\chi_p^2,\chi_p^2)}{\sqrt{\mathrm N_{K/\mathbb Q}(p)}} =-\alpha(p).

The Chinese remainder theorem extends these prime-modulus identities to (4.2) for squarefree nn, with one minus sign for each prime factor accounting for μ(n)\mu(n).

Proof of (4.6). For coprime squarefree primary a,ba,b prime to SS, the Chinese remainder theorem gives

γ2(ab)=γ2(a)γ2(b)χa(b)2χb(a)2,\gamma_2(ab)=\gamma_2(a)\gamma_2(b)\chi_a(b)^2\chi_b(a)^2,

and cubic reciprocity gives χa(b)2=χb(a)2\chi_a(b)^2=\chi_b(a)^2. Multiplying by α(ab)‾ξ(ab)\overline{\alpha(ab)}\xi(ab) therefore proves (4.6).

Proof of (4.1). We now prove (4.1) and the formula for G(n)G(n) in (4.3), including their dependence on a fixed ray class. We begin by evaluating the normalized quadratic Gauss sum. For c∈Oc \in\mathcal{O} coprime to 22, we will show that

Γquad(c):=∣c∣−1∑x mod ce(x2/c)=12∑y mod 2Oe(−cy2/4).\Gamma_{\rm quad}(c):=|c|^{-1}\sum_{x\bmod c}e(x^2/c) =\frac12\sum_{y\bmod2\mathcal O}e(-cy^2/4).

To justify (A.1), apply Poisson summation over z∈Oz \in\mathcal{O} to e(z2/c)e−πη∣z∣2e(z^2/c)e^{-\pi\eta|z|^2}, with η>0\eta>0. The Gaussian makes the sum convergent. The Fourier transform of the product is

∣c∣2rηe−πηNK/Q(c)∣y∣2/(4rη)e(−cy2/(4rη)),rη=1+3NK/Q(c)η2/16.\frac{|c|}{2\sqrt{r_\eta}} e^{-\pi\eta \mathrm N_{K/\mathbb Q}(c)|y|^2/(4r_\eta)}e(-cy^2/(4r_\eta)), \qquad r_\eta=1+3\mathrm N_{K/\mathbb Q}(c)\eta^2/16.

Divide both sides by the Gaussian mass

∑z∈Oe−πη∣z∣2≍η−1.\sum_{z\in\mathcal{O}}e^{-\pi\eta|z|^2}\asymp\eta^{-1}.

On the Fourier side, replacing rηr_\eta by 11 has total error

Oc ⁣(η2∑y∈O∣y∣2e−Ccη∣y∣2)=Oc(1),O_c\!\biggl(\eta^2\sum_{y\in\mathcal O}|y|^2e^{-C_c\eta|y|^2}\biggr) =O_c(1),

where Cc>0C_c>0 is fixed. The normalized error is therefore Oc(η)→0O_c(\eta)\to0. Grouping the original sum modulo cc and the Fourier sum modulo 2O2\mathcal{O}, then letting η↓0\eta\downarrow0, gives (A.1).

For c=a+bωc=a+b\omega, formula (A.1) becomes

Γquad(a+bω)=1+i−b+ia+ib−a2.\Gamma_{\mathrm{quad}}(a+b\omega)=\frac{1+i^{-b}+i^a+i^{b-a}}{2}.

Thus Γquad(c)\Gamma_{\mathrm{quad}}(c) depends only on cc mod 4O4\mathcal{O} and is invariant under multiplication by a square in (O/4O)×(\mathcal{O}/4\mathcal{O})^\times. Representatives for the four square classes and their values are

c1−1λ−λΓquad(c)11i−i.\begin{aligned}\begin{array}{c|rrrr} c&1&-1&\lambda&-\lambda\\ \hline \Gamma_{\rm quad}(c)&1&1&\mathrm i&-\mathrm i \end{array}\end{aligned}.

Consequently the function

R(c1,c2):=Γquad(c1c2)Γquad(c1)Γquad(c2)\mathcal R(c_1,c_2):=\frac{\Gamma_{\rm quad}(c_1c_2)}{\Gamma_{\rm quad}(c_1)\Gamma_{\rm quad}(c_2)}

on these square classes satisfies

R((−1)eλf,(−1)gλh)=(−1)eh+fg+fh,e,f,g,h∈{0,1},\mathcal{R}\left((-1)^e\lambda^f,(-1)^g\lambda^h\right)=(-1)^{eh+fg+fh},\qquad e,f,g,h\in\{0,1\},

and is a symmetric bicharacter. For a prime pp, each y∈O/(p)y \in\mathcal{O}/(p) has 1+χp3(y)1+\chi_p^3(y) square roots. Grouping by y=x2y=x^2 gives

Γquad(p)=1∣p∣∑y mod p(1+χp3(y))e(y/p)=γ3(p),\Gamma_{\rm quad}(p) =\frac1{|p|}\sum_{y\bmod p}\bigl(1+\chi_p^3(y)\bigr)e(y/p) =\gamma_3(p),

since ∑y mod pe(y/p)=0\sum_{y\bmod p}e(y/p)=0. The Chinese remainder theorem extends this equality to squarefree nn. Together with cubic reciprocity, it identifies the above R\mathcal{R} with the quotient χb(a)/χa(b)\chi_b(a)/\chi_a(b) for coprime primary a,ba,b. Furthermore χc(4)=(−2/c)3=(c/(−2))3\chi_c(4)=(-2/c)_3=(c/(-2))_3 is a character modulo 22. Consequently G(c)=χc(4)‾Γquad(c)G(c)=\overline{\chi_c(4)}\Gamma_{\rm quad}(c) factors through a fixed ray class group and satisfies (4.1), with R(c,c)=χc(−1)\mathcal{R}(c,c)=\chi_c(-1).

Proof of (4.7). In the fixed ray class group, the multiplicative relation for GG gives

G(a−1)=χa(−1)G(a)‾.G(a^{-1})=\chi_a(-1)\overline{G(a)}.

Using the symmetry and multiplicativity of R\mathcal{R}, we obtain

G(ba−1)=G(b)G(a−1)R(b,a−1)=χa(−1)G(a)‾G(b)R(a,b),G(ba^{-1})=G(b)G(a^{-1})\mathcal{R}(b,a^{-1})=\chi_a(-1)\overline{G(a)}G(b)\mathcal{R}(a,b),

which is (4.7).

Proof of (4.4)–(4.5). The second identity in (4.3) follows from the inverse-character Gauss identity: γ1(n)γ−1(n)=χn(−1)\gamma_1(n)\gamma_{-1}(n)=\chi_n(-1). It remains to prove (4.4)–(4.5). Multiplying the second identity in (4.2) by α(n)‾\overline{\alpha(n)} gives (4.4), and combining it with this inverse-character identity gives (4.5). □

The cubic theta transformation with fixed ray class twists

We prove the transformation formula of Proposition 6.2 for the completed sum T(X;Ψ)T(X;\Psi) in (5.3), together with the uniformity assertion of Lemma 6.3 and the coefficient and weight bounds of Lemma 6.4. The proof adapts the theta-transformation method of Dunn and Radziwiłł [7], which extends Patterson [38] and Yoshimoto [51]. The treatment of the fixed ray class twists and the uniformity assertions are supplied below. We use the setup preceding the proposition and the dependence of constants specified in these three statements: Ψ0\Psi_0, SS are fixed, while P\mathcal{P} is the varying finite set of primes outside SS and jpj_p are the exponents of their character factors in Ψ\Psi. Character powers follow the zero-extension convention preceding (6.5); in particular, χp0(x)=1p∤x\chi_p^0(x)=\mathbf1_{p\nmid x}. This convention also applies when we write χp(x)j\chi_p(x)^j.

Write w=(z,v)∈C×R>0w=(z,v)\in\mathbb{C}\times\mathbb{R}_{>0}, and use θ\theta from (6.1). The three sequences dσd_\sigma, σ∈{0,+,−}\sigma\in\{0,+,-\}, are defined by the Fourier expansions of θ(γσw)‾\overline{\theta(\gamma_\sigma w)} in (6.8), with the representatives (6.7). For the additive characters, write e˘(z)=exp⁡(2πi(z+zˉ))\breve e(z)=\exp(2\pi \mathrm i(z+\bar z)); thus e(z)=e˘(z/λ)e(z)=\breve{e}(z/\lambda), with ee as in (2.6).

Proof of Proposition 6.2 and Lemmas 6.3 and 6.4. Finite Fourier expansion. Define ϕ:O→C\phi:\mathcal{O}\to\mathbb{C} by

ϕ(n)={χn(λ)2Ψ0(n),n≡1(mod3), (n,S)=1,0,otherwise.\begin{aligned}\phi(n)= \begin{cases} \chi_n(\lambda)^2\Psi_0(n), & n\equiv1\pmod3,\ (n,S)=1,\\ 0, & \text{otherwise}. \end{cases}\end{aligned}

For Ψ(n)=Ψ0(n)∏p∈Pχpjp(n)\Psi(n)=\Psi_0(n)\prod_{p\in\mathcal{P}}\chi_p^{j_p}(n), define

ΘΨ(z,v)=∑ℓ∈λ−3Oℓ≠0τ(−ℓ)‾ϕ(λ3ℓ)(∏p∈Pχpjp(λ3ℓ))vK1/3(4π∣ℓ∣v)e˘(ℓz).\Theta_\Psi(z,v)= \sum_{\substack{\ell\in\lambda^{-3}\mathcal O\\\ell\ne0}} \overline{\tau(-\ell)}\phi(\lambda^3\ell) \Bigl(\prod_{p\in\mathcal P}\chi_p^{j_p}(\lambda^3\ell)\Bigr) vK_{1/3}(4\pi|\ell|v)\breve e(\ell z).

When Ψ=Ψk\Psi=\Psi_k, the multiplier of τ(−ℓ)‾\overline{\tau(-\ell)} is ϕk(λ3ℓ)\phi_k(\lambda^3\ell), so this definition recovers Θk\Theta_k in (6.2). We first write ΘΨ\Theta_\Psi as a finite sum of translates of θ‾\overline{\theta} and determine their reduced denominators. Choose once and for all a nonzero L∈OL\in\mathcal{O}, with prime divisors in SS, divisible by the conductor of Ψ0\Psi_0, every prime in SS, and sufficiently high powers of the primes above 22 and 33. The supplementary law of cubic reciprocity for λ\lambda [7], (1.5), which evaluates χn(λ)2=(λ/n)3\chi_n(\lambda)^2=(\lambda/n)_3, makes ϕ\phi periodic modulo LL. As before (6.4), define

ϕ^(h0)=1NK/Q(L)∑x mod Lϕ(x)e(−h0x/L),h0∈O/(L).\widehat\phi(h_0)=\frac1{\mathrm N_{K/\mathbb Q}(L)}\sum_{x\bmod L}\phi(x)e(-h_0x/L), \qquad h_0\in\mathcal O/(L).

Since ∣ϕ∣≤1|\phi|\le1, we have ∣ϕ^(h0)∣≤1|\widehat{\phi}(h_0)|\le1.

For p∈Pp\in\mathcal{P} and 0≤j≤50\le j\le5, we expand χpj\chi_p^j on O/(p)\mathcal{O}/(p) in additive characters. Its Fourier coefficients are defined, for hp∈O/(p)h_p\in\mathcal{O}/(p), by

Cp,j(hp):=1NK/Q(p)∑y mod pχpj(y)e(−hpy/p).C_{p,j}(h_p):=\frac1{\mathrm N_{K/\mathbb Q}(p)}\sum_{y\bmod p}\chi_p^j(y)e(-h_py/p).

Finite Fourier inversion then gives, for x∈Ox\in\mathcal{O},

χpj(x)=∑hp∈O/(p)Cp,j(hp)e(hpx/p).\chi_p^j(x)=\sum_{h_p\in\mathcal O/(p)}C_{p,j}(h_p)e(h_px/p).

Multiplying (A.3) over the primes in P\mathcal{P} gives, for x∈Ox\in\mathcal{O},

∏p∈Pχpjp(x)=∑hp∈O/(p)p∈P(∏p∈PCp,jp(hp))e ⁣(x∑p∈Phpp).\prod_{p\in\mathcal P}\chi_p^{j_p}(x) =\sum_{\substack{h_p\in\mathcal O/(p)\\p\in\mathcal P}} \Bigl(\prod_{p\in\mathcal P}C_{p,j_p}(h_p)\Bigr) e\!\Bigl(x\sum_{p\in\mathcal P}\frac{h_p}{p}\Bigr).

Each summand is indexed by a tuple (hp)p∈P(h_p)_{p\in\mathcal{P}}. At x=λ3ℓx=\lambda^3\ell, its additive character gives the shift λ2∑p∈Php/p\lambda^2\sum_{p\in\mathcal{P}}h_p/p of θ‾\overline\theta. Writing h=(h0,(hp)p∈P)\boldsymbol h=(h_0,(h_p)_{p\in\mathcal P}) with h0∈O/(L)h_0\in\mathcal{O}/(L), we claim that

ΘΨ(z,v)=∑hcF(h)θ(z+zh,v)‾,\Theta_\Psi(z,v)=\sum_{\boldsymbol h} c_{\mathrm F}(\boldsymbol h)\overline{\theta(z+z_{\boldsymbol h},v)},

where

zh=λ2(h0/L+∑p∈Php/p),cF(h)=ϕ^(h0)∏p∈PCp,jp(hp).z_{\boldsymbol h}=\lambda^2\Bigl(h_0/L+\sum_{p\in\mathcal P}h_p/p\Bigr), \qquad c_{\mathrm F}(\boldsymbol h)=\widehat\phi(h_0)\prod_{p\in\mathcal P} C_{p,j_p}(h_p).

To verify this identity, expand the fixed factor ϕ\phi as well. The finite Fourier expansions identify all nonzero Fourier coefficients, and the constant terms of the translates cancel because

∑hcF(h)=ϕ(0)∏p∈Pχpjp(0)=0.\sum_{\boldsymbol h} c_{\mathrm F}(\boldsymbol h)=\phi(0)\prod_{p\in\mathcal P}\chi_p^{j_p}(0)=0.

For j≢0(mod6)j\not\equiv0\pmod6, we use γj(p)\gamma_j(p) from (2.6). Changing yy to −y-y gives

1NK/Q(p)∑y mod pχp(y)je(−y/p)=χp(−1)jγj(p),∣γj(p)∣=1.\frac1{\sqrt{\mathrm N_{K/\mathbb Q}(p)}}\sum_{y\bmod p}\chi_p(y)^j e(-y/p) =\chi_p(-1)^j\gamma_j(p),\qquad |\gamma_j(p)|=1.

For hp≠0h_p\ne0, substitute y=hp−1uy=h_p^{-1}u in the definition of Cp,j(hp)C_{p,j}(h_p). For hp=0h_p=0 and j≠0j\ne0, use character orthogonality; for j=0j=0, sum the additive character over nonzero residues directly. These calculations give

Cp,j(hp)={NK/Q(p)−1/2χp(−1)jγj(p)χp(hp)−j,j≠0, hp≠0,0,j≠0, hp=0,−NK/Q(p)−1,j=0, hp≠0,1−NK/Q(p)−1,j=0, hp=0.\begin{aligned} C_{p,j}(h_p)= \begin{cases} \mathrm N_{K/\mathbb Q}(p)^{-1/2}\chi_p(-1)^j\gamma_j(p)\chi_p(h_p)^{-j}, &j\ne0,\ h_p\ne0,\\ 0, &j\ne0,\ h_p=0,\\ -\mathrm N_{K/\mathbb Q}(p)^{-1},&j=0,\ h_p\ne0,\\ 1-\mathrm N_{K/\mathbb Q}(p)^{-1},&j=0,\ h_p=0. \end{cases} \end{aligned}

For the tuple h\boldsymbol h, define its set of active primes by A(h):={p∈P:hp≠0}\mathcal A(\boldsymbol h):=\{p\in\mathcal P:h_p\ne0\}. By (A.5), cF(h)≠0c_{\mathrm F}(\boldsymbol h)\ne0 implies {p∈P:jp≠0}⊆A(h)\{p\in\mathcal P:j_p\ne0\}\subseteq\mathcal A(\boldsymbol h). Thus only primes with jp=0j_p = 0 can be inactive, as asserted in the proposition. Put r=∏p∈A(h)pr=\prod_{p\in\mathcal A(\boldsymbol h)}p, and let c0c_0 be the reduced denominator of λ2h0/L\lambda^2 h_0/L. At each active prime pp,

vp(λ2hp/p)=−1,vp(zh−λ2hp/p)≥0,v_p(\lambda^2h_p/p)=-1, \qquad v_p(z_{\boldsymbol h}-\lambda^2h_p/p)\ge0,

so the factor pp cannot cancel from the denominator. At primes dividing LL, all terms λ2hp/p\lambda^2 h_p/p are integral. Thus the reduced denominator is c=c0rc=c_0r up to a unit, including when (h0,L)≠1(h_0,L)\ne1. The finite Fourier expansion has therefore expressed ΘΨ\Theta_{\Psi} as OL(2∣P∣)O_L(2^{|\mathcal{P}|}) groups of translates θ(z+zh,v)‾\overline{\theta(z+z_{\boldsymbol h},v)}, with shifts zhz_{\boldsymbol h} and coefficients cF(h)c_{\mathrm F}(\boldsymbol h) given by (A.4). The groups are indexed by h0∈O/(L)h_0\in\mathcal{O}/(L) and the active set A\mathcal{A}. Each group has a common reduced denominator c=c0∏p∈Apc=c_0\prod_{p\in\mathcal{A}}p up to a unit. For each translate θ(z+zh,v)‾\overline{\theta(z+z_{\boldsymbol h},v)}, we next identify which of the three Fourier expansions in (6.8) will be used after changing coordinates at zhz_{\boldsymbol h}.

Reduced denominators and cusp coefficients. We next express each translate θ(z+zh,v)‾\overline{\theta(z+z_{\boldsymbol h},v)} using one of the three cusp expansions identified above. Our representatives γ0,γ+,γ−\gamma_0,\gamma_+,\gamma_- correspond to γ1,γ10,γ19\gamma_1,\gamma_{10},\gamma_{19}, respectively, in the numbering of [7], which follows [38]. We allow (h0,L)≠1(h_0,L)\ne1 and keep the reduced denominator c0c_0 of λ2h0/L\lambda^2h_0/L in the calculation; this is the extension beyond [7].

Set M=λ12L4M=\lambda^{12}L^4. Since (p,M)=1(p,M)=1 for p∈A(h)p\in\mathcal A(\boldsymbol h), the Chinese remainder theorem lets us choose representatives hp∈Oh_p\in\mathcal{O} with hp≡0(modM2)h_p\equiv0\pmod{M^2}. Changing representatives modulo pp changes zhz_{\boldsymbol h} by an element of λ2O=3O\lambda^2\mathcal{O}=3\mathcal{O}, under which θ\theta is periodic. To track dependence on the primes in rr, restrict rr to a residue class modulo M2M^2. For fixed h0h_0, active set, and this residue class, write zh=a/cz_{\boldsymbol h}=a/c in lowest terms, normalizing a≡1(mod3)a\equiv1\pmod3 if λ∣c\lambda\mid c, and c≡1(mod3)c\equiv1\pmod3 otherwise. Since

a=λ2c(h0/L+∑p∣rhp/p),a=\lambda^2c\left(h_0/L+\sum_{p\mid r}h_p/p\right),

these choices fix the residue of aa modulo Mc0Mc_0. Since (a,c)=1(a,c)=1 and (r,M)=1(r,M)=1, the Chinese remainder theorem gives δ′∈O\delta'\in\mathcal{O} satisfying the following congruences, where qq ranges over prime divisors of the indicated elements:

aδ′≡1(modqvq(Mc0))(q∣c0),a\delta'\equiv1 \pmod{q^{v_q(Mc_0)}} \qquad(q\mid c_0),
δ′≡0(modqvq(M))(q∣M, q∤c0),\delta'\equiv0 \pmod{q^{v_q(M)}} \qquad(q\mid M,\ q\nmid c_0),
aδ′≡1(modr).a\delta'\equiv1 \pmod r.

Put

bg=aδ′−1c,g=(abgcδ′)∈SL2(O).\begin{aligned}b_g=\frac{a\delta'-1}{c},\qquad g=\begin{pmatrix}a&b_g\\c&\delta'\end{pmatrix}\in\mathrm{SL}_2(\mathcal{O}).\end{aligned}

The congruences for δ′\delta' give

bg≡{0,q∣c0,−c−1,q∣M, q∤c0,(modqvq(M)).\begin{aligned} b_g\equiv \begin{cases} 0,&q\mid c_0,\\ -c^{-1},&q\mid M,\ q\nmid c_0, \end{cases} \pmod{q^{v_q(M)}}. \end{aligned}

To choose the cusp expansion at a/ca/c, put

H={I,3∣c,(10u01),vλ(c)=1,u0∈{λ,−λ},u0≡c(mod3),(u0−110),(c,λ)=1,u0≡a(mod3).\begin{aligned} H= \begin{cases} I,&3\mid c,\\[2pt] \bigl(\begin{smallmatrix}1&0\\u_0&1\end{smallmatrix}\bigr), &v_\lambda(c)=1,\quad u_0\in\{\lambda,-\lambda\},\quad u_0\equiv c\pmod3,\\[2pt] \bigl(\begin{smallmatrix}u_0&-1\\1&0\end{smallmatrix}\bigr), &(c,\lambda)=1,\quad u_0\equiv a\pmod3. \end{cases} \end{aligned}

In the last case choose u0∈Ou_0\in\mathcal{O} from a fixed set of representatives modulo 33. In each case, put g1=gH−1g_1=gH^{-1}. Then g=g1Hg=g_1H and g1≡I(mod3)g_1\equiv I\pmod3, as in [7]. The invariances

θ(γw)=θ(w)(γ∈SL2(Z)),θ(z+t,v)=θ(z,v)(t∈Z+3O)\theta(\gamma w)=\theta(w)\quad(\gamma\in\mathrm{SL}_2(\mathbb{Z})),\qquad \theta(z+t,v)=\theta(z,v)\quad(t\in\mathbb{Z}+3\mathcal{O})

reduce θ‾(Hw)\overline{\theta}(Hw) to one of the three functions θ‾(γσw)\overline{\theta}(\gamma_\sigma w), σ∈{0,−,+}\sigma\in\{0,-,+\}, with the matrices γσ\gamma_\sigma from (6.7). These representatives give the infinity expansion and the two additional cusp expansions needed here: translating by ω\omega or −ω-\omega, then applying inversion, gives the representatives labeled −- and ++, respectively [7], (5.9)–(5.15), Appendix A. Write tσ(ℓ)t_\sigma(\ell) for the coefficient of vK1/3(4π∣ℓ∣v)e˘(ℓz)v K_{1/3}(4\pi|\ell|v)\breve{e}(\ell z) in θ(γσw)\theta(\gamma_\sigma w). Let τ1,τ2:λ−4O∖{0}→C\tau_1,\tau_2:\lambda^{-4}\mathcal{O}\setminus\{0\}\to\mathbb{C} be the coefficient sequences defined in [7], (5.13), (5.14). Patterson’s cusp calculation [38] then gives

t0(ℓ)=τ(ℓ),t−(ℓ)=ω2τ1(ω2ℓ)e˘(ℓ),t+(ℓ)=ωτ2(ωℓ)e˘(ℓ).t_0(\ell)=\tau(\ell),\qquad t_-(\ell)=\omega^2\tau_1(\omega^2\ell)\breve{e}(\ell),\qquad t_+(\ell)=\omega\tau_2(\omega\ell)\breve{e}(\ell).

Complex conjugation changes the Fourier frequency from ℓ\ell to −ℓ-\ell. Thus the coefficients in the normalization (6.8) are

dσ(ℓ)=tσ(−ℓ)‾,σ∈{0,+,−}.d_\sigma(\ell)=\overline{t_\sigma(-\ell)},\qquad \sigma\in\{0,+,-\}.

For indices ℓ=uλmnb3\ell=u\lambda^mnb^3 as in (6.10), the coefficient magnitudes satisfy

∣t0(ℓ)∣≤{3k/2+2∣b∣,m=3k−4, k≥1,3k/2+5/2∣b∣,m=3k−3, k≥0.\begin{aligned}|t_0(\ell)|\le \begin{cases} 3^{k/2+2}|b|, & m=3k-4,\ k\ge1,\\ 3^{k/2+5/2}|b|, & m=3k-3,\ k\ge0. \end{cases}\end{aligned}

The m=−4m=-4 coefficients of t−t_- and t+t_+ have magnitude at most 9∣b∣9|b|. For the chosen HH, let σ∈{0,+,−}\sigma\in\{0,+,-\} index the corresponding expansion in (6.8). For each translate θ‾(z+a/c,v)\overline{\theta}(z+a/c,v), we have identified the Fourier expansion in zz of θ‾(H(z,v))\overline{\theta}(H(z,v)) as one of the three expansions in (6.8):

θ‾(H(z,v))=35/221σ=0v2/3+∑0≠ℓ∈λ−4Odσ(ℓ)vK1/3(4π∣ℓ∣v)e˘(ℓz).\overline{\theta}(H(z,v))=\frac{3^{5/2}}{2}\mathbf{1}_{\sigma=0}v^{2/3}+\sum_{0\ne\ell\in\lambda^{-4}\mathcal{O}}d_\sigma(\ell)vK_{1/3}(4\pi|\ell|v)\breve{e}(\ell z).

These coefficients satisfy the support restriction (6.10) and bound (6.11), proving the coefficient assertions of Lemma 6.4. The next step uses g=g1Hg=g_1H to express the original translate through this Fourier series evaluated at g−1(z+a/c,v)g^{-1}(z+a/c,v), and computes the accompanying multiplier and additive phase. The matrix gg maps ∞\infty to a/ca/c, which is why this is called an expansion at the cusp a/ca/c.

The local character transformation. The matrix factorization g=g1Hg=g_1H lets us apply the automorphy law

θ(g1w)=κ(g1)θ(w),κ(g1)=(c1/a1)3,g1=(a1b1c1d1)≡I(mod3),\begin{aligned} \theta(g_1w)=\kappa(g_1)\theta(w),\qquad \kappa(g_1)=(c_1/a_1)_3,\qquad g_1=\Bigl(\begin{matrix}a_1&b_1\\c_1&d_1\end{matrix}\Bigr)\equiv I\pmod3, \end{aligned}

where κ\kappa is Kubota’s cubic character [7]. We will combine its conjugate with the finite Fourier coefficients Cp,j(hp)C_{p,j}(h_p) in (A.5) to obtain the factors Bp,jB_{p,j} of (6.5). For the translated theta function, the automorphy law gives

θ(z+a/c,v)‾=κ(g1)‾ θ(Hg−1(z+a/c,v))‾,\overline{\theta(z+a/c,v)} =\overline{\kappa(g_1)}\,\overline{\theta\bigl(Hg^{-1}(z+a/c,v)\bigr)},

where

g−1(z+a/c,v)=(−δ′c−zˉc2(v2+∣z∣2),vNK/Q(c)(v2+∣z∣2)).g^{-1}(z+a/c,v) =\Bigl(-\frac{\delta'}c-\frac{\bar z}{c^2(v^2+|z|^2)}, \frac{v}{\mathrm N_{K/\mathbb Q}(c)(v^2+|z|^2)}\Bigr).

At z=0z=0, the term with Fourier index ℓ\ell in the expansion of θ‾(Hw)\overline{\theta}(Hw) therefore acquires the phase e˘(−δ′ℓ/c)\breve{e}(-\delta'\ell/c). We claim that the multiplier is given by

κ(g1)={(c0/a)3(a/r)3,3∣c,(−u0/(a−u0bg))3 ((c0/u0)/a)3(a/r)3,vλ(c)=1,(a/c0)3(a/r)3,(c,λ)=1.\begin{aligned} \kappa(g_1)= \begin{cases} (c_0/a)_3(a/r)_3,&3\mid c,\\ (-u_0/(a-u_0b_g))_3\,((c_0/u_0)/a)_3(a/r)_3, &v_\lambda(c)=1,\\ (a/c_0)_3(a/r)_3,&(c,\lambda)=1. \end{cases} \end{aligned}

To verify (A.10), use the determinant equation and cubic reciprocity. For the case vλ(c)=1v_\lambda(c)=1, the required congruences are

a(c−u0δ′)≡−u0(moda−u0bg),bgc≡−1(moda).a(c-u_0\delta')\equiv-u_0\pmod{a-u_0b_g},\qquad b_gc\equiv-1\pmod a.

For (c,λ)=1(c,\lambda)=1, the congruence

−bgc=1−aδ′≡1(mod9)-b_g c = 1-a\delta' \equiv1 \pmod{9}

removes the supplementary factors; reciprocity at the remaining primes then gives (δ′/(−bgc))3=1(\delta'/(-b_g c))_3=1. The factors other than (a/r)3(a/r)_3 depend only on the fixed residues at primes in SS. We may therefore write

κ(g1)=κ0(a/r)3=κ0∏p∣rχp(a)2,∣κ0∣=1,\kappa(g_1)=\kappa_0(a/r)_3=\kappa_0\prod_{p\mid r}\chi_p(a)^2,\qquad|\kappa_0|=1,

where κ0\kappa_0 is fixed once h0h_0, the active set, and r mod M2r\bmod M^2 are fixed.

Fix h0h_0 and the active set, so that the denominator c=c0rc=c_0r is fixed while the nonzero residues hph_p vary. Put D0=λ3c0D_0=\lambda^3c_0. For each active prime pp, define in O/(p)\mathcal{O}/(p)

σp=λ2c/p,ϵp=−((λ3c/p)σp)−1.\sigma_p=\lambda^2c/p,\qquad\epsilon_p=-((\lambda^3c/p)\sigma_p)^{-1}.

The inverse in ϵp\epsilon_p is taken in the field O/(p)\mathcal{O}/(p); it exists because rr is squarefree and p∤λc0p\nmid\lambda c_0. The expression for aa gives a≡σphp(modp)a\equiv\sigma_ph_p\pmod p. For a dual Fourier index ℓ∈λ−4O\ell\in\lambda^{-4}\mathcal{O}, the additive form of the Chinese remainder theorem yields

e˘(−δ′ℓ/c)=ψ(λ4ℓ)∏p∣re(ϵphp−1λ4ℓ/p),\breve e(-\delta'\ell/c)=\psi(\lambda^4\ell)\prod_{p\mid r}e(\epsilon_ph_p^{-1}\lambda^4\ell/p),

where

ψ(λ4ℓ):=e(−δ0′r−1λ4ℓ/D0),δ0′=δ′ mod D0.\psi(\lambda^4\ell):=e(-\delta'_0r^{-1}\lambda^4\ell/D_0),\qquad\delta'_0=\delta'\bmod D_0.

The residues δ0′\delta'_0 and r−1 mod D0r^{-1}\bmod D_0 are fixed by the choices above, so ψ\psi ranges over a fixed finite family of additive characters of O\mathcal{O}. We claim the local transformation identity

∑hp≠0Cp,j(hp)χp(a)−2e(ϵphp−1λ4ℓ/p)=χp(σp)−2ωp,jBp,j(λ4ℓ),\sum_{h_p\ne0}C_{p,j}(h_p)\chi_p(a)^{-2}e(\epsilon_ph_p^{-1}\lambda^4\ell/p)=\chi_p(\sigma_p)^{-2}\omega_{p,j}B_{p,j}(\lambda^4\ell),

where Bp,jB_{p,j} is defined in (6.5), and

ωp,j={χp(−1)jγj(p)γj+2(p)χp(ϵp)−j−2,j≠0,4,γ4(p),j=4,−γ2(p)χp(ϵp)−2,j=0 active.\begin{aligned} \omega_{p,j}= \begin{cases} \chi_p(-1)^j\gamma_j(p)\gamma_{j+2}(p)\chi_p(\epsilon_p)^{-j-2},&j\ne0,4,\\ \gamma_4(p),&j=4,\\ -\gamma_2(p)\chi_p(\epsilon_p)^{-2},&j=0\ {\rm active}. \end{cases} \end{aligned}

The indices of γj(p)\gamma_j(p) are read modulo six; in the second case, χp(−1)4=1\chi_p(-1)^4=1. In each case ∣ωp,j∣=1|\omega_{p,j}|=1. To prove (A.12), combine (A.5), the conjugate of (A.10), and (A.11). Setting y=hp−1y=h_p^{-1} changes the character exponent as follows:

χp(hp)−jχp(a)−2=χp(σp)−2χp(hp)−j−2=χp(σp)−2χp(y)j+2,y=hp−1.\chi_p(h_p)^{-j}\chi_p(a)^{-2}=\chi_p(\sigma_p)^{-2}\chi_p(h_p)^{-j-2}=\chi_p(\sigma_p)^{-2}\chi_p(y)^{j+2},\qquad y=h_p^{-1}.

For j≠0,4j\ne0,4, the character χpj+2\chi_p^{j+2} is nontrivial, and its Gauss sum gives the first case of Bp,jB_{p,j}. For j=4j=4, it is trivial on nonzero residues; after rescaling by ϵp≠0\epsilon_p\ne0, the sum is

∑y≠0e(λ4ℓy/p)=−1+NK/Q(p)1p∣λ4ℓ.\sum_{y\ne0}e(\lambda^4\ell y/p) =-1+\mathrm N_{K/\mathbb Q}(p)\mathbf1_{p\mid\lambda^4\ell}.

For active j=0j=0, the coefficient −NK/Q(p)−1-\mathrm N_{K/\mathbb Q}(p)^{-1} combines with a cubic Gauss sum to give

Bp,0(λ4ℓ)=NK/Q(p)−1/2χp(λ4ℓ)−2,B_{p,0}(\lambda^4\ell) =\mathrm N_{K/\mathbb Q}(p)^{-1/2}\chi_p(\lambda^4\ell)^{-2},

with the unit factors included in ωp,0\omega_{p,0}. This proves (A.12) in every case. An inactive prime has j=0j=0 and contributes only the scalar 1−NK/Q(p)−11-\mathrm N_{K/\mathbb Q}(p)^{-1}. Consequently, for each fixed h0h_0 and active set A\mathcal{A}, combining the Fourier coefficients dσ(ℓ)d_{\sigma}(\ell) of θ(H(z,v))‾\overline{\theta(H(z,v))} with the sums over hp≠0h_p \ne0 gives

dσ(ℓ)ψ(λ4ℓ)∏p∈ABp,jp(λ4ℓ),d_\sigma(\ell)\psi(\lambda^4\ell) \prod_{p\in\mathcal A}B_{p,j_p}(\lambda^4\ell),

up to a scalar independent of ℓ\ell. These are the arithmetic factors in the dual sum (6.9). In particular, Bp,1(λ4ℓ)=χp(λ4ℓ)3B_{p,1}(\lambda^4\ell)=\chi_p(\lambda^4\ell)^3 is the quadratic factor used in (6.18).

The archimedean transform. We now derive the transformed weight V∗♯V_*^\sharp in (6.6) and the scalar multiplying the dual sum. For Re⁡s>1\operatorname{Re}s>1, introduce the Dirichlet series associated with the completed sum (5.3):

T(s,Ψ)=(∑n∈On≡1 (3)(n,S)=1∗α(n)‾γ2(n)Ψ(n)NK/Q(n)−s)(∑b∈Ob≡1 (3)(b,S)=1α(b)‾ 3Ψ(b)3NK/Q(b)−3s+1/2).\mathcal T(s,\Psi)= \Bigl(\sum_{\substack{n\in\mathcal O\\n\equiv1\ (3)\\(n,S)=1}}^* \overline{\alpha(n)}\gamma_2(n)\Psi(n)\mathrm N_{K/\mathbb Q}(n)^{-s}\Bigr) \Bigl(\sum_{\substack{b\in\mathcal O\\b\equiv1\ (3)\\(b,S)=1}}\overline{\alpha(b)}^{\,3} \Psi(b)^3\mathrm N_{K/\mathbb Q}(b)^{-3s+1/2}\Bigr).

With V∗V_* from (5.3) and its Mellin transform defined before (6.6), Mellin inversion gives

T(X;Ψ)=12πi∫(σ)V^∗(s−12)T(s,Ψ)Xs−1/2 ds,σ>1.T(X;\Psi)=\frac1{2\pi \mathrm i}\int_{(\sigma)} \widehat V_*(s-\tfrac12)\mathcal T(s,\Psi)X^{s-1/2}\,d s, \qquad \sigma>1.

The identity (6.3), with Ψ\Psi in place of Ψk\Psi_k, identifies these coefficients with those of ΘΨ\Theta_\Psi after inserting the factor α(nb3)‾\overline{\alpha(nb^3)}. We now compute the normalization relating T(s,Ψ)\mathcal{T}(s,\Psi) to the Mellin transform of the derivative of ΘΨ\Theta_\Psi.

For z=x+iyz=x+\mathrm i y, use ∂zˉ=(∂x+i∂y)/2\partial_{\bar z}=(\partial_x+\mathrm i\partial_y)/2 and ∂z=(∂x−i∂y)/2\partial_z=(\partial_x-\mathrm i\partial_y)/2. Differentiation in z‾\overline{z} supplies a factor ℓ‾\overline{\ell} in each Fourier mode. Combined with the factor ∣ℓ∣−1|\ell|^{-1} from the Bessel integral below, this produces the required angular factor α(ℓ)‾\overline{\alpha(\ell)}. We therefore define, initially for Re⁡s>1\operatorname{Re}s>1, the Mellin transform

J(s)=∫0∞∂zˉΘΨ(z,v)∣z=0v2s−1 dv.\mathcal J(s)=\int_0^\infty \partial_{\bar z}\Theta_\Psi(z,v)\Bigr|_{z=0} v^{2s-1}\,d v.

For each nonzero Fourier mode, differentiation and the Mellin integral for K1/3K_{1/3} (see [35], (10.43.19)) give the following formulas:

∂zˉe˘(ℓz)=2πiℓˉ e˘(ℓz),\partial_{\bar z}\breve e(\ell z) =2\pi \mathrm i\bar\ell\,\breve e(\ell z),
∫0∞v2sK1/3(4π∣ℓ∣v) dv=22s−1Γ(s+1/3)Γ(s+2/3)(4π∣ℓ∣)2s+1Re⁡s>0.\int_0^\infty v^{2s}K_{1/3}(4\pi|\ell|v)\,d v =\frac{2^{2s-1}\Gamma(s+1/3)\Gamma(s+2/3)} {(4\pi|\ell|)^{2s+1}} \qquad \operatorname{Re} s>0.

For Re⁡s>1\operatorname{Re}s>1, the coefficient formula (2.10) and (A.15) show that the sum of the integrals of the absolute values is finite. We may therefore integrate the differentiated Fourier series term by term.

For ℓ=λ−3nb3\ell=\lambda^{-3}nb^3, the phase and scale simplify to

iα(ℓ)‾∣ℓ∣−2s=27sα(n)‾ α(b)‾ 3NK/Q(n)−sNK/Q(b)−3s.\mathrm i\overline{\alpha(\ell)}|\ell|^{-2s} =27^s\overline{\alpha(n)}\, \overline{\alpha(b)}^{\,3}\mathrm N_{K/\mathbb Q}(n)^{-s}\mathrm N_{K/\mathbb Q}(b)^{-3s}.

Combining this with (2.10) gives

J(s)=35/24(27(2π)2)sΓ(s+1/3)Γ(s+2/3)T(s,Ψ).\mathcal{J}(s)=\frac{3^{5/2}}{4}\left(\frac{27}{(2\pi)^2}\right)^s\Gamma(s+1/3)\Gamma(s+2/3)\mathcal{T}(s,\Psi).

We next show that J(s)\mathcal{J}(s) is entire and use (A.16) to continue T(s,Ψ)\mathcal{T}(s,\Psi).

Write (z′,v′)=g−1(z+a/c,v)(z',v')=g^{-1}(z+a/c,v). Differentiating (A.9) at z=0z=0 gives

∂z′∂zˉ∣z=0=−1c2v2,∂z′‾∂zˉ∣z=0=0,∂v′∂zˉ∣z=0=0.\frac{\partial z'}{\partial\bar z}\bigg|_{z=0} =-\frac1{c^2v^2}, \qquad \frac{\partial\overline{z'}}{\partial\bar z}\bigg|_{z=0}=0, \qquad \frac{\partial v'}{\partial\bar z}\bigg|_{z=0}=0.

Thus the derivative in cusp coordinates is −(cv)−2∂z′-(cv)^{-2}\partial_{z'}, with no height-derivative term. The defining Fourier series (A.2) for ΘΨ\Theta_\Psi controls v→∞v \to\infty. For v→0v \to0, use (A.8) and the cusp expansions (6.8). In each expansion the horizontal derivative removes the constant mode. The remaining modes decay exponentially as v→∞v \to\infty; at v→0v \to0, the transformed height 1/(NK/Q(c)v)1/(\mathrm N_{K/\mathbb Q}(c)v) tends to infinity, giving exponential decay in 1/v1/v. Thus the integral defining J(s)\mathcal{J}(s) converges for every ss and defines an entire function. Equation (A.16), after division by its gamma factors, also continues T(s,Ψ)\mathcal{T}(s,\Psi) to an entire function.

For the term indexed by h\boldsymbol h in (A.4), write ch,δh′,g1,h,Hhc_{\boldsymbol h},\delta'_{\boldsymbol h},g_{1,\boldsymbol h},H_{\boldsymbol h} for the corresponding choices above. Let σh\sigma_{\boldsymbol h} be the cusp index determined by HhH_{\boldsymbol h}. Define its cusp Mellin transform by

Jh∨(s)=∫0∞∂z{θ(Hh(z,v))‾}∣z=−δh′/chv2s−1 dv.\mathcal J_{\boldsymbol h}^\vee(s)=\int_0^\infty \partial_z\bigl\{\overline{\theta(H_{\boldsymbol h}(z,v))}\bigr\} \Bigr|_{z=-\delta'_{\boldsymbol h}/c_{\boldsymbol h}}v^{2s-1}\,d v.

Substituting v↦(NK/Q(ch)v)−1v\mapsto(\mathrm N_{K/\mathbb Q}(c_{\boldsymbol h})v)^{-1} in each translated term of J(s)\mathcal{J}(s), using (A.8), gives

J(s)=−∑hcF(h)κ(g1,h)‾ α(ch)‾ 2NK/Q(ch)1−2sJh∨(1−s).\mathcal J(s)=-\sum_{\boldsymbol h} c_{\mathrm F}(\boldsymbol h)\overline{\kappa(g_{1,\boldsymbol h})}\, \overline{\alpha(c_{\boldsymbol h})}^{\,2}\mathrm N_{K/\mathbb Q}(c_{\boldsymbol h})^{1-2s}\mathcal J_{\boldsymbol h}^\vee(1-s).

Here cF(h)c_{\mathrm F}(\boldsymbol h) is the coefficient in (A.4). The cusp coefficients dσh(ℓ)d_{\sigma_{\boldsymbol h}}(\ell) and (A.15) give, for Re⁡s>1\operatorname{Re}s>1,

Jh∨(s)=i Γ(s+1/3)Γ(s+2/3)4(2π)2s∑0≠ℓ∈λ−4Odσh(ℓ)α(ℓ)NK/Q(ℓ)se˘(−δh′ℓ/ch).\mathcal J_{\boldsymbol h}^\vee(s)=\frac{\mathrm i\,\Gamma(s+1/3)\Gamma(s+2/3)}{4(2\pi)^{2s}} \sum_{0\ne\ell\in\lambda^{-4}\mathcal O} \frac{d_{\sigma_{\boldsymbol h}}(\ell)\alpha(\ell)}{\mathrm N_{K/\mathbb Q}(\ell)^s} \breve e(-\delta'_{\boldsymbol h}\ell/c_{\boldsymbol h}).

The Dirichlet series T(s,Ψ)\mathcal{T}(s,\Psi) converges absolutely for Re⁡s>1\operatorname{Re}s>1, while the series in (A.19) at 1−s1-s converges absolutely for Re⁡s<0\operatorname{Re}s<0. The functional equation (A.18) and Stirling’s formula [22] therefore give polynomial bounds in ∣Im⁡s∣|\operatorname{Im}s| for T(s,Ψ)\mathcal{T}(s,\Psi) on both sides of the strip 0≤Re⁡s≤10\leq\operatorname{Re}s\leq1. Splitting the integral defining J(s)\mathcal{J}(s) at v=1v=1 gives finite order; Phragmén–Lindelöf [22] then gives the same type of bound inside the strip (compare the arguments of Dunn and Radziwiłł in [7]). Together with the rapid decay of V^∗\widehat V_* on vertical lines, these bounds justify moving the ss-contour in (A.13) to Re⁡s<0\operatorname{Re}s<0. No poles are crossed, since T(s,Ψ)\mathcal{T}(s,\Psi) is entire. Setting t=12−st=\frac{1}{2}-s then gives a line Re⁡t>1/2\operatorname{Re} t>1/2, where the dual coefficient series converges absolutely.

For fixed h0h_0 and active set, insert (A.19) into (A.18) and use (A.12) to sum over the nonzero hph_p. After the normalization in (A.16), the scalar CC in (6.9) for this group of translates is

C=−i81α(c)‾ 2ϕ^(h0)κ0‾∏p inactive(1−NK/Q(p)−1)∏p activeχp(σp)−2ωp,jp,∣C∣≤1/81.C=-\frac{\mathrm i}{81}\overline{\alpha(c)}^{\,2} \widehat\phi(h_0)\overline{\kappa_0} \prod_{p\ {\rm inactive}}(1-\mathrm N_{K/\mathbb Q}(p)^{-1}) \prod_{p\ {\rm active}}\chi_p(\sigma_p)^{-2}\omega_{p,j_p}, \qquad |C|\le1/81 .

Dividing (A.16) by its gamma factors and setting t=12−st=\frac{1}{2}-s produces the gamma quotient in (6.6). Its numerator gamma factors have their first pole at t=−5/6t=-5/6, and its reciprocal denominator gamma factors are entire. Thus no poles lie between the current contour Re⁡t>1/2\operatorname{Re}t>1/2 and Re⁡t=0\operatorname{Re}t=0. The rapid decay of V^∗\widehat V_*, together with Stirling’s formula, allows us to shift each kernel contour to Re⁡t=0\operatorname{Re}t=0. This gives the weight V∗♯V_*^\sharp defined in (6.6), evaluated at NK/Q(ℓ)X/NK/Q(c)2\mathrm N_{K/\mathbb Q}(\ell)X/\mathrm N_{K/\mathbb Q}(c)^2. Together with (A.12), this gives the dual sum (6.9).

For k0k_0 as in Lemma 6.3, fix h0h_0, the factors of Ψ\Psi other than χk0\chi_{k_0}, and the active/inactive choices at their primes. Every prime dividing k0k_0 has exponent jp=1j_p=1 and is therefore active, so

r=k0∏p∈Afixp.r=k_0\prod_{p\in\mathcal A_{\rm fix}}p.

Thus r/k0r/k_0 is fixed, and fixing k0k_0 mod M2M^2 fixes rr mod M2M^2. The choices of H,c0H,c_0 and the residues in (A.11) are therefore fixed in each such class. Thus d=dσd=d_\sigma, ψ\psi, and c0c_0 have precisely the asserted dependence on k0k_0.

It remains to bound the transformed weight. The first numerator pole of the gamma quotient in (6.6) is at t=−5/6t=-5/6. We may therefore shift the kernel contour to Re⁡t=−1/4\operatorname{Re} t=-1/4 for 0<x≤10<x\le1, obtaining the factor x1/4x^{1/4}, and to Re⁡t=A\operatorname{Re} t=A for x≥1x\ge1, obtaining x−Ax^{-A}. Each application of x∂xx\partial_x introduces a factor −t-t. Stirling’s formula on −1/4≤Re⁡t≤A-1/4\le\operatorname{Re} t\le A then gives, for every A>0A>0 and j≥0j\ge0,

∣(x∂x)jV∗♯(x)∣≪A,jmin⁡{x1/4,(1+x)−A}×sup⁡−A≤η≤1/4∫R(1+∣u∣)⌈4A⌉+j+2∣V^∗(η+iu)∣ du.\begin{aligned} |(x\partial_x)^jV_*^\sharp(x)| &\ll_{A,j}\min\{x^{1/4},(1+x)^{-A}\}\\ &\quad\times\sup_{-A\le\eta\le1/4}\int_{\mathbb R} (1+|u|)^{\lceil4A\rceil+j+2}|\widehat V_*(\eta+\mathrm i u)|\,d u. \end{aligned}

When V∗V_* is supported in a fixed compact interval I⊂(0,∞)I\subset(0,\infty), repeated integration by parts in its Mellin transform gives rapid decay in ∣u∣|u|, uniformly for −A≤η≤1/4-A\le\eta\le1/4. Thus, for a sufficiently large J=J(A,j)J=J(A,j), the supremum of integrals in (A.20) is ≪A,j,I∥V∗∥CJ(I)\ll_{A,j,I}\|V_*\|_{C^J(I)}. This proves (6.12).

Thus (6.9) expresses T(X;Ψ)T(X;\Psi) as OΨ0,S(2∣P∣)O_{\Psi_0,S}(2^{|\mathcal{P}|}) dual sums with ∣C∣≤1/81|C|\le1/81, the asserted ray class dependence of (d,ψ,c0)(d,\psi,c_0), and the weight decay just established. Together with the support and coefficient bounds proved above, this completes the proof of Proposition 6.2 and Lemmas 6.3 and 6.4. □

Separating variables in smooth weights

This appendix justifies the separation of smooth weights in the completed mean-square estimate of Section 6, and the treatment of kernels depending on the row in the Poisson reductions of Sections 4 and 7.

Separating the variables

We use Mellin inversion to separate the variables of a smooth weight, with coefficient bounds controlled by finitely many derivatives.

Lemma B.1. Fix a box I=I1×⋯×Id\mathcal I=I_1\times\cdots\times I_d, where each IjI_j is a compact interval in (0,∞)(0,\infty). Every K∈Cc∞((0,∞)d)\mathcal K\in C_c^\infty((0,\infty)^d) supported in I\mathcal I has a representation

K(x)=∫Rdb(t)∏j=1dxjitj dt,xj>0.\mathcal K(\mathbf x)=\int_{\mathbb R^d}b(\mathbf t) \prod_{j=1}^d x_j^{\mathrm i t_j}\,d\mathbf t, \qquad x_j>0.

For J≥0J\ge0 and every even integer q>J+dq>J+d, the coefficient satisfies

∫Rd∣b(t)∣(1+∣t∣)J dt≪I,J,q,d∥K∥Cq(I).\int_{\mathbb R^d}|b(\mathbf t)|(1+|\mathbf t|)^J\,d\mathbf t \ll_{\mathcal I,J,q,d}\|\mathcal K\|_{C^q(\mathcal I)}.

Proof. Take the Fourier transform in logarithmic coordinates:

b(t)=1(2π)d∫RdK(ey1,…,eyd)e−it⋅y dy.b(\mathbf t)=\frac1{(2\pi)^d}\int_{\mathbb R^d} \mathcal K(e^{y_1},\ldots,e^{y_d}) e^{-\mathrm i\mathbf t\cdot\mathbf y}\,d\mathbf y.

Applying Fourier inversion to this gives (B.1). Since the intervals IjI_j are fixed and bounded away from zero, the CqC^q norm in logarithmic coordinates is bounded by a constant times ∥K∥Cq(I)\|\mathcal K\|_{C^q(\mathcal I)}. Applying (1−Δy)q/2(1-\Delta_{\mathbf{y}})^{q/2} under the integral therefore gives

∣b(t)∣≪I,q,d(1+∣t∣)−q∥K∥Cq(I)(q≥0 even).|b(\mathbf t)|\ll_{\mathcal I,q,d} (1+|\mathbf t|)^{-q}\|\mathcal K\|_{C^q(\mathcal I)} \qquad(q\ge0\text{ even}).

Multiplication by (1+∣t∣)J(1+|\mathbf{t}|)^J and integration proves the bound when q>J+dq>J+d. See also the multivariable Mellin formulation in [39], §10.1, (130)–(131).

In our applications, the weight has the more specific form

KR(x)=∏j=1dWj(xj)F(R∏j=1dxjaj),R>0,\mathcal{K}_R(\mathbf{x})=\prod_{j=1}^{d}W_j(x_j)F\left(R\prod_{j=1}^{d}x_j^{a_j}\right),\qquad R>0,

where the exponents aja_j are fixed real numbers and Wj∈Cc∞((0,∞))W_j\in C_c^\infty((0,\infty)) is supported in IjI_j. Assume that F∈C∞((0,∞))F\in C^\infty((0,\infty)) satisfies

∥F∥A,q:=max⁡0≤m≤qsup⁡u>0(1+u)A∣(u∂u)mF(u)∣<∞(A≥0, q∈Z≥0).\|F\|_{A,q}:=\max_{0\le m\le q}\sup_{u>0}(1+u)^A\left|(u\partial_u)^mF(u)\right|<\infty\qquad(A\ge0,\ q\in\mathbb{Z}_{\ge0}).

On the fixed box I\mathcal{I}, the argument of FF is comparable to RR. The chain rule and the pointwise bound in the proof give coefficients bRb_R satisfying, for A≥0A\ge0 and even q≥0q\ge0,

∣bR(t)∣≪(1+R)−A(1+∣t∣)−q∥F∥A,q∏j=1d∥Wj∥Cq(Ij).|b_R(\mathbf{t})|\ll(1+R)^{-A}(1+|\mathbf{t}|)^{-q}\|F\|_{A,q}\prod_{j=1}^{d}\|W_j\|_{C^q(I_j)}.

Consequently, for J≥0J\ge0 and even q>J+dq>J+d,

∫Rd∣bR(t)∣(1+∣t∣)J dt≪(1+R)−A∥F∥A,q∏j=1d∥Wj∥Cq(Ij).\int_{\mathbb{R}^d}|b_R(\mathbf{t})|(1+|\mathbf{t}|)^J\,d\mathbf{t}\ll(1+R)^{-A}\|F\|_{A,q}\prod_{j=1}^{d}\|W_j\|_{C^q(I_j)}.

The constants depend only on AA, JJ, qq, dd, the intervals IjI_j, and the exponents aja_j. Thus separation preserves the decay in RR, and its cost is controlled by finitely many derivatives of the original weights.

Weights depending on the row

The next lemma extends mean-square bounds for a common test function to smooth kernels depending on the row, with losses controlled by uniform bounds on their derivatives.

Fix compact intervals I⊂int⁡I∗I\subset\operatorname{int}I_* in (0,∞)(0,\infty) and an integer m≥0m\ge0.

Lemma B.2. Consider two families of finite sums

Sj,r(U)=∑naj,r(n)U(xj,r,n),j=1,2,xj,r,n>0,S_{j,r}(U)=\sum_n a_{j,r}(n)U(x_{j,r,n}),\qquad j=1,2,\qquad x_{j,r,n}>0,

with a finite set of rows rr and nonnegative weights wrw_r. Suppose that, for every U∈Cc∞(I∗)U\in C_c^\infty(I_*),

∑rwr∣Sj,r(U)∣2≤Mj∥U∥Cm(I∗)2,j=1,2.\sum_r w_r|S_{j,r}(U)|^2\le M_j\|U\|_{C^m(I_*)}^2,\qquad j=1,2.

Then for any ∣cr∣≤wr|c_r|\le w_r and smooth kernels Kr\mathcal{K}_r supported in I2I^2,

∣∑rcr∑n1,n2a1,r(n1)a2,r(n2)‾Kr(x1,r,n1,x2,r,n2)∣≪I,I∗,mM1M2sup⁡r∥Kr∥C2m+4(I2).\left|\sum_r c_r\sum_{n_1,n_2}a_{1,r}(n_1)\overline{a_{2,r}(n_2)}\mathcal{K}_r(x_{1,r,n_1},x_{2,r,n_2})\right|\ll_{I,I_*,m}\sqrt{M_1M_2}\sup_r\|\mathcal{K}_r\|_{C^{2m+4}(I^2)}.

Proof. Choose a real V∈Cc∞(I∗)V\in C_c^\infty(I_*) equal to one on II, and set Ut(x)=V(x)xitU_t(x)=V(x)x^{\mathrm i t}. Apply (B.1) of Lemma B.1 to each Kr\mathcal{K}_r with d=2d=2 and I=I2\mathcal{I}=I^2. Replacing the second Mellin variable by its negative and multiplying by V(x)V(y)V(x)V(y), which equals one on the support of Kr\mathcal{K}_r, gives

Kr(x,y)=∫R2br(s,t)Us(x)Ut(y)‾ ds dt,\mathcal{K}_r(x,y)=\int_{\mathbb{R}^2}b_r(s,t)U_s(x)\overline{U_t(y)}\,ds\,dt,

with the following uniform bound, obtained from the pointwise coefficient estimate in the proof of Lemma B.1 with d=2d=2 and I=I2\mathcal{I}=I^2:

sup⁡r∣br(s,t)∣≪I,qsup⁡r∥Kr∥Cq(I2)(1+∣s∣+∣t∣)q(q≥0 even).\sup_r|b_r(s,t)|\ll_{I,q}\frac{\sup_r\|\mathcal{K}_r\|_{C^q(I^2)}}{(1+|s|+|t|)^q}\qquad(q\ge0\text{ even}).

The hypothesis applies to each UtU_t, and ∥Ut∥Cm(I∗)≪I,I∗,m(1+∣t∣)m\|U_t\|_{C^m(I_*)} \ll_{I,I_*,m} (1+|t|)^m. Taking the common bound for brb_r before integrating and applying Cauchy–Schwarz in rr therefore bounds the absolute value in (B.4) by a constant depending on II, I∗I_*, mm, qq times

M1M2sup⁡r∥Kr∥Cq(I2)∫R2(1+∣s∣)m(1+∣t∣)m(1+∣s∣+∣t∣)q ds dt.\sqrt{M_1M_2}\sup_r\|\mathcal{K}_r\|_{C^q(I^2)}\int_{\mathbb{R}^2}\frac{(1+|s|)^m(1+|t|)^m}{(1+|s|+|t|)^q}\,ds\,dt.

Taking q=2m+4q=2m+4 makes the integral converge and proves the claim.

Recombining the separated sums

We use weighted Cauchy–Schwarz to pass from mean-square bounds for the separated sums to a bound for their weighted integral.

Lemma B.3. Let rr, ι\iota range over finite sets, let d≥1d\ge1, and let wr≥0w_r\ge0. Suppose the measurable coefficients cι,rc_{\iota,r} satisfy

∣cι,r(t)∣≤bι(t),M:=∑ι∫Rdbι(t) dt<∞,|c_{\iota,r}(\mathbf t)|\le b_\iota(\mathbf t),\qquad M:=\sum_\iota\int_{\mathbb R^d}b_\iota(\mathbf t)\,d\mathbf t<\infty,

where bι:Rd→[0,∞)b_\iota:\mathbb{R}^d\to[0,\infty). For measurable Fι,t(r)F_{\iota,\mathbf t}(r), whenever the right-hand side is finite, one has

∑rwr∣∑ι∫cι,r(t)Fι,t(r) dt∣2≤M∑ι∫bι(t)∑rwr∣Fι,t(r)∣2 dt.\sum_r w_r\Bigl|\sum_\iota\int c_{\iota,r}(\mathbf t) F_{\iota,\mathbf t}(r)\,d\mathbf t\Bigr|^2 \le M\sum_\iota\int b_\iota(\mathbf t) \sum_r w_r|F_{\iota,\mathbf t}(r)|^2\,d\mathbf t.

The same assertion holds for finite sums with the integrals omitted.

Proof. For each fixed rr, weighted Cauchy–Schwarz gives

∣∑ι∫cι,r(t)Fι,t(r) dt∣2≤M∑ι∫bι(t)∣Fι,t(r)∣2 dt.\Bigl|\sum_\iota\int c_{\iota,r}(\mathbf t) F_{\iota,\mathbf t}(r)\,d\mathbf t\Bigr|^2 \le M\sum_\iota\int b_\iota(\mathbf t) |F_{\iota,\mathbf t}(r)|^2\,d\mathbf t.

Multiply by wrw_r and sum over rr to obtain (B.6). The same argument with sums in place of integrals proves the finite version.

References

  1. [1]M. Agrawal, N. Kayal, and N. Saxena, PRIMES is in P, Ann. of Math. (2) 160 (2004), no. 2, 781–793.DOI
  2. [2]S. Bettin and S. M. Gonek, The θ=∞ conjecture implies the Riemann hypothesis, Mathematika 63 (2017), no. 1, 29–33. doi:10.1112/S0025579316000139.DOI
  3. [3]V. Bhargava, G. Ivanyos, R. Mittal, & N. Saxena. (2017). Irreducibility and deterministic r-th root finding over finite fields. In Proceedings of ISSAC 2017 (pp. 37–44). ACM.
  4. [4]G. Bhowmik and I. Z. Ruzsa, Average Goldbach and the quasi-Riemann hypothesis, Anal. Math. 44 (2018), no. 1, 51–56.arxiv.org/abs/1711.06442
  5. [5]V. Blomer, L. Goldmakher, & B. Louvel. (2014). L-functions with n-th-order twists. In Int. Math. Res. Not. (Number 7, pp. 1925–1955). https://doi.org/10.1093/imrn/rns257
  6. [6]H. Davenport. (2000). Multiplicative Number Theory (Vol. 74). Springer.
  7. [7]A. Dunn, & M. Radziwiłł. (2024). Bias in cubic Gauss sums: Patterson's conjecture. In Ann. of Math. (Vol. 200, Number 3, pp. 967–1057). https://doi.org/10.4007/annals.2024.200.3.3
  8. [8]A.-S. Elsenhans, J. Klüners, and F. Nicolae, Imaginary quadratic number fields with class groups of small exponent, Acta Arith. 193 (2020), no. 3, 217–233.arxiv.org/abs/1803.02056
  9. [9]K. Ford. (2002). Zero-free regions for the Riemann zeta function. In Number Theory for the Millennium, II (pp. 25–56). A K Peters. https://arxiv.org/abs/1910.08205
  10. [10]E. Frenkel, R. Langlands, and B. C. Ngô, Formule des traces et fonctorialité: le début d'un programme, Ann. Sci. Math. Québec 34 (2010), no. 2, 199–243. arXiv:1003.4578.arxiv.org/abs/1003.4578
  11. [11]S. D. Galbraith, Mathematics of Public Key Cryptography, Cambridge University Press, 2012.DOI
  12. [13]L. Goldmakher and B. Louvel, A quadratic large sieve inequality over number fields, Math. Proc. Cambridge Philos. Soc. 154 (2013), no. 2, 193–212. doi:10.1017/S0305004112000370.DOI
  13. [14]T. H. Grönwall, Sur les séries de Dirichlet correspondant à des caractères complexes, Rend. Circ. Mat. Palermo 35 (1913), 145–159. doi:10.1007/BF03015596.DOI
  14. [15]F. Grube. (1874). Ueber einige Euler'sche Sätze aus der Theorie der quadratischen Formen. In Z. Math. Phys. (Vol. 19, pp. 492–519).
  15. [16]J. Hadamard, Sur la distribution des zéros de la fonction ζ(s) et ses conséquences arithmétiques, Bull. Soc. Math. France 24 (1896), 199–220. doi:10.24033/bsmf.545.DOI
  16. [17]H. Hasse, Vorlesungen über Zahlentheorie, Grundlehren der mathematischen Wissenschaften 59, Springer, Berlin, 1950. doi:10.1007/978-3-642-52795-1.DOI
  17. [18]D. R. Heath-Brown. (1995). A mean value estimate for real character sums. In Acta Arith. (Vol. 72, Number 3, pp. 235–275).
  18. [19]D. R. Heath-Brown. (2000). Kummer's conjecture for cubic Gauss sums. In Israel J. Math. (Vol. 120, pp. 97–124). https://doi.org/10.1007/s11856-000-1273-y
  19. [20]D. R. Heath-Brown, & S. J. Patterson. (1979). The distribution of Kummer sums at prime arguments. In J. Reine Angew. Math. (Vol. 310, pp. 111–130). https://doi.org/10.1515/crll.1979.310.111
  20. [21]E. Hecke. (1918). Eine neue Art von Zetafunktionen und ihre Beziehungen zur Verteilung der Primzahlen, I, II. In Math. Z.: 1; 6 (pp. 357–376; 11–51).
  21. [22]H. Iwaniec and E. Kowalski, Analytic Number Theory, American Mathematical Society Colloquium Publications 53, American Mathematical Society, Providence, RI, 2004. doi:10.1090/coll/053.DOI
  22. [23]E. Kani. (2011). Idoneal numbers and some generalizations. In Ann. Sci. Math. Québec (Vol. 35, Number 2, pp. 197–227).
  23. [24]N. M. Korobov. (1958). Estimates of trigonometric sums and their applications. In Uspekhi Mat. Nauk (Vol. 13, Numbers 4(82), pp. 185–192). https://www.mathnet.ru/eng/rm7458mathnet.ru/eng/rm7458
  24. [25]E. Kowalski, Y. Lin, Ph. Michel, & W. Sawin. (2020). Periodic twists of GL₃-automorphic forms. In Forum Math. Sigma (Vol. 8, p. e15). https://doi.org/10.1017/fms.2020.7
  25. [26]T. Kubota. (1969). On automorphic functions and the reciprocity law in a number field. Kinokuniya Book-Store. https://hdl.handle.net/2433/84907hdl.handle.net/2433/84907
  26. [27]E. Landau, Neuer Beweis des Primzahlsatzes und Beweis des Primidealsatzes, Math. Ann. 56 (1903), 645–670. doi:10.1007/BF01444310.DOI
  27. [28]J. E. Littlewood. (1912). Quelques conséquences de l’hypothèse que la fonction ζ(s) de Riemann n’a pas de zéros dans le demi-plan R(s)>1/2. In C. R. Acad. Sci. Paris (Vol. 154, pp. 263–266).
  28. [29]J. E. Littlewood, On the class-number of the corpus P(√-k), Proc. London Math. Soc. (2) 27 (1928), no. 1, 358–372.DOI
  29. [31]H. Maier and H. L. Montgomery, The sum of the Möbius function, Bull. London Math. Soc. 41 (2009), no. 2, 213–226. doi:10.1112/blms/bdn119.DOI
  30. [32]G. L. Miller. (1976). Riemann's hypothesis and tests for primality. In J. Comput. System Sci. (Vol. 13, Number 3, pp. 300–317).
  31. [33]H. L. Montgomery and R. C. Vaughan, Multiplicative Number Theory I: Classical Theory, Cambridge Studies in Advanced Mathematics 97, Cambridge University Press, 2007.
  32. [34]M. Ram Murty and A. Sankaranarayanan, Averages of exponential twists of the Liouville function, Forum Math. 14 (2002), no. 2, 273–291.DOI
  33. [35]NIST Digital Library of Mathematical Functions. (n.d.). NIST Digital Library of Mathematical Functions. Retrieved https://dlmf.nist.gov/10.43.E19dlmf.nist.gov/10.43.E19
  34. [36]OpenAI. (2026). The Quasi-Riemann Hypothesis: A Zero-Free Half-Plane Re(s)>7/8. https://github.com/openai/math/blob/main/preprints/The-Quasi-Riemann-Hypothesis-September-30-2026/paper.pdfgithub.com/openai/math/blob/main/preprints/The-Quasi-Riemann-Hypothesis-September-30-2026/paper.pdf
  35. [38]S. J. Patterson, A cubic analogue of the theta series, J. Reine Angew. Math. 296 (1977), 125–161. doi:10.1515/crll.1977.296.125.DOI
  36. [39]I. Petrow and M. P. Young, A generalized cubic moment and the Petersson formula for newforms, Math. Ann. 373 (2019), nos. 1–2, 287–353.DOI
  37. [40]I. Petrow and M. P. Young, The Weyl bound for Dirichlet L-functions of cube-free conductor, Ann. of Math. (2) 192 (2020), no. 2, 437–486. doi:10.4007/annals.2020.192.2.3.DOI
  38. [41]B. Riemann. (1859). Ueber die Anzahl der Primzahlen unter einer gegebenen Grösse. In Monatsber. Königl. Preuss. Akad. Wiss. Berlin (pp. 671–680). https://www.maths.tcd.ie/pub/HistMath/People/Riemann/Zeta/maths.tcd.ie/pub/HistMath/People/Riemann/Zeta
  39. [42]K. A. Rodosskiĭ. (1956). On non-residues and zeros of L-functions. In Izv. Akad. Nauk SSSR Ser. Mat. (Vol. 20, Number 3, pp. 303–306).
  40. [45]E. C. Titchmarsh, A divisor problem, Rend. Circ. Mat. Palermo 54 (1930), 414–429; correction, 57 (1933), 478–479. doi:10.1007/BF03021203; correction.DOI
  41. [46]C.-J. de la Vallée Poussin. (1896). Recherches analytiques sur la théorie des nombres premiers. Première partie: La fonction ζ(s) de Riemann et les nombres premiers en général. In Ann. Soc. Sci. Bruxelles (Vol. 20, pp. 183–256).
  42. [47]C.-J. de la Vallée Poussin, Sur la fonction ζ(s) de Riemann et le nombre des nombres premiers inférieurs à une limite donnée, Mém. Couronnés Autres Mém. Acad. Roy. Sci. Lettres Beaux-Arts Belgique, collection in-8^°, 59 (1899), 1–74. doi:10.3406/marb.1899.2449.DOI
  43. [48]I. M. Vinogradov. (1958). A new estimate of the function ζ(1+it). In Izv. Akad. Nauk SSSR Ser. Mat. (Vol. 22, Number 2, pp. 161–164). https://www.mathnet.ru/eng/im3962mathnet.ru/eng/im3962
  44. [49]I. M. Vinogradov. (1985). Selected Works. Springer.
  45. [50]P. J. Weinberger, Exponents of the class groups of complex quadratic fields, Acta Arith. 22 (1973), 117–124. doi:10.4064/aa-22-2-117-124.DOI
  46. [51]A. Yoshimoto, On the cubic theta function, Nagoya Math. J. 105 (1987), 153–167.DOI

Paper details

Contents