Almost commuting matrices

It is shown that if A and B are n × n complex matrices with $\text{A = A}{}^1\text{and}\text{∥AB ? BA∥</}\text{2?}{}^2\text{(n ? 1)}$, then there exist n × n matrices A′ and B′ with $\text{A′ = A′}{}^1\text{such that}\text{A′B′ = B′A′}\text{and}\text{∥A ? A′∥? ?, ∥B ? B′∥? ?}$.     

Asymptotic properties of general symmetric hyperbolic systems

Results on partition of energy and on energy decay are derived for solutions of the Cauchy problem $\text{?u}\text{?t}\text{+ ∑}{}_{\mathrm{j\; =\; 1}}{}^{\mathrm{n}}\text{A}{}_{\mathrm{j}}\text{?u}\text{?x}{}_{\mathrm{j}}\text{= 0, u(0, x) = ?(x)}$. Here the A_j's are constant, k × k Hermitian matrices, x = (x₁,…, x_n), t represents time, and u = u(t, x) is a k-vector. It is shown that the energy of Mu approaches a limit $\text{E}{}_{\mathrm{M}}\text{(?)}\text{as}\text{¦ t ¦ → ∞}$, where M is an arbitrary matrix; that there exists a sufficiently large subspace of data ?, which is invariant under the solution group U₀(t) and such that $\text{U}{}_0\text{(t)? = 0}\text{for}\text{¦ x ¦ ? a ¦ t ¦ ? R, a}\text{and}\text{R}$ depending on ? and that the local energy of nonstatic solutions decays as $\text{¦ t ¦ → ∞}$. More refined results on energy decay are also given and the existence of wave operators is established, considering a perturbed equation $\text{E(x)}\text{?u}\text{?t}\text{+ ∑}{}_{\mathrm{j\; =\; 1}}{}^{\mathrm{n}}\text{A}{}_{\mathrm{j}}\text{?u}\text{?x}{}_{\mathrm{j}}\text{= 0,}\text{where}\text{¦ E(x) ? I ¦ = O(¦ x ¦}{}^{\mathrm{?1\; ?\; ?}}\text{)}$ at infinity.     

Eigenvalue inequalities for products of matrix exponentials

Motivated by models from stochastic population biology and statistical mechanics, we proved new inequalities of the form $\text{(1) ?(e}{}^{\mathrm{A}}\text{e}{}^{\mathrm{B}}\text{)??(e}{}^{\mathrm{A+B}}\text{)}$, where A and B are n × n complex matrices, 1<n<∞, and ? is a real-valued continuous function of the eigenvalues of its matrix argument. For example, if A is essentially nonnegative, B is diagonal real, and ? is the spectral radius, then (1) holds; if in addition A is irreducible and B has at least two different diagonal elements, then the inequality (1) is strict. The proof uses Kingman's theorem on the log-convexity of the spectral radius, Lie's product formula, and perturbation theory. We conclude with conjectures.     

Problèmes de Neumann non linéaires sur les variétés riemanniennes

Nonlinear Neumann problems on riemannian manifolds. Let ($\text{M}\text{, g}$) be a C^∞ compact riemannian manifold of dimension n ? 2 whose boundary B is an (n ? 1)-dimensional submanifold and let $\text{M =}\text{M}\text{?B}$ be the interior of $\text{M}$. Study of Neumann problems of the form: $\text{Δφ +?(φ, x) = 0}\text{in}\text{M, (}\text{dφ}\text{dn}\text{) + g(φ, y) = 0}\text{on}\text{B}$, where, for every $\text{(t, x, y) ? R × M × B, ¦?(t, x)¦}\text{and}\text{¦g(t, y)¦}$ are bounded by $\text{C(1 + ¦t¦}{}^{\mathrm{a}}\text{)}\text{or}\text{C}\text{exp}\text{(¦t¦}{}^{\mathrm{a}}\text{)}$. Application to the determination of a conformal metric for which the scalar curvature of M and the mean curvature of B take prescribed values.     

The asymptotic eigenvalue-distribution for a certain class of random matrices

For an n × n Hermitean matrix A with eigenvalues λ₁, …, λ_n the eigenvalue-distribution is defined by $\text{G(x, A) :=}\text{1}\text{n}$ · number {λ_i: λ_i ? x} for all real x. Let A_n for n = 1, 2, … be an n × n matrix, whose entries a_ik are for i, k = 1, …, n independent complex random variables on a probability space (Ω, $\text{R}$, p) with the same distribution F_a. Suppose that all moments $\text{E}$ | a | ^k, k = 1, 2, … are finite, $\text{E}$a=0 and $\text{E}$ | a | ². Let

$\text{M(A)=}\text{∑}\text{σ=1}\text{s}\text{θ}{}_{\mathrm{\sigma }}\text{P}{}_{\mathrm{\sigma }}\text{(A,A}{}^1\text{)}$

