Pith. sign in

REVIEW 1 major objections 8 minor 34 references

Random Multiplicative Functions and Making Squares from Polynomial Values

T0 review · 1 major / 8 minor · reviewed 2026-07-08 · glm-5.2

Pith's one-line read Gaussian limits proven for random multiplicative functions on polynomial values

desk verdict New CLTs for Rademacher and extended Rademacher RMFs along polynomial values, with sharp Diophantine estimates for the quadratic case read the letter →

arxiv 2607.06398 v1 pith:SWNPJI4A submitted 2026-07-07 math.NT math.PR

classification math.NTmath.PR
keywords functionsmultiplicativerademacherachievecentralcountingequationsestablish
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper proves that when you evaluate a polynomial P(n) at integers n = 1, 2, ..., N and feed those values into a random multiplicative function f (where f assigns random ±1 signs to each prime and extends multiplicatively), the partial sum Σ f(P(n)) behaves like a Gaussian random variable after normalization. The authors establish this central limit theorem for two types of random multiplicative functions — Rademacher (where f vanishes on non-square-free integers, modeling the Möbius function) and extended Rademacher (where f is defined completely multiplicatively). For the extended Rademacher case, the result is unconditional for all separable polynomials whose irreducible factors have degree at most 3. For the Rademacher case, the result holds under the assumption that a positive proportion of P(n) are square-free, which is known unconditionally in the same degree-at-most-3 setting. The core difficulty is not probabilistic but Diophantine: proving the central limit theorem reduces, via martingale methods, to showing that the equation P(n₁)P(n₂)P(n₃)P(n₄) = square has surprisingly few solutions where the nᵢ are all distinct — a paucity phenomenon. The paper establishes this paucity by combining tools from Diophantine geometry (bounds on integral points due to Evertse–Silverman and Reuss), results on Pell–Fermat equations (for the quadratic case, where the estimates are sharpest), and new combinatorial decompositions based on the complete graph K₄ that organize information about pairwise greatest common divisors of the P(nᵢ).

What carries the argument

The proof reduces the probabilistic central limit theorem to Diophantine counting via the McLeish martingale CLT. The key Diophantine tools are: (1) Hooley's results on the distribution of fundamental solutions to Pell–Fermat equations, used to bound off-diagonal solutions when deg P = 2; (2) Evertse–Silverman uniform bounds on integral points of superelliptic curves y² = d·P(x), used when deg P ≥ 3; (3) Reuss's determinant-method estimates for counting power-free values of cubic polynomials, used to control the set of n where P(n) has large square divisors (condition H_δ); (4) a combinatorial decomposition of the solution set based on labeling the six edges of the complete graph K₄ with gcd

What would settle it

A counterexample would be a separable polynomial P with all irreducible factors of degree ≤ 3 for which the normalized sum (1/√N) Σ g(P(n)) does not converge to N(0,1), or for which the off-diagonal solutions to P(n₁)P(n₂)P(n₃)P(n₄) = square are not o(N²).

Watch

Extended reading notes

Core claim

The central discovery is a paucity phenomenon for the Diophantine equation P(n₁)P(n₂)P(n₃)P(n₄) = square: among all quadruples (n₁, n₂, n₃, n₄) with 1 ≤ nᵢ ≤ N, the number of solutions where all four nᵢ are distinct is o(N²), which is small enough to confirm Gaussian behavior of the partial sums. For degree-2 polynomials, the bound is O(N^{3/2+ε}), proved using the multiplicative structure of quadratic norm forms and growth estimates for solutions of Pell–Fermat equations. For degree ≥ 3 (with irreducible factors of degree ≤ 3), the bound is O(N^{14/9+ε}), proved using Evertse–Silverman uniform bounds on integral points, Reuss's determinant method for power-free values, and a graph-theoretic

Load-bearing premise

The Rademacher case (Theorem 1.1) assumes that a positive proportion of polynomial values P(n) are square-free. This is known unconditionally only when every irreducible factor of P has degree ≤ 3; for general polynomials it follows from the ABC Conjecture or the Square-Free Sieve Conjecture.

