REVIEW 3 major objections 4 minor
Breaking Universality in the Lower Order Terms in the 1-level and 2-level Density of Holomorphic Cusp Newforms
T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read This paper shows that the lower-order terms in the 1- and 2-level densities of holomorphic cusp newforms are not universal: they depend on how the prime factors of the level N grow, with explicit new corrections when one prime factor stays
desk verdict A substantial extension of Miller's lower-order-term program, but the 2-level density claim is incomplete because Eq. (1.19)'s root-number average is never evaluated or bounded. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the averaged explicit formula, in which n-level densities are rewritten as sums over primes and Hecke eigenvalues, together with weighted moment sums A'_r, A_r, B''_{r1,r2}, B'_{r1,r2}, B_{r1,r2}, evaluated by the trace formula and its extension to arbitrary level. Rational-function identities convert shifted moment sums M_{c,k}(p) into closed forms in p and lambda_f(p); the surviving universal terms are controlled by Euler-type constants and by Catalan-number coefficients. Whether the level-divisor moments A'_r and B'' survive determines the regime dichotomy.
What would settle it
For N=2q with q a large prime, fixed weight k, and an even Schwartz function phi whose Fourier transform is supported in [-0.1,0.1], compute the weighted 1-level density via the explicit formula and the trace formula to precision 1/log^5 R. After subtracting the universal term from Theorem 7.2, the remaining coefficient of hat phi(0)/log R should equal -2 log 2 / (2^2 - 1) = -(2/3) log 2; any other value would refute the fixed-prime correction.
Extended reading notes
Core claim
The central discovery is that the factorization of the level controls lower-order terms in the weighted density of low-lying zeros. In four regimes—N prime, N=q1q2 with q1 fixed and q2 growing, N=q1q2 with both primes growing as powers of N, and N=p^2—the main terms coincide, but the lower-order terms agree with the prime-level case only when no prime factor is small. If q1 is fixed, the 1-level density gains explicit corrections such as SA'(F) = -2 (log q1)/(log R) * hat phi(0)/(q1^2-1) - (log q1)/(log R) * hat phi''(0)/(q1^2-1) + O(1/log^5 R), and analogous corrections appear in the 2-level density. The paper identifies the oldform space as the mechanism: a small prime factor leaves oldfor
Load-bearing premise
The argument assumes that the error terms in the trace-formula estimates for the weighted harmonic sums are small enough to survive the repeated sums over primes at precision 1/log^4 R; if a hidden dependence on k or N invalidates this, the regime classification collapses.
Editorial extensions
If this is right
- If correct, the lower-order terms can serve as arithmetic fingerprints: two families with the same symmetry type can be distinguished by their 1- and 2-level densities at scale 1/log^4 R.
- The explicit fixed-prime formula gives testable predictions, e.g. SA'(F) = -2 log(q1)/log(R) * hat phi(0)/(q1^2-1) - log(q1)/log(R) * hat phi''(0)/(q1^2-1) + O(1/log^5 R).
- The N=p^2 case joins the prime-level and two-growing-factors cases, so square levels do not break universality at this precision.
- The conjectured generalization implies that only prime factors growing slower than the reciprocal of the target error create genuinely new lower-order terms.
- Because the 2-level density error is dictated by the 1-level density through inclusion-exclusion, sharper 1-level control automatically improves n-level computations.
Reading between the lines
- The dichotomy likely persists for higher n-level densities: the inclusion-exclusion argument makes the n-level error inherit the 1-level error, so the same factor-growth regimes should produce family-dependent lower-order terms at every level.
- The fixed-q1 corrections are proportional to log(q1) and to 1/(q1^2-1), so comparing levels with different small fixed primes (say q1=2 and q1=3) would provide a clean numerical test of the formula's shape.
- A sharper threshold between 'small prime at most a given size' and 'largest factor at least a fixed power smaller than N' may depend on the error exponent; intermediate regimes where the small prime grows like log N, log^2 N, or log^3 N could map the boundary.
- The same framework could be adapted to other GL(2) families, or to varying weight with fixed level, where analogous lower-order terms may show a complementary dichotomy.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the weighted 1- and 2-level densities for families of holomorphic cusp newforms of level N, with the harmonic weight w_R(f)=Z_N(1,f)/Z(1,f). Four level regimes are considered: N prime, N=q_1q_2 with q_1 fixed and q_2→∞, N=q_1q_2 with both primes growing at complementary power rates, and N=p^2. The main technical content is a computation of the terms S_1(F,φ) and S_2(F,φ_1,φ_2) up to O(log^{-4}R), followed by explicit lower-order terms in Theorems 7.1–7.6. The advertised conclusion is that the lower-order terms agree with the prime-level case whenever both prime factors go to infinity, but differ when one prime factor is fixed, with an explicit formula such as SA'(F) = -2 log(q_1)/log(R) · φ̂(0)/(q_1^2-1) - log(q_1)/log(R) · φ̂''(0)/(q_1^2-1) + O(log^{-5}R).
Significance. If the derivation were complete, the paper would give a genuine sharpening of Miller's 2009 lower-order-term analysis and an explicit family-dependent breaking of universality in the 1- and 2-level densities. The manuscript has real strengths: the skeleton via the Petersson trace formula is standard but carefully executed, the support restriction σ<0.22 is explicit and used in the tail estimates, and the appendices contain detailed proofs of several auxiliary lemmas. However, the central 2-level density claim is not established as written because a term in Eq. (1.19) involving the harmonic average of the root number ε_f is never evaluated or bounded. In addition, there is a concrete error in the principal fixed-q_1 formula, Eq. (7.1)/(7.3). These are load-bearing issues, not presentation problems.
major comments (3)
- [Eq. (1.19); Theorems 7.3–7.6] The claimed O(log^{-4}R) computation of the 2-level density is incomplete. Equation (1.19) expresses D_2(F,φ_1,φ_2) as the computed S_1 and S_2 terms plus φ_1(0)φ_2(0) · (Σ_{f∈F} w_R(f)(1-ε_f))/(2W_R(F)). This last factor is the harmonic average of the root number ε_f. Sections 4–7 compute SA′, SA, SB″, SB′, SBf, and SB∞, but never evaluate or bound Σ_f w_R(f)ε_f / W_R(F). Since ε_f=±1, this term can be a nonzero constant independent of R unless a separate argument is given. Even granting Propositions 2.3–2.5 exactly, D_2 is not determined to the claimed precision. The authors must prove that the root-number average is 1+O(log^{-4}R), or else compute its family-dependent lower-order terms and include them in Theorems 7.3–7.6.
- [Theorem 7.1, Eq. (7.1)/(7.3)] The second term in Eq. (7.1) does not follow from the expansion in Eq. (7.3). Set a=log(q_1)/log(R). The contribution of φ̂''(0) to SA′(F) is -2a Σ_{r≥1} q_1^{-2r} · (1/2)φ̂''(0)(2ra)^2 = -4a^3 φ̂''(0) Σ_{r≥1} r^2 q_1^{-2r}, which is O((log q_1)^3/log^3 R). The paper instead writes -a φ̂''(0)/(q_1^2-1), which is O(log q_1/log R). The two differ by a factor of order (log q_1/log R)^2, and the claimed O(log^{-5}R) error cannot absorb the discrepancy. This is a central formula for the fixed-q_1 case; it must be corrected or the notation clarified.
- [Theorems 7.3–7.6, proof sketches] The four main theorems on S_2 are only given proof sketches. The generic tail bounds in Lemmas 7.7–7.8 are useful, but the substitutions from Lemmas 6.4–6.7 into SB′, SBf, and SB∞ involve many mixed sums over p_1,p_2, including fixed-q_1 exceptional terms with implied constants O_{q_1,k} in Lemma 6.5 (e.g., B′_{r_1,r_2}(q_1,p_2) for odd r_1). The manuscript does not verify that every such error is O(log^{-4}R). Since these theorems are the central results, the proof needs to be expanded or a precise case-by-case table supplied. As written, the reader cannot check the claimed uniformity in the q_1-fixed regime.
minor comments (4)
- [Section 1, Eq. (1.11)] Eq. (1.11) states α_f(n)β_f(n)=1 for all n, which is not the standard multiplicative relation for composite n. For prime p, α_f(p)β_f(p)=1 is correct; for general n the Satake parameters are indexed by prime powers and one uses multiplicativity. Please rephrase.
- [Section 6, opening paragraph] Typographical error: 'Peterson trace formula' should be 'Petersson trace formula'.
- [Theorem 7.2, Eq. (7.7)] The constant γPNT3 is listed as 1+∫_1^∞ E(t)/t^2 dt ≈ -1.33258. Since this same constant is used again in Theorems 7.4 and 7.5, it would help to state once in Section 7 that E(t)=θ(t)-t and that all integrals over E(t) are convergent by the prime number theorem.
- [General notation] The functions φ_1,φ_2 are consistently assumed to be even Schwartz with Fourier transforms supported in [-σ,σ], but the notation sometimes writes φ̂''(0) for the second derivative of the Fourier transform and sometimes φ̂(0); this is standard but the dependence on the original φ should be stated more explicitly in the theorem statements.
Circularity Check
No circularity found: the lower-order terms are obtained by direct Petersson trace formula computation; the unevaluated root-number term in (1.19) is a completeness gap, not a circular input.
full rationale
The derivation chain is: explicit formula (1.12) plus Petersson trace formula estimates (Props 2.3–2.5 from ILS and Barrett et al.) and Hecke eigenvalue identities. Sections 4–5 express S1 and S2 as sums of moments A', A, B'', B', B, which Section 6 computes from the trace formula and Section 7 substitutes. The level-dependence of the lower-order terms is not assumed at the start: it emerges only after substituting the computed A' and B'' moments in the fixed-prime-factor case. The self-citations to Miller 2009 and Barrett et al. 2016 are to prior work whose cited ingredients are independent and externally checkable (the exact M3,0 expansion and the generalized trace formula with stated coprimality hypotheses); they do not assume the paper's target formulas. The numerical constants are integrals and convergents involving the prime-number-theorem error, not parameters fitted to the target densities. The real caveat is Eq. (1.19): the final term phi1(0)phi2(0) * sum_f wR(f)(1-eps_f)/(2 WR(F)) is never evaluated or bounded in Sections 4–7, so the assertion that computing Ak,N, S1, and S2 determines D2 up to O(log^-4 R) lacks a proof. That is a correctness/completeness problem, not circularity: the missing root-number average is an extra contribution, not a rearrangement of what was already assumed. Therefore the circularity score is 0.
Assumptions & free parameters
assumptions (6)
- domain assumption Generalized Riemann Hypothesis
- standard math Petersson trace formula estimates of ILS (Props 2.3 and 2.4) for squarefree level
- standard math Barrett et al. trace formula (Prop 2.5) for general level, used for N=p^2
- standard math Deligne bound |alpha_f(p)^m + beta_f(p)^m| <= 2
- standard math Prime Number Theorem with E(t) = theta(t)-t
- domain assumption Test functions with Fourier support in [-sigma,sigma], sigma < 0.22
Cite this review
Pith. "Pith review of Breaking Universality in the Lower Order Terms in the 1-level and 2-level Density of Holomorphic Cusp Newforms." pith.science (2026). https://pith.science/paper/SITEDAUY
@misc{pith2026250821691,
author = {Pith},
title = {Pith review of: Breaking Universality in the Lower Order Terms in the 1-level and 2-level Density of Holomorphic Cusp Newforms},
year = {2026},
howpublished = {\url{https://pith.science/paper/SITEDAUY}},
note = {Machine review of arXiv:2508.21691}
}
abstract
The Katz-Sarnak density conjecture states that, as the analytic conductor $R \to \infty$, the distribution of the normalized low-lying zeros (those near the central point $s = 1/2$) converges to the scaling limits of eigenvalues clustered near 1 of subgroups of $U(N)$. There is extensive evidence supporting this conjecture for many families, including the family of holomorphic cusp newforms. Interestingly, there are very few choices for the main term of the limiting behavior. In 2009, S. J. Miller computed lower-order terms for the 1-level density of families of elliptic curve $L$-functions and compared to cuspidal newforms of prime level; while the main terms agreed, the lower order terms depended on the arithmetic of the family. We extend his work by identifying family-dependent lower-order correction terms in the weighted 1-level and 2-level densities of holomorphic cusp newforms up to $O\left(1/\log^4 R\right)$ error, sharpening Miller's $O\left(1/\log^3 R\right)$ error. We consider cases where the level is prime or when the level is a product of two, not necessarily distinct, primes. We show that the rates at which the prime factors of the level tend to infinity lead to different lower-order terms, breaking the universality of the main behavior.
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.