I initially asked on math.stackoverflow but have since come to understand this forum may be more appropriate, as this is indeed a question that arose in writing a research article.
Denote $g$ a symmetric positive definite matrix field (defined over, say, $\mathbb R^n$), and $g_{ij}$ its components in some basis (the same for all $g(x), x \in \mathbb R^n$). So $M$ is a smooth function that, to $x\in \mathbb R^n$, associates a positive definite matrix $g(x)$.
I am carrying out some computations, and they turn out to be valid only if:
$\sum_{s=1}^n g_{si}\partial_s g_{jk} - g_{sk}\partial_s g_{ji} = 0$
where $\partial_t$ denotes differentiation wrt to the $t$-th variable (components in the same basis).
I first thought this condition to be unreasonable. I now have the suspicion this may be some kind of metric compatibility condition, but I don't manage to make sense of the results I find on this topic. Understanding this constraint better could help me know if it is realistic to expect such an $g$ be produced in my circumstances, and perhaps how.
Another lead may be to look at this as a PDE. This resembles the divergence of some field, but I am unfortunately not an expert of this either. Perhaps this corresponds to a well-known PDE.
Per suggestion, for $n=2$. For $j=1$:
$g_{11} \partial_1 g_{11} - g_{21} \partial_2 g_{11} = 0$
$g_{11} \partial_1 g_{12} - g_{22} \partial_2 g_{11} = 0$
$g_{12} \partial_1 g_{12} - g_{22} \partial_2 g_{12} = 0$
For $j=2$:
$g_{11} \partial_1 g_{21} - g_{21} \partial_2 g_{21} = 0$
$g_{11} \partial_1 g_{22} - g_{22} \partial_2 g_{21} = 0$
$g_{12} \partial_1 g_{22} - g_{22} \partial_2 g_{22} = 0$
Despite that there are fewer terms than possible combinations, the following repeat:
$g_{11} \partial_1 g_{12}$ in eqs. (2) and (4).
$g_{22} \partial_2 g_{12}$ in eqs. (3) and (5).
I thought I might add some context. Originally, I am interested in a smooth function $F : \mathbb R^n \to \mathbb R^n$. If I want to impose $J_F(x) = g(F(x))M$ ($M$ a matrix allowing this to hold with $g$ spd) then, for all $i,k$,
$\partial_i F_k = \sum_s g_{ks}(F) M_{si}$
Now if we differentiate:
$\partial_{ij} F_k = \sum_{st} \partial_j F_t \partial_t g_{ks}(F) M_{si}$
$\partial_j F_t$ can be substituted using the previous relationship:
$\partial_{ij} F_k = \sum_{stu} g_{tu}(F) M_{uj} \partial_t g_{ks}(F) M_{si}$
rearranging the sums:
$\partial_{ij} F_k = \sum_{su} M_{uj}M_{si} \sum_t g_{tu}(F) \partial_t g_{ks}(F) $
Fixing $k$, the left-hand side is symmetric wrt $(i,j)$ thus so is the right-hand side. Thus:
$\sum_{su} M_{uj}M_{si} \sum_t g_{tu}(F) \partial_t g_{ks}(F) = \sum_{su} M_{uj}M_{si} \sum_t g_{ts}(F) \partial_t g_{ku}(F)$
i.e.:
$\sum_{su} M_{uj}M_{si} \sum_t (g_{tu}(F) \partial_t g_{ks}(F) - g_{ts}(F) \partial_t g_{ku}(F)) = 0$
In my application case, I can choose $M$ "arbitrary enough" (but perhaps there is some pattern I didn't see that makes some of these relations trivial). Thus I need that, for all $u,s$,
$\sum_t (g_{tu} \partial_t g_{ks} - g_{ts} \partial_t g_{ku})$
Maybe my trouble is I am writing $J_F(x) = g(F(x))M$ in a sense too restricted, and a more general framework in which $\partial_{ij} F \ne \partial_{ji} F$ would be more suited.
Also, following a slightly different path (imposing a different condition than $J_F = ...$), this condition does not appear as for all $x$, but is only necessary at a given $x_0$.
Update 3. If we take the condition punctually at some $x_0$$x$. Posing $G : x \mapsto g(x_0) x$$G_{x} : y \mapsto g(x) y$, the Jacobian is $g(x_0)$$g(x)$, then
$\partial_i (g_{jk} \circ G) = \sum_s g_{si} \partial_s g_{jk}(G)$$\partial_i (g_{jk} \circ G_x) = \sum_s g_{si} \partial_s g_{jk}(G_x)$
The condition then writes, in these terms, $\partial_i(g_{jk} \circ G) = \partial_k(g_{ji} \circ G)$$\partial_i(g_{jk} \circ G) = \partial_k(g_{ji} \circ G_x)$. Moreover as $g$ is symmetric, all index permutations are allowed.
In summary, if we consider the condition only at some fixed $x_0$, this amounts to saying the $3$-tensor $\nabla (g\circ G)$$\nabla (g\circ G_x)$ is symmetric.
I am still not clear on the significance of this (g eating itself has a symmetric "Jacobian tensor") but this seems like a step forward.