Fact-checked by Grok 2 weeks ago

Neumann series

In mathematics, particularly in functional analysis and linear algebra, the Neumann series is an infinite power series expansion that provides the inverse of a bounded linear operator I - A on a Banach space, where I is the identity operator and \|A\| < 1 in the operator norm, given by (I - A)^{-1} = \sum_{k=0}^{\infty} A^k.^[1] This series converges absolutely under the norm condition, analogous to the geometric series \sum_{k=0}^{\infty} z^k = \frac{1}{1-z} for |z| < 1, and it serves as a foundational tool for solving linear equations involving perturbations of the identity.^[2] Introduced by Carl Neumann in 1877 within the context of potential theory, the series was originally developed to address problems in integral equations and has since become integral to operator theory.^[3] Key properties include its guaranteed convergence when the spectral radius of A is less than 1, ensuring the invertibility of I - A in the space of bounded operators, and its extension to more general settings like higher-norm topologies under suitable conditions on A.^[1] The series is particularly notable for its role in iterative methods, such as the Neumann iteration x_{k} = A x_{k-1} + b for solving (I - A)x = b, where convergence rates depend on the eigenvalues of A.^[3] Applications of the Neumann series span diverse fields, including the solution of Fredholm integral equations of the second kind, where it expands the resolvent kernel, and computational techniques like thermoacoustic tomography, where efficient algorithms leverage the series for inverse problems.^[4] In wave scattering problems, such as the Helmholtz equation, it reformulates the Lippmann-Schwinger equation as u_s = (I - \hat{S}(q \cdot))^{-1} \hat{S}(q u_i), converging for weak scatterers when \|q\|_{L^\infty} < 1/(2 C_1(k, \Omega)).^[2] More advanced variants, like the Liouville-Neumann series, extend its utility to nonlinear and probabilistic analyses, highlighting its enduring influence in modern mathematics and physics.^[3]

Introduction

Historical background

The Neumann series was introduced by Carl Gottfried Neumann in 1877 as a method for solving integral equations of the second kind arising in potential theory, particularly in the context of the Dirichlet problem for elliptic partial differential equations.^[5] In his work Untersuchungen über das logarithmische und Newton'sche Potential, Neumann developed this iterative series expansion to address boundary value problems, building on earlier ideas in electrostatics and gravitation while providing a rigorous approach to convergence under suitable conditions on the kernel.^[6] This innovation marked a significant advancement in handling integral equations that model physical potentials, offering a perturbative solution technique that could be applied to logarithmic and Newtonian potentials in two and three dimensions.^[7] Carl Gottfried Neumann (1832–1925), a prominent German mathematician, was born in Königsberg, Prussia and studied at the University of Königsberg under his father Franz Ernst Neumann, Ernst Kummer, and Carl Jacobi, earning his doctorate in 1855.^[8] He held professorships at Halle, Tübingen, and Leipzig, where he contributed to the formalization of potential theory and its applications to physics, including electrodynamics and analytical mechanics.^[8] Neumann's broader mathematical legacy includes a proof of the Dirichlet principle for certain elliptic boundary value problems in 1868, which addressed existence questions for solutions to Laplace's equation, and his role as a co-founder and editor of Mathematische Annalen starting in 1871, which became a leading journal for pure and applied mathematics.^[8] His work bridged classical analysis and emerging areas of partial differential equations, influencing the development of modern boundary value problem techniques. In the late 19th and early 20th centuries, the Neumann series gained traction among mathematicians studying integral equations, notably through David Hilbert's foundational contributions to the field around 1904–1910.^[9] Hilbert incorporated and extended Neumann's iterative method in his theory of Fredholm integral equations of the second kind, using it to expand solutions in series form and establish spectral properties that laid the groundwork for functional analysis.^[10] This adoption helped integrate the series into the broader framework of operator theory, facilitating solutions to problems in physics and engineering that resisted direct analytical methods.^[9]

Overview

