[guided]The purpose of the Bessel potentials is to convert an $H^s \to H^t$ estimate into an $L^2 \to L^2$ estimate. For $r \in \mathbb{R}$, the Bessel potential is the map
\begin{align*}
J^r: \mathcal{S}(\mathbb{R}^n) \to \mathcal{S}'(\mathbb{R}^n)
\end{align*}
defined on Schwartz functions by the Fourier multiplier
\begin{align*}
\widehat{J^r f}(\xi) = (1+|\xi|^2)^{r/2}\hat{f}(\xi).
\end{align*}
It extends by completion to an isometric isomorphism $J^r:H^r(\mathbb{R}^n)\to L^2(\mathbb{R}^n)$.
Thus controlling $\widetilde{S}f$ in $H^t(\mathbb{R}^n)$ is the same as controlling $J^t\widetilde{S}f$ in $L^2(\mathbb{R}^n)$, while the input norm is $\|J^s f\|_{L^2(\mathbb{R}^n)}$.
We now compute the kernel of the conjugated operator in Fourier variables. The mixed Fourier transform of the compactly supported smooth kernel $\widetilde{k}$ is the function $a: \mathbb{R}^n \times \mathbb{R}^n \to \mathbb{C}$ given by
\begin{align*}
a(\xi,\eta) = \frac{1}{(2\pi)^n}\int_{\mathbb{R}^n}\int_{\mathbb{R}^n}\widetilde{k}(x,y)e^{-ix\cdot \xi}e^{iy\cdot \eta}\,d\mathcal{L}^n(y)\,d\mathcal{L}^n(x).
\end{align*}
Because $\widetilde{k}$ is smooth and compactly supported in both variables, [integration by parts](/theorems/2098) in $x$ and $y$ gives rapid decay in both $\xi$ and $\eta$. Hence $a$ is a Schwartz function on $\mathbb{R}^n \times \mathbb{R}^n$.
For $f \in C_c^\infty(\mathbb{R}^n)$, Fourier inversion and [Fubini's theorem](/theorems/2961) apply because all functions involved are smooth and compactly supported or rapidly decreasing. They give
\begin{align*}
\widehat{\widetilde{S}f}(\xi) = \int_{\mathbb{R}^n} a(\xi,\eta)\hat{f}(\eta)\,d\mathcal{L}^n(\eta).
\end{align*}
After applying $J^t$ to the output and writing $\hat{f}(\eta)=(1+|\eta|^2)^{-s/2}\widehat{J^s f}(\eta)$, the relevant kernel is the function $b: \mathbb{R}^n \times \mathbb{R}^n \to \mathbb{C}$ defined by
\begin{align*}
b(\xi,\eta) = (1+|\xi|^2)^{t/2}a(\xi,\eta)(1+|\eta|^2)^{-s/2}.
\end{align*}
The factor $(1+|\xi|^2)^{t/2}(1+|\eta|^2)^{-s/2}$ grows at most polynomially in $(\xi,\eta)$, while $a$ decays faster than every polynomial. Therefore
\begin{align*}
b \in L^2(\mathbb{R}^n \times \mathbb{R}^n).
\end{align*}
Define the linear map
\begin{align*}
A: L^2(\mathbb{R}^n) \to L^2(\mathbb{R}^n)
\end{align*}
initially for $g \in C_c^\infty(\mathbb{R}^n)$ by
\begin{align*}
(Ag)(\xi) = \int_{\mathbb{R}^n} b(\xi,\eta)g(\eta)\,d\mathcal{L}^n(\eta).
\end{align*}
The $L^2$ hypothesis on $b$ is exactly the Hilbert-Schmidt condition. Indeed, Cauchy-Schwarz in the $\eta$ variable gives
\begin{align*}
|(Ag)(\xi)|^2 \leq \left(\int_{\mathbb{R}^n}|b(\xi,\eta)|^2\,d\mathcal{L}^n(\eta)\right)\left(\int_{\mathbb{R}^n}|g(\eta)|^2\,d\mathcal{L}^n(\eta)\right).
\end{align*}
Now integrate both sides over $\xi \in \mathbb{R}^n$ with respect to $\mathcal{L}^n$. This gives
\begin{align*}
\|Ag\|_{L^2(\mathbb{R}^n)}^2 \leq \|b\|_{L^2(\mathbb{R}^n \times \mathbb{R}^n)}^2\|g\|_{L^2(\mathbb{R}^n)}^2.
\end{align*}
Taking square roots proves that $A$ is bounded on $L^2(\mathbb{R}^n)$.
Finally set $g=J^s f$. The previous estimate becomes
\begin{align*}
\|\widetilde{S}f\|_{H^t(\mathbb{R}^n)} = \|J^t\widetilde{S}f\|_{L^2(\mathbb{R}^n)} \leq \|b\|_{L^2(\mathbb{R}^n \times \mathbb{R}^n)}\|J^s f\|_{L^2(\mathbb{R}^n)}.
\end{align*}
Thus
\begin{align*}
\|\widetilde{S}f\|_{H^t(\mathbb{R}^n)} \leq \|b\|_{L^2(\mathbb{R}^n \times \mathbb{R}^n)}\|f\|_{H^s(\mathbb{R}^n)}.
\end{align*}
Since $C_c^\infty(\mathbb{R}^n)$ is dense in $H^s(\mathbb{R}^n)$ for every $s \in \mathbb{R}$, the estimate extends $\widetilde{S}$ uniquely to a bounded linear map from $H^s(\mathbb{R}^n)$ to $H^t(\mathbb{R}^n)$.[/guided]