Appendix 1
The time derivatives of \(A\left( {T_{1} ,T_{2} } \right)\) and \(B\left( {T_{1} ,T_{2} } \right)\) w.r.t \(T_{2}\) are given by
$$ \begin{gathered} \frac{{\partial A\left( {T_{1} ,T_{2} } \right)}}{{\partial T_{2} }} = \frac{I}{{8\left( {\omega_{2} + 2\omega_{1} } \right)\left( {2\omega_{1} - \omega_{2} } \right)\varepsilon \omega_{1}^{2} \omega_{2}^{2} }}( - 24\omega_{1} B\left( {T_{1} ,T_{2} } \right)^{2} \lambda^{2} \overline{A}\left( {T_{1} ,T_{2} } \right)\varepsilon \omega_{2}^{2} \hfill \\ e^{{ - 2iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} + 16\omega_{1}^{3} \lambda B\left( {T_{1} ,T_{2} } \right)^{2} \overline{A}\left( {T_{1} ,T_{2} } \right)\omega_{2}^{2} \varepsilon e^{{ - 2iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 32\omega_{1}^{2} B\left( {T_{1} ,T_{2} } \right)^{2} \lambda^{2} \hfill \\ \overline{A}\left( {T_{1} ,T_{2} } \right)\omega_{2} \varepsilon e^{{ - 2iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 4\lambda \omega_{1} B\left( {T_{1} ,T_{2} } \right)^{2} \overline{A}\left( {T_{1} ,T_{2} } \right)\omega_{2}^{4} \varepsilon e^{{ - 2iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} + 32\omega_{1}^{3} \hfill \\ B\left( {T_{1} ,T_{2} } \right)^{2} \lambda^{2} \overline{A}\left( {T_{1} ,T_{2} } \right)\varepsilon e^{{ - 2iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} + 24\lambda^{2} A\left( {T_{1} ,T_{2} } \right)^{2} \overline{A}\left( {T_{1} ,T_{2} } \right)\omega_{2}^{2} \varepsilon \omega_{1} - 192\omega_{1}^{3} \hfill \\ B\left( {T_{1} ,T_{2} } \right)\lambda^{2} \overline{B}\left( {T_{1} ,T_{2} } \right)A\left( {T_{1} ,T_{2} } \right)\varepsilon + 48\omega_{1}^{3} \lambda A\left( {T_{1} ,T_{2} } \right)^{2} \overline{A}\left( {T_{1} ,T_{2} } \right)\varepsilon \omega_{2}^{2} - 12\lambda \omega_{1} A\left( {T_{1} ,T_{2} } \right)^{2} \hfill \\ \overline{A}\left( {T_{1} ,T_{2} } \right)\omega_{2}^{3} \varepsilon - 12\lambda \omega_{1} A\left( {T_{1} ,T_{2} } \right)^{2} \overline{A}\left( {T_{1} ,T_{2} } \right)\omega_{2}^{4} \varepsilon + 16I\omega_{1}^{3} G_{h} B\left( {T_{1} ,T_{2} } \right)\omega_{2}^{3} \varepsilon \hfill \\ e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 4i\omega_{1}^{2} \mu_{h1} A\left( {T_{1} ,T_{2} } \right)\omega_{2}^{4} \varepsilon - 4i\omega_{1} G_{h} B\left( {T_{1} ,T_{2} } \right)\omega_{2}^{5} \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} \hfill \\ + 16i\omega_{1}^{4} H_{h} A\left( {T_{1} ,T_{2} } \right)\omega_{2}^{2} \varepsilon + 16i\omega_{1}^{4} \mu_{h1} A\left( {T_{1} ,T_{2} } \right)\omega_{2}^{2} \varepsilon - 4i\omega_{1}^{2} H_{h} A\left( {T_{1} ,T_{2} } \right)\omega_{2}^{4} \varepsilon \hfill \\ - 64\omega_{1}^{3} \lambda^{2} A\left( {T_{1} ,T_{2} } \right)^{2} \overline{A}\left( {T_{1} ,T_{2} } \right)\varepsilon - 4\omega_{1}^{3} f_{h} {\Omega }^{2} \omega_{2}^{2} e^{{\frac{{iT_{1} \left( { - \omega_{1} + {\Omega }} \right)}}{\varepsilon }}} + 4\omega_{1}^{2} f_{h} {\Omega }^{3} \omega_{2}^{2} e^{{\frac{{iT_{1} \left( { - \omega_{1} + {\Omega }} \right)}}{\varepsilon }}} \hfill \\ + \omega_{1} f_{h} {\Omega }^{2} \omega_{2}^{4} e^{{\frac{{iT_{1} \left( { - \omega_{1} + {\Omega }} \right)}}{\varepsilon }}} - f_{h} {\Omega }^{3} \omega_{2}^{4} e^{{\frac{{iT_{1} \left( { - \omega_{1} + {\Omega }} \right)}}{\varepsilon }}} + 32\omega_{1}^{3} \lambda B\left( {T_{1} ,T_{2} } \right)\overline{B}\left( {T_{1} ,T_{2} } \right)A\left( {T_{1} ,T_{2} } \right) \hfill \\ \omega_{2}^{2} \varepsilon + 16\omega_{1} B\left( {T_{1} ,T_{2} } \right)\lambda^{2} \overline{B}\left( {T_{1} ,T_{2} } \right)A\left( {T_{1} ,T_{2} } \right)\omega_{2}^{2} \varepsilon - 8\lambda \omega_{1} B\left( {T_{1} ,T_{2} } \right)\overline{B}\left( {T_{1} ,T_{2} } \right)A\left( {T_{1} ,T_{2} } \right) \hfill \\ \omega_{2}^{4} \varepsilon ), \hfill \\ \end{gathered} $$
$$ \begin{gathered} \frac{{\partial B\left( {T_{1} ,T_{2} } \right)}}{{\partial T_{2} }} = \frac{1}{{8e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} \left( {2\omega_{1} - \omega_{2} } \right)\left( {2\omega_{1} + \omega_{2} } \right)\varepsilon \omega_{2}^{3} }}( - 4\omega_{1}^{2} f_{h} {\Omega }^{2} \omega_{2}^{2} \hfill \\ e^{{\frac{{i\left( { - \omega_{1} T_{0} \varepsilon - T_{1} \omega_{2} + \omega_{2} \varepsilon T_{0} + T_{1} {\Omega }} \right)}}{\varepsilon }}} + 16G_{h} A\left( {T_{1} ,T_{2} } \right)\varepsilon \omega_{2}^{2} \omega_{1}^{3} + 4\mu_{h2} B\left( {T_{1} ,T_{2} } \right)\omega_{2}^{5} \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} \hfill \\ + 8i\lambda^{2} A\left( {T_{1} ,T_{2} } \right)^{2} \overline{B}\left( {T_{1} ,T_{2} } \right)\omega_{2}^{2} \varepsilon e^{{iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 4i\lambda A\left( {T_{1} ,T_{2} } \right)^{2} \overline{B}\left( {T_{1} ,T_{2} } \right)\omega_{2}^{4} \varepsilon e^{{iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} \hfill \\ + 32i\omega_{1}^{2} \lambda A\left( {T_{1} ,T_{2} } \right)\overline{A}\left( {T_{1} ,T_{2} } \right)B\left( {T_{1} ,T_{2} } \right)\omega_{2}^{2} \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - f_{h} {\Omega }^{3} \omega_{2}^{3} \hfill \\ e^{{\frac{{i\left( { - \omega_{1} T_{0} \varepsilon - T_{1} \omega_{2} + \omega_{2} \varepsilon T_{0} + T_{1} {\Omega }} \right)}}{\varepsilon }}} - 12i\lambda B\left( {T_{1} ,T_{2} } \right)^{2} \overline{B}\left( {T_{1} ,T_{2} } \right)\varepsilon \omega_{2}^{4} e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} + 4H_{h} B\left( {T_{1} ,T_{2} } \right) \hfill \\ \omega_{2}^{4} e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} + 4H_{h} B\left( {T_{1} ,T_{2} } \right)\omega_{2}^{5} \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 480i\omega_{1}^{2} B\left( {T_{1} ,T_{2} } \right)^{2} \lambda^{2} \overline{B}\left( {T_{1} ,T_{2} } \right) \hfill \\ \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 16\omega_{1}^{2} H_{h} B\left( {T_{1} ,T_{2} } \right)\omega_{2}^{3} \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 32i\omega_{1} A\left( {T_{1} ,T_{2} } \right)^{2} \lambda^{2} \overline{B}\left( {T_{1} ,T_{2} } \right) \hfill \\ \omega_{2} \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 4\omega_{1} G_{h} A\left( {T_{1} ,T_{2} } \right)\omega_{2}^{4} \varepsilon + 16iB\left( {T_{1} ,T_{2} } \right)\lambda^{2} A\left( {T_{1} ,T_{2} } \right)\overline{A}\left( {T_{1} ,T_{2} } \right)\varepsilon \omega_{2}^{2} \hfill \\ e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} + 120iB\left( {T_{1} ,T_{2} } \right)^{2} \lambda^{2} \overline{B}\left( {T_{1} ,T_{2} } \right)\varepsilon \omega_{2}^{2} e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 192i\omega_{1}^{2} B\left( {T_{1} ,T_{2} } \right)\lambda^{2} \hfill \\ A\left( {T_{1} ,T_{2} } \right)\overline{A}\left( {T_{1} ,T_{2} } \right)\varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} + 16i\omega_{1}^{2} \lambda A\left( {T_{1} ,T_{2} } \right)^{2} \overline{B}\left( {T_{1} ,T_{2} } \right)\omega_{2}^{2} \varepsilon e^{{iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} \hfill \\ + 4\omega_{1}^{2} f_{h} {\Omega }^{3} \omega_{2} e^{{\frac{{i\left( { - \omega_{1} T_{0} \varepsilon - T_{1} \omega_{2} + \omega_{2} \varepsilon T_{0} + T_{1} {\Omega }} \right)}}{\varepsilon }}} + 48i\omega_{1}^{2} \lambda B\left( {T_{1} ,T_{2} } \right)^{2} \overline{B}\left( {T_{1} ,T_{2} } \right)\varepsilon \omega_{2}^{2} \hfill \\ e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 16\omega_{1}^{2} \mu_{h2} B\left( {T_{1} ,T_{2} } \right)\omega_{2}^{3} \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} - 8i\lambda A\left( {T_{1} ,T_{2} } \right)\overline{A}\left( {T_{1} ,T_{2} } \right)B\left( {T_{1} ,T_{2} } \right) \hfill \\ \omega_{2}^{4} \varepsilon e^{{ - iT_{0} \left( { - \omega_{2} + \omega_{1} } \right)}} + f_{h} {\Omega }^{2} \omega_{2}^{4} e^{{\frac{{i\left( { - \omega_{1} T_{0} \varepsilon - T_{1} \omega_{2} + \omega_{2} \varepsilon T_{0} + T_{1} {\Omega }} \right)}}{\varepsilon }}} . \hfill \\ \end{gathered} $$
Appendix 2
The coefficients of the Jacobian matrix [\({\text{J}}\)] as shown in Eq. (27) are expressed as
$$ J_{11} = - \frac{1}{2}\mu_{1} - \frac{1}{2}H - \frac{\lambda }{{8\omega_{1} }}a_{20}^{2} \sin \left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{a_{20}^{2} \sin \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} }}{{4\omega_{1} \left( {2\omega_{1} - \omega_{2} } \right)\omega_{2}^{2} }}\left( {2\omega_{1} - 3\omega_{2} } \right), $$
$$ J_{12} = \frac{1}{{\omega_{1} }}Ga_{20} \omega_{2} \sin \left( {\phi_{10} - \phi_{20} } \right) - \frac{\lambda }{{4\omega_{1} }}a_{20}^{2} a_{10} \cos \left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{a_{20}^{2} \cos \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} a_{10} }}{{2\omega_{1} \left( {2\omega_{1} - \omega_{2} } \right)\omega_{2}^{2} }}\left( {2\omega_{1} - 3\omega_{2} } \right) - \frac{{f{\Omega }^{2} \left( { - 3\omega_{1} + {\Omega }} \right)}}{{4\omega_{1}^{2} }}\cos \left( {\phi_{10} } \right), $$
$$ J_{13} = - \frac{1}{{2\omega_{1} }}G\omega_{2} \cos \left( {\phi_{10} - \phi_{20} } \right) - \frac{\lambda }{{4\omega_{1} }}a_{20} a_{10} \sin \left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{a_{20} \sin \left( {4\phi_{10} - 4\phi_{20} } \right)\lambda^{2} a_{10} }}{{2\omega_{1} \left( {2\omega_{1} - \omega_{2} } \right)\omega_{2}^{2} }}\left( {2\omega_{1} - 3\omega_{2} } \right), $$
$$ J_{14} = - \frac{1}{{2\omega_{1} }}Ga_{20} \omega_{2} \sin \left( {\phi_{10} - \phi_{20} } \right) + \frac{\lambda }{{4\omega_{1} }}a_{20}^{2} a_{10} \cos \left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{4a_{20}^{2} \cos \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} a_{10} }}{{\omega_{1} \left( {2\omega_{1} - \omega_{2} } \right)\omega_{2}^{2} }}\left( {2\omega_{1} - 3\omega_{2} } \right), $$
$$ J_{21} = \frac{{f{\Omega }^{2} \left( {{\Omega } - 3\omega_{1} } \right)}}{{4a_{10}^{2} \omega_{1}^{2} }}\cos \left( {\phi_{10} } \right) - \frac{{3\lambda a_{10} }}{{4\omega_{1} }} + \frac{{\lambda^{2} }}{{4\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{1} \omega_{2}^{2} }}\left( { - 6a_{10} \omega_{2}^{2} + 16\omega_{1}^{2} a_{10} } \right) - \frac{1}{{2a_{10}^{2} \omega_{1} }}Ga_{20} \omega_{2} \sin \left( {\phi_{10} - \phi_{20} } \right), $$
$$ J_{22} = \frac{{f{\Omega }^{2} \left( {{\Omega } - 3\omega_{1} } \right)}}{{4a_{10} \omega_{1}^{2} }}\sin \left( {2\phi_{10} } \right) + \frac{{\lambda a_{20}^{2} }}{{4\omega_{1} }}\sin \left( {2\phi_{10} - 2\phi_{10} } \right) + \frac{{a_{20}^{2} \sin \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} }}{{2\omega_{1} \left( {2\omega_{1} - \omega_{2} } \right)\omega_{2}^{2} }}\left( {2\omega_{1} - 3\omega_{2} } \right) + \frac{1}{{2a_{10} \omega_{1} }}Ga_{20} \omega_{2} \cos \left( {\phi_{10} - \phi_{20} } \right), $$
$$ J_{23} = - \frac{\lambda }{{8\omega_{1} }}\left( {4a_{20} + 2a_{20} \cos \left( {2\phi_{10} - 2\phi_{20} } \right)} \right) + \frac{{\lambda^{2} }}{{4\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{1} \omega_{2}^{2} }}(48\omega_{1}^{2} a_{20} - 4\omega_{2}^{2} a_{20} ) - \frac{{a_{20} \cos \left( {2\phi_{10} - 2\phi_{10} } \right)\lambda^{2} }}{{2\omega_{1} \left( {2\omega_{1} - \omega_{2} } \right)\omega_{2}^{2} }}\left( {2\omega_{1} - 3\omega_{2} } \right) + \frac{1}{{2a_{10} \omega_{1} }}G\omega_{2} \sin \left( {\phi_{10} - \phi_{20} } \right), $$
$$ J_{24} = - \frac{{\lambda a_{20}^{2} }}{{4\omega_{1} }}\sin \left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{a_{20}^{2} \sin \left( {2\phi_{10} - 2\phi_{10} } \right)\lambda^{2} }}{{2\omega_{1} \left( {2\omega_{1} - \omega_{2} } \right)\omega_{2}^{2} }}\left( {2\omega_{1} - 3\omega_{2} } \right) - \frac{1}{{2a_{10} \omega_{1} }}Ga_{20} \omega_{2} \cos \left( {\phi_{10} - \phi_{20} } \right), $$
$$ J_{31} = \frac{1}{{2\omega_{2} }}G\omega_{1} \cos \left( {\phi_{10} - \phi_{20} } \right) + \frac{\lambda }{{4\omega_{2} }}a_{10} a_{20} \sin \left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{a_{10} \sin \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} a_{20} }}{{2\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{2}^{2} }}\left( { - \omega_{2} + 4\omega_{1} } \right), $$
$$ J_{32} = - \frac{1}{{2\omega_{2} }}G\omega_{1} a_{10} \sin \left( {\phi_{10} - \phi_{20} } \right) + \frac{\lambda }{{4\omega_{2} }}a_{10}^{2} a_{20} {\text{cos}}\left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{a_{10}^{2} \cos \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} a_{20} }}{{2\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{2}^{2} }}\left( { - \omega_{2} + 4\omega_{1} } \right), $$
$$ J_{33} = - \frac{1}{2}\mu_{2} - \frac{1}{2}H + \frac{\lambda }{{8\omega_{2} }}a_{10}^{2} \sin \left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{a_{10}^{2} \sin \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} }}{{4\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{2}^{2} }}\left( { - \omega_{2} + 4\omega_{1} } \right), $$
$$ J_{34} = \frac{1}{{\omega_{2} }}Ga_{10} \omega_{1} \sin \left( {\phi_{10} - \phi_{20} } \right) - \frac{\lambda }{{4\omega_{2} }}a_{10}^{2} a_{20} {\text{cos}}\left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{a_{10}^{2} {\text{cos}}\left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} a_{20} }}{{2\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{2}^{2} }}\left( { - \omega_{2} + 4\omega_{1} } \right) - \frac{{f{\Omega }^{2} \left( {{\Omega } - 3\omega_{2} } \right)}}{{4\omega_{2}^{2} }}\sin \left( {2\phi_{20} } \right), $$
$$ J_{41} = - \frac{\lambda }{{8\omega_{2} }}\left( {4a_{10} + 2a_{10} \cos \left( {2\phi_{10} - 2\phi_{20} } \right)} \right) + \frac{{\lambda^{2} }}{{4\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{2}^{3} }}( - 4a_{10} \omega_{2}^{2} + 48\omega_{1}^{2} a_{10} ) + \frac{{a_{10} \cos \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} }}{{2\omega_{2}^{2} \left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)}}\left( {4\omega_{1} - \omega_{2} } \right) + \frac{1}{{2a_{20} \omega_{2} }}G\omega_{1} \sin \left( {2\phi_{10} - 2\phi_{20} } \right), $$
$$ J_{42} = \frac{{\lambda a_{10}^{2} }}{{4\omega_{2} }}\sin \left( {2\phi_{10} - 2\phi_{20} } \right) - \frac{{\lambda^{2} a_{10}^{2} \sin \left( {2\phi_{10} - 2\phi_{20} } \right)}}{{2\omega_{2}^{2} \left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)}}\left( {4\omega_{1} - \omega_{2} } \right) + \frac{1}{{2a_{20} \omega_{2} }}Ga_{10} \omega_{1} \cos \left( {\phi_{10} - \phi_{20} } \right), $$
$$ J_{43} = \frac{{f{\Omega }^{2} \left( {{\Omega } - 3\omega_{2} } \right)}}{{4a_{20}^{2} \omega_{2}^{2} }}{\text{sin}}\left( {\phi_{20} } \right) - \frac{{3\lambda a_{20} }}{{4\omega_{2} }} + \frac{{\lambda^{2} }}{{4\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{2}^{3} }}\left( {120\omega_{1}^{2} a_{20} - 30\omega_{2}^{2} a_{20} } \right) - \frac{1}{{2a_{20}^{2} \omega_{2} }}Ga_{10} \omega_{1} \sin \left( {\phi_{10} - \phi_{20} } \right), $$
$$ J_{44} = - \frac{{f{\Omega }^{2} \left( {{\Omega } - 3\omega_{2} } \right)}}{{4a_{20} \omega_{2}^{2} }}{\text{cos}}\left( {\phi_{20} } \right) - \frac{\lambda }{{4\omega_{2} }}a_{10}^{2} \sin \left( {2\phi_{10} - 2\phi_{20} } \right) + \frac{{a_{10}^{2} \sin \left( {2\phi_{10} - 2\phi_{20} } \right)\lambda^{2} }}{{2\left( {4\omega_{1}^{2} - \omega_{2}^{2} } \right)\omega_{2}^{2} }}\left( {4\omega_{1} - \omega_{2} } \right) - \frac{1}{{2a_{20} \omega_{2} }}Ga_{10} \omega_{1} \cos \left( {\phi_{10} - \phi_{20} } \right). $$
Appendix 3
The coefficients of characteristic Eq. (29) are \(\alpha_{1}\), \(\alpha_{2}\), \(\alpha_{3}\) and \(\alpha_{4}\) can be written as
$$ \alpha_{1} = - J_{44} - J_{22} - J_{11} - J_{33} , $$
$$ \alpha_{2} = - J_{14} J_{41} - J_{12} J_{21} - J_{13} J_{31} + J_{11} J_{44} + J_{22} J_{44} - J_{42} J_{24} + J_{33} J_{44} - J_{34} J_{43} + J_{11} J_{33} + J_{22} J_{33} - J_{32} J_{23} + J_{11} J_{22} , $$
$$ \begin{gathered} \alpha_{3} = J_{14} J_{41} J_{22} + J_{42} J_{24} J_{33} - J_{22} J_{33} J_{44} + J_{34} J_{43} J_{11} - J_{12} J_{31} J_{23} + J_{34} J_{43} J_{22} \hfill \\ \quad - J_{42} J_{14} J_{21} + J_{13} J_{31} J_{44} - J_{22} J_{11} J_{44} + J_{12} J_{21} J_{33} - J_{11} J_{33} J_{44} - J_{14} J_{31} J_{43} \hfill \\ \quad + J_{32} J_{23} J_{11} - J_{12} J_{41} J_{24} + J_{32} J_{23} J_{44} - J_{32} J_{13} J_{21} - J_{22} J_{33} J_{11} - J_{42} J_{23} J_{34} \hfill \\ \quad + J_{34} J_{13} J_{41} + J_{14} J_{41} J_{33} + J_{42} J_{24} J_{11} + J_{12} J_{21} J_{44} - J_{32} J_{43} J_{24} + J_{13} J_{31} J_{22} , \hfill \\ \end{gathered} $$
$$ \begin{gathered} \alpha_{4} = J_{22} J_{11} J_{33} J_{44} - J_{32} J_{13} J_{41} J_{24} + J_{32} J_{23} J_{14} J_{41} - J_{12} J_{31} J_{43} J_{24} - J_{42} J_{23} J_{14} J_{31} \hfill \\ \quad - J_{34} J_{43} J_{22} J_{11} - J_{34} J_{43} J_{22} J_{11} + J_{42} J_{13} J_{31} J_{24} - J_{12} J_{41} J_{23} J_{34} + J_{32} J_{43} J_{24} J_{11} \hfill \\ \quad + J_{42} J_{14} J_{21} J_{33} - J_{32} J_{23} J_{11} J_{44} + J_{34} J_{13} J_{41} J_{22} + J_{12} J_{31} J_{23} J_{44} + J_{42} J_{23} J_{34} J_{11} \hfill \\ \quad - J_{32} J_{43} J_{14} J_{21} - J_{42} J_{24} J_{11} J_{33} + J_{12} J_{41} J_{24} J_{33} - J_{13} J_{31} J_{22} J_{44} - J_{14} J_{41} J_{22} J_{33} \hfill \\ \quad - J_{42} J_{13} J_{21} J_{34} - J_{12} J_{21} J_{33} J_{44} + J_{14} J_{31} J_{43} J_{22} + J_{32} J_{13} J_{21} J_{44} + J_{12} J_{21} J_{34} J_{43} . \hfill \\ \end{gathered} $$
Appendix 4
The coefficients of Eq. (32) are written as
$$ \begin{gathered} q_{6} = 256\omega_{1}^{6} \lambda^{4} - 36\omega_{1}^{2} \lambda^{3} \omega_{2}^{6} - 72\omega_{1}^{4} \lambda^{2} \omega_{2}^{6} + 144\omega_{1}^{6} \lambda^{2} \omega_{2}^{4} - 384\omega_{1}^{6} \lambda^{3} \omega_{2}^{2} - 192\omega_{1}^{4} \lambda^{4} \omega_{2}^{2} \hfill \\ \quad + 240\omega_{1}^{4} \lambda^{3} \omega_{2}^{4} + 36\omega_{1}^{2} \lambda^{4} \omega_{2}^{4} + 9\omega_{1}^{2} \lambda^{2} \omega_{2}^{8} . \hfill \\ \end{gathered} $$
$$ p_{6} = 900\lambda^{4} + 9\lambda^{2} \omega_{2}^{4} - 180\lambda^{3} \omega_{2}^{2} . $$
$$ q_{41} = 1024\omega_{1}^{7} \lambda^{2} \omega_{2}^{2} + 384\omega_{1}^{5} \lambda \omega_{2}^{6} + 96\omega_{1}^{3} \lambda^{2} \omega_{2}^{6} - 640\omega_{1}^{5} \lambda^{2} \omega_{2}^{4} - 48\omega_{1}^{3} \lambda \omega_{2}^{8} - 768\omega_{1}^{7} \lambda \omega_{2}^{4} . $$
$$ p_{41} = - 48\omega_{2}^{5} \lambda {\Omega } + 480\lambda^{2} \omega_{2}^{3} . $$
$$ q_{40} = - 1024\omega_{1}^{8} \lambda^{2} \omega_{2}^{2} - 384\omega_{1}^{6} \lambda \omega_{2}^{6} - 96\omega_{1}^{4} \lambda^{2} \omega_{2}^{6} + 640\omega_{1}^{6} \lambda^{2} \omega_{2}^{4} + 48\omega_{1}^{4} \lambda \omega_{2}^{8} + 768\omega_{1}^{8} \lambda \omega_{2}^{4} . $$
$$ p_{40} = 48\lambda \omega_{2}^{6} - 480\lambda^{2} \omega_{2}^{4} . $$
$$ q_{22} = 1024\omega_{1}^{8} \omega_{2}^{4} - 512\omega_{1}^{6} \omega_{2}^{6} + 64\omega_{1}^{4} \omega_{2}^{8} . $$
$$ p_{22} = 64\omega_{2}^{6} . $$
$$ q_{21} = - 2048\omega_{1}^{9} \omega_{2}^{4} - 128\omega_{1}^{5} \omega_{2}^{8} + 1024\omega_{1}^{7} \omega_{2}^{6} . $$
$$ p_{21} = - 128\omega_{2}^{7} . $$
$$ q_{20} = - 256\omega_{1}^{6} H\omega_{2}^{6} \mu_{1} + 16\omega_{1}^{4} \omega_{2}^{8} \mu_{1}^{2} + 64\omega_{1}^{6} \omega_{2}^{8} - 128\omega_{1}^{6} H^{2} \omega_{2}^{6} - 512\omega_{1}^{8} \omega_{2}^{6} + 16\omega_{1}^{4} H^{2} \omega_{2}^{8} + 512\omega_{1}^{8} H\mu_{1} \omega_{2}^{4} + 1024\omega_{1}^{10} \omega_{2}^{4} + 32\omega_{1}^{4} H\mu_{1} \omega_{2}^{8} + 256\omega_{1}^{8} H^{2} \omega_{2}^{4} - 128\omega_{1}^{6} \mu_{1}^{2} \omega_{2}^{6} + 256\omega_{1}^{8} \omega_{2}^{4} \mu_{1}^{2} . $$
$$ p_{20} = 16\omega_{2}^{6} \mu_{2}^{2} + 32\omega_{2}^{6} H\mu_{2} + 64\omega_{2}^{8} + 16\omega_{2}^{6} H^{2} . $$
$$ q_{06} = - 64\omega_{1}^{4} \omega_{2}^{4} + 32\omega_{1}^{2} \omega_{2}^{6} - 4\omega_{2}^{8} . $$
$$ p_{06} = - 4\omega_{2}^{2} . $$
$$ q_{05} = - 192\omega_{2}^{6} \omega_{1}^{3} + 24\omega_{1} \omega_{2}^{8} + 384\omega_{2}^{4} \omega_{1}^{5} . $$
$$ p_{05} = 24\omega_{2}^{3} . $$
$$ q_{04} = - 576\omega_{1}^{6} \omega_{2}^{4} + 288\omega_{1}^{4} \omega_{2}^{6} - 36\omega_{2}^{8} \omega_{1}^{2} . $$
$$ p_{04} = - 36\omega_{2}^{4} . $$