The Neumann series provides a method to express the inverse of an operator of the form I - T, where I is the identity operator and T is a bounded linear operator on a Banach space, as an infinite sum \sum_{k=0}^\infty T^k, provided that the operator norm satisfies \|T\| < 1. This series converges in the operator norm, yielding (I - T)^{-1} = \sum_{k=0}^\infty T^k, which allows for the explicit construction of the inverse without direct inversion techniques.^[11]^[12] In the scalar case, this construction reduces to the familiar geometric series \sum_{k=0}^\infty r^k = (1 - r)^{-1} for |r| < 1, serving as a motivating analogy that extends naturally to operators in infinite-dimensional spaces. The series plays a crucial role in solving linear equations of the form (I - T)x = y by expressing the solution as x = \sum_{k=0}^\infty T^k y, facilitating iterative approximations with controlled error in settings where direct methods are infeasible.^[11]^[12]^[13] By generalizing finite-dimensional matrix inversion techniques to abstract operator theory, the Neumann series bridges classical linear algebra with functional analysis, enabling the analysis of problems in partial differential equations and integral equations within Hilbert or Banach spaces. Named after the mathematician Carl Neumann who introduced it in the 19th century, it remains a foundational tool for establishing invertibility and resolvent properties in operator algebras.^[12]^[12]

Formulation

Definition

In a normed vector space X over the real or complex numbers, equipped with a norm \|\cdot\|, the Neumann series of a bounded linear operator T: X \to X is the formal infinite sum

S = \sum_{k=0}^{\infty} T^k,

where T^0 = I denotes the identity operator on X, and T^k for k \geq 1 is the k-fold composition of T with itself.^[14] This series expands explicitly as

S = I + T + T^2 + T^3 + \cdots.

The partial sums are defined as

S_n = \sum_{k=0}^{n} T^k = I + T + \cdots + T^n,

with the full series converging to S in the operator norm if \lim_{n \to \infty} S_n = S exists.^[14] When X is a Banach space—a complete normed vector space—the operator norm \|T\| = \sup_{\|x\| \leq 1} \|Tx\| provides a measure of the size of T, facilitating analysis of the series in the space B(X) of bounded linear operators on X, which itself forms a Banach space under this norm.^[15]^[14]

Convergence criteria

The Neumann series \sum_{k=0}^{\infty} T^k for a bounded linear operator T on a Banach space converges absolutely in the operator norm if \sum_{k=0}^{\infty} \|T^k\| < \infty, which guarantees convergence to a bounded operator.^[16] This absolute convergence implies pointwise convergence on the space and norm convergence of the series due to the completeness of the Banach algebra of bounded operators.^[17] A sufficient condition for convergence in Banach spaces is \|T\| < 1, under which \|T^k\| \leq \|T\|^k for all k, ensuring the series behaves like a geometric series with ratio less than 1.^[16] In this case, the partial sums S_n = \sum_{k=0}^{n} T^k satisfy the remainder estimate \|S - S_n\| \leq \|T\|^{n+1} / (1 - \|T\|), where S is the sum of the series, providing a quantitative bound on the approximation error.^[16] More generally, the series converges in the operator norm if the spectral radius \rho(T) < 1, where \rho(T) = \limsup_{n \to \infty} \|T^n\|^{1/n} by Gelfand's formula.^[17] This condition extends the norm-based criterion, as \rho(T) \leq \|T\| always holds, but \rho(T) < 1 allows convergence even if \|T\| \geq 1, provided the powers decay sufficiently fast.^[17]

Properties

Inversion formula

The Neumann series yields the inverse of the operator I - T whenever the series converges. Let S = \sum_{k=0}^\infty T^k denote the sum of the Neumann series, where T is a bounded linear operator on a Banach space and the series converges in the operator norm. Then S = (I - T)^{-1}.^[16] To establish this, compute the product (I - T)S. Because the series converges absolutely, term-by-term multiplication is valid, yielding

(I - T) \sum_{k=0}^\infty T^k = \sum_{k=0}^\infty T^k - \sum_{k=0}^\infty T^{k+1} = \sum_{k=0}^\infty (T^k - T^{k+1}).

The right-hand side forms a telescoping sum: all terms T^{k+1} for k \geq 0 cancel with subsequent -T^{k+1}, leaving the identity operator I. An analogous computation shows S(I - T) = I.^[16] A standard sufficient condition ensuring convergence (and thus invertibility) is \|T\| < 1 with respect to the operator norm, in which case

