Grassmann variation of a chiral Majorana mass term ($i/\partial \psi -m^{\ast }{\psi }^c=0,i/\partial {\psi }^c-m\psi =0$)

For anticommuting components and an antisymmetric numerical matrix $A$, $\delta ({\psi }^{T}A\psi )=2\delta {\psi }^{T}A\psi$: move the second varied field past the first and use $A^T=-A$. Applying this to $C^{-1}$ and to $C$ cancels the factor $1/2$ in $\mathcal{L}=\bar{\psi }i/\partial \psi +\frac{1}{2}(m{\psi }^T{C}^{-1}\psi -m^{\ast }\bar{\psi }C{\bar{\psi }}^T)$. Varying the Dirac adjoint gives the first displayed equation. Integration by parts in the spinor variation gives $i{\gamma }^{\mu T}{\partial }_{\mu }{\bar{\psi }}^T+m{C}^{-1}\psi =0$; multiplication by $C$ gives the second. Squaring the coupled equations yields $(\square +|m{|}^2)\psi =0$, so the physical Majorana fermion mass is $|m|$.