Editorial extensions

If this is right

  • The unconditional extended Rademacher CLT for polynomials with irreducible factors of degree ≤ 3 confirms the Gaussian prediction of Chinis–Shala for a broad class of polynomials, extending their degree-2 result.
  • The paucity estimates for P(n₁)P(n₂)P(n₃)P(n₄) = square are of independent Diophantine interest and may find applications in analytic number theory beyond random multiplicative functions.
  • The function-field analog of the square-free sieve condition (H_δ) is known unconditionally for arbitrary polynomials, suggesting that a polynomial Rademacher CLT over function fields F_q(t) may now be obtainable.
  • For specific polynomial families like P(X) = b²Q(X)⁴ − 1, the variance computation simplifies dramatically (the off-diagonal count vanishes identically by Bennett–Walsh), pointing to special polynomial structures where the Diophantine obstruction is minimal.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

1 major / 8 minor

Summary. The paper proves central limit theorems for sums of Rademacher and extended Rademacher random multiplicative functions evaluated at polynomial values P(n). The main results (Theorems 1.1 and 1.2) confirm Gaussian limiting distributions for separable P with deg P >= 2, positive leading coefficient, and all irreducible factors of degree <= 3. The Rademacher case (Theorem 1.1) additionally assumes a positive density of square-free values (condition (2)), while the extended Rademacher case (Theorem 1.2) is unconditional. The proofs reduce, via the McLeish martingale CLT (following Soundararajan-Xu), to Diophantine counting problems: paucity of solutions to P(n_1)P(n_2) = square and P(n_1)P(n_2)P(n_3)P(n_4) = square. For deg P = 2, the authors develop a new method using Pell-Fermat equations (Propositions 5.1, 6.1). For deg P >= 3, they use Evertse-Silverman, Reuss's determinant method, and Bombieri-Pila (Propositions 8.1, 8.2, 3.6). A key optimization in Lemma 6.3 is verified using the Z3 SMT solver with publicly available code.

Significance. The paper makes substantial contributions to the study of random multiplicative functions along polynomial phases. The extended Rademacher CLT (Theorem 1.2) is unconditional and confirms the Gaussian prediction of Chinis-Shala for a broad class of polynomials, going beyond the Steinhaus case previously settled by Klurman-Shkredov-Xu. The Diophantine estimates are of independent interest: Proposition 5.1 gives a sharp O(sqrt(N)) bound for the variance problem in the quadratic case via Pell-Fermat theory, and Proposition 6.1 gives the key fourth-moment bound O(N^{3/2+epsilon}). The Z3-verified optimization in Lemma 6.3 (showing E <= 3/2) is independently checkable, and the code is publicly available. The K_4-based combinatorial decomposition of the fourth-moment counting problem is a novel and interesting technical contribution. The restriction to irreducible factors of degree <= 3 is clearly tied to the available Diophantine tools (Reuss, Evertse-Silverman) and is honestly disclosed.

major comments (1)
  1. Lemma 3.1 (p. 7): The hypothesis states that the coefficients of P have absolute value at most epsilon^{-1} N^{epsilon - 1}, and similarly for a, b. For fixed epsilon > 0 and large N, the bound epsilon^{-1} N^{epsilon - 1} tends to 0, which is unsatisfiable for integers a, b >= 1. This is almost certainly a typo for epsilon^{-1} N^{1 + epsilon} (or similar). Since Lemma 3.1 is used in Proposition 8.1 (p. 19) and Proposition 8.2 (p. 20) to bound solutions to aP(n_1) = bP(n_2), the correct hypothesis needs to be stated. The proof itself (via Bombieri-Pila) is standard and the bound is correct once the hypothesis is fixed, but the statement as written is technically vacuous.
