The source survey poses two questions. (1) The approximate functional equation computes a central L-value L(π,1/2)L(\pi,1/2) in time C(π)1/2+ε\mathrm{C}(\pi)^{1/2+\varepsilon}, where C\mathrm{C} is the analytic conductor; can the ideas entering the proofs of subconvexity be reinterpreted as faster, “subconvex-time” algorithms with exponent below 1/21/2? (2) Let N(X)N(X) count cuspidal representations of GLn(Q)\mathrm{GL}_n(\mathbb{Q}) with analytic conductor at most XX; determine the asymptotic behavior of N(X)N(X), in particular whether log⁡N(X)/log⁡X→n+1\log N(X)/\log X\to n+1. Question (2) is answered affirmatively by Brumley–Milićević for n≤2n\le 2 without restriction and for n≥3n\ge 3 within the archimedean spherical spectrum — pure power growth ∣F(Q)∣∼C(F)Qn+1|\mathfrak{F}(Q)|\sim\mathcal{C}(\mathfrak{F})Q^{n+1} with no logarithmic factors — leaving open only the case of non-spherical archimedean components for n≥3n\ge 3; no follow-up removing this restriction has appeared. Question (1) remains open in the general conductor aspect: aspect-specific subconvex-time algorithms due to Hiary (Riemann zeta; Dirichlet LL-functions at smooth moduli), Vishe (GL2\mathrm{GL}_2 spectral aspect; twists by smooth moduli), and amortized regimes (Odlyzko–Schönhage, Booker–Then) settle only special aspects, and no uniform C1/2−δC^{1/2-\delta} algorithm for GL1/GL2\mathrm{GL}_1/\mathrm{GL}_2 central values is known. The open core is the algorithmic question, together with the residual counting case.

Background

An automorphic representation of GLn\mathrm{GL}_n over Q\mathbb{Q} is an irreducible constituent of the space of complex-valued functions on the adèle group GLn(AQ)\mathrm{GL}_n(\mathbb{A}_\mathbb{Q}) (the restricted product of GLn(R)\mathrm{GL}_n(\mathbb{R}) with the compact pp-adic groups GLn(Zp)\mathrm{GL}_n(\mathbb{Z}_p)) transforming under the left action of GLn(Q)\mathrm{GL}_n(\mathbb{Q}); it is cuspidal when the functions are orthogonal to the images of all proper parabolic subgroups. For n=1n=1 the cuspidal representations are Dirichlet characters (multiplicative characters of the invertible residue classes modulo an integer, necessarily primitive when viewed as automorphic forms), and for n=2n=2 they include classical holomorphic cusp forms and Maass cusp forms (non-holomorphic eigenfunctions of the hyperbolic Laplacian) on the upper half-plane. To each unitary cuspidal representation π\pi is attached an LL-function

L(π,s)=∏vLv(π,s),L(\pi,s)=\prod_v L_v(\pi,s),

an Euler product over the places vv of Q\mathbb{Q} that converges in a right half-plane, extends to an entire function, and satisfies a functional equation relating ss and 1−s1-s. For unitary π\pi the functional equation is symmetric about the central point s=1/2s=1/2, and L(π,1/2)L(\pi,1/2) is the central value.

The complexity of π\pi is measured by the analytic conductor C(π)\mathrm{C}(\pi), a positive real number packaging the finite ramification (the product of the local conductors q(πv)q(\pi_v) over the finite primes) and the archimedean spectral parameters (archimedean = at the real place; roughly, a polynomial in the eigenvalues of the Lie algebra of GLn(R)\mathrm{GL}_n(\mathbb{R})): for ζ(1/2+it)\zeta(1/2+it) it is comparable to 1+∣t∣1+|t|, and for a Dirichlet character of modulus qq with trivial archimedean data it is qq. The functional equation combined with the Phragmén–Lindelöf principle gives the convexity bound ∣L(π,1/2)∣≪C(π)1/4+ε|L(\pi,1/2)|\ll \mathrm{C}(\pi)^{1/4+\varepsilon}, and the subconvexity problem asks for a bound C(π)1/4−δ\mathrm{C}(\pi)^{1/4-\delta} for some fixed δ>0\delta>0, uniformly in the aspect in which π\pi varies (spectral parameter, conductor, or character twists). The source monograph-length survey of Michel–Venkatesh established subconvexity for all LL-functions of GL1\mathrm{GL}_1 and GL2\mathrm{GL}_2 over number fields, uniformly in all aspects, and organized the open questions of the area.

