Vorono\"{i} summation via switching cusps

We consider the Fourier expansion of a Hecke (resp.\ Hecke--Maa\ss) cusp form of general level $N$ at the various cusps of $\Gamma_{0}(N)\bs\Hb$. We explain how to compute these coefficients via the local theory of $p$-adic Whittaker functions and establish a classical Vorono\"i summation formula allowing an arbitrary additive twist. Our discussion has applications to bounding sums of Fourier coefficients and understanding the (generalised) Atkin--Lehner relations.


Introduction
The purpose of this note is to prove two formulae. The first, a Voronoï-type summation formula, relating additively twisted Fourier coefficients of a cuspidal Hecke eigenform f to a dual sum of coefficients at another cusp, related to the original twist. Secondly, we give an explicit formula for the Fourier coefficients f at an arbitrary cusp. As applications, we prove upper bounds for sums of Fourier coefficients and we revise the classical Atkin-Lehner relations in a p-adic setting.
Whilst similar formulae are, in certain settings, well known, we provide a perspective here that links the two fundamentally and lays bare the mechanics behind additive twists and switching cusps. That perspective is derived from noting that classical Fourier coefficients are special values of p-adic Whittaker functions on GL 2 (Q p ). In that vein, we intend the following exposition to serve pedagogically as a reference for the analysis of classical Fourier coefficients with adelic machinery.
For the remainder of this introduction we detail our main results 1 (Proposition 4.1 and Corollary 6.3 below) for a holomorphic cusp form f (z) = n≥1 a f (n)e(nz) of level N ∈ N, weight k ∈ 2N and nebentypus χ. Fix a modulus q ∈ N for the additive character n → e(an/q) where (a, q) = 1. Then for a Bruhat-Schwartz function F : R >0 → C we prove that Date: 3 rd April 2019. 1 Our results explicitly handle the case of holomorphic and of Maaß forms simultaneously. 1 where H f F is the Hankel transform of F (see (38) and (39) below) and a f (n; a) denotes the n-th Fourier coefficient of f at the cusp a = a/q of (extended) cusp width δ(a), as defined in (8). Formula (1) constitutes one case of Corollary 6.3 given later.
Many of the technicalities in evaluating such formulae explicitly, in terms of the coefficients a f (n), are buried inside the coefficients a f (n; a), bringing us to our second result, formally given in Proposition 4.1. Write n 1 = p|N p np for some n p ≥ 0 and n 0 ∈ N with (n 0 , N) = 1. Then at the cusp a = a/q for some q | N we have (2) a f (n 1 n 0 ; a) = Ω χ,q,δ(a) a f (n 0 , ∞)e na δ(a)L where W p is the p-th component of the global Whittaker function associated to f (see (24) below), its matrix argument is given in Proposition 4.1; ω χ,p is the p-th component of the Hecke-character associated to χ; and Ω χ,q,δ(a) ∈ C |z|=1 is defined in (33) and in particular satisfies Ω χ,q,δ(a) = 1 whenever χ = 1 or q = 1. By the results of [1], one may evaluate the local factors at p | n 1 explicitly. Furthermore, combining (2) with (1) we thus obtain a general Voronoï summation formula for GL 2 in the N-q aspect.
Remark 1.1. Formula (2) implies the well known phenomenon that the Fourier coefficients a f (n; a) in general are not multiplicative objects. Moreover, due to the possible large values of W p , derived in [1], the coefficients a f (n; a) can be quite large in terms of the level N.
We are able to give some novelty results as immediate corollaries to (2). For example, in Corollary 4.4 we compute the bound n≤X |a f (n; a)| 2 ≪ k 2 N ǫ δ(a) k X k+ǫ + X k− 1 2 +ǫ (L, N L ) k .
Another application of (2) is to derive an adelic approach to the classical Atkin-Lehner relations. We will see in §5 that, via (2), the Atkin-Lehner involution arises naturally through local matrix identities. The corresponding identities relating Fourier coefficients at different cusps arise from several functional equations of the local Whittaker functions W p .
Remark 1.2. Such formulae as (1) and (2) are naturally sensitive to the inducing information that determines the local representation attached to f at p. Hence a Voronoï formula is inherently non-uniform in the N-q aspect outside the extreme cases (N, q) = 1 and N | q, where the local representations do not interact. We remark that our formula trivially agrees with the 'go-to' reference [13, §A] in such extreme cases (see Corollaries 6.5 and 6.6). Uniformity may nonetheless be recovered should one concede to looking at families of fixed (or restricted) representation type at p, in which case such a general Voronoï formula is quickly derived from the formulae given in [1], in conjunction with (1)  (π, C), for an automorphic representation π of GL 2 , in terms of well chosen model periods in the intermediate Hom-spaces; one wall of the Gelfand formation of course gives rise to the Whittaker model. The equality arising from equating the expansions coming from each side determines the functional equation for a general Voronoï-type formula in the representation theoretic setting. Classically, this can be thought of as swapping from the Fourier expansion of a cusp form at ∞ to the Fourier expansion at 0. This is because the group N (resp. wNw −1 ) corresponds to the stabiliser of ∞ (resp. 0). We make this remark precise in Theorem 6.2.
Common notation. Put e(z) := e 2πiz for z ∈ C. If g = ( a b c d ) ∈ GL 2 (C) we write gz = az+b cz+d and j(g, z) = cz + d. If moreover g ∈ SL 2 (Z) then, by convention, a Dirichlet character χ defines a function χ(g) := χ(d). Let (a, b) and [a, b] denote the greatest common divisor and the least common multiple of a, b ∈ N, respectively.
We write a | b ∞ to denote that there exists a k ∈ N such that a | b k ; similarly we let (a, b ∞ ) := max k≥1 (a, b k ). We let ⌊y⌋ and ⌈y⌉ denote the floor and ceiling of y ∈ R. The p-adic valuation of x ∈ Q p is denoted by v p (x). If A is a logical assertion, we use Dirac's symbol δ(A) to denote a 1 if A is true and a 0 otherwise.