(I - T)^{-1} = \sum_{k=0}^\infty T^k.

This inverse is unique within the algebra of bounded linear operators, since the invertible elements form a group under operator composition.^[16]

Resolvent properties

In spectral theory, the resolvent operator associated with a bounded linear operator T on a Banach space is defined as R(\lambda, T) = (\lambda I - T)^{-1} for \lambda in the resolvent set \rho(T), which is the set of complex numbers \lambda for which \lambda I - T is invertible.^[18] When \lambda = 1, this reduces to the case of the Neumann series expansion for (I - T)^{-1} = \sum_{k=0}^{\infty} T^k, provided that the spectral radius \rho(T) < 1. This parameterization allows the Neumann series to be generalized to the resolvent form R(\lambda, T) = \lambda^{-1} \sum_{k=0}^{\infty} (T/\lambda)^k, connecting iterative approximations to the broader spectral structure of T.^[18] The resolvent R(\lambda, T) is holomorphic as a function of \lambda on the complement of the spectrum \sigma(T), the set of \lambda where \lambda I - T fails to be invertible. This analyticity follows from the fact that small perturbations in \lambda away from \sigma(T) preserve invertibility, allowing local power series expansions of the resolvent via the Neumann series in suitable disks.^[18] Such holomorphy enables the application of complex analysis tools, like contour integrals, to extract spectral projections and eigenvalues from the resolvent.^[19] The Neumann series for the resolvent converges precisely when |\lambda| > \rho(T), where \rho(T) = \sup \{ |\mu| : \mu \in \sigma(T) \} is the spectral radius of T. In this region, the operator norm satisfies \|(T/\lambda)^k\| \to 0 as k \to \infty, ensuring the series sums to R(\lambda, T).^[18] This convergence domain highlights the role of the spectrum in delimiting where the iterative expansion is valid, with the boundary |\lambda| = \rho(T) often marking the edge of analytic continuation.^[20] Around an isolated eigenvalue \mu \in \sigma(T), the resolvent admits a Laurent series expansion R(\lambda, T) = \sum_{n=-m}^{\infty} A_n (\lambda - \mu)^n, where m is the order of the pole (typically 1 for simple eigenvalues), and the coefficients A_{-1}, \dots, A_{-m} form the principal part related to the generalized eigenspace projection.^[18] The regular (holomorphic) part \sum_{n=0}^{\infty} A_n (\lambda - \mu)^n can be expressed using a Neumann series adapted to the reduced resolvent on the complement of the eigenspace, facilitating the decomposition of T into its spectral components.^[21] This expansion is central to perturbation analyses of isolated eigenvalues under operator perturbations.^[18]

Examples

Matrix example

To illustrate the Neumann series in the finite-dimensional case, consider the 2×2 matrix

T = \begin{pmatrix} 0.4 & 0.3 \\ 0.2 & 0.5 \end{pmatrix},

which satisfies \|T\|_\infty = 0.7 < 1, where \|\cdot\|_\infty denotes the induced infinity norm (maximum absolute row sum).^[22] The Neumann series \sum_{k=0}^\infty T^k then converges to (I - T)^{-1}.^[23] The exact inverse is

(I - T)^{-1} = \begin{pmatrix} 2.0833 & 1.25 \\ 0.8333 & 2.5 \end{pmatrix},

computed via the formula for the inverse of a 2×2 matrix.^[23] The partial sums S_n = \sum_{k=0}^n T^k approximate this inverse, with explicit computations as follows:

S_0 = I = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix},

T^1 = T = \begin{pmatrix} 0.4 & 0.3 \\ 0.2 & 0.5 \end{pmatrix}, \quad S_1 = \begin{pmatrix} 1.4 & 0.3 \\ 0.2 & 1.5 \end{pmatrix},

T^2 = \begin{pmatrix} 0.22 & 0.27 \\ 0.18 & 0.31 \end{pmatrix}, \quad S_2 = \begin{pmatrix} 1.62 & 0.57 \\ 0.38 & 1.81 \end{pmatrix},

