Interesting relationships between Cholesky decomposition and diagonalization

First, let me rephrase your remark. Let $L=HU$ be the polar factorization of $L$ ($H$ hermitian positive definite, $U$ unitary). Then $\Sigma=LL^\ast=H^2$ tells you that the Hermitian part of $L$ is $\sqrt\Sigma$. Then $U=L\Sigma^{-1/2}$ is its unitary part.

On the other hand, you have $\sqrt\Sigma=LQ^\ast=QL^\ast$. This is exactly the $QR$-factorization of $\sqrt\Sigma$: $L^\ast$ is the upper triangular part with positive real diagonal in the $QR$-factorization of $\sqrt\Sigma$.


Computing a Cholesky factor is a much easier task than computing the eigendecomposition. For instance, you can do Cholesky while staying in the original field (modulo some square roots on the diagonal, which can be avoided if you go for the $LDL^*$ form), while arbitrarily complicated algebraic field extensions can be required for the eigendecomposition. So there is no hope to obtain $(P,\Lambda)$ from $L$ through simple algebraic formulas.

The other direction is typically much less useful, at least from a computational point of view --- on the top of my head, everything you can do with Cholesky, you can do it also with the eigenvalue decomposition, and it's more stable.