Resolvent estimates for the Schrödinger operator with $L^\infty$ electric and magnetic potentials and applications to the local energy decay

1. Problem Statement and Context

The paper addresses the problem of local energy decay for solutions to the wave equation associated with the magnetic Schrödinger operator:
$$P = (i\nabla + b(x))^2 + V(x)$$
defined on a domain $\Omega = \mathbb{R}^d \setminus \mathcal{O}$, where $\mathcal{O}$ is a bounded obstacle (possibly empty). The potentials $V$ (electric) and $b$ (magnetic) are assumed to be in $L^\infty$.

The Core Challenge:
Previous work (specifically reference [13]) established exponential local energy decay for potentials decaying exponentially at infinity (e.g., $|V(x)| + |b(x)| \leq Ce^{-c|x|}$). This relied on the analytic continuation of the weighted resolvent. However, many physical potentials decay slower than exponential (e.g., sub-exponential decay like $e^{-|x|^s}$ with $0 < s < 1$). The authors aim to:

1. Extend resolvent estimates to a larger class of potentials that decay sub-exponentially but faster than any polynomial.
2. Derive the corresponding decay rates for local energy for these slower-decaying potentials.
3. Handle two specific geometric settings:
   • Case (a): $\Omega = \mathbb{R}^d$ ($d \geq 3$) with a non-trivial magnetic potential ($b \not\equiv 0$).
   • Case (b): $\Omega = \mathbb{R}^d \setminus \mathcal{O}$ with a non-trapping obstacle $\mathcal{O}$, but zero magnetic potential ($b \equiv 0$).

2. Methodology

The authors employ a multi-step strategy combining spectral theory, complex analysis, and microlocal analysis techniques.

A. Generalized Weight Functions and Gevrey Classes

Instead of exponential weights $e^{-c|x|}$, the authors introduce a general weight function $\mu(x) = \sqrt{\Theta(c\langle x \rangle)}$, where $\Theta$ satisfies specific decay and regularity conditions:

• Decay: $\Theta(r) = O((r+1)^{-k})$ for all integers $k$ (super-polynomial decay).
• Regularity: Derivatives of $\Theta$ are controlled by $\Theta$ itself.
• Convexity-like condition: $\Theta(r_1)\Theta(r_2) \leq C_1 \Theta(C_2(r_1+r_2))$.
• Sequence Control: The decay rate is characterized by a sequence $m_k$ (related to Gevrey classes). For example, if $\Theta(r) = e^{-(r+1)^s}$ with $0 < s < 1$, then $m_k \sim (k!)^{1/s}$.

B. Resolvent Estimates

The central technical achievement is proving uniform estimates for the derivatives of the weighted resolvent $\mu(P - \lambda^2)^{-1}\mu$ with respect to the spectral parameter $\lambda$.

• Free Resolvent (Section 2): They establish bounds for the free Laplacian $P_0 = -\Delta$. Using the representation of the wave propagator via the Fourier transform and Huygens' principle (for odd dimensions) or analytic properties of the kernel (for even dimensions), they prove that the $k$-th derivative of the weighted free resolvent grows like $C^{k+1} m_k$.
• Perturbed Resolvent (Section 4): Using resolvent identities (specifically the second resolvent equation), they relate the perturbed resolvent $(P - \lambda^2)^{-1}$ to the free one.
  • They utilize Lemma 3.1 and 3.2 to show that the class of operator-valued functions satisfying the $C^{k+1}m_k$ bound is stable under products and inverses.
  • This allows them to transfer the estimates from the free operator to the perturbed operator $P$, even in the presence of the magnetic potential (Case a) or the obstacle (Case b).

C. Limiting Absorption Principle

The authors rely on the limiting absorption principle to ensure the resolvent exists on the real axis (excluding zero). They provide a concise proof in Appendix B showing that for short-range potentials, the weighted resolvent extends continuously from the lower half-plane to $\mathbb{R} \setminus \{0\}$.

D. Time Decay via Tauberian Arguments

In Section 5, they convert the resolvent estimates into time decay rates.

• They construct a Gevrey cutoff function $\zeta$ (Appendix A) compatible with the sequence $m_k$.
• By analyzing the solution formula $u(t) = \cos(t\sqrt{P})f_1 + P^{-1/2}\sin(t\sqrt{P})f_2$ and using the resolvent bounds, they derive an integral inequality for the local energy.
• By optimizing the order of differentiation $k$ relative to time $t$ (specifically choosing $k$ such that $m_k^{1/k} \sim t$), they extract the decay rate.

3. Key Contributions and Results

Main Theorem 1.1 (Resolvent Estimates)

Under the assumption that the potentials satisfy $|V(x)| + |b(x)| \leq C \Theta(c\langle x \rangle)$ with $\Theta$ satisfying the Gevrey-type conditions, the authors prove:
$$\left\| \frac{d^k}{d\lambda^k} \mu \nabla^\ell (P - \lambda^2)^{-1} \mu \right\| \leq C^{k+1} m_k (|\lambda|+1)^{\ell-1}$$
for all $k \geq 0$ and $\lambda \in \mathbb{R}$ (with minor caveats for low frequencies in even dimensions).

• Significance: This generalizes previous results by allowing potentials that decay sub-exponentially (e.g., $e^{-|x|^s}$ for $0 < s < 1$) rather than strictly exponentially. The bounds are uniform in the spectral parameter and the order of derivation.

Main Theorem 1.2 (Local Energy Decay)

As a consequence of the resolvent estimates, the local energy of the wave equation solution decays at a specific rate determined by the decay of the potential.

• If $|V(x)| + |b(x)| \leq C e^{-c|x|^s}$ with $0 < s < 1$, the local energy decays as:
  $$E_{loc}(t) \leq C_0 e^{-c_0 t^s}$$
• Significance: This establishes a precise link between the decay rate of the potential and the decay rate of the energy.
  • Exponential potential ($s=1$) $\to$ Exponential energy decay.
  • Sub-exponential potential ($0 < s < 1$) $\to$ Sub-exponential (stretched exponential) energy decay.
  • The result holds for both the whole space with magnetic fields and exterior domains with obstacles (provided $b=0$ in the latter).

4. Significance and Impact

1. Broadening the Scope of Potentials: The paper moves beyond the restrictive assumption of exponential decay. It demonstrates that even potentials decaying slower than exponential (but faster than polynomial) still yield strong local energy decay, albeit at a slower, stretched-exponential rate.
2. Unified Framework for Derivatives: The introduction of the sequence $m_k$ and the Gevrey-class framework provides a robust tool for handling derivatives of operator-valued functions, which is crucial for high-frequency analysis and resolvent estimates.
3. Handling Magnetic Fields and Obstacles: The work successfully extends these estimates to the magnetic Schrödinger operator in $\mathbb{R}^d$ ($d \geq 3$) and to the Dirichlet realization in exterior domains (with $b=0$). The separation of cases highlights the technical difficulties introduced by magnetic potentials in the presence of obstacles.
4. Sharpness of Decay Rates: The results confirm that the decay rate of the local energy is intrinsically linked to the decay rate of the coefficients. This resolves questions about whether "slow" decay of potentials leads to "slow" decay of energy, providing a quantitative answer ($e^{-t^s}$).

5. Conclusion

This paper represents a significant advancement in the spectral theory of Schrödinger operators. By developing a refined resolvent estimate technique based on Gevrey-class weights, the authors successfully characterize the local energy decay for a wide class of sub-exponentially decaying potentials. The results bridge the gap between the classical exponential decay regime and polynomial decay regimes, offering a precise mathematical description of how the "tail" of the potential influences the long-time behavior of wave propagation.