Consider the Dirac delta distribution $\delta$ in $\mathbf{R}^d$. It is quite standard to approximate it by functions $g_n$ with $\|g_n\|_{L^1} = 1$.
Is it possible to choose a sequence of test functions $\{g_n\}$ converging to $\delta$ as a distribution such that their derivatives are $L^1$-bounded? I mean, such that for a given a natural number $k$ there is a constant $M>0$ so that
$\|\partial^\mu g_n \|_{L^1} \le M$
for all $n$'s and all multi-indices $|\mu| \le k$?