In M.I. Friedlin's famous paper "On the Factorization of Non-Negative Definite Matrices", he shows that if a non-negative definite symmetric matrix $a(x)=\{a^{ij}(x)\}_{i,j=1}^n$ is in $C^2(\mathbb{R}^n)$, then there exists a symmetric Lipschitz continuous matrix $\sigma(x)$ such that $\sigma^2(x)=a(x)$.
Let $m>0$ be an integer. Are there any reasonable conditions on $a(x)$ that ensure there exists a symmetric matrix $\sigma(x)$ such that $\sigma^2(x)=a(x)$ and $\sigma\in C^m(\mathbb{R}^d)$?