Matthew J. Colbrook,Andrew Horning,Tianyiwa Xie
Abstract:We introduce a simple, general, and convergent scheme to compute generalized eigenfunctions of self-adjoint operators with continuous spectra on rigged Hilbert spaces. Our approach does not require prior knowledge about the eigenfunctions, such as asymptotics or other analytic properties. Instead, we carefully sample the range of the resolvent operator to construct smooth and accurate wave packet approximations to generalized eigenfunctions. We prove high-order convergence in key topologies, including weak-star convergence for distributional eigenfunctions, uniform convergence on compact sets for locally smooth generalized eigenfunctions, and convergence in seminorms for separable Frechet spaces, covering the majority of physical scenarios. The method's performance is illustrated with applications to both differential and integral operators, culminating in the computation of spectral measures and generalized eigenfunctions for an operator associated with Poincare's internal waves problem. These computations corroborate experimental results and highlight the method's utility for a broad range of spectral problems in physics.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to calculate the generalized eigenfunctions of self - adjoint operators with continuous spectra in rigged Hilbert spaces. Specifically, the authors propose a simple, general, and convergent scheme to calculate these generalized eigenfunctions without prior knowledge of the specific properties of the eigenfunctions, such as asymptotic behavior or other analytic characteristics.
### Main problems
1. **Calculation of generalized eigenfunctions**: For self - adjoint operators with continuous spectra, traditional eigenvectors cannot form a basis for Hilbert spaces. Therefore, a new method is required to calculate generalized eigenfunctions.
2. **Challenges in numerical calculation**: Generalized eigenfunctions are linear functionals on infinite - dimensional spaces. How to perform effective numerical approximation with limited resources on a computer is a key issue.
3. **Theoretical basis**: Traditional spectral measure calculation theory is based on Hilbert spaces, but non - trivial generalized eigenfunctions are not elements in Hilbert spaces. How to provide a strict theoretical basis for such calculations is also a challenge.
4. **Capturing fine properties**: In many practical cases, generalized eigenfunctions have important additional regularity (e.g., local smoothness). How to ensure that general approximation schemes can capture these fine properties?
### Solutions
The authors solve the above problems in the following ways:
- **Wave - packet approximation**: By carefully sampling the resolvent operator on the complex plane, construct smooth wave - packet approximations of generalized eigenfunctions consisting of narrow - band spectral content.
- **High - order convergence**: Prove the weak - \* topology, uniform convergence on compact sets, and seminorm convergence in separable Fréchet spaces for these wave - packets, which are applicable to most physical scenarios.
- **Application examples**: Demonstrate the effectiveness of this method through the application of differential and integral operators, and finally apply it to calculate the spectral measures and generalized eigenfunctions of operators related to the Poincaré interior wave problem.
### Mathematical formula representation
Some important formulas involved in the paper include:
1. **Resolvent operator**:
\[
R_A(z)=(A - z)^{-1}
\]
where \( A \) is a self - adjoint operator, \( z\notin\Lambda(A) \), and its operator norm is:
\[
\|R_A(z)\|=\frac{1}{\text{dist}(z,\Lambda(A))}\leq\frac{1}{|\text{Im}(z)|}
\]
2. **Wave - packet approximation**:
\[
u_\lambda^{(\epsilon)}=\int_{\mathbb{R}}K_\epsilon(\lambda - \tilde{\lambda})dE(\tilde{\lambda})f,\quad f\in H
\]
where \( K_\epsilon \) is a family of approximate identity operators or smoothing kernels.
3. **Convolution form of Poisson kernel**:
\[
u_\lambda^{(\epsilon)}=\frac{1}{2\pi i}[R_A(\lambda - i\epsilon)f - R_A(\lambda + i\epsilon)f]
\]
4. **High - order rational convolution kernel**:
\[
K(x)=\frac{1}{2\pi i}\sum_{j = 1}^m\left(\frac{\alpha_j}{x - a_j}-\frac{\alpha_j}{x - \bar{a}_j}\right)
\]
where \( a_1,\ldots,a_m \) are distinct points in the upper half - plane, and \( \alpha_1,\ldots,\alpha_m \) are their corresponding complex residues.
Through these methods, the authors provide a theoretically rigorous and easy - to - implement tool for calculating generalized eigenfunctions in rigged Hilbert spaces, thus opening a new door for mathematicians and scientists in computational spectral analysis.