It is well known that for any inner function $\theta $ defined in the unit disk $\mathbb {D}$, the following two conditions: (i) there exists a sequence of polynomials $\{p_n\}_n$ such that $\lim _{n \to \infty } \theta (z) p_n(z) = 1$ for all $z \in \mathbb {D}$ and (ii) $\sup _n \| \theta p_n \|_\infty < \infty $, are incompatible, i.e., cannot be satisfied simultaneously. However, it is also known that if we relax the second condition to allow for arbitrarily slow growth of the sequence $\{ \theta (z) p_n(z)\}_n$ as $|z| \to 1$, then condition (i) can be met for some singular inner function. We discuss certain consequences of this fact which are related to the rate of decay of Taylor coefficients and moduli of continuity of functions in model spaces $K_\theta $. In particular, we establish a variant of a result of Khavinson and Dyakonov on nonexistence of functions with certain smoothness properties in $K_\theta $, and we show that the classical Aleksandrov theorem on density of continuous functions in $K_\theta $ is essentially optimal. We consider also the same questions in the context of de Branges–Rovnyak spaces $\mathcal {H}(b)$ and show that the corresponding approximation result also is optimal.

Given a holomorphic self-map $\varphi $ of $\mathbb {D}$ (the open unit disc in $\mathbb {C}$ ), the composition operator $C_{\varphi } f = f \circ \varphi $ , $f \in H^2(\mathbb {\mathbb {D}})$ , defines a bounded linear operator on the Hardy space $H^2(\mathbb {\mathbb {D}})$ . The model spaces are the backward shift-invariant closed subspaces of $H^2(\mathbb {\mathbb {D}})$ , which are canonically associated with inner functions. In this paper, we study model spaces that are invariant under composition operators. Emphasis is put on finite-dimensional model spaces, affine transformations, and linear fractional transformations.