Rapidly convergent series expansions for a class of resolvents

Graeme W. Milton
2024-08-01
Abstract:Following advances in the abstract theory of composites, we develop rapidly converging series expansions about $z=\infty$ for the resolvent ${\bf R}(z)=[z{\bf I}-{\bf P}^\dagger{\bf Q}{\bf P}]^{-1}$ where ${\bf Q}$ is an orthogonal projection and ${\bf P}$ is such that ${\bf P}{\bf P}^\dagger$ is an orthogonal projection. It is assumed that the spectrum of ${\bf P}^\dagger{\bf Q}{\bf P}$ lies within the interval $[z^-,z^+]$ for some known $z^+\leq 1$ and $z^-\geq 0$ and that the actions of the projections ${\bf Q}$ and ${\bf P}{\bf P}^\dagger$ are easy to compute. The series converges in the entire $z$-plane excluding the cut $[z^-,z^+]$. It is obtained using subspace substitution, where the desired resolvent is tied to a resolvent in a larger space and ${\bf Q}$ gets replaced by a projection $\underline{\bf Q}$ that is no longer orthogonal. When $z$ is real the rate of convergence of the new method matches that of the conjugate gradient method.
Numerical Analysis,Mathematical Physics,Analysis of PDEs,Classical Analysis and ODEs
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to develop a rapidly convergent series expansion method for the inverse matrix of a class of analytic functions (i.e., the solutions of analytic functions on the complex plane). Specifically, the paper focuses on analytic functions of the form \( R(z) = [zI - P^\dagger QP]^{-1} \), where \( Q \) is an orthogonal projection operator and \( P \) is an operator such that \( PP^\dagger \) becomes an orthogonal projection. ### Core Problems of the Paper 1. **Rapidly Convergent Series Expansion**: - The author aims to develop a series expansion method that converges rapidly near \( z = \infty \) for calculating \( R(z) \). - This series expansion is very important for applications in numerical analysis and composite material theory, especially when dealing with large - scale systems. 2. **Spectral Constraints**: - Assume that the spectrum of \( P^\dagger QP \) lies in the interval \([z_-, z_+]\), and \( z_+ \leq 1 \) and \( z_- \geq 0 \). - These assumptions simplify the problem and ensure that the series converges in the entire region of the \( z \)-plane except for the cut line \([z_-, z_+]\). 3. **Subspace Replacement Technique**: - By introducing the technique of non - orthogonal subspace replacement, the required analytic function is related to the analytic function in a larger space. - This method can significantly improve the convergence rate of the series, especially when \( z \) is real, and its convergence rate matches that of the conjugate gradient method. ### Markdown Representation of Formulas - \( R(z) = [zI - P^\dagger QP]^{-1} \) - \( I \geq A \geq 0 \) - \( z_+ I \geq P^\dagger AP \geq z_- I \) - \( \sigma=\frac{z - 1}{z} \) ### Key Steps 1. **Rescaling and Transformation**: - By appropriately rescaling and transforming \( A \) and \( z \), make \( I \geq A \geq 0 \). - After using these transformations, \( R(z) = [zI - P^\dagger AP]^{-1} \) is still satisfied. 2. **Equivalent Equations**: - Transform the original equation into the equivalent form \([Γ(I - A/z)]E = h \), where \( E = Py=\Gamma E \), \( h = Pb/z=\Gamma h \). 3. **Series Expansion**: - Using the subspace replacement technique, transform the problem into solving a new series expansion form. - The new series expansion has a faster convergence rate, especially when dealing with specific types of spectral ranges. ### Conclusion By introducing the subspace replacement technique and making reasonable assumptions about the spectral range, this paper successfully develops a rapidly convergent series expansion method. This method is not only of great theoretical significance but also shows superior performance in practical applications, especially when dealing with complex problems in composite materials and numerical analysis.