[solution]
**Step 1: Smoothness via differentiation under the integral.** For any multi-index $\beta$, we differentiate under the integral sign: $\partial^\beta(f * g)(x) = (\partial^\beta f * g)(x)$. To justify this, we verify the hypotheses of the Leibniz integral rule. The integrand $f(x - y) g(y)$ is smooth in $x$ for each $y$, and the derivative $\partial_x^\beta f(x - y) g(y)$ is dominated by $\|\partial^\beta f\|_{0,0} \cdot |g(y)|$, which is integrable since $g \in \mathcal{S} \subseteq L^1$. By induction on $|\beta|$, the Leibniz rule applies at each stage, giving $\partial^\beta(f * g) = (\partial^\beta f) * g$, and since $\partial^\beta f \in L^1$ and $g \in L^1$, this convolution is a well-defined continuous function. Hence $f * g \in C^\infty(\mathbb{R}^n)$.
**Step 2: Semi-norm estimates.** Fix $\alpha, \beta \in \mathbb{N}_0^n$. By Step 1,
\begin{align*}
|x^\alpha \, \partial^\beta(f * g)(x)| &= \left|x^\alpha \int_{\mathbb{R}^n} (\partial^\beta f)(x - y) \, g(y) \, d\mathcal{L}^n(y)\right|.
\end{align*}
We use the binomial expansion for multi-indices: $x^\alpha = ((x-y) + y)^\alpha = \sum_{\gamma \leq \alpha} \binom{\alpha}{\gamma} (x-y)^\gamma \, y^{\alpha - \gamma}$. Substituting:
\begin{align*}
|x^\alpha \, \partial^\beta(f*g)(x)| &\leq \sum_{\gamma \leq \alpha} \binom{\alpha}{\gamma} \int_{\mathbb{R}^n} |(x-y)^\gamma \, (\partial^\beta f)(x-y)| \cdot |y^{\alpha-\gamma} \, g(y)| \, d\mathcal{L}^n(y).
\end{align*}
**Step 3: Bound each factor.** The first factor satisfies $|(x-y)^\gamma (\partial^\beta f)(x-y)| \leq \|f\|_{\gamma, \beta}$ by the definition of the Schwartz semi-norm, since $f \in \mathcal{S}(\mathbb{R}^n)$. For the second factor, we need $|y^{\alpha-\gamma} g(y)|$ to be integrable. Choose $N > n$. Since $g \in \mathcal{S}(\mathbb{R}^n)$, we have $|y^{\alpha-\gamma} g(y)| \leq \|g\|_{\alpha - \gamma + N\mathbf{e}, 0} \cdot (1 + |y|)^{-N}$ where $N\mathbf{e}$ is a multi-index of sufficient size, and $(1+|y|)^{-N}$ is integrable over $\mathbb{R}^n$ since $N > n$. Therefore
\begin{align*}
\int_{\mathbb{R}^n} |y^{\alpha-\gamma} g(y)| \, d\mathcal{L}^n(y) &\leq C_N \, \|g\|_{\alpha-\gamma + N\mathbf{e}, 0} < \infty.
\end{align*}
**Step 4: Conclude.** Combining:
\begin{align*}
\|f * g\|_{\alpha,\beta} &\leq \sum_{\gamma \leq \alpha} \binom{\alpha}{\gamma} \|f\|_{\gamma,\beta} \cdot C_N \, \|g\|_{\alpha - \gamma + N\mathbf{e}, 0} < \infty.
\end{align*}
Since $\alpha$ and $\beta$ were arbitrary, $f * g \in \mathcal{S}(\mathbb{R}^n)$.
[/solution]