T^3 = \begin{pmatrix} 0.142 & 0.201 \\ 0.134 & 0.209 \end{pmatrix}, \quad S_3 = \begin{pmatrix} 1.762 & 0.771 \\ 0.514 & 2.019 \end{pmatrix}.

These partial sums approach the exact inverse as n increases.^[23] The error \| (I - T)^{-1} - S_n \|_\infty decreases geometrically with n. For instance, the errors are approximately 2.333 for n=0, 1.633 for n=1, 1.143 for n=2, and 0.800 for n=3. This aligns with the bound \| (I - T)^{-1} - S_n \|_\infty \leq \|T\|_\infty^{n+1} / (1 - \|T\|_\infty) = 0.7^{n+1} / 0.3, which yields exactly 2.333, 1.633, 1.143, and 0.800 for the respective n.^[22]

Operator example

A prominent example of the Neumann series in infinite-dimensional spaces arises with the Volterra integral operator, defined on functions in the space L^2[0,1] by

(Tf)(x) = \int_0^x K(x,y) f(y) \, dy,

where the kernel K(x,y) is continuous and satisfies |K(x,y)| \leq M for some M > 0, ensuring the operator is compact.^[24]^[25] To guarantee the norm condition \|T\| < 1, select a kernel scaled such that M \cdot 1 < 1, as the operator norm satisfies \|T\| \leq M on the unit interval due to the bounded kernel and finite integration length; this holds particularly for compact Volterra operators, where the triangular structure limits the norm relative to the interval size.^[24] Consider solving the equation (I - T)f = g for f \in L^2[0,1], where g is given. The Neumann series provides the solution as f = \sum_{n=0}^\infty T^n g. This series is generated via successive approximations: start with f_0 = g, and iterate f_{n+1} = g + T f_n, yielding f_n = \sum_{k=0}^n T^k g, with the limit f = \lim_{n \to \infty} f_n as the full sum.^[25]^[24] For a concrete illustration, take K(x,y) = \frac{1}{2} (so M = \frac{1}{2} < 1) and g(x) = 1. The first few iterates are:

f_0(x) = 1,

f_1(x) = 1 + \frac{1}{2} \int_0^x 1 \, dy = 1 + \frac{x}{2},

f_2(x) = 1 + \frac{1}{2} \int_0^x \left(1 + \frac{y}{2}\right) dy = 1 + \frac{x}{2} + \frac{1}{8} x^2.

The series converges in the L^2[0,1] norm to the exact solution f(x) = e^{x/2}, which satisfies the original equation f(x) - \frac{1}{2} \int_0^x f(y) \, dy = 1, with the error bounded by \|f - f_n\| \leq \frac{(M)^{n+1}}{1 - M} \|g\|, reflecting the geometric convergence under \|T\| < 1. This convergence follows from the condition \|T\| < 1, which ensures geometric decay of the norms \|T^n\| \leq \|T\|^n. While the compactness of Volterra operators provides additional analytical benefits, such as approximation by finite-rank operators, it is not essential for the Neumann series convergence here.^[25]^[24] This process mirrors Picard iteration, serving as its continuous counterpart in function spaces, whereas the matrix example provides the finite-dimensional discrete analog for intuition.^[26]^[25]

Applications

Approximate inversion