minor comments (8)
  1. p. 7, Lemma 3.3: The remark mentions 'Gemini 3.1 Pro suggested to us in an incorrectly stated form' and 'GPT-5.5 Pro simplified the proof.' While acknowledging AI assistance is fine, the specific naming of AI models and their roles is unusual for a mathematics journal and could be streamlined.
  2. p. 8, Proposition 3.6: The constraint delta_0 < 1/13 arises from the interplay between the three ranges I_1, I_2, I_3. The choice of theta requires delta_0 < theta < min(1/11, (1 - 11*delta_0)/2). It would help the reader to explicitly verify that (1 - 11*delta_0)/2 > delta_0 when delta_0 < 1/13, ensuring the interval for theta is non-empty.
  3. p. 11, proof of Proposition 5.1: The factorization (14) and the subsequent case split are intricate. In the line after (15), 'we deduced << |Delta| N' should read 'we deduce << |Delta| N' for grammatical consistency.
  4. p. 14, Lemma 6.2: The bound (18) involves a factor d^epsilon. It would be helpful to clarify whether this is d^epsilon or (log d)^{O(1)}, as the proof uses 2^{omega(d)} which is bounded by d^epsilon but also by tau(d).
  5. p. 16, Lemma 6.3: The three bounds R, S_{ijk}, T_{ij} are defined in (23)-(27). The third bound T_{ij} uses the notation '2 - R + (beta_i + beta_j)/2'. In the proof, this is explained for (i,j) = (4,3) as '2 - (beta_1 + beta_2)/2 = 2 - R + (beta_4 + beta_3)/2'. This is correct but the general formula could be stated more explicitly for arbitrary (i,j).
  6. p. 19, Proposition 8.1: The bound M_4(N, A) << A^{11/9} N^{14/9 + epsilon} is stated for A <= N^5. The proof takes N^delta = A^{1/9} N^{4/9}, which requires A <= N^5 for delta <= 1. This is consistent but the role of the exponent 5 could be flagged more prominently.
  7. p. 22, Appendix A: The Z3 code is clearly presented. The sharpness check output gives specific values of alpha variables. It would be useful to include a brief human-readable verification that these values indeed give E = 3/2, to complement the machine check.
  8. References: The arXiv preprint [5] (Bhargava) and [2] (Alpoge) are cited with arXiv numbers. If published versions exist, they should be updated.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the positive assessment. The referee raises one major comment, which we address below.

read point-by-point responses
  1. Referee: Lemma 3.1 (p. 7): The hypothesis states that the coefficients of P have absolute value at most ε^{-1} N^{ε - 1}, and similarly for a, b. For fixed ε > 0 and large N, the bound ε^{-1} N^{ε - 1} tends to 0, which is unsatisfiable for integers a, b >= 1. This is almost certainly a typo for ε^{-1} N^{1 + ε} (or similar). Since Lemma 3.1 is used in Proposition 8.1 (p. 19) and Proposition 8.2 (p. 20) to bound solutions to aP(n_1) = bP(n_2), the correct hypothesis needs to be stated. The proof itself (via Bombieri-Pila) is standard and the bound is correct once the hypothesis is fixed, but the statement as written is technically vacuous.

    Authors: The referee is entirely correct. The bound ε^{-1} N^{ε-1} in the hypothesis of Lemma 3.1 is a typo; it should read ε^{-1} N^{1+ε}. With the corrected hypothesis, the condition is satisfiable for integers a, b ≥ 1 when N is large, and the proof via Bombieri–Pila goes through unchanged. We have verified that in all applications of Lemma 3.1—specifically in Propositions 8.1 and 8.2, and in the third bound of Lemma 6.3—the parameters a, b arise as products of square-free parts and their divisors, which are indeed bounded by N^{1+ε} (up to N^ε factors absorbed by the ε in the exponent). The polynomial P is fixed throughout, so the condition on its coefficients is automatically satisfied for large N. We will correct the statement in the revised manuscript. revision: yes

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity found; the derivation chain is self-contained against external benchmarks.

full rationale

