ABSTRACT We extend the range of validity of the artis 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae (SNe Ia) are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium population and ionization solver, a new multifrequency radiation field model, and a new atomic data set with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionization, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf (WD) in which the possible effects of gravitational settling of 22Ne prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M⊙ WD model obtained when 5.5 Gyr of very efficient settling is assumed to a similar model without settling. We find that this degree of 22Ne settling has only a modest effect on the resulting nebular spectra due to increased 58Ni abundance. Due to the high ionization in sub-Chandrasekhar models, the nebular [Ni ii] emission remains negligible, while the [Ni iii] line strengths are increased and the overall ionization balance is slightly lowered in the model with 22Ne settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe iii] emission relative to [Fe ii] in comparison to observations of normal SNe Ia.