with complex numbers θ_σ and finite products P_σ of factors A and $\text{A}{}^1$ (= Hermitean conjugate) be a function which assigns to each matrix A an Hermitean matrix M(A). The following limit theorem is proved: There exists a distribution function G₀(x) = G₁x) + G₂(x), where G₁ is a step function and G₂ is absolutely continuous, such that with probability $\text{1 G(x, M(}\text{A}{}_{\mathrm{n}}\text{n}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{))}$ converges to G₀(x) as n → ∞ for all continuity points x of G₀. The density g of G₂ vanishes outside a finite interval. There are only finitely many jumps of G₁. Both, G₁ and G₂, can explicitly be expressed by means of a certain algebraic function f, which is determined by equations, which can easily be derived from the special form of M(A). This result is analogous to Wigner's semicircle theorem for symmetric random matrices (E. P. Wigner, Random matrices in physics, SIAM Review9 (1967), 1–23). The examples $\text{A}{}_{\mathrm{r}}\text{A}{}^{\mathrm{1r}}$, $\text{A}{}^{\mathrm{r}}\text{+ A}{}^{\mathrm{1r}}$, $\text{A}{}^{\mathrm{r}}\text{A}{}^{\mathrm{1r}}\text{± A}{}^{\mathrm{1r}}\text{A}{}^{\mathrm{r}}$, r = 1, 2, …, are discussed in more detail. Some inequalities for random matrices are derived. It turns out that with probability 1 the sharpened form

$\text{lim sup}\text{n→∞}\text{∑}\text{i=1}{}^{\mathrm{n}}\text{|γ}{}_{\mathrm{i}}{}^{\mathrm{(n)}}\text{|}{}^2\text{?6A}{}_{\mathrm{n}}\text{6}{}^2\text{? 0.8228…}$

of Schur's inequality for the eigenvalues λ_i⁽ⁿ⁾ of A_n holds. Consequently random matrices do not tend to be normal matrices for large n.     

Arnold's canonical matrices and the asymptotic simplification of ordinary differential equations

Let A(x,ε) be an n×n matrix function holomorphic for |x|?x₀, 0<ε?ε₀, and possessing, uniformly in x, an asymptotic expansion $\text{A(x,ε)?Σ}{}^{\mathrm{\infty }}{}_{\mathrm{r=0}}\text{A}{}_{\mathrm{r}}\text{(x) ε}{}^{\mathrm{r}}$, as ε→0⁺. An invertible, holomorphic matrix function P(x,ε) with an asymptotic expansion $\text{P(x,ε)?Σ}{}^{\mathrm{\infty }}{}_{\mathrm{r=0}}\text{P}{}_{\mathrm{r}}\text{(x)ε}{}^{\mathrm{r}}$, as ε→0⁺, is constructed, such that the transformation y = P(x,ε)z takes the differential equation $\text{ε}{}^{\mathrm{h}}\text{dy}\text{dx}\text{= A(x,ε)y,h}$ a positive integer, into $\text{ε}{}^{\mathrm{h}}\text{dz}\text{dx}\text{= B(x,ε)z}$, where B(x,ε) is asymptotically equal, to all orders, to a matrix in a canonical form for holomorphic matrices due to V.I. Arnold.     

On Lp estimates for the wave equation

New and more elementary proofs are given of two results due to W. Littman: (1) Let $\text{n ? 2, p ?}\text{2n}\text{(n ? 1)}$. The estimate $\text{∫∫ (¦▽u¦}{}^{\mathrm{p}}\text{+ ¦u}{}_{\mathrm{t}}\text{¦}{}^{\mathrm{p}}\text{) dx dt ? C ∫∫ ¦□u¦}{}^{\mathrm{p}}\text{dx dt}$ cannot hold for all u?C₀^∞(Q), Q a cube in $\text{R}{}^{\mathrm{<mi\; mathvariant="italic"\; is="true">n</mi>}}\text{× R}$, some constant C. (2) Let n ? 2, p ≠ 2. The estimate $\text{∫ (¦▽(t)¦}{}^{\mathrm{p}}\text{+ ¦u}{}_{\mathrm{t}}\text{(t)¦}{}^{\mathrm{p}}\text{) dx ? C(t) ∫ (¦▽u(0)¦}{}^{\mathrm{p}}\text{+ ¦u}{}_{\mathrm{t}}\text{(0)¦}{}^{\mathrm{p}}\text{) dx}$ cannot hold for all C^∞ solutions of the wave equation □u = 0 in $\text{R}{}^{\mathrm{n}}\text{x R}$; all t ?$\text{R}$; some function C: $\text{R}$ → $\text{R}$.     

A generalization of the fast LUP matrix decomposition algorithm and applications