The paper's central claims (Theorems 1.1 and 1.2) are derived from a chain of Diophantine estimates and a standard CLT framework (McLeish martingale CLT via Theorem 2.1/2.2). The key Diophantine inputs—Hooley's results on Pell equations (Proposition 5.1), Evertse–Silverman (Lemma 3.3), Reuss's determinant method (Proposition 3.6), and Bombieri–Pila (Lemma 3.1)—are all externally sourced results, not self-citations. The self-citations present ([34] Wang–Xu, [32] Soundararajan–Xu, [24] Klurman–Shkredov–Xu) provide the Steinhaus-case framework and paucity results that the present work extends, but the present results are not equivalent to those by construction. The Z3-verified optimization in Lemma 6.3 is independently checkable via the provided code repository. The CLT criteria in Theorem 2.1 are attributed to Soundararajan–Xu [32] but are a routine modification of a standard McLeish CLT application, not a self-serving definition. No step in the derivation chain reduces to its own inputs by construction.

Assumptions & free parameters 0 free parameters · 8 assumptions · 0 invented entities

No free parameters are fitted to data. No new mathematical entities (particles, forces, dimensions) are postulated. The axioms are all standard or domain-specific results from the prior literature, clearly cited. The main conditional assumption (H_δ / condition (2)) is explicitly stated as an open problem, not postulated as true.

assumptions (8)
  • standard math McLeish Martingale Central Limit Theorem (Theorem 2.1, §2): the standard probabilistic framework reducing CLT to moment conditions.
    Invoked in §2 to reduce the probabilistic claim to Diophantine counting. Standard tool in the field.
  • domain assumption Huxley's bound on polynomial congruence roots (eq. 8, §3): ρ_P(q) ≤ |disc(P)|^{1/2} (deg P)^{ω(q)}.
    Used throughout §3–8 to bound the number of solutions to polynomial congruences. External result [23].
  • domain assumption Evertse–Silverman uniform bound on integral points (Lemma 3.3, §3): #{(x,y) ∈ Z² : y² = dP(x)} ≤ 2^{O_P(1+ω(d))}.
    Used in Lemma 3.3 and Proposition 8.1 to bound solutions to P(x) = dy². External result [12].
  • domain assumption Reuss determinant-method estimate (§3, Proposition 3.6): bounds on N(N;A,B) for cubic factors.
    Used in Proposition 3.6 to establish H_δ for polynomials with irreducible factors of degree ≤ 3. External result [31].
  • domain assumption Condition (2): lim inf #{n ≤ N : μ²(P(n)) = 1}/N > 0 (positive density of square-free values).
    Load-bearing assumption for Theorem 1.1 (Rademacher case). Known for deg ≤ 3 unconditionally; requires ABC or Square-Free Sieve Conjecture in general.
  • domain assumption H_δ (condition (9), §3): B(N, N^{1/2}) ≪ N^{1−δ} — power-saving Square-Free Sieve Conjecture.
    Assumed in Theorem 8.3 for degree ≥ 4 extended Rademacher case. Proposition 3.6 proves H_δ for factors of degree ≤ 3.
  • standard math Bombieri–Pila bound on integral points on curves (Lemma 3.1, §3).
    Used to bound solutions to aP(n₁) = bP(n₂). External result [8].
  • domain assumption Hooley's result on Pell equations (§5, cited as [21, Theorem 1]).
    Used in Proposition 5.1 to bound S(D) and establish the variance estimate for quadratic P.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Random Multiplicative Functions and Making Squares from Polynomial Values." pith.science (2026). https://pith.science/paper/SWNPJI4A

@misc{pith2026260706398,
  author       = {Pith},
  title        = {Pith review of: Random Multiplicative Functions and Making Squares from Polynomial Values},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SWNPJI4A}},
  note         = {Machine review of arXiv:2607.06398}
}
abstract

For a large family of polynomials $P(X)\in \mathbb{Z}[X]$, we prove central limit theorems for $\sum_{n\le N} f(P(n))$ for both Rademacher and extended Rademacher multiplicative functions $f$. To achieve this, we establish a paucity phenomenon in counting solutions to \[P(n_1)P(n_2)P(n_3)P(n_4) = \square, \quad 1\le n_1, n_2, n_3, n_4 \le N.\] Results of Hooley, Evertse--Silverman, and Reuss play an important role in the proof. Our estimates are sharpest for $\deg P = 2$, thanks to the rich theory of Pell--Fermat equations.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

