Let $M = \mathbb{R}^3 \setminus B_1$ where $B_1$ is the unit ball (equip $M$ with the euclidean metric for simplicity, but it will be replaced by an arbitrary asymptotically flat metric). Let $N: L^2(\partial M) \to L^2(\partial M)$ be the Dirichlet to Neumann map defined by $N(f) = \frac{\partial u}{\partial r}$ where $u$ is the unique harmonic extension of $f$ on $M$ going to $0$ at infinity.
Consider the operator $H: H^1(\partial M) \to L^2(\partial M)$ defined by $H(f) = |\nabla f|^2 - N(f)^2 $ where $\nabla$ is the covariant derivative with respect to the round sphere $\partial M$.
What do we know about the solvability of the PDE: $$H(f) = g$$ where $g \in L^2 (\partial M)$? (or any other function space that makes the question easy to answer). How do I approach this problem? Also $H$ should be a pseudo-differential operator. How do I go about finding its order?
What are the different methods/tools for showing an operator like $H$ is surjective?
What if $M$ is equipped with an arbitrary asymptotically flat metric? Does that affect the solvability? It is known that the first and 0^th order of $N$ only depend on the metric and second fundamental form on $\partial M$; that might suggest that that the solvability is independent of the metric on $M$ but only depends on intrinsic and extrensic geometry on $\partial M$.
Any suggestions or references are appreciated.