Pointwise bounds on Dirichlet Green's functions for a singular drift term

1. Problem Statement

The paper investigates the pointwise upper and lower bounds for the Dirichlet Green's function, $G(x, y)$, associated with a specific class of linear second-order elliptic operators in the unit ball $\Omega = B(0, 1) \subset \mathbb{R}^n$ ($n \ge 3$).

The operator is defined as:
$$Lu = -\Delta u + B \cdot \nabla u = 0$$
The core challenge addressed is the presence of a singular drift term $B(X)$ that diverges near the boundary $\partial \Omega$. Specifically, the drift satisfies:
$$\delta(X)^{1-\beta} |B(X)| \le M$$
where $\delta(X) = \text{dist}(X, \partial \Omega)$, $0 < \beta < 1$, and $M < \infty$.

Key Constraints and Conditions:

• Non-Coercivity: The bilinear form associated with the operator is not coercive. Standard energy estimates (which rely on coercivity) fail in this setting, necessitating new techniques for pointwise bounds.
• Divergence Condition: The drift satisfies a pointwise divergence condition:
  $$-\frac{M}{\delta(x)^{2-\beta}} \le \nabla \cdot B \le 0$$
  This implies $\nabla \cdot B$ is negative in the distributional sense, ensuring the existence of solutions to the adjoint equation.
• Regularity: The Green's function is assumed to be locally in $C^{3,\alpha}$ for some $\alpha > 0$.

The paper aims to establish bounds of the form:
$$G(x, y) \asymp \frac{1}{|x-y|^{n-2}}$$
despite the singularity of the drift, which does not belong to standard $L^p$ or Stummel-Kato spaces typically required for such estimates.

2. Methodology

The author introduces a novel technique that combines modified Lorentz norms, level set analysis, and spherical geometry to overcome the lack of coercivity.

A. Lower Bounds (Section 3)

The lower bound proof follows an extension of arguments found in [GW82].

• Test Functions: Uses smooth cut-off functions localized around the pole $y$.
• Energy Estimates: By testing the adjoint equation with $G(\cdot, y)\eta$, the author derives integral inequalities involving $|\nabla G|^2$ and $G^2$.
• Harnack's Inequality: Combines these integral bounds with Harnack's inequality to establish the lower bound:
  $$G(y, z) \ge \frac{K(M, \lambda)}{|z-y|^{n-2}}$$
  for $|z-y| \le \frac{1}{2}\delta(y)$.

B. Upper Bounds (Section 4)

This is the primary technical contribution, involving a sophisticated iterative argument on the level sets of the Green's function.

1. Modified Lorentz Norms: The author defines a functional based on the measure of super-level sets $\Omega_s = \{x : G(x, 0) \ge s\}$. Let $f(s) = |\Omega_s|^{\frac{n-2}{n}}$. The goal is to show $f(s) \le K/s$.
2. Differential Inequalities on Level Sets:
   • The proof analyzes the behavior of $G$ at points of minimum and maximum distance from the origin on the level sets $\partial \Omega_s$.
   • Case (a): Gradient is Large. When the gradient is large, the author uses Taylor expansions and the specific structure of the Laplacian in spherical coordinates to derive differential inequalities.
   • Case (b): Gradient is Small. When the gradient is small, the decay is easier to estimate.
   • Key Equation (62): The proof utilizes the radial form of the Laplacian at the point of tangency between a level set and a sphere centered at the origin:
     $$\frac{\partial^2 G}{\partial r^2} + \frac{n-1}{r}\frac{\partial G}{\partial r} + B \cdot \hat{r} \frac{\partial G}{\partial r} \le 0$$
     This inequality is crucial because it relates the drift term $B$ directly to the radial decay, allowing the author to control the singularity.
3. Contradiction Argument:
   • The author assumes the upper bound fails (i.e., the level sets become too small relative to the value $s$).
   • By constructing sequences of maximum and minimum points converging to the pole, they derive a contradiction with Harnack's inequality. Specifically, the ratio of the Green's function values at these points would violate the uniform bounds required by Harnack's principle if the drift were too singular relative to the geometry.
