Let $L = \sum\_{i,j=1}^n -\frac{\partial}{\partial x^i} (a^{ij}(x)\frac{\partial}{\partial x^j}) + \sum\_{i=1}^n b^i(x) \frac{\partial}{\partial x^i} + c(x)$ be a second order elliptic operator with smooth coefficients, $\Omega$ a bounded open domain with smooth boundary in $\mathbb{R}^n$, and $f$ be a function in $L^p(\Omega)$. We say that $u \in W\_0^{1,p}(\Omega)$ (one weak derivative in $L\_p$ and vanishing boundary values) is a weak solution of $Lu = f$ if for all $g \in W\_0^{1,q}(\Omega)$ $(q = p^\*)$ we have 


$\int\_{\Omega} \sum\_{i,j=1}^n a^{ij}(x)\frac{\partial u}{\partial x^i}\frac{\partial g}{\partial x^j} + \sum\_{i=1}^n b^i(x)\frac{\partial u}{\partial x^i}g(x) + c(x)u(x)g(x) = f(x)$.


The standard result is of course that all such weak solutions $u$ actually belong to $W\_0^{2,p}(\Omega)$.


I am trying to complete the following proof of this statement:


(1) First we establish an a priori estimate for strong solutions $v \in W\_0^{2,p}(\Omega)$ of $Lv = f$:


$\vert\vert v\vert \vert\_{W\_0^{2,p}(\Omega)} \leq C(\vert\vert f\vert\vert\_{L^p(\Omega)} + \vert\vert v\vert\vert\_{L^p(\Omega)})$


This is non-trivial but can be established by proving the relevant estimate for the Laplacian with a Newton Potential argument and then using the freezing coefficients technique.


(2) Next we observe that if $L$ is injective on $W\_0^{1,p}$, then we are done. This is because $L$ injective implies


$\vert\vert v\vert\vert\_{L^p(\Omega)} \leq C\vert\vert Lv\vert\vert\_{L^p(\Omega)}$


One proves this by assuming it was false and then using Rellich compactness to produce a non-zero solution to $Lv = 0$.


Having established this estimate, we consider the smooth mollifications $f\_{\epsilon}$ of $f$. By $L^2$ theory we can find smooth $v\_{\epsilon}$ strong solutions of $Lv\_{\epsilon} = f\_{\epsilon}$. Since we have


$\vert\vert v\_{\epsilon} - v\_{\epsilon'}\vert\vert\_{W\_0^{2,p}} \leq C\vert\vert f\_{\epsilon}-f\_{\epsilon'}\vert\vert\_{L^p(\Omega)}$


The $v\_{\epsilon}$ converge to some $v \in W\_0^{2,p}(\Omega)$ which will solve $Lv = f$ strongly. Since strong solutions are clearly weak solutions, by the injectivity of $L$ on $W\_0^{1,p}(\Omega)$ we conclude that $u = v \in W\_0^{2,p}(\Omega)$ and we are done. 


(3) We have no way to guarantee that $L$ is injective, for example $0$ might be an $L^2$ eigenvalue of $L$. However, if $p=2$ then we could guarantee that $L\_{\lambda} = L + \lambda I$ is injective for some large $\lambda$. If we could establish this fact in the general case we would be done since $L\_{\lambda}u = f + \lambda u \in L^p$ and $L\_{\lambda}$ injective imply that (2) is applicable.


**Question**: What is the simplest way to prove that $L\_{\lambda} = L + \lambda I$ is injective on $W\_0^{1,p}(\Omega)$ for large $\lambda$? Do there exist **weak** $L^P$ maximum principles?


Of course, I would prefer that the proof not use $L^p$ regularity. 



Proof of L^p Elliptic Regularity
I stumbled across the book *Second Oder Elliptic Equations and Elliptic Systems* by Yah-Ze Chen which appears to contain an answer to my question. It is available on google books here 