2.
Classical Fourier expansions at the cusps of Γ 0 (N)\H 2.1. Classical newforms. The results of this work apply specifically to newforms on Γ 0 (N)\H which vanish in the cusps. These may come from either the world of Maaß forms of holomorphic modular forms. 2 We introduce notation here to deal with both of these cases. From a pedagogical perspective, much of the theory we introduce in this section and the next applies more generally, without the assumptions that the form is new or indeed cuspidal [4,8]. We shall uniformly refer to a function f : H → C as a cuspidal newform of weight k ∈ Z ≥0 if f falls into one of the following two categories: • Cuspidal holomorphic Hecke newforms, in which case k > 0.
We are able to speak of both genres of newform in the same breath by introducing the notation Since ( 1 1 0 1 ) ∈ Γ 0 (N) for any N ≥ 1, we have in particular that f (z + 1) = f (z) by (3). Then with the specialist notation in (4), we deduce the usual Fourier expansion for f (at the cusp ∞) 2.2. The cusps of Γ 0 (N)\H. A detailed account of the construction of the cusps of the Riemann surface Γ 0 (N)\H is given in [14, §3.4.1], to which we refer. We denote this set of cusps by C(N); that is, the set of boundary points of H (modulo the left Möbius action of Γ 0 (N)) that are stabilised by a non-scalar element of Γ 0 (N), or simply C(N) = Γ 0 (N)\P 1 (Q). Indeed, SL 2 (Z) acts transitively on P 1 (Q) and the stabiliser of the point ∞ = [1 : 0] ∈ P 1 (Q) is Γ ∞ = {± ( 1 n 0 1 ) : n ∈ Z}. Thus C(N) is identified with Γ 0 (N)\ SL 2 (Z)/Γ ∞ as a left Γ 0 (N)-set. We write its elements as a = σ −1 ∞ for a choice of representatives {σ} of C(N).
On the other hand, there is a (transitive) right action of SL 2 (Z) on P 1 (Z/NZ) which, as in [14, §3.4.1], may be used to parameterise the set of cusps by classes [q : d] ∈ P 1 (Z/NZ) such that q | N and d ∈ (Z/(q, N/q)Z) × . From this perspective one refers to [q : d] as the "fraction" d/q and calls q the denominator of the cusp. It follows that the number of cusps of denominator q is φ((q, N/q)) and furthermore #C(N) = q|N φ((q, N/q)). In this description, the cusp a = ∞ corresponds to [0 : 1] and is the unique cusp of denominator N.
To summarise these constructions, any cusp a ∈ C(N) may be identified by an element such that σa = ∞. Then the cusp a has denominator q; it corresponds to the class [q : d] ∈ P 1 (Z/NZ) (as above); and if σ ≡ 1 ∈ Γ 0 (N)\ SL 2 (Z)/Γ ∞ then a is equal to the rational cusp a/q ∈ P 1 (Q) where (a, N) = 1 and ad ≡ 1 (mod (q, N/q)).
Equivalently, a more tactile definition may be given by noting that w(a) is the least integer n ≥ 1 such that the stabiliser of a ∈ P 1 (Q), under the left action of SL 2 (Z), contains ( 1 n 0 1 ). In fact, one may compute the width precisely as w(a) = N/(q 2 , N), or equivalently w(a) = [q 2 , N]/q 2 (see [14, §3.4

.1]).
Later on it shall be handy to take any fraction a/q with (a, q) = 1 and consider the equivalence class defining the cusp a. This relates to our notion above as follows: with σ −1 as in (7), let r, s ∈ Z such that qr + Nds = (q, N) and rescale σ so that ar Replacing the coefficient ar by its unique minimal representative a ′ ≥ 0 modulo (q, N), another representative of the cusp a is then a ′ /(q, N), the upshot being that the latter has denominator dividing N.
2.3. The classical Fourier expansion at an arbitrary cusp. For a cusp a ∈ C(N), pick a representative σ ∈ SL 2 (Z) such that σa = ∞. The expansion of f about a is given by the expansion of f | k σ −1 at the cusp ∞. If the nebentypus were trivial, χ = 1, then f | k σ −1 would be periodic on vertical strips of width w(a) in H. However, to account for general characters of conductor M | N, we introduce the extended width As with w(a), the extended width only depends on the equivalence class of a. Suppose that σ = d −b −q a ; that is, a = σ −1 ∞ corresponds to the cusp a/q as described in §2.2. Then we have for each t ∈ Z. A few applications of the modular identity (3) imply that As the conductor of χ is M we see that χ(1−aqw(a)t) = 1 if and only if M | qw(a)t, or equivalently, Hence, our definition of δ(a) is justified in that it is the minimally chosen positive integer satisfying the periodicity relation The resulting Fourier expansion f at the cusp a is of the form a)κ f (ny/δ(a))e(nx/δ(a)).
(One may refer to [9, Theorems 3.4.1 and 3.7.4] for both the holomorphic and Maaß cases, respectively.) When a = ∞ we recover the familiar coefficients a f (n; a) = a f (n); however, unlike a f (n), the coefficients a f (n; a) are not multiplicative in general (the reason for this can be found in Section 4.2 below). We refer to f as a cusp form whenever its constant terms vanish, a f (0; a) = 0, with respect to each a ∈ C(N). By definition, if f is holomorphic then its Laurent expansion contains no negatively indexed terms: a f (n; a) = 0 for all n < 0.
Remark 2.1. Note that δ(a) = w(a) whenever the conductor M is sufficiently small; for example, if M | N 1 where N 1 is the smallest integer such that N | N 2 1 . On the other hand, if q is sufficiently square-full above its prime divisors then δ(a) = w(a), likewise. As with w(a), the extended width depends only on (q, N) as can be seen from δ(a) = w(a) M (qw(a),M ) and the fact that (M | N). Finally, suppose we wish to choose a different representative τ satisfying τ a = ∞. This includes representing a by a different fraction without necessarily asking for the numerator and denominator to be coprime. Then This shows explicitly that replacing the matrix σ by an equivalent one τ only skews the Fourier coefficients by the root of unity χ(γ)e(nm/δ(a)).

