I am trying to expand the flat-space action
$$ S_{BI} = -T_p \int{d^{p+1}} \sigma \ \mathrm{Tr}\left( e^{-\phi} \sqrt{ -\det(\eta_{ab} + 4\pi^2\alpha^2 \partial_a\Phi^i\partial_b\Phi^i + 2\pi \alpha F_{ab}) \det(Q^{i}_{j}) } \right).\tag{1} $$
After some manipulation, I want to use the power series expansion of the natural log to expand the term
$$ \text{ln}[\delta^{c}_{b} + \lambda^2\eta^{cd} \partial_d\Phi^i\partial_b\Phi^i ]\tag{2} $$ to the fourth power in $\eta^{cd} \partial_d\Phi^i\partial_b\Phi^i.$
My problem is, I don't understand how to compute powers of $\eta^{cd} \partial_d\Phi^i\partial_b\Phi^i $, i.e. $(\eta^{cd} \partial_d\Phi^i\partial_b\Phi^i)^2$. How do I treat the indices?