Let . In the first two components of the vector field, replacing by their negatives changes the sign of the whole component; in the third, the product is unchanged. Thus , establishing the symmetry, with solutions mapped to solutions by .
At an equilibrium point, positivity of makes . The remaining equations give and . Hence
where the two nonzero equilibrium points exist when . They merge with the origin at , in a pitchfork bifurcation.
At the origin the Jacobian matrix and its eigenvalues are
At these are . To see the local pitchfork bifurcation dynamics more explicitly put , , and use on the center manifold. The equations become
At , invariance at cubic order gives , hence the reduced equation is
It has the supercritical pitchfork bifurcation form as increases, with the two branches attracting near onset while the origin becomes unstable.
For , the determinant of the planar block is , so the origin has one positive and two negative eigenvalues and is a saddle equilibrium. At , the planar eigenvalues are explicitly . For , the planar trace is negative and its determinant positive, so the origin is asymptotically stable. To distinguish its type put
For both planar eigenvalues are real negative, giving a stable node. At they coincide at ; the planar block is not a scalar matrix, so it is a defective stable node. For they are nonreal with negative real part, giving a stable spiral in the planar block together with the stable third direction. This refines the requested node-versus-saddle classification.
The origin cannot have a Hopf bifurcation: a planar imaginary pair would require trace zero, namely , but the determinant there is , so the pair is real and of opposite signs.
For either nonzero equilibrium point put . The Jacobian matrix is
Its characteristic polynomial is
A nonzero imaginary root of an equation requires, by separating real and imaginary parts, and . Equivalently
so the candidate Hopf bifurcation values are
At such an admissible value the factorization is , giving the pair and a negative third eigenvalue , since on the nonzero branches.
The mathematical model parameter restrictions matter. A real candidate requires , and the nonzero equilibrium points also require . The upper-sign value is admissible for ; the lower-sign value additionally requires . For fixed , differentiation of the characteristic equation at the imaginary root of an equation gives the eigenvalue crossing speed
Thus the transverse eigenvalue crossing condition holds except at the coincident-root case , . At a transverse eigenvalue crossing, a generic nonzero first Lyapunov coefficient gives a Hopf bifurcation. The displayed radical is the spectral threshold; it does not justify an unconditional generic-Hopf claim for all positive mathematical model parameters, nor does it by itself exclude a degenerate Hopf bifurcation point. For example yields no real candidate at all. The cubic Routh-Hurwitz stability criterion gives local Lyapunov stability on the side with and , where
For the modified equation, the planar block at the origin is , with determinant for the stipulated positive mathematical model parameters. Thus the origin is always a saddle equilibrium for nonzero in that mathematical model parameter range. The nonzero equilibrium points now satisfy , and , so they are
For these branches stay a positive distance from the origin and do not merge with it as varies. The origin also has no zero eigenvalue for . Therefore the pitchfork bifurcation is absent in the modified system.
A transverse eigenvalue crossing occurs when the real part of an eigenvalue passes through zero with nonzero parameter derivative. This spectral condition is separate from the nonlinear nondegeneracy conditions in the Hopf theorem.