Whittaker-Fourier expansions and adele groups
In this section we take an expository route to give explicit realisations of classical forms recast in the adelic theory of representations for G := GL 2 . We hope to exhibit the advantage gained by strong approximation when working with adelic constructions -this is the key insight in converting global ramification problems into p-adic analytic ones.
3.1.1. Notation for adele groups. Let A denote the adele ring of Q and A f its ring of finite adeles; these are given respectively by the restricted direct products A = ′ p≤∞ Q p and A f = ′ p<∞ Q p with respect to the additive subgroups Z p ≤ Q p . We reserve the letter G to denote GL 2 and introduce shorthand for the p-adic 3.1.3. Strong approximation. We diagonally embed Q ֒→ A as a discrete subgroup. The quotient Q\A is compact and we have the strong approximation theorem (11) A [15,). We additionally have strong approximation theorems for the adele groups A × = GL 1 (A) and G(A) = GL 2 (A) which state the following: for each N ≥ 1 (see [8,Ch. 3]). The adelic congruence subgroups are then related to the classical ones via (14) G(Q) ∩ G + ∞ × K 0 (N) = Γ 0 (N). In the spirit of clarity, we specify the place-value of matrices γ ∈ G(Q) by defining the inclusions ι p : 3.1.4. Dirichlet charactersà la Hecke. A Dirichlet character χ (of conductor M) may be associated to a Hecke character ω χ : where pr denotes the projection of R >0 × p<∞ Z × p after applying strong approximation (12) to A × ; and Ch M denotes the cannonical surjection induced by the Chinese remainder theorem, As with all Hecke characters, we have the tensor factorisation ω χ = ⊗ p≤∞ ω χ,p . Observe in particular that ω χ,∞ | R + = 1. See [12, §12.1] for extended details. Abusing notation, we extend the character ω χ = ⊗ p≤∞ ω χ,p to a character of K 0 (N) by setting 3.1.5. Adelisation of classical modular forms. We refer the unfamiliar reader to [8] for a friendly introduction and to the following construction; see also [12, §12]. For our classical newform f , the 'adelisation process' is to associate an automorphic form ϕ : G(Q)\G(A) → C to f by defining (13) with γ ∈ G(Q), g ∞ ∈ G + ∞ and k 0 ∈ K 0 (N). By choice of f , we consider ϕ as fixed throughout our exposition. The central character of ϕ is given by ω χ , which is understood as follows. Via (12) as ω χ (r ∞ ) = 1 by construction. In the adelic setting, the central character also finds itself as the one dimensional representation of K 0 (N) generated by ϕ since for all k 0 ∈ K 0 (N) and g ∈ G(A). In particular, ϕ is right K 1 (N)-invariant.
3.1.6. Automorphic representations. The right regular action 3 , ω χ ), containing only forms with central character ω χ . Let π denote the fixed automorphic representation representation of G(A) generated by ϕ (or, by abuse of language, f ). The property that f is a cusp form is equivalent to π being contained in the cuspidal part of the spectrum of L 2 (Z(A)G(Q)\G(A), ω χ ); that f is moreover a newform directly implies that π is irreducible. Letπ denote the contragredient representation attached to π. Theñ π has central character ω −1 χ and satisfies the isomorphismπ ∼ = (ω −1 χ • det) ⊗ π. The tensor product theorem (see [4]) states that an automorphic representation π may be realised as a restricted tensor product π ∼ = ⊗ ′ p π p where each π p is some (complex) representation of G p . This algebraic construction is tied fundamentally to the theory of Euler products of L-functions in that the coefficients λ π (n) of (17) L(s, π) = p≤∞ L(s, π p ) = L(s, π ∞ ) n≥1 λ π (n) n s are local objects. It is the G(Q)-invariance of π that implies the global functional equation about s → 1 − s involving the epsilon factor ε(s, π) = p≤∞ ε(s, π p ). Note that in our particular example that λ π = λ f and the completed L-function attached to f is given by Λ(s, f ) = L(s, π).
3.2. The Whittaker-Fourier Expansion. Let ψ : Q\A → C × be the standard additive character on A defined as follows: is satisfied for each x ∈ A and g ∈ G(A). By the uniqueness of Whittaker models for G(A), there exists a subspace of such functions which, under the right regular action of G(A), is isomorphic to π; this subspace itself is called the Whittaker model and is denoted by W(π, ψ). The implied G(A)-isomorphism may be explicated via the map where we take the invariant measure dx on A, which naturally descents to a probability measure on Q\A. It is then apparent that W(π, ψ) contains the Fourier transforms of the function x → ϕ(n(x)g) on Q\A whose dual group is identified as Q by considering the characters x → ψ(ξx) for each ξ ∈ Q. By Fourier inversion, we thus derive an expansion of ϕ in terms of W ϕ so that where we have implicitly used the cuspidal property of ϕ to deduce the vanishing of the constant term (at ξ = 0). This adelic Fourier expansion is rather marvellous in that it encodes the classical Fourier expansion (9) of f at all cusps a of Γ 0 (N)\H simultaneously. We now conclude this section with some remarks on using the vantage point of (19) to compute the coefficients a f (n; a).
Shifting to an arbitrary cusp a = σ −1 ∞, we combine (19) with our observation above in (20) to find Note the conceptual part played by strong approximation 13 for this identity to hold. The corollary to (21) is that we may determine the values W ϕ (a(ξ)g z ι f (σ)) in two ways: firstly by returning full circle in recovering the classical Fourier expansion (9), and secondly by explicating their values locally via the associated representation theory. Beginning with the former, a special case of the following result is given (for amusement) in [8,Lem. 3.6]. We give full generality here.
Proof. Using the left invariance of ϕ by G(Q) and of the additive Haar measure where we have defined Firstly, let ξ = n/δ(a) with n ∈ Z. Consider the fundamental domain where we associate R/δ(a)Z = [0, δ(a)). This is chosen using strong approximation (11). For any x p ∈ δ(a)Z p we have σ −1 n(x p )σ ∈ K p (N), considering the inclusion Γ 0 (N) ֒→ K p (N). Then (16) implies that (22). By construction (of δ(a)), we have (15) we find

