Coalescence of binary neutron star give rise to electromagnetic emission, kilonova, powered by radioactive decays of r-process nuclei. Observations of kilonova associated with GW170817 provided unique opportunity to study the heavy element synthesis in the Universe. However, atomic data of r-process elements to decipher the light curves and spectral features of kilonova are not fully constructed yet. In this paper, we perform extended atomic calculations of neodymium (Nd, Z=60) to study the impact of accuracies in atomic calculations to the astrophysical opacities. By employing multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods, we calculate energy levels and transition data of electric dipole transitions for Nd II, Nd III, and Nd IV ions. Compared with previous calculations, our new results provide better agreement with the experimental data. The accuracy of energy levels was achieved in the present work 10 %, 3 % and 11 % for Nd II, Nd III and Nd IV, respectively, comparing with the NIST database. We confirm that the overall properties of the opacity are not significantly affected by the accuracies of the atomic calculations. The impact to the Planck mean opacity is up to a factor of 1.5, which affects the timescale of kilonova at most 20 %. However, we find that the wavelength dependent features in the opacity are affected by the accuracies of the calculations. We emphasize that accurate atomic calculations, in particular for low-lying energy levels, are important to provide predictions of kilonova light curves and spectra.