We characterize the model spaces KΘ in which functions with smooth boundary extensions are dense. It is shown that such approximations are possible if and only if the singular measure associated to the singular inner factor of Θ is concentrated on a countable union of Beurling–Carleson sets. In fact, we use a duality argument to show that if there exists a restriction of the associated singular measure which does not assign positive measure to any Beurling–Carleson set, then even larger classes of functions, such as Hölder classes and large collections of analytic Sobolev spaces, fail to be dense. In contrast to earlier results on density of functions with continuous extensions to the boundary in KΘ and related spaces, the existence of a smooth approximant is obtained through a constructive method.