A Tauberian approach to the orthorecursive expansion of unity
Benoit Cloitre
TL;DR
This paper establishes a Tauberian framework that links the unknown asymptotics of unweighted partial sums $A(x)=\sum_{n\le x} a_n$ to known asymptotics of weighted sums $A_g(x)=\sum_{n\le x} a_n g(n/x)$ via a modified Mellin transform $g^{\star}$ and its zeros. The key device is the identity $L(a,s)=s\frac{A_g^{*}(-s)}{g^{\star}(-s)}$, enabling a Perron-type contour analysis to translate weighted behavior into unweighted asymptotics, with the asymptotics governed by $\alpha_g$ (the leftmost zero of $g^{\star}$) and the main term parameters. The authors apply this to the orthorecursive expansion of unity studied by Kalmynin and Kosenko, transforming the recurrence into a weighted-sum setting with $g(t)=\frac{2}{1+t}$ and showing $A_g(x)=\mathcal{O}(x^{-3/2})$. They identify the principal zero of $g^{\star}$ as $\alpha_1=\operatorname{Re}(\rho_1)\approx1.3465$, and deduce the optimal decay $C_N=\sum_{n\le N} c_n=\mathcal{O}(N^{-\alpha_1+\varepsilon})$ for any $\varepsilon>0$, thereby resolving the open problem. The results highlight a versatile Tauberian method that combines Mellin–Flajolet techniques with zero analysis of digamma-structured transforms, with potential applicability to other $L^2$-approximation problems and related Tauberian questions.
Abstract
We establish a Tauberian theorem connecting the unknown asymptotic behavior of the partial sums $\sum_{n\le x}a_{n}$ to the known asymptotics of weighted sums $\sum_{n\le x}a_{n}g(n/x)$, as $x\rightarrow\infty$, where $g:(0,1]\to\mathbb{R}$ is a given function. Our approach relies on an identity relating a modified Mellin transform of $g$ to the Dirichlet series $\sum_{n\ge1}a_{n}n^{-s}$. As an application, we solve an open problem posed by Kalmynin and Kosenko regarding the "orthorecursive expansion of unity" associated with a sequence $(c_{n})_{n\geq0}$. Specifically, we improve their partial-sum bound $C_{N}=\sum_{0\leq n\leq N}c_{n}=\mathcal{O}(N^{-1/2})$, by obtaining the optimal estimate $C_{N}=\mathcal{O}(N^{-α_{1}+ε})$, where $α_{1}\approx1.3465165$ is the smallest real part among the zeros of a transcendental function related to the digamma function.