Computing Fourier coefficients via Whittaker functionals
The key application of Proposition 3.1 is that we can now give an explicit product formula for the classical coefficients a f (n; a) in terms of the factorisation of W ϕ into local Whittaker functionals; this is stated in Proposition 4.1. We then exhibit several immediate corollaries to this product formula.
and thus on all of G(R) by the Iwasawa decomposition. In particular, we find W ϕ (1) = κ f (1). At the finite primes p ∤ N, right invariance by K p determines the normalised spherical vector W p uniquely: by the Iwasawa decomposition it is enough to determine W p on A(Q p ) which is described by the well-known formula (26) W p (a(n)) = |n| 1 2 p λ f (p vp(n) ) for any n ∈ N (see [4,Th. 4.6.5]). In fact (26) this holds more generally for all p < ∞ (see [19,Section 2.4]).
The theme of the present article is to understand the new vectors W p when p | N. By a generalisation of Atkin-Lehner theory [5] the new vector is uniquely determined by the property that W p | K 1,p (vp(N )) = 1. However, the interesting twist in such ramification problems is that a complete description of W p is not obtained by its values (26) on A(Q p ), alone.
By Proposition 3.1, (25) and (26) we recover the well known identity for n > 0. Implicit in (27)  where n 0 ≥ 1 with (n 0 , N) = 1. We extend Proposition 3.1 here, alongside (25) and (26), to allow for n < 0. If f is holomorphic, we intend the archimedean component κ f to account for the vanishing on n < 0. We obtain (28) a f (n; a) = W ∞ (a(n/δ(a))) κ f (n/δ(a)) p<∞ W p (a(n/δ(a))σ) We evaluate the p | n 0 sum via (26) as before so that To handle the terms W p (a(n/δ(a))σ) for p | N -and to suggest how one might go about evaluating them -we invoke the Bruhat decomposition of σ which, with notation as in §2.3, is given by We now consider the expression prime-by-prime; to this end write q = p|N p qp , M = p|N p Mp , N = p|N p Np and recall that 4 Then, for some p | N, the right-K p,1 (N) invariance of W p implies (31) W p a n q 2 δ(a) w n −a q = ω χ,p n(q 2 δ(a), p ∞ ) (n, p ∞ )q 2 δ(a) W p a(p np−dπ p (qp) )wn u p p −qp where u p = −a (n, p ∞ )δ(a)q n(δ(a)q, p ∞ ) ∈ Z × p . A convenient notation for us now shall be for t, l ∈ Z and v ∈ Z × p . This notation is in accordance with [1,2,7,17,18] and we use it here to summarise our computations with the following proposition.  Proof. After the deduction of (31), the proof now follows from the observation that, Remark 4.2. By Proposition 4.1, we reduce the problem of understanding the coefficients a f (n; a) to the local problem of evaluating the terms W p (g t,l,v ) for p | N.
There is a quite beautiful method to do so via taking the Fourier expansion of the function v → W p (g t,l,v ) on Z × p /(1 + p Np Z p ) and computing the Fourier coefficients via the Jacquet-Langlands' local functional equation for GL 2 (see [18, §2]). These computations were performed to some extent in [1] and then in greater detail in the forthcoming (University of Bristol) PhD thesis of the first named author. Whilst ξf may not be a newform itself (even if f is), there exists a unique newform f ξ which generates the automorphic representation ω ξ π. From the identification of the associated local Whittaker newforms W p as in (24), the Fourier expansion of f ξ at the cusp infinity is easily determined by means of (27), say. However, at other cusps a = ∞ this is not the case. Proposition 4.1 shows that a f ξ (n; a) can be expressed in terms of the Whittaker new vectors in W(ω ξ,p π p , ψ p ), however at least when ω ξ ramifies at p, there is no uniform formula relating W p and ξf .