34 extracted references · 34 canonical work pages

  1. [1]

    Akhtari, The Diophantine equationaX 4 −bY 2 = 1,J

    S. Akhtari, The Diophantine equationaX 4 −bY 2 = 1,J. Reine Angew. Math.630(2009), 33–57. 10

  2. [2]

    Note on a theorem of Professor X

    L. Alp¨ oge, Note on a theorem of Professor X, arXiv preprint arXiv:2109.14328 (2021). 7

  3. [3]

    Basquin, Sommes friables de fonctions multiplicatives al´ eatoires,Acta Arith.152(2012), no

    J. Basquin, Sommes friables de fonctions multiplicatives al´ eatoires,Acta Arith.152(2012), no. 3, 243–

  4. [4]

    M. A. Bennett and P. G. Walsh, The Diophantine equationb 2X 4 −dY 2 = 1,Proc. Amer. Math. Soc. 127(1999), no. 12, 3481–3491. 22

  5. [5]

    The geometric sieve and the density of squarefree values of invariant polynomials

    M. Bhargava, The geometric sieve and the density of squarefree values of invariant polynomials, arXiv preprint arXiv:1402.0031 (2014). 17

  6. [6]

    Bhargava, A

    M. Bhargava, A. Shankar, T. Taniguchi, F. Thorne, J. Tsimerman, and Y. Zhao, Bounds on 2-torsion in class groups of number fields and integral points on elliptic curves,J. Amer. Math. Soc.33(2020), no. 4, 1087–1099. 7

  7. [7]

    Bombieri and W

    E. Bombieri and W. Gubler,Heights in Diophantine geometry, New Mathematical Monographs, 4, Cambridge Univ. Press, Cambridge, 2006. 7

  8. [8]

    Bombieri and J

    E. Bombieri and J. Pila, The number of integral points on arcs and ovals,Duke Math. J.59(1989), 337–357. 3, 7