We show that any m × n matrix A, over any field, can be written as a product, LSP, of three matrices, where L is a lower triangular matrix with l's on the main diagonal, S is an m × n matrix which reduces to an upper triangular matrix with nonzero diagonal elements when the zero rows are deleted, and P is an n × n permutation matrix. Moreover, L, S, and P can be found in O(m^α?1n) time, where the complexity of matrix multiplication is O(m^α). We use the LSP decomposition to construct fast algorithms for some important matrix problems. In particular, we develop O(m^α?1n) algorithms for the following problems, where A is any m × n matrix: (1) Determine if the system of equations $\text{A}\text{x}\text{=}\text{b}$ (where $\text{b}$ is a column vector) has a solution, and if so, find one such solution. (2) Find a generalized inverse, $\text{A}{}^1$, of A (i.e., $\text{AA}{}^1\text{A = A}$). (3) Find simultaneously a maximal independent set of rows and a maximal independent set of columns of A.     

Maximum zero strings of Bell numbers modulo primes

For any prime p, the sequence of Bell exponential numbers B_n is shown to have p ? 1 consecutive values congruent to zero (mod p), beginning with B_m, where $\text{m ≡ 1 ?}\text{(p}{}^{\mathrm{p}}\text{? 1)}\text{(p ? 1)}{}^2$ ($\text{mod}\text{(p}{}^{\mathrm{p}}\text{? 1)}\text{(p ? 1)}$). This is an improvement over previous results on the maximal strings of zero residues of the Bell numbers. Similar results are obtained for the sequence of generalized Bell numbers A_n generated by .     

On the convergence of the Neumann series in interval analysis

Let $\text{A}$ be a real or complex n × n interval matrix. Then it is shown that the Neumann series $\text{Σ}{}^{\mathrm{\infty }}{}_{\mathrm{k=0}}\text{A}{}^{\mathrm{k}}$ is convergent iff the sequence {$\text{A}$^k} converges to the null matrix $\text{O}$, i.e., iff the spectral radius of the real comparison matrix $\text{B}$ constructed in [2] is less than one.     

Continuity of positive functionals on topological 1-algebras

It is shown that the set $\text{C}$_{m × n} of complex m × n matrices forms a lower semilattice under the partial ordering A ? B defined by $\text{A}{}^1\text{A = A}{}^1\text{B,}{}^1\text{AA}{}^1\text{= BA}{}^1\text{,}\text{where}\text{A}{}^1$ denotes the conjugate transpose of A. As a special case of a result for division rings, it is further shown that, over any field F, form = n = 2 and any proper involution 1 of F_{2 × 2}, the corresponding intersections A ∩ B all exist.     

The resolvent of a singularly perturbed elliptic operator

Let A and B be uniformly elliptic operators of orders 2m and 2n, respectively, m > n. We consider the Dirichlet problems for the equations (?^{2(m ? n)}A + B + λ²ⁿI)u_? = f and (B + λ²ⁿI)u = f in a bounded domain Ω in R^k with a smooth boundary ?Ω. The estimate is derived. This result extends the results of [7, 9, 10, 12, 14, 15, 18]by giving estimates up to the boundary, improving the rate of convergence in ?, using lower norms, and considering operators of higher order with variable coefficients. An application to a parabolic boundary value problem is given.     

Some existence and nonexistence theorems for solutions of degenerate parabolic equations

The initial and boundary value problem for the degenerate parabolic equation v_t = Δ(?(v)) + F(v) in the cylinder $\text{Ω × ¦0, ∞), Ω ? R}{}^{\mathrm{n}}$ bounded, for a certain class of point functions ? satisfying ?′(v) ? 0 (e.g., $\text{?(v) = ¦v¦}{}^{\mathrm{m}}\text{sign}\text{v}$) is considered. In the case that F(v) sign $\text{v ? C(1 + ¦?(v)¦}{}^{\mathrm{\alpha }}\text{), α < 1}$, the equation has a global time solution. The same is true for α = 1 provided the measure of Ω is sufficiently small. In the case that $\text{F(v)}\text{?(v)}$ is nondecreasing a condition is given on the initial state v(x, 0) which implies that the solution must blow up in finite time. The existence of such initial states is discussed.     

On Lp-spectrum of elliptic operators

For elliptic operators $\text{A = ∑}{}_{\mathrm{¦\alpha ¦\; ?\; m}}\text{a}{}_{\mathrm{\alpha }}\text{(x) D}{}^{\mathrm{\alpha }}$ on Rⁿ and certain of their singular perturbations $\text{B = ∑}{}_{\mathrm{¦\alpha ¦\; ?\; m}}\text{b}{}_{\mathrm{\alpha }}\text{(x)D}{}^{\mathrm{\alpha }}$ relative compactness of B with respect to A is established. This result applies to the study of L^p-spectra of elliptic operators for different p.     

