Construction of elliptic equation with Neumann boundary condition from a minimization problem My question mainly concern about how to construct a elliptic equation with Neumann boundary condition from a minimization problem.
Let $B=B_1 \subset \mathbb{R}^3$ and $E : H^1(B) \to \mathbb{R}$
$$E(u)= \int_{B}|\nabla u|^2+(u^2-1)^2 dx - \int_{\partial B}Q(u)d\mathcal{H}^2$$
We assume that $u_0 \in W^{1,2}$ to be the minimizer of the functional $E$ in the configuration space
$$K=\{u\in W^{1,2}(B:\mathbb{R})\}.$$
Since $u_0$ is the critical point of the functional, we let $\xi \in K$, we obtain the equation 
$$\int_B \nabla u \cdot \nabla \xi + 4(u^2-1)u \xi dx - \int_{\partial B }Q'(u)\xi d\mathcal{H}^2 = 0. $$
If we further require that $\xi$ vanishes on the boundary, we have the EL equation
$$\int_B \nabla u \cdot \nabla \xi + 4(u^2-1)u \,\xi dx = 0. $$
Suppose we also have that $u \in H^2(B)$, we have
$$\int_B -\Delta u \, \xi + 4(u^2-1)u \xi dx + \int_{\partial B } \dfrac{\partial u}{\partial n}\xi- Q'(u)\xi d\mathcal{H}^2 = 0. $$
We finally obtain the equation
$$\Delta u = 4(u^2-1)u \,\text{  in } B \,\text{ and }\, \dfrac{\partial u}{\partial n}=Q'(u) \,\text{ on }\, \partial B.$$
My main goal is to prove the minimizer $u_0$ solve the above equation weakly with the desried Neumann boundary condition. However, my question is how to obtain the $H^2$ bound of $u$? I think we can apply standard estimate to obtain $H^2_{loc}$. If we do not have the fact that $u \in H^2(B)$, we may hard to have the existence of $\dfrac{\partial u}{\partial n}$ on the boundary by trace theorem.
 A: Let $u\in H^1(B)$ be a minimising value of $E:H^1(B)\to \mathbb{R}_+$. Then $E'(u)=0$ in the sense that 
$$
E'(u)v = 
\int_B [\nabla u \cdot \nabla v + 4(u^2-1)uv ] dx - \int_{\partial B} Q'(u)v \hspace{.5pc}d\mathcal{H}^2 = 0
$$
for all $v\in H^1(B)$. This implies, in particular, that 
$$
\int_B [-u \Delta v + 4(u^2-1)uv] \hspace{.5pc}dx = 0
$$
for $v\in \mathscr{D}(B)$ (this uses the fact that $v$ and all its derivatives vanish on $\partial B$ if it has compact support in $B$). Having this last equation hold for all $v\in \mathscr{D}(B)$ is exactly the statement that
$$
\Delta u = 4(u^2-1)u \hspace{1pc}\mbox{in $B$  }
$$
in the sense of distributions. The interior equation in the distributional sense is therefore just what you get from regarding $E'(u)$ itself as a distribution (i.e. by restricting it to $\mathscr{D}(B))$. 
For each compact subset $K$ of $B$, take $\psi_K\in \mathscr{D}(B)$ with $\psi_{|K}=1$. Then considering the action of $E'(u)$ on $\{v(1-\psi_K):v\in C^\infty(\bar B),\mbox{$K\subset B$ compact}\}$ should - I think - also give you a weak form of your boundary condition. 
I might be wrong, but I think that it's only once you have your interior equation in a (weak or distributional) form which does a priori involve the assumption that $u\in H^2(B)$ that it's time time to start worrying about regularity. In particular, you may be able to establish that $u\in H^2(B)$ using the fact that it is a $H^1$ solution to your PDE. 