4.3.
Bounds for Fourier coefficients on average. We now explore some amusing corollaries to Proposition 4.1.
Corollary 4.4. Let 0 ≤ θ ≤ 7 64 + ǫ be an exponent towards the Ramanujan conjecture for Hecke eigenvalues (see [3]). Then we have Proof. According to Proposition 4.1 we have to bound 0<n≤X |a f (n; a)| 2 = Observe that the function is (q, N/q) periodic, see also [17,Lemma 3.12]. Thus, by the Chinese reminder theorem and [17, Proposition 2.10] we have The result follows after estimating the remaining sums, in particular the n 1 -sum using the Rankin trick.

4.4.
A worked example: principle series. Let f be a holomorphic modular form of nebentypus χ. Furthermore, assume that M = N = p h for even h. Thus, π p is the irreducible principle series ω χ,p | · | s ⊞ | · | −s in the sense of [11] (see also [7, §2.1.6] for notation). In this specific case one can extract a precise expression for W p from [1, Lemma 5.8] for the sake of completeness we will give a proof here.
Lemma 4.6. We have Proof. For the proof of this lemma we will temporarily introduce some notation from [1,18]. Indeed, let W πp be the unique K 1,p (N) ′ -fixed vector normalised by W πp (1) = 1. Using [18, Lemma 2.18, Propositon 2.28] we find that . Furthermore according to [18, (11)] we expand Fortunately the finite Fourier coefficients have been computed in [1,Lemma 2.3]. With this at hand we obtain We now apply [18, Lemma 2.37] to find b χ ∈ Z × p such that The claim lemma now follows directly from orthogonality of characters.
Note that individually |a f (n, a/p h 2 )| 2 ≪ ǫ n k−1+ǫ N − k 2 + 1 2 . Thus for large X the average bound above exploits the narrow support of the coefficients a f (n, a/p h 2 ).

