The interaction (oblique collision) of two ion acoustic solitons (IASs) in a magnetized relativistic degenerate plasma with relativistic degenerate electrons and non-degenerate cold ions is studied. The extended Poincaré–Lighthill–Kuo (PLK) method is used to obtain two Korteweg deVries (KdV) wave equations that describe the interacting IASs, then the phase shifts due to interaction are calculated. We studied influence of the fluid number density on the interaction process, interacting solitons phase shifts and also phase velocities. The introduced model is valid for astrophysical objects with high density matter such as white dwarfs, neutron stars, degenerate electrons gas in metals and laboratory degenerate plasma. An inverse proportionality between the phase shifts, phase velocity and the equilibrium electron fluid number density $$n_{eo}$$
n
eo
was established in the range $$10^{35}\,{\text {m}}^{-3}>n_{eo}>10^{38}\,{\text {m}}^{-3}$$
10
35
m
-
3
>
n
eo
>
10
38
m
-
3
. We found that the soliton waves get sharper (narrower) and higher with increasing the electrons fluid number density $$n_{eo}$$
n
eo
, and hence less spacial occupying. The phase shifts and the phase velocity remain approximately unchanged in the range of $$10^{35}\,{\text {m}}^{-3}<n_{eo}<10^{38}\,{\text {m}}^{-3}$$
10
35
m
-
3
<
n
eo
<
10
38
m
-
3
. The impact of the obliqueness angle $$\theta $$
θ
on the soliton interaction process is also studied.