The energy level schemes of the neodymium oxyhalides (NdOX, X = F, Cl, and Br) were studied and simulated with a phenomenological model accounting simultaneously for both the free ion interactions and the crystal fleld effect. The former included the electrostatic and interconfigurational interactions as well as the spin-orbit coupling. The simulations were carried out by using the data from the optical absorption and luminescence as well as the inelastic neutron scattering measured at low temperatures between 2.5 and 77 K. The values of the Slater integral F 2 describing the eJectrostatic interactions decrease while F 4 and F6 increase as a function of the ionic radius of the halide anion. The strength of the spin-orbit coupling is quite the same in all three matrices. The crystal field effect — measured as the crystal fleld strength parameter S — is almost twice as strong in the hexagonal NdOF matrix (650 cm -1 ) than in the tetragonal NdOCI or NdOBr (367 and 378 cm -1 , respectively). Similar evolution was obtained for the shortand mid-range crystal field strengths related to the spatial extension of the interaction.