Supersymmetry (SUSY) with bilinearly broken $$R$$ R parity (bRPV) offers an attractive possibility to explain the origin of neutrino masses and mixings. In such scenarios, the study of neutralino decays at colliders gives access to neutrino sector parameters. The ILC offers a very clean environment to study the neutralino properties as well as its subsequent decays, which typically involve a $$W$$ W or $$Z$$ Z boson and a lepton. This study is based on ILC beam parameters according to the Technical Design Report for a center of mass energy of $$500~\hbox {GeV}$$ 500 GeV . A full detector simulation of the International Large Detector (ILD) has been performed for all Standard Model backgrounds and for neutralino pair production within a simplified model. The bRPV parameters are fixed according to current neutrino data. In this scenario, the $$\tilde{\chi }^0_1$$ χ ~ 1 0 mass can be reconstructed with an uncertainty of $$\delta m^{{\mathrm {fit}}}_{\widetilde{\chi }^{0}_{1}}=(40\text {(stat.)} \oplus 50\text {(syst.)})~\hbox {MeV}$$ δ m χ ~ 1 0 fit = ( 40 (stat.) ⊕ 50 (syst.) ) MeV for an integrated luminosity of $$500\,\hbox {fb}^{-1}$$ 500 fb - 1 from direct $$\tilde{\chi }^0_1$$ χ ~ 1 0 pair production, thus, to a large extent independently of the rest of the SUSY spectrum. The achievable precision on the atmospheric neutrino mixing angle $$\sin ^2 \theta _{23}$$ sin 2 θ 23 from measuring the neutralino branching fractions BR( $$\tilde{\chi }_1^0\rightarrow W\mu $$ χ ~ 1 0 → W μ ) and BR( $$\tilde{\chi }_1^0\rightarrow W \tau $$ χ ~ 1 0 → W τ ) at the ILC is in the same range than current uncertainties from neutrino experiments. Thus, the ILC could have the opportunity to unveil the mechanism of neutrino mass generation.