Computationally, the standard tool is the approximate functional equation: a representation of L(π,1/2)L(\pi,1/2) as a smoothly truncated sum over roughly C(π)1/2\mathrm{C}(\pi)^{1/2} Dirichlet coefficients (Hecke eigenvalues) with explicit test functions. Since each term is computable in polylogarithmic time given the coefficients, a single central value can be computed to any fixed precision in about C(π)1/2+ε\mathrm{C}(\pi)^{1/2+\varepsilon} arithmetic operations. The exponent 1/21/2 is thus the algorithmic twin of the analytic convexity exponent 1/41/4, and one asks whether the extra cancellation exploited in subconvexity proofs — which replace the soft convexity argument by genuine structure of the family — can be reinterpreted as an algorithm running in time C1/2−δ\mathrm{C}^{1/2-\delta}. There was a precedent: for ζ(1/2+it)\zeta(1/2+it), Riemann–Siegel-type methods already achieve essentially t1/2t^{1/2}, and refined summation technologies — the amortized method of Odlyzko–Schönhage for many evaluation points, and Hiary's algorithms with exponents 4/134/13 and, in variants, 1/31/3 — show that the barrier can be beaten in special settings. Throughout, run time is measured either in unit-cost arithmetic operations on numbers of O(log⁡C)O(\log\mathrm{C}) bits, or in bit complexity, with a demand of pp bits of precision contributing polynomially in pp.

