ABSTRACT Tellurium is a primary candidate for the identification of the 2.1 $\, \mu$m emission line in kilonovae (KNe) spectra AT2017gfo and GRB230307A. Despite this, there is currently an insufficient amount of atomic data available for this species. We calculate the required atomic structure and collisional data, particularly the data required for accurate non-local-thermodynamic-equilibrium (NLTE) modelling of the low temperatures and densities in KNe. We use a multiconfigurational Dirac–Hartree–Fock method to produce optimized one-electron orbitals for Te i-iii. As a result energy levels and Einstein A-coefficients for Te i-iii have been calculated. These orbitals are then employed within Dirac R-matrix collision calculations to provide electron-impact-excitation collision strengths that were subsequently averaged according to a thermal Maxwellian distribution. Subsequent tardis simulations using this new atomic data reveal no significant changes to the synthetic spectra due to the very minor contribution of Te at early epochs. NLTE simulations with the colradpy package reveal optically thin spectra consistent with the increasing prominence of the Te iii 2.1 $\, \mu$m line as the KNe ejecta cools. This is reinforced by the estimation of luminosities at nebular KNe conditions. New line ratios for both observation and laboratory benchmarks of the atomic data are proposed.