5
$\begingroup$

I asked this question some time ago in MSE but I didn't recieved any feedback.

https://math.stackexchange.com/questions/4672664/diagonalization-of-symmetric-matrices-of-functions

This problem arised to me when I was trying to find an analog to orthogonal reference frames for singular metric tensors.

Let $U\subseteq\mathbb{R}^m$ be an open subset such that $0\in U$ and $n\leq m$. Let $G(p)=\begin{pmatrix} g_{11}(p) & g_{12}(p) & \cdots & g_{1n}(p)\\ g_{12}(p) & g_{22}(p) & \cdots & g_{2n}(p)\\ \vdots & \vdots & \ddots & \vdots \\ g_{1n}(p) & g_{2n}(p) & \cdots & g_{nn}(p)\\ \end{pmatrix} $ be a symmetric matrix of infinitely differentiable functions $g_{ij}:U\rightarrow\mathbb{R}$ such that $G(0)=0\in\mathcal{M}_{n\times n}(\mathbb{R})$, that is $g_{ij}(0)=0$ for all $i,j$.

To diagonalize the matrix near $0\in\mathbb{R}^m$, the Gram-Schmidt process doesn't work because it will divide by $0\in\mathbb{R}$ in some steps. You can not even do an adaptation of this process beacause there is no initial orthonormal basis to start with or more acurately any basis will be orthogonal but not normal since the "norms" will be zero. I write norms between quotation marks because $G(p)$ isn't positive semidefinite so it doesn't define any norm and I don't want to put any condition about the "definitness" of the matrix.

Said that, the question is: is it possible to find an open subset $V\subseteq U$ such that $0\in V$ where such a matrix is diagonalizable? Maybe with a different process than Gram-Schmidt or even leaving the ring of infinitely differentiable functions. I just need to prove the existence theoretically, in other words, I don't need it to be computable.

Thank you in advance!

$\endgroup$
3
  • $\begingroup$ Your matrix is symmetric, so is diagonalizable everywhere. (A different question that is more interesting was answered by @RobertBryant below.) $\endgroup$ May 16, 2023 at 16:50
  • 1
    $\begingroup$ @ChristianRemling no it's not. That's the case for matrices over a field, but not over a general ring, and that's what I'm asking. $\endgroup$ May 16, 2023 at 17:38
  • $\begingroup$ This old answer by Peter Michor links to a survey of results that could help if you relax the premise of your question. Namely, you could potentially find a continuous mapping $V\to U$, which "blows up" the points where $G(p)$ is singular such that the diagonalization of $G$ is continuous over $V$. $\endgroup$ May 18, 2023 at 11:27

1 Answer 1

8
$\begingroup$

In general, this cannot be done. For example, in dimension $2$ in coordinates $(x,y)$, let $$ G(x,y) = \left[\begin{matrix}x&y\\y&-x\end{matrix}\right]. $$ If $G$ could be diagonalized by a differentiable invertible matrix $A(x,y)$, i.e., if $$ A^T G A = \left[\begin{matrix}\lambda_1&0\\ 0&\lambda_2\end{matrix}\right] $$ where $\lambda_1$ and $\lambda_2$ were differentiable, then the $\lambda_i$ would have to vanish at $x=y=0$. Taking determinants yields $$ -(x^2+y^2)(\det A)^2 = \lambda_1\lambda_2\,. $$ Then, looking at the lowest order terms on each side (the terms of order $2$), you'd have $x^2+y^2$ written as a product of two factors linear in $x$ and $y$, which is impossible.

For similar reasons, you cannot achieve $$ G = A^T\left[\begin{matrix}\lambda_1&0\\ 0&\lambda_2\end{matrix}\right]A $$ for a differentiable $A$ and $\lambda_i$. The above argument shows that $A$ could not be invertible, so we would have to have $\det A$ vanishing at $x=y=0$. Then $-(x^2+y^2) = (\det A)^2\lambda_1\lambda_2$ would imply that $\det A$ vanishes at most to order 1 at $x=y=0$ and that $\lambda_1$ and $\lambda_2$ do not vanish at $x=y=0$, which again gives a contradiction, since $x^2+y^2$ is not the square of a linear term.

In fact, one cannot have $$ G = A^T\left[\begin{matrix}\lambda_1&0\\ 0&\lambda_2\end{matrix}\right]A $$ with $A$ and $\lambda_i$ being merely continuous on some disk $ x^2+y^2\le \epsilon^2$ for some $\epsilon>0$.

Here is why: The relation $-(x^2+y^2) = (\det A)^2\lambda_1\lambda_2$, shows that $\det A$, $\lambda_1$ and $\lambda_2$ must be nonzero away from $(x,y)=(0,0)$, so each $\lambda_i$ cannot change sign and we must have $\lambda_1\lambda_2<0$ away from $(x,y)=(0,0)$. Without loss of generality, we can assume that $\lambda_1 = \mu_1^2$ and $\lambda_2=-\mu_2^2$, where the $\mu_i$ are continuous, so by modifying $A$ in the obvious way, we can reduce to the case that $\lambda_1 = -\lambda_2 = 1$ and, moreover, that $\det A = \sqrt{x^2+y^2}>0$.

The mapping $f:S^1\to\mathrm{SL}(2,\mathbb{R})$ defined by $$ f(\theta) = \frac1{\sqrt{\epsilon}} A(\epsilon\cos\theta,\epsilon\sin\theta) $$ then satisfies $$ f(\theta)^T\left[\begin{matrix}1&0\\ 0&-1\end{matrix}\right]f(\theta) = \left[\begin{matrix}\cos\theta&\sin\theta\\ \sin\theta&-\cos\theta\end{matrix}\right]. $$ Now, consider the map $s:\mathrm{SL}(2,\mathbb{R})\to H$ (where $H$, a hyperboloid of one sheet, is the quadric surface in the symmetric $2$-by-$2$ matrices defined by setting the determinant equal to $-1$) defined by $$ s(A) = A^T\left[\begin{matrix}1&0\\ 0&-1\end{matrix}\right]A. $$ Both $\mathrm{SL}(2,\mathbb{R})$ and $H$ are homotopic to the circle and hence have $\pi_1\simeq \mathbb{Z}$. The map $s$ carries the generator of $\pi_1(\mathrm{SL}(2,\mathbb{R}))$ to twice a generator of $\pi_1(H)$. However, the above formula for $f$ shows that $s\circ f:S^1\to H$ carries a generator of $\pi_1(S^1)$ to a generator of $\pi_1(H)$, which is impossible.

$\endgroup$
2
  • $\begingroup$ Thank you Robert. I'm surprised that this can not be done even continuously. Do you know if it has been studied under which conditions this can be done continuously? $\endgroup$ May 17, 2023 at 5:29
  • $\begingroup$ @user1234567890: In the case that $\det G$ vanishes at an isolated point, this can be reduced to a topological question, because then the $(p,q)$-type of $G$ away from the isolated singularity must be constant, and you are essentially asking when a map from a punctured $m$-disk (which is homotopic to $S^{m-1}$) into the homogeneous space $\mathrm{GL}(n,\mathbb{R})/\mathrm{O}(p,q)$ can be lifted continuously to $\mathrm{GL}(n,\mathbb{R})$. $\endgroup$ May 17, 2023 at 9:04

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.