The second question concerns the size of the family of all such objects. Let F\mathfrak{F} denote the universal family of unitary cuspidal representations of GLn(AQ)\mathrm{GL}_n(\mathbb{A}_\mathbb{Q}), taken up to twisting by the unitary characters ∣det⁡∣it|\det|^{it} (without this quotient — or without fixing the central character, as in the source's finiteness corollary — every unitary twist orbit would contribute a continuum of members), and let N(X)N(X) be the number of its members with analytic conductor at most XX. Michel–Venkatesh proved that N(X)N(X) is bounded by a power of XX and asked for the correct growth rate, suggesting log⁡N(X)/log⁡X→n+1\log N(X)/\log X\to n+1; the answer controls, for example, the number of twists required in converse-theorem characterizations of automorphic representations. A sanity check at n=1n=1: the number of Dirichlet characters of conductor at most XX is asymptotic to (18/π4)X2(18/\pi^4)X^2, consistent with the exponent n+1=2n+1=2. Brumley–Milićević resolved the question in the main cases: denoting by F(Q)\mathfrak{F}(Q) the members of conductor at most QQ, they proved the Weyl–Schanuel law

∣F(Q)∣∼C(F) Qn+1,|\mathfrak{F}(Q)|\sim\mathcal{C}(\mathfrak{F})\,Q^{n+1},

pure power growth with no logarithmic factors, for n≤2n\le 2 unconditionally and for n≥3n\ge 3 restricted to the archimedean spherical spectrum, i.e. to representations π\pi whose archimedean component π∞\pi_\infty contains a nonzero vector fixed by the maximal compact subgroup On(R)\mathrm{O}_n(\mathbb{R}) (the analogue of weight zero at infinity). The restriction is technical, stemming from the explicit spectral inversion of archimedean test functions in the Arthur trace formula (a noncommutative generalization of the Selberg trace formula, relating the spectrum of GLn\mathrm{GL}_n to orbital integrals), and the authors state their belief that it can be removed in a subsequent work. Their method rests on an Effective Limit Multiplicity property for the trace formula and yields as by-products uniform Weyl laws with logarithmic savings in the level and quantitative bounds on the non-tempered spectrum (representations whose local components fail the expected Ramanujan-type size bounds).

Problem Statement

The source poses two questions; they are logically independent but originate from the same survey, and the problem is completed only when each surviving core below is resolved.

(A) Fast computation of central L-values — open; the primary surviving core. By the approximate functional equation, L(π,1/2)L(\pi,1/2) for a unitary cuspidal representation π\pi of GLn(AQ)\mathrm{GL}_n(\mathbb{A}_\mathbb{Q}), n≤2n\le 2, is computable to any prescribed precision in time C(π)1/2+ε\mathrm{C}(\pi)^{1/2+\varepsilon}, where C\mathrm{C} is the analytic conductor. Does there exist an algorithm computing L(π,1/2)L(\pi,1/2) in time C(π)1/2−δ+o(1)\mathrm{C}(\pi)^{1/2-\delta+o(1)} for some fixed δ>0\delta>0 — in the source's phrasing, an algorithmic interpretation of the ideas entering the proofs of subconvexity, though an answer by any method resolves the mathematical target? The following aspects are settled and form the boundary of the open core: ζ(1/2+it)\zeta(1/2+it) (Hiary: exponent 4/134/13, with variants at 2/52/5 and 1/31/3, improving the classical 1/21/2 of Riemann–Siegel, 3/83/8 of Schönhage and 1/31/3 of Heath-Brown) and Dirichlet LL-functions at power-full (e.g. smooth) moduli (Hiary: q1/3+εq^{1/3+\varepsilon}); central values of GL2\mathrm{GL}_2 cusp forms in the spectral aspect (Vishe: T7/8+o(1)T^{7/8+o(1)}); twists L(1/2,f×χq)L(1/2,f\times\chi_q) by Dirichlet characters of smooth or highly composite modulus qq (Vishe: q5/6+o(1)q^{5/6+o(1)}); and amortized regimes for many evaluation points (Odlyzko–Schönhage; Booker–Then). The open core is the general conductor aspect — for example Dirichlet LL-functions at prime or arbitrary moduli, GL2\mathrm{GL}_2 central values varying in the level aspect, or ideally a uniform statement over all aspects — and, beyond that, the extension to n≥3n\ge 3, where even the analytic subconvexity bounds remain unknown. A proof that no such algorithm exists in a stated computational model would equally resolve this part.

(B) The asymptotic count of the universal family — largely resolved; residual core. Let N(X)N(X) be the number of cuspidal representations of GLn(Q)\mathrm{GL}_n(\mathbb{Q}) with analytic conductor at most XX, counted up to unitary twist — that is, the universal family F\mathfrak{F} of the Background. Determine the asymptotic behavior of N(X)N(X), in particular whether log⁡N(X)/log⁡X→n+1\log N(X)/\log X\to n+1. Brumley–Milićević prove ∣F(Q)∣∼C(F)Qn+1|\mathfrak{F}(Q)|\sim\mathcal{C}(\mathfrak{F})Q^{n+1} for n≤2n\le 2 without restriction and for n≥3n\ge 3 under the archimedean spherical hypothesis, which answers the source question affirmatively in those cases. The surviving core is the case n≥3n\ge 3 with arbitrary — in particular non-spherical — archimedean component: establishing the Effective Limit Multiplicity property, or the Weyl–Schanuel law directly, for non-spherical inducing data, thereby removing the restriction that Brumley–Milićević identify as purely technical but do not remove; or a disproof of the n+1n+1 prediction for the unrestricted family.

Completion boundary: the problem is answered when (A) is resolved by a proven subconvex-time algorithm for a conductor aspect outside the settled list above (the uniform all-aspects statement over GL1/GL2\mathrm{GL}_1/\mathrm{GL}_2 being the natural full target), or by a proven impossibility result in a stated computational model; and (B) is resolved by the asymptotic for the unrestricted GLn\mathrm{GL}_n universal family, n≥3n\ge 3, or by a refutation. Partial progress on either core alone does not complete the problem, but is meaningful.

References

  1. Philippe Michel and Akshay Venkatesh, “The subconvexity problem for GL2”, Publications Mathématiques de l'IHÉS 111 (2010), 171–271. DOI 10.1007/s10240-010-0025-8, https://doi.org/10.1007/s10240-010-0025-8
  2. Farrell Brumley and Djordje Milićević, “Counting cusp forms by analytic conductor”, Ann. Sci. Éc. Norm. Supér. (4) 57, no. 5 (2024), 1473–1597. DOI 10.24033/asens.2595, https://doi.org/10.24033/asens.2595
  3. Farrell Brumley, Didier Lesesvre, and Djordje Milićević, “Conductor zeta function for the GL(2) universal family”, arXiv:2105.02068 (2021), https://arxiv.org/abs/2105.02068
  4. Ghaith A. Hiary, “Fast methods to compute the Riemann zeta function”, Annals of Mathematics (2) 174, no. 2 (2011), 891–946. DOI 10.4007/annals.2011.174.2.4, https://doi.org/10.4007/annals.2011.174.2.4
  5. Ghaith A. Hiary, “Computing Dirichlet character sums to a power-full modulus”, Journal of Number Theory 140 (2014), 122–146. DOI 10.1016/j.jnt.2013.12.005, https://doi.org/10.1016/j.jnt.2013.12.005
  6. Pankaj Vishe, “Rapid computation of L-functions for modular forms”, International Mathematics Research Notices 2013, no. 11 (2013), 2624–2656. DOI 10.1093/imrn/rns112, https://doi.org/10.1093/imrn/rns112
  7. Pankaj Vishe, “A fast algorithm to compute L(1/2, f×χ_q)”, Journal of Number Theory 133, no. 5 (2013), 1502–1524. DOI 10.1016/j.jnt.2012.10.005, https://doi.org/10.1016/j.jnt.2012.10.005
  8. Andrew R. Booker and Holger Then, “Rapid computation of L-functions attached to Maass forms”, International Journal of Number Theory 14, no. 5 (2018), 1459–1485. DOI 10.1142/S1793042118500896, https://doi.org/10.1142/S1793042118500896
  9. Fredrik Johansson, “Rapid computation of special values of Dirichlet L-functions”, arXiv:2110.10583 (2021), https://arxiv.org/abs/2110.10583
  10. Akshay Venkatesh, “Sparse equidistribution problems, period bounds and subconvexity”, Annals of Mathematics (2) 172, no. 2 (2010), 989–1094. DOI 10.4007/annals.2010.172.989, https://doi.org/10.4007/annals.2010.172.989
  11. Andrew M. Odlyzko and Arnold Schönhage, “Fast algorithms for multiple evaluations of the Riemann zeta function”, Transactions of the American Mathematical Society 309, no. 2 (1988), 797–809. DOI 10.1090/S0002-9947-1988-0961614-2, https://doi.org/10.1090/S0002-9947-1988-0961614-2
  12. Ian Petrow, “The Weyl law for algebraic tori”, Journal of the European Mathematical Society 26 (2024), 2441–2532. DOI 10.4171/JEMS/1465, https://doi.org/10.4171/JEMS/1465
  13. Ernest Hunter Brooks and Ian Petrow, “Counting automorphic forms on norm one tori”, Acta Arithmetica 183 (2018), 117–143. DOI 10.4064/aa161214-22-1, https://doi.org/10.4064/aa161214-22-1
  14. Manami Roy, Ralf Schmidt, and Shaoyun Yi, “On counting cuspidal automorphic representations for GSp(4)”, Forum Mathematicum 33, no. 3 (2021), 821–843. DOI 10.1515/forum-2020-0313, https://doi.org/10.1515/forum-2020-0313
  15. Philippe Michel, “Exposé Bourbaki 1190: Recent progress on the subconvexity problem”, Astérisque, Société Mathématique de France (2023). DOI 10.24033/ast.1191, https://doi.org/10.24033/ast.1191
  16. Valentin Blomer, Farrell Brumley, and Ilya Khayutin, “The mixing conjecture under GRH”, arXiv:2212.06280 (2022), https://arxiv.org/abs/2212.06280
  17. Valentin Blomer, Farrell Brumley, and Maksym Radiwiłł, “Joint Linnik problems”, arXiv:2603.05609 (2026), https://arxiv.org/abs/2603.05609