The truncated Neumann series provides a practical method for approximating the inverse of an invertible matrix A when a scaling can be applied such that \|I - A\| < 1 in a suitable matrix norm. In this setting, the inverse is given by the infinite series A^{-1} = \sum_{k=0}^\infty (I - A)^k, and a finite truncation yields the approximation A^{-1} \approx S_n = \sum_{k=0}^n (I - A)^k. This approach is particularly useful in numerical linear algebra for matrices where direct inversion via methods like Gaussian elimination would be computationally prohibitive. The primary computational advantage of the truncated Neumann series lies in its iterative nature, which transforms the problem of explicit matrix inversion—typically requiring O(n^3) operations for an n \times n dense matrix—into a sequence of matrix-vector multiplications or similar operations costing O(n^2) per iteration for dense cases, and O(\mathrm{nnz}) (where \mathrm{nnz} is the number of nonzeros) for sparse matrices. This makes it especially effective for large-scale sparse systems arising in scientific computing, where only a modest number of terms n is needed to achieve sufficient accuracy under the convergence condition. For instance, in applications involving diagonally dominant or well-conditioned matrices after scaling, convergence can occur rapidly, avoiding the full factorization overhead of direct solvers. Error analysis for the approximation is grounded in the remainder of the geometric series. Specifically, the truncation error satisfies \|A^{-1} - S_n\| \leq \frac{\|I - A\|^{n+1}}{1 - \|I - A\|} \cdot \|A^{-1}\|, providing a posteriori or a priori bounds that guide the choice of n based on the spectral properties or norm estimates of A. This bound ensures that the approximation error decreases geometrically with n, contingent on the initial scaling to enforce \|I - A\| < 1. In practice, monitoring the residual or using adaptive truncation can further refine the estimate. Implementation of the truncated Neumann series is well-suited for preconditioned iterative solvers in numerical linear algebra, where it serves as an approximate inverse preconditioner for systems like Ax = b. By approximating the preconditioner matrix as a low-order polynomial in A via the series, it facilitates efficient iterations in methods such as GMRES or conjugate gradients, particularly for ill-conditioned or sparse systems where exact inverses are infeasible. This integration enhances convergence rates without requiring sparse LU factorizations, though care must be taken in scaling A (e.g., via diagonal shifts) to satisfy the norm condition.

Integral equations

The Neumann series provides a method for solving Fredholm integral equations of the second kind, which take the form \phi(x) = f(x) + \lambda \int_a^b K(x,y) \phi(y) \, dy, where \lambda is a parameter, K(x,y) is the kernel, and f(x) is a given function.^[25] This equation can be rewritten as \phi = f + T \phi, with the integral operator T = \lambda K defined by (T \phi)(x) = \lambda \int_a^b K(x,y) \phi(y) \, dy.^[25] To solve it, the Neumann series expands \phi as the infinite sum \phi = \sum_{n=0}^\infty T^n f, where T^0 f = f and higher powers involve iterated kernels, such as T^2 f (x) = \lambda^2 \int_a^b \int_a^b K(x,y_1) K(y_1,y_2) f(y_2) \, dy_2 \, dy_1.^[25] The successive approximations \phi_n = \sum_{k=0}^n T^k f converge to the unique solution in appropriate function spaces if the operator norm satisfies \|T\| < 1, often ensured by |\lambda| \max |K(x,y)| (b-a) < 1.^[25] This condition guarantees the series sums to (I - T)^{-1} f.^[26] For Volterra integral equations of the second kind, \phi(x) = f(x) + \lambda \int_a^x K(x,y) \phi(y) \, dy, the Neumann series \phi = \sum_{n=0}^\infty T^n f converges uniformly for any \lambda, without requiring a norm condition on T.^[25] The triangular structure of the Volterra operator, where integrations are over y \leq x, introduces factorial denominators in the iterated terms, bounding the norm by \|T^n\| \leq M^n (x-a)^n / n! for kernel bound M, yielding exponential convergence akin to the series for e^{|\lambda| M (x-a)}.^[25] In both cases, the Neumann series sums to an expression involving the resolvent kernel R(x,y;\lambda) = \sum_{n=1}^\infty \lambda^{n-1} K_n(x,y), where K_n are the iterated kernels, allowing the solution \phi(x) = f(x) + \int_a^b R(x,y;\lambda) f(y) \, dy.^[25] This resolvent integral operator formalizes the inverse of I - T.^[26]

Advanced topics

Openness of invertible operators

In the context of bounded linear operators on a Banach space X, the set \mathrm{GL}(X) of invertible operators in the Banach algebra \mathcal{B}(X) is open with respect to the operator norm topology. This topological property arises from the Neumann series, which provides a perturbation argument demonstrating that small deviations from an invertible operator remain invertible. Specifically, the openness ensures that there exists a neighborhood around any invertible operator consisting entirely of invertible operators, induced by the metric d(S, T) = \|S - T\| for S, T \in \mathcal{B}(X). Consider an invertible operator S \in \mathcal{B}(X) and another operator T \in \mathcal{B}(X) such that \|S^{-1}(S - T)\| < 1. Then T is invertible, with the inverse given by the Neumann series

