Let
$$ \hat{\mathbf{x}}_1^2=1\\ \hat{\mathbf{x}}_2^2=1\\ \hat{\mathbf{x}}_3^2=1\\ \hat{\mathbf{x}}_1\hat{\mathbf{x}}_2+\hat{\mathbf{x}}_2\hat{\mathbf{x}}_1=0\\ \hat{\mathbf{x}}_1\hat{\mathbf{x}}_3+\hat{\mathbf{x}}_3\hat{\mathbf{x}}_1=0\\ \hat{\mathbf{x}}_2\hat{\mathbf{x}}_3+\hat{\mathbf{x}}_3\hat{\mathbf{x}}_2=0 $$
Then $$ \mathbf{e}_1\wedge \mathbf{e}_2 \wedge \mathbf{e}_3=\sqrt{|\det g|}\hat{\mathbf{x}}_1\wedge \hat{\mathbf{x}}_2\wedge \hat{\mathbf{x}}_3 $$
for reference see my previous question (How can I show $\mathbf{e}_0\mathbf{e}_1\mathbf{e}_2\mathbf{e}_3=\sqrt{|g|}\gamma_0\gamma_1\gamma_2\gamma_3$)
My question is what is the following expression in terms of $g$?
$$ (\mathbf{e}_1\wedge \mathbf{e}_2+\mathbf{e}_1\wedge \mathbf{e}_3)^2 $$
Let
$$ g=\pmatrix{ g_{11}&g_{12}&g_{13}\\ g_{21}&g_{22}&g_{23}\\ g_{31}&g_{32}&g_{33} } $$
This is what I got so far:
$$ \begin{align} (\mathbf{e}_1\wedge \mathbf{e}_2+\mathbf{e}_1\wedge \mathbf{e}_3)^2&=(\mathbf{e}_1\wedge \mathbf{e}_2+\mathbf{e}_1\wedge \mathbf{e}_3)(\mathbf{e}_1\wedge \mathbf{e}_2+\mathbf{e}_1\wedge \mathbf{e}_3)\\ &=(\mathbf{e}_1\wedge \mathbf{e}_2)^2+(\mathbf{e}_1\wedge \mathbf{e}_2)(\mathbf{e}_1\wedge \mathbf{e}_3)+(\mathbf{e}_1\wedge \mathbf{e}_3)(\mathbf{e}_1\wedge \mathbf{e}_2)+(\mathbf{e}_1\wedge \mathbf{e}_3)^2 \end{align} $$
Let us treat each term individually:
- The part $(\mathbf{e}_1\wedge \mathbf{e}_2)^2$ is:
$$ (\mathbf{e}_1\wedge \mathbf{e}_2)^2=((g_{11}g_{22}-g_{12}g_{21})\hat{\mathbf{x}}_1 \wedge \hat{\mathbf{x}}_2)^2=-(g_{11}g_{22}-g_{12}g_{21})^2 $$
- The part $(\mathbf{e}_1\wedge \mathbf{e}_3)^2$ is:
$$ (\mathbf{e}_1\wedge \mathbf{e}_3)^2=((g_{11}g_{33}-g_{13}g_{31})\hat{\mathbf{x}}_1 \wedge \hat{\mathbf{x}}_3)^2=-(g_{11}g_{33}-g_{13}g_{31})^2 $$
- The part $(\mathbf{e}_1\wedge \mathbf{e}_2)(\mathbf{e}_1\wedge \mathbf{e}_3)$ is:
$$ \begin{align} (\mathbf{e}_1\wedge \mathbf{e}_2)(\mathbf{e}_1\wedge \mathbf{e}_3)&=((g_{11}g_{22}-g_{12}g_{21}) \hat{\mathbf{x}}_1 \wedge \hat{\mathbf{x}}_2)((g_{11}g_{33}-g_{13}g_{31}) \hat{\mathbf{x}}_1\wedge \hat{\mathbf{x}}_3)\\ &=-(g_{11}g_{22}-g_{12}g_{21})(g_{11}g_{33}-g_{13}g_{31}) (\hat{\mathbf{x}}_2 \wedge \hat{\mathbf{x}}_1)( \hat{\mathbf{x}}_1\wedge \hat{\mathbf{x}}_3)\\ &=-(g_{11}g_{22}-g_{12}g_{21})(g_{11}g_{33}-g_{13}g_{31}) ( \hat{\mathbf{x}}_2\wedge \hat{\mathbf{x}}_3) \end{align} $$
- The part $(\mathbf{e}_1\wedge \mathbf{e}_3)(\mathbf{e}_1\wedge \mathbf{e}_2)$ is:
$$ \begin{align} (\mathbf{e}_1\wedge \mathbf{e}_3)(\mathbf{e}_1\wedge \mathbf{e}_2)&=((g_{11}g_{33}-g_{13}g_{31}) \hat{\mathbf{x}}_1 \wedge \hat{\mathbf{x}}_3)((g_{11}g_{22}-g_{12}g_{21}) \hat{\mathbf{x}}_1\wedge \hat{\mathbf{x}}_2)\\ &=-(g_{11}g_{33}-g_{13}g_{31})(g_{11}g_{22}-g_{12}g_{21}) (\hat{\mathbf{x}}_3 \wedge \hat{\mathbf{x}}_1)( \hat{\mathbf{x}}_1\wedge \hat{\mathbf{x}}_2)\\ &=-(g_{11}g_{33}-g_{13}g_{31})(g_{11}g_{22}-g_{12}g_{21}) ( \hat{\mathbf{x}}_3\wedge \hat{\mathbf{x}}_2)\\ &=(g_{11}g_{33}-g_{13}g_{31})(g_{11}g_{22}-g_{12}g_{21}) ( \hat{\mathbf{x}}_2\wedge \hat{\mathbf{x}}_3) \end{align} $$
So the cross-terms cancel, and we get a sum of areas:
$$ (\mathbf{e}_1\wedge \mathbf{e}_2+\mathbf{e}_1\wedge \mathbf{e}_3)^2=-(g_{11}g_{22}-g_{12}g_{21})^2-(g_{11}g_{33}-g_{13}g_{31})^2 $$
What is the geometric interpretation of what I have done? Do the cross-terms really cancel for areas even in the midst of arbitrary curved space -- this result so so surprising to me.