The Sagdeev pseudopotential (SP) method is used to study ion acoustic solitary waves (IASWs) in a warm, magnetized plasma with relativistic electrons. Employing the pseudopotential approach allows for the investigation of solitary wave (SW) structures across arbitrary amplitudes. The study highlights the simultaneous occurrence of compressive \(\left( {N > 1} \right)\) subsonic \(\left( {M < 1} \right)\) solitons, as well as rarefactive \(\left( {N < 1} \right)\) subsonic and supersonic \(\left( {M > 1} \right)\) solitons, under specific parametric conditions. Notably, it is seen that as the direction cosine of wave propagation \(k_{z}\) increases, both the amplitude of SWs and the depth of the potential well decrease. The reduction in amplitude indicates a closer alignment between the magnetic field lines and the direction of wave propagation. The coexistence of compressive subsonic, rarefactive subsonic, and supersonic solitons in this plasma model is a rich and complex phenomenon that has both fundamental and practical implications in plasma physics. It reflects the intricate interplay of nonlinear effects, particle dynamics, and wave propagation in plasmas, with potential applications in both laboratory and astrophysical contexts.