T^{-1} = \sum_{k=0}^{\infty} \left[S^{-1}(S - T)\right]^k S^{-1},

which converges in the operator norm due to the contraction condition on S^{-1}(S - T). This follows from rewriting T = S(I - S^{-1}(S - T)), where I - S^{-1}(S - T) is invertible via the geometric series expansion \sum_{k=0}^{\infty} [S^{-1}(S - T)]^k, as the spectral radius of S^{-1}(S - T) is at most its norm, which is less than 1. Equivalently, the condition \|S - T\| < \|S^{-1}\|^{-1} guarantees the series converges and yields the explicit inverse.^[27] The proof proceeds by verifying that the partial sums of the series approximate the inverse: the zeroth partial sum is S^{-1}, and higher terms correct for the perturbation, with the remainder bounded by a geometric series tail estimate \|r_n\| \leq r^{n+1}/(1 - r) where r = \|S^{-1}(S - T)\| < 1. This perturbation argument shows that the ball of radius \|S^{-1}\|^{-1} around S lies entirely in \mathrm{GL}(X), confirming openness. Since \|S^{-1}\| is finite for each invertible S, such a ball always exists, and the union over all S \in \mathrm{GL}(X) covers the set. Moreover, the inversion map \mathrm{GL}(X) \to \mathrm{GL}(X), T \mapsto T^{-1}, is continuous on this open set, as the series expansion bounds the difference \|T^{-1} - S^{-1}\| \leq \|S - T\| \|S^{-1}\|^2 / (1 - \|S^{-1}(S - T)\|).^[27] A key corollary is that the spectrum \sigma(S) of any S \in \mathcal{B}(X) is closed in \mathbb{C}. The resolvent set \rho(S) = \{\lambda \in \mathbb{C} : \lambda I - S \text{ is invertible}\} is open, as the same Neumann series argument applies to perturbations around points \lambda \in \rho(S): if |\mu - \lambda| < \|(\lambda I - S)^{-1}\|^{-1}, then \mu I - S is invertible with inverse \sum_{k=0}^{\infty} [(\lambda I - S)^{-1} (\lambda - \mu) I]^k (\lambda I - S)^{-1}. Thus, the complement \sigma(S) is closed, establishing a fundamental property of spectral theory via the openness of invertibles.^[27]

Generalizations

The Neumann series has connections to the theory of strongly continuous semigroups generated by unbounded operators on Banach spaces. For an analytic semigroup \{T(t)\}_{t \geq 0} with generator A, the semigroup can be expressed as T(t) = e^{tA} using the exponential power series \sum_{k=0}^\infty \frac{(tA)^k}{k!}, which converges in the strong operator topology for small t > 0 in a sector of the complex plane. This exponential series provides a way to handle evolution equations like the heat equation, where the generator A is unbounded, extending the scope beyond bounded perturbations of the identity.^[28] In Hilbert spaces, for bounded self-adjoint operators H, the resolvent (I - zH)^{-1} can be expanded using the Neumann series \sum_{n=0}^\infty z^n H^n, which converges in the operator norm if |z| < 1 / \|H\|. In weaker topologies, such as the strong or weak operator topology, convergence can hold under broader spectral conditions, which is useful in the context of von Neumann algebras for spectral projections and functional calculus. Analogous Neumann series arise in p-adic analysis within number theory, where convergence is governed by the non-Archimedean p-adic norm rather than the real absolute value, allowing formal power series expansions in p-adic Banach algebras or operator algebras over \mathbb{Q}_p. In p-adic functional analysis, for a closed operator A in a p-adic Hilbert space analogue, the geometric series \sum_{n=0}^\infty A^n inverts I - A provided the spectral radius condition holds with respect to the p-adic valuation, mirroring the classical case but with ultrametric convergence properties that enable applications to p-adic L-functions and algebraic closures.^[29] These series, often realized via Mal'cev-Neumann formal expansions with well-ordered supports, facilitate the construction of invertible elements in p-adic division rings and support theorems like the p-adic spectral theorem for normal operators.^[29] Developments in quantum mechanics have leveraged perturbation theory using series expansions analogous to the Neumann series for unbounded Hamiltonians and open systems, particularly in extensions to von Neumann entropy perturbations and quantum information measures since 2000. For instance, expansions of the von Neumann entropy S(\rho) = -\operatorname{Tr}(\rho \log \rho) for perturbed density operators \rho = \rho_0 + \epsilon V employ perturbation series for the logarithm, converging in trace norm for small perturbations \epsilon, yielding explicit formulas for entropy shifts in non-degenerate quantum states.^[30] Similarly, in quantum information, perturbation series for functions of linear operators on density matrices provide bounds on fidelity and relative entropy, with applications to error correction in noisy quantum channels where the series converges on finite-dimensional subspaces of infinite-dimensional Hilbert spaces.^[31] These extensions emphasize mathematical rigor in handling degeneracy and unbounded perturbations, distinct from classical bounded cases.^[31]

