my question is related to evaluating the derivatives of products of tensors using indices. I have gone through the questions already asked on math.stackexchange, (search on tensor derivative yields 902 results, only this one is closer but still unanswered) but after going through them for one week, I couldn't find the answer. So unfortunately here comes another question related to tensor derivatives.
I have a tensor ${\bf F}$ (deformation gradient in solid mechanics) and I have dozens of relations where I need derivatives of tensorial functions
$$\mathbb G(\bf F) = \frac{\partial {\boldsymbol \sigma}(\bf F)}{\partial \bf F}$$
I know that derivative of $\bf F$ with itself is the fourth rank identity tensor $\mathbb I$:
$$
\frac{\partial \bf F}{\partial \bf F} = {\bf I}\overline {\otimes} {\bf I} = \mathbb I
$$
or in indical notation:
$$
\frac{\partial F_{ij}}{\partial F_{pq}} = \delta_{ip} \delta_{jq}
$$
given that components of $\bf F$ are not inter-related through further dependency.
The problem comes in when evaluating product of two tensors, example: $$ \frac{\partial (\bf F^{-1}\cdot \bf F)}{\partial \bf F} = \frac{\partial (F^{-1}_{ij}\cdot F_{jk})}{\partial F_{pq}} {\bf e_i}\otimes{\bf e_k}\otimes{\bf e_p}\otimes{\bf e_q} = {\bf 0} $$ Question: is it correct to write above equation, using chain rule as: $$ \begin{split} \frac{\partial (\bf F^{-1}\cdot \bf F)}{\partial \bf F} &= \frac {\partial\bf F^{-1}}{\partial \bf F} \cdot {\bf F} + {\bf F^{-1}} \cdot \frac {\partial\bf F}{\partial \bf F} \\ &= \left[ \frac{\partial F^{-1}_{ij}}{\partial F_{pq}}\,F_{jk} + F^{-1}_{ij}\,\frac{\partial F_{jk}}{\partial F_{pq}} \right] {\bf e_i}\otimes{\bf e_k}\otimes{\bf e_p}\otimes{\bf e_q} = {\bf 0} \end{split} $$ The terms in square brackets now become: $$ \begin{split} \, & \frac{\partial F^{-1}_{ij}}{\partial F_{pq}}\,F_{jk} + F^{-1}_{ij}\,\frac{\partial F_{jk}}{\partial F_{pq}} = 0 \\ \Longrightarrow & \frac{\partial F^{-1}_{ij}}{\partial F_{pq}} = - F^{-1}_{ij}\,\frac{\partial F_{jk}}{\partial F_{pq}} F^{-1}_{jk} \\ &\quad \qquad = - F^{-1}_{ij}\,\delta_{jp} \delta_{kq} F^{-1}_{jk} \end{split} $$ So the final result is: $$ \frac{\partial F^{-1}_{ij}}{\partial F_{pq}} = - F^{-1}_{ip}\, F^{-1}_{jq} $$ However, in Matrix Cook Book (page 10 eq. 60) (MCB), and this Wikipedia article, the result is: $$ \frac{\partial F^{-1}_{ij}}{\partial F_{pq}} = - F^{-1}_{ip}\, F^{-1}_{qj} $$ The interchanged indices mean that my formulation has transposed the "actual" result (as per MCB and Wikipedia).
What is wrong with the derivation method using indices, which I have used.
Thanks.
Note: please note that my question is not about the derivative of inverse of a tensor, but the method itself, i.e. derivative of an expressions consisting of products of tensors, with respect to a tensor!