[guided]We now compute $d\beta_F$ for $\beta_F = F_1\, dx_2 \wedge dx_3 + F_2\, dx_3 \wedge dx_1 + F_3\, dx_1 \wedge dx_2$. Since $d$ annihilates the constant-coefficient $2$-forms $dx_i \wedge dx_j$, the Leibniz rule gives
\begin{align*}
d\beta_F = dF_1 \wedge dx_2 \wedge dx_3 + dF_2 \wedge dx_3 \wedge dx_1 + dF_3 \wedge dx_1 \wedge dx_2.
\end{align*}
For each $i$, write $dF_i = \sum_j \partial_{x_j} F_i\, dx_j$. The key observation is that in each term $dF_i \wedge (\text{the }2\text{-form not involving } dx_i)$, all summands $\partial_{x_j} F_i\, dx_j \wedge \cdots$ with $j \ne i$ contain a repeated factor and therefore vanish. Only the summand $j = i$ survives.
For the first term, the $2$-form is $dx_2 \wedge dx_3$, so only $j = 1$ survives:
\begin{align*}
dF_1 \wedge dx_2 \wedge dx_3 = \partial_{x_1} F_1\, dx_1 \wedge dx_2 \wedge dx_3.
\end{align*}
For the second term, the $2$-form is $dx_3 \wedge dx_1$, so only $j = 2$ survives:
\begin{align*}
dF_2 \wedge dx_3 \wedge dx_1 = \partial_{x_2} F_2\, dx_2 \wedge dx_3 \wedge dx_1.
\end{align*}
We want this in the standard order $dx_1 \wedge dx_2 \wedge dx_3$. The cyclic permutation $(2,3,1) \to (1,2,3)$ is an even permutation (it is the square of a transposition), so $dx_2 \wedge dx_3 \wedge dx_1 = dx_1 \wedge dx_2 \wedge dx_3$.
Similarly for the third term, with $(3,1,2) \to (1,2,3)$ also even:
\begin{align*}
dF_3 \wedge dx_1 \wedge dx_2 = \partial_{x_3} F_3\, dx_3 \wedge dx_1 \wedge dx_2 = \partial_{x_3} F_3\, dx_1 \wedge dx_2 \wedge dx_3.
\end{align*}
This is exactly why the $2$-form basis was chosen to be the cyclically ordered $(dx_2 \wedge dx_3,\, dx_3 \wedge dx_1,\, dx_1 \wedge dx_2)$ rather than the lexicographically ordered $(dx_1 \wedge dx_2,\, dx_1 \wedge dx_3,\, dx_2 \wedge dx_3)$: cyclic ordering ensures all three terms contribute with the same sign to the volume form.
Summing:
\begin{align*}
d\beta_F = \left( \frac{\partial F_1}{\partial x_1} + \frac{\partial F_2}{\partial x_2} + \frac{\partial F_3}{\partial x_3} \right) dx_1 \wedge dx_2 \wedge dx_3.
\end{align*}
The coefficient is the divergence $\nabla \cdot F = \sum_i \partial_{x_i} F_i$, so $d\beta_F = (\nabla \cdot F)\, dx_1 \wedge dx_2 \wedge dx_3$.[/guided]