Let $ A(x)=(a_{ij}(x)):\mathbb{R}^d\to\mathbb{R}^{d\times d} $ be a matrix satisfying ellipticity condition, that is
\begin{align}
\sum_{i,j=1}^da_{ij}(x)\xi_i\xi_j\geq\mu|\xi|^2
\end{align}\begin{align}
\mu^{-1}|\xi|^2\geq \sum_{i,j=1}^da_{ij}(x)\xi_i\xi_j\geq\mu|\xi|^2
\end{align}
with $ \mu>0 $ being a positive constant and $ \xi\in\mathbb{R}^d $. Since the spectrum of $ L=-\operatorname{div}(A(x)\nabla) $, for $ \zeta\in\mathbb{C}\backslash\mathbb{R}_+ $, $ (L-\zeta I)^{-1} $ exists. I read some paper and know that for some special $ A $ and $ p $, we can derive that
\begin{align}
\left\|(L-\zeta I)^{-1}\right\|_{L^p\to L^p}\leq \frac{C}{1+|\zeta|}.
\end{align}
where $ C $ may depends on $ \arg\zeta $. I want to know what is the motivation to study such estimates and where can I use such estimates.