Let $u\in H^1(\mathbb{R}^3)$ and $a\in C^\infty(\mathbb{R}^3)$ such that $|\nabla a|,\Delta a \in L^\infty(\mathbb{R}^3)$, where the norm in the Sobolev space $H^1(\mathbb{R}^3)$ is given by $\|u\|_{H^1}=\|u\|_{L^2}+\|\nabla u\|_{L^2}$. I would like to show that $$\int_{\mathbb{R}^3}\nabla u\cdot\nabla a\, dx=-\int_{\mathbb{R}^3}u\Delta a\, dx.$$ My thoughts are to use integration by parts and show that the integral over the boundary is zero, which intuitively makes sense because $|\nabla a |$ is bounded and $u$ decays. In the reference given in this answer the authors states, in particular, that if $B_r\subset \mathbb{R}^3$ is an open ball with radius $r$ and $f,g\in C^1(\mathbb{R}^3)$, then $$\int_{B_r} D_j f(x)g(x)dx=\int_{\partial {B_r}}(fg)(y)\nu_j(y)d_{2}y-\int_{B_r} f(x)D_jg(x)dx, $$ where $d_{2}y$ is the euclidean $2$-dimensional density on $\partial B_r$ and $\nu (y)$ the outer normal at $y\in\partial B_r$, which in this case i believe is $y/r$. Since $u$ is not necessarily $C^1$, i thought that i would need to approximate $u$ by smooth functions and take the limit $r\rightarrow \infty$, does this work? Also, i currently don't know how to evaluate the integral $$\int_{\partial {B_r}}(uD_ja)(y)\nu_j(y)d_{2}y.$$
Hints or references will be appreciated.
As $S(\mathbb{R}^3)\subset L^1(\mathbb{R}^3)$, taking $u_n\rightarrow u$ a.e., $u_n\in S(\mathbb{R}^3)$ Schwartz functions, for all $n\in\mathbb{N}$, we have $$ \left|\int_{\partial {B_r}}(u_nD_ja)(y)\frac{y_j}{r}dy\right|\leq \|{\nabla a}\|_{\infty}\int_{\partial {B_r}}|u_n|dy\stackrel{r\rightarrow\infty}{\longrightarrow} 0. $$ it follows that $$\int_{\mathbb{R}^3} D_j u_n(x)D_ja(x)dx=-\int_{\mathbb{R}^3} u_n(x)D^2_ja(x)dx.$$
but now i don't know what ensures that i can take the limit $n\rightarrow\infty$ inside the integrals. The dominated convergence theorem seems to require $u\in L^1(\mathbb{R}^3)$, i cannot think of any $L^1$ function to bound the sequence.
On the other hand, i could take $a_n \rightarrow a$ with $a_n$ smooth and compactly supported. As $u\in H^1(\mathbb{R}^3)$, integration by parts holds by the definition of weak derivate, but i think i cannot take $n\rightarrow \infty$.