<http://books.google.com/books?id=eQcbiPQPweQC&pg=PA49&dq=strong+solution+dirichlet+problem+Lp&hl=en&ei=byKiTeXXO6GG0QGG7dGgBQ&sa=X&oi=book_result&ct=result&resnum=1&ved=0CCcQ6AEwADgK#v=onepage&q=strong%20solution%20dirichlet%20problem%20Lp&f=false>


To save time for those who are interested, here is the relevant argument:


For large $\lambda > 0$ we want to show that $L\_{\lambda} = L - \lambda I$ is injective on $W\_0^{1,p}$.


Claim: Let $L^T\_{\lambda}$ be the transpose of $L^{\lambda}$ with respect to the paring that defines weak solutions. Then we claim that $L^T\_{\lambda}$ inective on $W\_0^{2,p}$ implies that $L\_{\lambda}$ is injective on $W\_0^{1,p}$


Proof: Suppose that $L^T\_{\lambda}$ is injective on $W\_0^{2,p}$. Then, by an argument contained in the original post above, for every $f \in L^p(\Omega)$ we can find $u \in W\_0^{2,p}(\Omega)$ such that $L^T\_{\lambda}u = f$. Now, suppose that $L\_{\lambda}v = 0$ for some $v \in W\_0^{1,p}$. After an integration by parts and the definition of weak solution, we see that $\varphi \in W\_0^{2,q}$ implies that 


$\int\_{\Omega}uL^T\_{\lambda}\varphi = 0$. 


Now choose $\Omega'' \subset\subset \Omega' \subset\subset \Omega$ and a bump function $\rho$ identically one in $\Omega''$ with support in $\Omega'$. $\rho\text{sgn}(u)$ is in $L^q$, and we can find $g \in W\_0^{2,q} $ such that $L^T\_{\lambda}g = \rho\text{sgn}(u)$. Plugging this $g$ into the above equality gives


$\int\_{\Omega''}|u| = -\int\_{\Omega\setminus\Omega''}\rho |u|$


Due to the arbitrariness of $\Omega''$, this implies that $\int\_{\Omega} |u| = 0$ and hence $u$ is $0$ a.e.


Claim: For $\lambda$ large enough, $L\_{\lambda}$ is injective on $W\_0^{2,p}$. 


Proof: Suppose $L\_{\lambda}u = 0$ for $u \in W\_0^{2,p}$. Let $\tilde{\Omega} = \Omega \times (-1,1)$, and $\tilde{\Omega'} = \Omega \times (-1/2,1/2)$. Let $(x,t)$ be the coordinates on $\Omega \times (-1,1)$. Then define $v(x,t) = \cos(\sqrt{\lambda}t)u(x)$. Let $\hat{L\_{\lambda}} = L\_{\lambda} + \partial\_t^2$. We have $\hat{L\_{\lambda}}v = 0$. The strong solution estimates give


$\vert\vert v\vert\vert\_{W^{2,p}(\tilde{\Omega'})} \leq C\vert\vert v\vert\vert\_{L^p(\tilde{\Omega})} \leq C\vert\vert u\vert\vert\_{L^p(\Omega)} \Rightarrow $


$\vert\vert \partial\_t^2v\vert\vert\_{L^p(\tilde{\Omega'})} \leq C\vert\vert u\vert\vert\_{L^p(\Omega)} \Rightarrow $


$\lambda\vert\vert u\vert\vert\_{L^p(\Omega)}(\int\_{-1/2}^{1/2}|\cos(\sqrt{\lambda}t|^p)^{1/p} \leq C\vert\vert u \vert\vert\_{L^p(\Omega)} \Rightarrow$


$\lambda^{1 - \frac{1}{2p}}\vert\vert u\vert\vert\_{L^p(\Omega)}(\int\_{-1/2}^{1/2}|\cos t|^p)^{1/p} \leq C\vert\vert u\vert\vert\_{L^p(\Omega)}$


Now taking $\lambda$ large enough implies that $u = 0$ almost everywhere.


