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$?
As paul garrett commented, we can employ integration by parts. Let $g_n$ be a sequence of functions tending, as distributions, to the dirac $\delta$. Integration by parts yields $$ \lim_{n\to\infty}\int g_n^\prime(x)\,f(x/\epsilon)\,\mathrm{d}x=-f'(0)/\epsilon\tag{1} $$ Note that $\|f(x/\epsilon)\|_{L^\infty}=\|f(x)\|_{L^\infty}$. Therefore, Hölder's Inequality says $$ \lim_{n\to\infty}\|g_n'\|_{L^1}\|f\|_{L^\infty}\ge\left|\,\lim_{n\to\infty}\int g_n^\prime(x)\,f(x/\epsilon)\,\mathrm{d}x\,\right|=|f'(0)/\epsilon|\tag{2} $$ Since $\epsilon$ is arbitrary, we get that $$ \lim_{n\to\infty}\|g_n'\|_{L^1}=\infty\tag{3} $$