Generalised Atkin-Lehner relations
In this section we will show how to relate the Fourier coefficients at different cusps. Classically, this corresponds to the Atkin-Lehner involution. From our point of view it will appear in terms of certain 'functional equations' of the local Whittaker functions W p . The upshot is that we get a large variety of relations between the cusps. Even though these relations become combinatorially quite involved and require a bit of notation we will maintain a description of the complete picture.
We start by providing the key identities. Indeed, this section is based on the following three lemmata. The first one is a local incarnation of the fact that cusps of the form a/q with N | q are equivalent to the cusp ∞.
p and all t ∈ Z. This is a reformulation of [2, Lemma 2.3]. We include a shorter more instructive proof here.
Proof. Note that Since l ≥ N p , the final matrix above is in K 1,p (N). Thus, exploiting the transformation behaviour of W p and (26) yields the result at once.
The next result is a local take on switching the cusps 0 and ∞ using the classical Atkin-Lehner involution. Proof. We put W ′ = π p (w)W p and define The local functional equation for this zeta integral is given by .
We can rewrite this as The statement follows after evaluating the remaining archimedean integral.
The final lemma lies at the heart of the generalised Atkin-Lehner involution.
HereW p is the normalised new vector inπ p .
This is exactly [18, Proposition 2.28] after accounting for an unramified twist to remove the assumption ω χ,p (p) = 1. To avoid shuffling around the unramified twist we give a self-contained proof.
In particular, after checking the right transformation behaviour, we find that π p (w n )W p = ω χ,p (p Np )ε(1/2, π p )W ′ p . For 0 ≤ l ≤ N p we have the matrix identity We compute . To conclude the proof we must identify W ′ p in terms ofW p . But given the isomorphism and since Ψ(W ′ p )| K 1,p = 1, the uniqueness of the new vector inπ p implies that We now piece together these local results. Given a cusp a/q, we can flip the sign of v p (q) according to the lemmata above at any place p | N. Let us moreover assume the cusp to be in standard form; that is, q | N. χ is the adelisation of χ (p) .) With this at hand we define The idea above is to replace π p (up to unramified twist) withπ p at each place p | M S . Note that the conductor of ω S π is M S , while the conductor of π S remains N. Finally, let f S be the newform associated to π S . This turns out to be the correct (generalised) Atkin-Lehner partner of f for given S. The generalised (pseudo)-Atkin-Lehner eigenvalue will turn out to be related to η a (f, S) := p∈S ε(1/2, π p )ω χ,p (−a).
We are now ready to establish the following proposition.
Proposition 5.4. Let S be as above. Then we have The cusp a S = a S /q S , where a S is uniquely determined modulo δ(a)q = δ(a S )q S defined by Thus, using Proposition 4.1 for a f S (n; a S ), yields the desired result.