References

[1]
[PDF] [6pt] @let@token Linear Analysis Lecture 8
Neumann Series. Consider the power series. 1. 1 − z. = ∞. X k=0 zk with radius of convergence 1 at z = 0. Lemma.(C. Neumann). If L ∈ B(V) and kLk < 1, then (I ...
[2]
None
### Definition and Summary of Neumann Series from the Paper
[3]
[PDF] A probabilistic analysis of the Neumann series iteration
Neumann series was introduced by Carl Neumann in 1877 in the context of potential theory [12]. Neumann series, or the more advanced Liouville-Neumann series ...
[4]
[PDF] arXiv:1909.07506v1 [math.PR] 16 Sep 2019
Sep 16, 2019 · Introduction. Neumann series was introduced by Carl Neumann in 1877 in the context of potential theory [NN77]. Neumann series, or the more ...
[5]
Catalog Record: Untersuchungen über das logarithmische und...
Untersuchungen über das logarithmische und Newton'sche Potential, von Dr. C. Neumann. ; Language(s): English ; Published: Leipzig : B. G. Teubner, 1877. ; Subjects ...Missing: Theorie der Potentiale
[6]
Lebesgue's criticism of Carl Neumann's method in potential theory
Jul 9, 2019 · Abstract. In the 1870s, Carl Neumann proposed the so-called method of the arithmetic mean for solving the Dirichlet problem on convex ...
[7]
Carl Neumann (1832 - 1925) - Biography - MacTutor
Carl Neumann was a German mathematician who worked on the Dirichlet principle and on integral equations. He was one of the founders of Mathematische Annalen.
[8]
David Hilbert (1862 - 1943) - Biography - MacTutor
Hilbert's work in integral equations in about 1909 led directly to 20th -century research in functional analysis (the branch of mathematics in which functions ...Missing: adoption | Show results with:adoption
[9]
[PDF] FREDHOLM, HILBERT, SCHMIDT Three Fundamental Papers on ...
Dec 15, 2011 · From this work emerged four general forms of integral equations now called Volterra and Fredholm equations of the first and second kinds (a ...Missing: adoption | Show results with:adoption
[10]
[PDF] Functional Analysis Lecture notes for 18.102, Spring 2020
Here we have used the convergence of the Neumann series, so (Id −B)−1 does exist. ... definition of the coefficients of the Fourier series on an expanding ...
[11]
[PDF] AMATH 731: Applied Functional Analysis Lecture Notes
... Neumann series. So the limit of the iteration sequence is the function to which the Neumann series converges. Note that this result holds regardless of the ...
[12]
[PDF] Functional Analysis Lecture Notes - MSU math
The Neumann series shows that R(ζ) has a convergent power series at each point: R(ζ + h) = ((ζ + h)I − A). −1. = ∞. ∑ n=0 h n. (ζI − Z) n+1 for small ...
[13]
[PDF] Lecture 12: Infinite Series. Operator Functions. Neumann Series.
Proposition 18. A normed linear space X is a Banach space if and only if every absolutely convergent infinite series with terms in X is convergent. Proof ...
[14]
[PDF] Neumann Series and Green's Functions
Jun 24, 2009 · In the above definition of the Green's function, K(x,r,(t − τ)) is the quantum free propagator.
[15]
[PDF] 7. Banach algebras Definition 7.1. A is called a Banach ... - OU Math
Neumann series converges, too (see Exercise 2.22). By the continuity. Page 4. 70. Christian Remling of the multiplication in A,. (e − x). ∞. X n=0 xn = (e − x) ...
[16]
[PDF] Math 261y: von Neumann Algebras (Lecture 2)
Sep 2, 2011 · If ρ(x) ≤ 1, then complex analysis implies that this series converges absolutely |λ| > 1. ... Then the spectral radius ρ(x) coincides with ...<|control11|><|separator|>
[17]
[PDF] arXiv:math/0004049v1 [math.FA] 8 Apr 2000
We show that the Gelfand formula for spectral radius and Neumann series can still be naturally interpreted for operators on topological vector spaces. Of course ...
[18]
Perturbation Theory for Linear Operators - SpringerLink
Free delivery 14-day returnsReviews. "The monograph by T. Kato is an excellent textbook in the theory of linear operators in Banach and Hilbert spaces. · Authors and Affiliations ...
[19]
[PDF] Spectrum (functional analysis)
Mar 12, 2013 · The space of bounded linear operators B(X) on a Banach space X is an example of a unital Banach algebra. Since the definition of the spectrum ...
[20]
[PDF] 1 A Note on Spectral Theory
The Neumann series applied to an element of. R(A) yields a polynomial, but to make sure that the series converges we need to consider the series in a Banach ...
[21]
[PDF] Linear Algebra via Complex Analysis
We compute the Laurent expansion of the resolvent about the eigenvalues of A. Using the Laurent expansion, we prove the Jordan decomposition theorem, prove ...
[22]
[PDF] 1 Matrix calculus - Cornell: Computer Science
Aug 31, 2016 · Numerical linear algebra is not just about algebra, but also about analysis, ... useful example: the matrix version of the geometric series.
[23]
[PDF] Matrix Analysis & Applied Linear Algebra - VIK Wiki
... linear algebra is usually en- countered between basic problem solving ... example, traditional one-semester undergraduate courses can be taught from ...<|control11|><|separator|>
[24]
[PDF] Integral equations - Galileo
Integral equations whose upper limit isx are known as Volterra equations. ... constructed by Neumann series expansion. The preceding equation can then be ...<|separator|>
[25]
[PDF] INTEGRAL EQUATIONS
7.1 Normed Vector Spaces. We will start with some definitions and results from the theory of normed vector ... Neumann series can be convergent. Example 7.1.
[26]
[PDF] 3 Neumann Series and Separable Kernels
The Neumann series technique is usually applied to Fredholm equations with fixed integration limits, though it can be applied to Volterra equations. Example 5: ...
[27]
[PDF] FUNCTIONAL ANALYSIS - ETH Zürich
Jun 8, 2017 · It then intro- duces Banach algebras and shows that the group of invertible elements is an open set. ... invertible operators on a Banach space is.
[28]
PERTURBATION THEORY FOR SEMI-GROUPS OF LINEAR ...
This paper studies the effect of adding a bounded linear operator to the infinitesimal generator of a semi-group of linear bounded operators.
[29]
[PDF] arXiv:2212.07137v3 [math.FA] 1 May 2023
May 1, 2023 · The problem of connecting the operator parameters that label the same self-adjoint extension of a given symmetric operator, respectively, within.
[30]
[PDF] Aspects of p-adic operator algebras - arXiv
Jul 16, 2019 · In this article, we propose a p-adic analogue of complex Hilbert space and consider generalizations of some well-known theorems from functional.
[31]
[0902.4733] Perturbation theory of von Neumann Entropy - arXiv
Feb 27, 2009 · We develop the perturbation theory systematically for calculating von Neumann entropy of non-degenerate systems as well as degenerate systems.Missing: series 2000
[32]
[2106.05533] Perturbation Theory for Quantum Information - arXiv
Jun 10, 2021 · We report lowest-order series expansions for primary matrix functions of quantum states based on a perturbation theory for functions of linear operators.