The local limit theorem for complex valued sequences: the parabolic case

Jean-François Coulombel,Grégory Faye
2024-08-23
Abstract:We give a complete expansion, at any accuracy order, for the iterated convolution of a complex valued integrable sequence in one space dimension. The remainders are estimated sharply with generalized Gaussian bounds. The result applies in probability theory for random walks as well as in numerical analysis for studying the large time behavior of numerical schemes.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the local limit theorem of complex - valued sequences, especially the iterative convolution expansion in the parabolic case. Specifically, the author aims to provide a complete asymptotic expansion of the iterative convolution of complex - valued integrable sequences in one - dimensional space and make an accurate estimate of the remainder term, using the generalized Gaussian bound. This result is not only of great significance for the study of random walks in probability theory, but also helpful for studying the large - time behavior of numerical schemes in numerical analysis. ### Main problem summary: 1. **Asymptotic expansion**: The goal of the paper is to give a complete asymptotic expansion of the one - dimensional iterative convolution of complex - valued sequences, with an accuracy that can reach any order. 2. **Remainder term estimate**: Make an accurate estimate of the remainder term, using the generalized Gaussian bound to ensure the optimality of the estimate. 3. **Application fields**: This result is applicable to the study of random walks in probability theory and the large - time behavior in numerical analysis. ### Specific problem description: - **Probability theory background**: In probability theory, the local limit theorem gives the asymptotic expansion of the probability of the sum of discrete - time random variables. Usually, these expansion formulas correspond to the increasing powers of \(n^{- 1/2}\), and the main term is a properly scaled Gaussian function. - **Convolution representation**: When all random variables are independent and identically distributed, the probability \(P(X_1+\cdots+X_n = j)\) corresponds to the value at index \(j\) of the \((n - 1)\)-th iterative convolution of the sequence \((P(X_1=\ell))_{\ell\in\mathbb{Z}}\). - **Beyond the probability framework**: This paper attempts to give a complete asymptotic expansion of such iterative convolution without assuming positivity, thus expanding the framework of traditional probability theory. - **Parabolic and dispersive behaviors**: The paper distinguishes between parabolic and dispersive behaviors, and this paper focuses on the parabolic case, that is, the stable case where the iterative convolution is bounded in the \(\ell^1\) norm. ### Mathematical formulas: - **Fourier series**: \[ F_a:\kappa\in\mathbb{C}\mapsto\sum_{\ell\in\mathbb{Z}}a_\ell\kappa^\ell \] - **Convolution definition**: \[ a\star b:=\left(\sum_{\ell'\in\mathbb{Z}}a_{\ell-\ell'}b_{\ell'}\right)_{\ell\in\mathbb{Z}} \] - **Main result of asymptotic expansion**: \[ a_\ell^{(n)}\sim\frac{1}{\sqrt{4\pi\beta n}}\exp\left(-\frac{(\ell-\alpha n)^2}{4\beta n}\right)+\sum_{m\geq1}\frac{1}{n^{(m + 1)/2}}\left(P_m(-d/dx)H_\beta^2\right)\left(\frac{\ell-\alpha n}{\sqrt{n}}\right) \] Through these formulas and methods, the author is able to provide a broader theoretical framework to deal with the iterative convolution problem of complex numerical sequences and provide a basis for error estimation in numerical analysis.