Let $\Delta$ denote the Laplace operator in $\mathbb R^2$, and let $u_1$ and $u_2$ be two orthogonal eigenfunction of $-\Delta$, i.e., $-\Delta u_1=\alpha_1 u_1$ and $-\Delta u_2=\alpha_2 u_2$, where $\alpha_1$ and $\alpha_2$ are the corresponding eigenvalue. I am wondering do we have $$ \Delta^{-1}(u_1+u_2)=\Delta^{-1}u_1+\Delta^{-1}u_2\tag 1 $$ and if yes, how we proof it?
I know $(1)$ does not hold for general function, but I think it might work for orthogonal functions... PS: by orthogonal I mean $$ \int u_1 u_2dx=0. $$ PPS: by $\Delta^{-1}$ I mean the inverse of Laplace operator, i.e., $\Delta^{-1}\Delta=I$
Note that if we have $\Delta \phi(x,y)=u(x,y)$, then
$$\phi(x,y)=\int_{\mathbb{R}^2} u(x',y')G(x,y|x',y')\,dx'\,dy' \tag 1$$
where $G(x,y|x',y')=\frac{1}{2\pi}\log\left(\sqrt{(x-x')^2+(y-y')^2}\right)$ is the "free-space" Green (or Green's) Function (SEE HERE).
Therefore, the inverse operator of the Laplacian is a linear operator and we assert that $\Delta^{-1}(u_1+u_2)=\Delta^{-1}(u_1)+\Delta^{-1}(u_2)$.