4. Extension to General Points (Theorem 3):
   • Once the bound is established for a pole at the origin ($G(x, 0)$), the result is extended to an arbitrary pole $y$ using the symmetry $G(x, y) = G^T(y, x)$ (where $G^T$ is the Green's function for the adjoint operator) and repeated applications of Harnack's inequality along a chain of balls connecting $x$ and $y$.

3. Key Contributions

1. First Estimates for Non-Coercive Singular Drifts: This paper provides the first pointwise upper and lower bounds for Green's functions where the drift term diverges like $\delta(X)^{-(1-\beta)}$ with $\beta < 1$. Previous results typically required drifts to be in $L^p$ ($p \ge n$) or Stummel-Kato spaces, which this singular drift does not satisfy.
2. Novel Technique for Singular Terms: The method of using modified Lorentz norms combined with a detailed analysis of the radial behavior of the Green's function on level sets (specifically exploiting the spherical form of the Laplacian) is a new approach for handling singular drifts.
3. Uniformity in Subdomains: The constants in the upper estimates are shown to be uniform in any ball $B(0, r)$ for $r < 1$, with explicit dependence on $r$.
4. Generalization to Morrey Spaces: The results suggest that the drift need only be majorized by a radially integrable function, placing the drift in generalized Morrey spaces $M^{n}_{2-2\beta}$ and $M^{n}_{2-\beta}$, which are more general than standard $L^p$ spaces.

4. Main Results

The paper establishes three main theorems:

• Theorem 1 (Lower Bound): For the operator $L$ with coefficients satisfying the ellipticity and drift conditions, there exists a constant $K(M, \lambda)$ such that for any $z, y \in \Omega$ with $|z-y| \le \frac{1}{2}\delta(y)$:
  $$G(y, z) \ge \frac{K(M, \lambda)}{|z-y|^{n-2}}$$

• Theorem 2 (Upper Bound at Origin): For the operator $L$ with drift satisfying the divergence condition, there exists a constant $K'(M, \lambda)$ such that for any $x \in \Omega$ with $|x| \le 1/2$:
  $$G(x, 0) \le \frac{K'(M, \lambda)}{|x|^{n-2}}$$

• Theorem 3 (General Upper Bound): Under the same conditions, for any $y \in B(0, r)$ ($r < 1$) and $x \in \Omega$ with $|x-y| \le \frac{1}{2}\delta(y)$:
  $$G(x, y) \le \frac{K'(M, \lambda, r)}{|x-y|^{n-2}}$$
  Here, the constant depends on the distance of the subdomain from the boundary.

5. Significance and Future Directions

• Significance: The results fill a critical gap in the theory of elliptic PDEs with singular coefficients. By proving that Green's functions retain the standard fundamental solution decay ($|x-y|^{-(n-2)}$) even in the presence of non-coercive, singular drifts, the paper validates the robustness of these estimates. This is particularly important for applications in stochastic processes (where drifts represent singular forces) and mathematical physics.
• Future Directions:
  • General Domains: Extending the results from the unit ball to general bounded Lipschitz domains.
  • Schrödinger Operators: Applying the technique to time-independent Schrödinger equations with singular potentials.
  • Parabolic Equations: Investigating Green's functions for the corresponding parabolic operators ($\partial_t u + B \cdot \nabla u - \Delta u = 0$).
  • Principal Term Perturbations: Studying operators where the principal part is a general elliptic operator $-\nabla \cdot (A \nabla u)$ rather than the Laplacian, though the author notes the current proof relies heavily on the spherical symmetry of the Laplacian.

In summary, Aritro Pathak's work introduces a powerful analytical framework to handle singular drifts that were previously intractable using standard energy methods, establishing that the fundamental singularity structure of the Green's function remains intact under these challenging conditions.