  Starting from the Vlasov-Maxwell equations describing the dynamics of various
species in a quasi-neutral plasma immersed in an external solenoidal magnetic
field and utilizing a technique known as the hydrodynamic substitution, a
relativistic hydrodynamic system of equations governing the dynamics of various
species has been obtained. Based on the method of multiple scales, a system
comprising three nonlinear Schrodinger equation for the transverse envelopes of
the three basic whistler modes, has been derived. Using the method of formal
series of Dubois-Violette, a traveling wave solution of the derived set of
coupled nonlinear Schrodinger equations in both the relativistic and the non
relativistic case has been obtained. An intriguing feature of our description
is that whistler waves do not perturb the initial uniform density of plasma
electrons. The plasma response to the induced whistler waves consists in
transverse velocity redistribution, which follows exactly the behaviour of the
electromagnetic whistlers. This property may have an important application for
transverse focusing of charged particle beams in future laser plasma
accelerators. Yet another interesting peculiarity are the selection rules
governing the nonlinear mode coupling. According to these rules self coupling
between modes in the non relativistic regime is absent, which is a direct
consequence of the vector character of the interaction governed by the Lorentz
force.