The Voronoï summation formulae
We now discuss some variants of the Voronoï summation formula for GL 2 . Our starting point is the following pre-Voronoï formula. Lemma 6.1. Let ζ ∈ A and let φ be a cuspidal automorphic form in the representation space of π. Then we have The trick used alongside the pre-Voronï summation formula is to choose the cuspidal function in various ways φ to determine a variety of Voronoï summation formulae. The following one is our take on a commonly used classical formula. This previously appeared in [2, Lemma 2.1] and is based on [20, Theorem 3.1].
Theorem 6.2. Let f be a modular form and let F : R >0 → R be a compactly supported smooth function and a = p/q be a cusp given by scaling matrix Then n∈N e n a b a f (n; a) n δ(a) where the cusp b := (apδ(a) + bs)/(qaδ(a) + rb) is not written in lowest terms, and the Hankel transform H f F of F is defined by if f is holomorphic and by (39) if f is Maaß form.
Proof. We will start by setting up the left hand side of (37) accordingly. We choose Furthermore we pick φ as follows. Write φ = (v φ,f , v φ,∞ ) under the isomorphism π = π f ⊗ π ∞ . Then φ will be uniquely determined by assuming that v φ,f = π f (a(δ(a) −1 )σ a )v ϕ,f and that in the Kirillov model of With this at hand we find that Furthermore φ is chosen such that W ∞,φ (a(ξ)) = |ξδ(a) −1 | We now evaluate the right-hand side of (37). We do so by treating archimedean and finite places separately. This is possible aŝ This holds since the action of w in the Kirillov model is precisely described by this Hankel transform. See [6] for details. In order to deal with the finite places we artificially write The last equality follows from j(σ, a b + iy) = −iby and σ( a b + iy) = − a b + i yb 2 .
Remark 6.4. This example shows that Voronoï summation is precisely the transition between the Fourier expansions of f at different cusps. While classically it is a very subtle matter to include an archimedean test function, it is straightforward adelically.
A reasonable generalisation of Corollary 6.3 would be to replace the cusp a = ∞ on the left-hand side with an arbitrary cusp. One obtains the following extension of our theorem given in Corollary 6.5. We end this note by reproducing two of the most classical Voronoï formula from Corollary 6.3. These agree with the formula given in [13], for instance. The statement then follows by using Lemma 5.2 and using the additional characters to twist a f (n 0 ) to af (n 0 ).
Of course this corollary exploits that if (a, N) = 1 then the cusp b = a/b is equivalent to 0, which switches to ∞ under the classical Atkin-Lehner involution. Lastly we consider the other extreme case N | b. In this case b = a/b is in fact equivalent to ∞. Once again, this result would follow from the application of Corollary 6.3 together with Proposition 4.1 and Lemma 5.1, but this would be like shooting pigeons with canons. We give a much simpler proof here.
Proof. In this situation δ(a) = 1 so that by Corollary 6.3 we have that Furthermore, because N | b the scaling matrix for a is in Γ 0 (N). Thus according to (10) we find that a f (n, a) = χ(a) −1 a f (n). This concludes the proof.