Geometry of the numerical range of matrices

Some techniques for the study of the algebraic curve C(A) which generates the numerical range W(A) of an n×n matrix A as its convex hull are developed. These enable one to give an explicit point equation of C(A) and a formula for the curvature of C(A) at a boundary point of W(A). Applied to the case of a nonnegative matrix A, a simple relation is found between the curvature of the function Φ(A)=p((1?α)A+ αA^T) (pbeingthePerronroot) at $\text{α=}\text{1}\text{2}$ and the curvature of W(A) at the Perron root of $\text{1}\text{2}\text{(A+A}{}^{\mathrm{T}}\text{)}$. A connection with 2-dimensional pencils of Hermitian matrices is mentioned and a conjecture formulated.     

Oscillation theory at a finite singularity

It is well known that every weak solution (with boundary values 0) of a semilinear equation $\text{Au + ?(x, u) = g}$ is a regular solution if ? fulfils the growth condition $\text{(1) ¦?(x, u)¦? c ¦u¦}{}^{\mathrm{<mtext>(n\; +\; 2m)</mtext><mtext>(n\; ?\; 2m)</mtext><mtext>\; ?\; ?</mtext>}}$. Here 2m is the order of A. In this paper we weaken this condition to $\text{c ¦u ¦}{}^{\mathrm{<mtext>(n\; +\; 2m)</mtext><mtext>(n\; ?\; 2m)</mtext><mtext>\; +\; 1</mtext>}}\text{? ?(x, u)u ? ?c ¦u ?}{}^{\mathrm{<mtext>(n\; +\; 2m)</mtext><mtext>(n\; ?\; 2m)</mtext><mtext>\; +\; 1\; ?\; ?</mtext>}}$. This requires a technique completely different from that which may be applied in case (1).     

Lp estimates for the wave equation

Let u(x, t) be the solution of u_tt ? Δ_xu = 0 with initial conditions $\text{u(x, 0) = g(x)}\text{and}\text{u}{}_{\mathrm{t}}\text{(x, 0) = ?;(x)}$. Consider the linear operator $\text{T: ?; → u(x, t)}$. (Here g = 0.) We prove for t fixed the following result. Theorem 1: T is bounded in L^p if and only if . Theorem 2: If the coefficients are variables in C and constant outside of some compact set we get: (a) If n = 2k the result holds for $\text{¦ p}{}^{\mathrm{?1}}\text{? 2}{}^{\mathrm{?1}}\text{¦ < (n ? 1)}{}^{\mathrm{?1}}$. (b) If n = 2k ? 1, the result is valid for $\text{¦ p}{}^{\mathrm{?1}}\text{? 2}{}^{\mathrm{?1}}\text{¦ ? (n ? 1)}$. This result are sharp in the sense that for p such that $\text{¦ p}{}^{\mathrm{?1}}\text{? 2}{}^{\mathrm{?1}}\text{¦ > (n ? 1)}{}^{\mathrm{?1}}$ we prove the existence of $\text{?; ? L}{}^{\mathrm{P}}$ in such a way that $\text{T?; ? L}{}^{\mathrm{P}}$. Several applications are given, one of them is to the study of the Klein-Gordon equation, the other to the completion of the study of the family of multipliers $\text{m(ξ) = ψ(ξ) e}{}^{\mathrm{i¦\xi ¦}}\text{¦ ξ ¦}{}^{\mathrm{?b}}$ and finally we get that the convolution against the kernel $\text{K(x) = ?(x)(1 ? ¦ x ¦)}{}^{\mathrm{?1}}$ is bounded in H¹.     

Nonnegative factorization of completely positive matrices

Let A be a real symmetric n × n matrix of rank k, and suppose that A = BB′ for some real n × m matrix B with nonnegative entries (for some m). (Such an A is called completely positive.) It is shown that such a B exists with $\text{m?}\text{1}\text{2}\text{k(k+1)?N}$, where 2N is the maximal number of (off-diagonal) entries which equal zero in a nonsingular principal submatrix of A. An example is given where the least m which works is $\text{(k+1)}{}^2\text{4}$ (k odd),$\text{k(k+2)}\text{4}$ (k even).     

A note on polynomial matrix functions over a finite field

Let F_n denote the ring of n×n matrices over the finite field F=GF(q) and let A(x)=A_Nx^N+ ?+ A₁x+A₀?F_n[x]. A function $\text{?:F}{}_{\mathrm{n}}\text{→F}{}_{\mathrm{n}}$ is called a right polynomial function iff there exists an A(x)?F_n[x] such that $\text{?(B)=A}{}_{\mathrm{N}}\text{B}{}^{\mathrm{N}}\text{+?+A}{}_1\text{B+ A}{}_0$ for every B?F_n. This paper obtains unique representations for and determines the number of right polynomial functions.