I asked the following question on math.SE (<https://math.stackexchange.com/questions/2420298/bvps-for-elliptic-pdos-when-do-green-functions-l2-inverses-define-pseudo-d>) just over two months ago, and it only received one, rather unsatisfactory to me, answer there. I'm wondering if people here can have a look. Some related questions have been posed [here](https://mathoverflow.net/questions/188815/greens-operator-of-elliptic-differential-operator) and [here](https://mathoverflow.net/questions/69521/estimates-on-the-green-function-of-an-elliptic-second-order-differential-operato) (and I noted in particular the literature link provided in the answer to the latter question). However, neither of those discussions cover the possibility of boundary conditions being present, or possible lack of compactness.


Repost of the original question
===============================


Let me illustrate my question by starting with the simplest possible example: Let us consider $P := - \mathrm{d}^2/\mathrm{d}x^2$, an elliptic partial differential operator on $\mathbb{R}$; let us also consider the following boundary-value problem on the interval $\overline{\Omega} = [0,1]$:
\begin{equation}
P u = f, \qquad u(0)=u(1)=0.
\end{equation}
As is (I think) well-known, when seen as an operator $L^2(\Omega) \to L^2(\Omega)$, $P$ is unbounded. However, it is closed on the dense domain $D(P) := H^2(\Omega) \cap H\_0^1(\Omega)$ where $H\_0^1(\Omega)$ is the closure of $C\_{\mathrm{c}}^\infty(\Omega)$ in the $H^1$ norm (so that any element of this space has vanishing trace on $\partial \Omega = \{0,1\}$, i.e. it satisfies the Dirichlet boundary condition above in a weak sense). Furthermore, $0$ is in the resolvent of $(P,D(P))$, i.e. there exists a bounded inverse $P^{-1} : L^2(\Omega) \to L^2(\Omega)$. In fact, in this example the inverse is easily computed: it is the integral operator defined by the (continuous, as it happens) kernel
\begin{equation}
G(x,y) = \begin{cases} x(1-y) & x \leq y \\ y(1-x) &x > y \end{cases}, \quad (x,y) \in \Omega \times \Omega.
\end{equation}
Of course, when viewed as a distribution in $\mathscr{D}'(\Omega \times \Omega)$, $G$ is the Schwartz kernel of $P^{-1}$ which we know on abstract grounds must exist since $P^{-1} : C\_{\mathrm{c}}^{\infty}(\Omega) \to \mathscr{D}'(\Omega)$ is continuous.


My question is the following: in this example and in more general examples where $P$ is a second-order elliptic differential operator on, say, an open (and not necessarily compact) region $\Omega$ with smooth boundary in $\mathbb{R}^n$, and assuming that we can find a suitable dense domain $D(P)$ for $P$ as above so that $(P,D(P))$ has a bounded inverse $P^{-1} : L^2(\Omega) \to L^2(\Omega)$, **does the Schwartz kernel $G$ of $P^{-1}$ always define a pseudodifferential operator on $\Omega$?**


Addendum for MO
===============


User mcd on math.SE points out that the Boutet de Monvel calculus ought to be relevant here. Aside from wishing to see exactly how this is, I wonder whether the possible lack of compactness (of $\Omega$) might cause problems in such an approach.


UPDATE
======


I have reduced my question to the following subproblem: [Composition of a smoothing operator with an $L^2$-bounded operator, non-compact Riemannian manifold](https://mathoverflow.net/questions/287167/composition-of-a-smoothing-operator-with-an-l2-bounded-operator-non-compact), as explained in a comment there.



BVPs for elliptic PDOs: When do Green functions ($L^2$ inverses) define pseudo-differential operators in the interior?
By complementing [Deane Yang](https://mathoverflow.net/users/613/deane-yang)'s strategy with [Jochen Wengenroth](https://mathoverflow.net/users/21051/jochen-wengenroth)'s observations in the related question [Composition of a smoothing operator with an $L^2$-bounded operator, non-compact Riemannian manifold](https://mathoverflow.net/questions/287167/composition-of-a-smoothing-operator-with-an-l2-bounded-operator-non-compact), I can now write a fully detailed answer to my question.


**Theorem.** Let $M$ be a smooth manifold equipped with a smooth measure $\mathrm{d}\mu$ (arising from example from a metric on $M$). Let also $P$ be an elliptic partial differential operator on $M$, and suppose there exists a domain $D(P)$ such that $C\_\mathrm{c}^\infty(M) \subseteq D(P) \subseteq L^2(M, \mathrm{d} \mu)$ and $P|\_{D(P)} : D(P) \to L^2(M, \mathrm{d} \mu)$ has a bounded left inverse $P^{-1}$. Then the restriction of $P^{-1}$ to $C\_\mathrm{c}^\infty(M)$ defines a pseudo-differential operator.


*Proof*. In two steps:


1. $P^{-1}$ maps $C\_\mathrm{c}^\infty(M)$ to $C^\infty(M)$ by elliptic regularity. Now we use a Lemma, whose proof is given below:


**Lemma.** Let $T : C\_\mathrm{c}^\infty(M) \to \mathscr{D}'(M)$ be linear, continuous relative to the standard inductive limit topology of $C\_\mathrm{c}^\infty(M)$ and the weak topology of $\mathscr{D}'(M)$, and with image contained in $C^\infty(M)$. Then $T : C\_\mathrm{c}^\infty(M) \to C^\infty(M)$ is continuous relative to the standard Fréchet space topology of the codomain.


Now, the composition $C\_\mathrm{c}^\infty(M) \hookrightarrow L^2(M, \mathrm{d} \mu) \xrightarrow{P^{-1}} L^2(M, \mathrm{d} \mu) \hookrightarrow \mathscr{D}'(M)$ is continuous and has image contained in $C^\infty(M)$; hence, by the Lemma, $P^{-1}$ maps $C\_\mathrm{c}^\infty(M)$ continuously to the Fréchet space $C^\infty(M)$.
2. Also since $P$ is elliptic, there is a pair $(Q,R)$ of properly supported pseudo-differential operators such that $R$ is smoothing and $PQ = I + R$. This equality is valid when both sides are applied to $\mathscr{D}'(M)$, and when we restrict it to $C\_\mathrm{c}^\infty(M)$ we can compose both sides with $P^{-1}$ to obtain
$$(\*) \qquad P^{-1} = Q - P^{-1}R \quad \text{on } C\_\mathrm{c}^\infty(M),$$
and in fact the equality is clearly also valid on the larger set $Q^{-1}[D(P)]$. Now, since $R$ is properly supported and smoothing, it is continuous from $\mathscr{E}'(M)$ to $C\_\mathrm{c}^\infty(M)$, and hence $\mathscr{E}'(M) \xrightarrow{R} C\_\mathrm{c}^\infty(M) \xrightarrow{P^{-1}} C^\infty(M)$ is continuous by part 1 of the proof. It follows from known results enunciated in the [related MO question](https://mathoverflow.net/questions/287167/composition-of-a-smoothing-operator-with-an-l2-bounded-operator-non-compact) that the Schwartz kernel of $P^{-1}R$ is smooth, i.e. that $P^{-1}R$ is a smoothing integral operator. Hence, the right-hand side of Equation (\*) is a pseudo-differential operator extending $P^{-1}$. This completes the proof.


*Proof of the Lemma.* $T$ is continuous as a map $C\_\mathrm{c}^\infty(M) \to C^\infty(M)$ when the codomain is equipped with the subspace topology from $\mathscr{D}'(M)$. Hence its graph is closed in the appropriate product topology, and it is also closed in the finer product topology obtained by equipping $C^\infty(M)$ with its Fréchet space topology. The result then follows from a version of the closed graph theorem applicable to linear maps from LF-spaces to Fréchet spaces.