Show all 34 references
  1. [9]

    A. R. Booker and T. D. Browning, Square-free values of reducible polynomials,Discrete Anal.2016, Paper No. 8, 18 pp. 4, 8, 9

  2. [10]

    Chinis and B

    J. Chinis and B. Shala, Random Chowla’s conjecture for Rademacher multiplicative functions,Trans. Amer. Math. Soc.378(2025), 8025–8053. 3, 4, 7, 11, 17

  3. [11]

    Drappeau, Remarques sur les moyennes des fonctions de Piltz sur les entiers friables,Quart

    S. Drappeau, Remarques sur les moyennes des fonctions de Piltz sur les entiers friables,Quart. J. Math. 67(2016), no. 4, 507–517. 18

  4. [12]

    Evertse and J

    J.-H. Evertse and J. H. Silverman, Uniform bounds for the number of solutions toY n =f(X),Math. Proc. Cambridge Philos. Soc.100(1986), no. 2, 237–248. 3, 7 24 R ´EGIS DE LA BRET `ECHE, VICTOR Y. WANG, AND MAX WENQIANG XU

  5. [13]

    Fouvry, On the size of the fundamental solution of the Pell equation,J

    ´E. Fouvry, On the size of the fundamental solution of the Pell equation,J. Reine Angew. Math.717 (2016), 1–33. 11

  6. [14]

    Fouvry and F

    ´E. Fouvry and F. Jouve, Fundamental solutions to Pell equation with prescribed size,Tr. Mat. Inst. Steklova276(2012), Teoriya Chisel, Algebra i Analiz, 46–56; translation inProc. Steklov Inst. Math. 276(2012), no. 1, 40–50. 11

  7. [15]

    Fouvry and F

    ´E. Fouvry and F. Jouve, Size of regulators and consecutive square-free numbers,Math. Z.273(2013), 869–882. 11, 14

  8. [16]

    A. J. Granville,ABCallows us to count squarefrees,Internat. Math. Res. Notices1998, no. 19, 991–

  9. [17]

    A. J. Granville, Rational and integral points on quadratic twists of a given hyperelliptic curve,Internat. Math. Res. Notices2007, no. 8, Art. ID 027, 25 pp. 8

  10. [18]

    A. J. Granville, Smooth numbers: computational number theory and beyond, inAlgorithmic num- ber theory: lattices, number fields, curves and cryptography, 267–323, Math. Sci. Res. Inst. Publ., 44, Cambridge Univ. Press, Cambridge, 2008. 18

  11. [19]

    A. J. Harper, On the limit distributions of some sums of a random multiplicative function.J. Reine Angew. Math.678(2013), 95–124. 5

  12. [20]

    Hooley, On the power free values of polynomials,Mathematika14(1967), 21–26

    C. Hooley, On the power free values of polynomials,Mathematika14(1967), 21–26. 4

  13. [21]

    Hooley, On the Pellian Equation and the Class Number of Indefinite Binary Quadratic Forms,J

    C. Hooley, On the Pellian Equation and the Class Number of Indefinite Binary Quadratic Forms,J. Reine Angew. Math.353(1984), 98–131. 3, 11, 13

  14. [22]

    Hua,Introduction to number theory, translated from the Chinese by Peter Shiu, Springer, Berlin- New York, 1982

    L.-K. Hua,Introduction to number theory, translated from the Chinese by Peter Shiu, Springer, Berlin- New York, 1982. 7

  15. [23]

    M. N. Huxley, A note on polynomial congruences. InRecent progress in analytic number theory, Vol. 1 (Durham, 1979), pages 193–196. Academic Press, London-New York, 1981. 3, 8

  16. [24]

    Klurman, I

    O. Klurman, I. D. Shkredov, and M. W. Xu, On the random Chowla conjecture,Geom. Funct. Anal. 33(2023), no. 3, 749–777. 2, 17

  17. [25]

    Y.-K. Lau, G. Tenenbaum and J. Wu, On mean values of random multiplicative functions,Proc. Amer. Math. Soc.141(2013), no. 2, 409–420. 5

  18. [26]

    Lemmermeyer, The ambiguous class number formula revisited,J

    F. Lemmermeyer, The ambiguous class number formula revisited,J. Ramanujan Math. Soc.28(2013), no. 4, 415–421. 7

  19. [27]

    S. J. Miller, One- and two-level densities for rational families of elliptic curves: evidence for the under- lying group symmetries,Compos. Math.140(2004), no. 4, 952–992. 4

  20. [28]

    Nagell,Introduction to number theory

    T. Nagell,Introduction to number theory. 2nd ed., Chelsea, 1964. 8

  21. [29]

    Poonen, Squarefree values of multivariable polynomials,Duke Math

    B. Poonen, Squarefree values of multivariable polynomials,Duke Math. J.118(2003), no. 2, 353–373. 4, 21

  22. [30]

    Ramsay, Square-free values of polynomials in one variable over function fields,Internat

    K. Ramsay, Square-free values of polynomials in one variable over function fields,Internat. Math. Res. Notices1992, no. 4, 97–102. 4, 21

  23. [31]

    Reuss, Power-free values of polynomials,Bull

    T. Reuss, Power-free values of polynomials,Bull. Lond. Math. Soc.47(2015), no. 2, 270–284. 3, 4, 9

  24. [32]

    Soundararajan and M

    K. Soundararajan and M. W. Xu, Central limit theorems for random multiplicative functions,J. Anal. Math.151(2023), no. 1, 343–374. 5, 6

  25. [33]

    R. C. Vaughan and T. D. Wooley, Further improvements in Waring’s problem,Acta Math.174(1995), no. 2, 147–240. 7

  26. [34]

    V. Y. Wang and M. W. Xu, Paucity phenomena for polynomial products,Bull. Lond. Math. Soc.56 (2024), no. 8, 2718–2726. 3, 7 Universit´e Paris Cit´e, Sorbonne Universit ´e, CNRS UMR 7586, Institut Universitaire de France, Institut de Math ´ematiques de Jussieu-Paris Rive Gauche,...

Pith tools

Reviewed July 8, 2026 · model on record in the stance chip above.