Neumann gradient estimate for nonlinear heat equation under integral Ricci curvature bounds

: In this paper, we consider a Li-Yau gradient estimate on the positive solution to the following nonlinear parabolic equation


Introduction
We investigate the Li-Yau gradient estimate for positive solutions to the following nonlinear parabolic equation: This equation is subject to Neumann boundary conditions on a compact Riemannian manifold that fulfills the integral Ricci curvature assumption.Here, p ≥ 0 is a real constant.It is noteworthy that, when a = 0, the aforementioned equation aligns with the universally recognized heat equation.Empirical studies have underscored the extensive applicability of scrutinizing solutions to parabolic equations on a Riemannian manifold in various fields, encompassing biology, physics, dynamical systems, and image processing.Furthermore, it serves as a crucial instrument in the domain of geometric analysis.In the pursuit of understanding such correlated issues, it becomes essential to contemplate initial estimates of the solution.In this regard, the Li-Yau type gradient estimation emerges as a significant subject.
The exploration of Li-Yau gradient estimates for the heat equation was initiated in the seminal paper by Li and Yau [1].In this work, they derived gradient estimates for positive solutions of the heat equation on a complete Riemannian manifold with bounded Ricci curvature, employing the maximum principle.Subsequently, Hamilton utilized a similar technique to establish Harnack inequalities for the Ricci flow.Moreover, Hamilton extended the Li-Yau gradient estimates of the heat equation to a full matrix version in [2].From this point forward, differential Harnack inequalities have played a pivotal role in the research of geometric flows.In 2002, Perelman [3] applied this method to formulate the Harnack inequality for the fundamental solution of the conjugate heat equation.Leveraging the Harnack inequality, Perelman demonstrated the non-collapse theorem for the Ricci flow.This ground breaking achievement laid the essential foundation for proving the feasibility of the Poincaré conjecture.Numerous papers have been published on the study of parabolic equations based on the Ricci curvature condition, including works by Hamilton [2], Li-Xu [4], and others [5][6][7][8][9][10].Additionally, the authors [11,12] examined a heat equation problem within the framework of control theory.
It is a logical approach to investigate the gradient estimation of positive solutions to equations under integral curvature conditions.It is worth mentioning that this research direction has seen substantial contributions.Prior to presenting the results, we establish some notations.For each x∈M n , let g(x) represent the smallest eigenvalue for the Ricci tensor Ric: T x M → T x M. For any given number K, the amount of n-dimensional Ricci curvature lying below (n − 1)K.Now consider, for any constants q, r > 0, Note that k(q, r) = 0 if and only if Ric ≥ (n − 1)K.It is crucial to recognize that the integral curvature bound is a natural condition, which is substantially weaker than a lower bound Ricci curvature condition.Consequently, the relationship between the integral Ricci curvature condition and the topology and geometry of manifolds necessitates comprehensive examination.Specifically, by formulating a global version of gradient estimates, Wang [13] achieved the Li-Yau type gradient estimates for the positive solution of the heat equation with Neumann boundary condition.In this scenario, the boundary of M is not convex, but it satisfies the interior rolling R-ball condition, with constants It should be noted that our research also encompasses the second fundamental form of ∂M and the interior rolling R-ball condition.The interior rolling ball condition was employed by Chen in [14] to derive an estimate for the first nonzero Neumann eigenvalue of M.
The central focus in gradient estimation under integral curvature conditions, tracing back to the paper by Petersen and Wei [15], revolves around an extension of the Bishop-Gromov relative volume comparison estimate to scenarios where there is only an integral bound for the portion of the Ricci curvature that falls below a specified number.Utilizing Yang's estimates on the Sobolev constants from [16] and capitalizing on the relative volume comparison for integral Ricci curvature, Petersen and Wei [17] broadened several geometric results for manifolds with lower Ricci curvature bounds to integral lower bounds.Building on this groundwork, Zhang and Zhu [18] proposed a Li-Yau gradient bound for positive solutions of the heat equation on compact manifolds.This result was predicated on the assumption that either |Ric − | ∈ L p (M) for some p > n 2 and the manifold is noncollapsed, or a certain Kato type of norm of |Ric − | is finite and the heat kernel possesses a Gaussian upper bound.Subsequently, Zhang and Zhu [19] expanded the result in [18] with the aid of the Sobolev inequality demonstrated by Dai, Wei, and Zhang [20].
In a similar vein, the application of analogous methods to the study of general equations has yielded significant results and findings.The utilization of such methodologies facilitates valuable insights and advancements in our comprehension of diverse techniques.Wang [21] extended the gradient estimates for the heat equation where a is a constant.These outcomes represented substantial generalizations of the Li-Yau [1], Hamilton [2], and Li-Xu [4] type gradient estimates under the integral Ricci curvature bounds.Wang [22] and Wang [23] conducted research on the elliptic gradient estimates on a Riemannian manifold with control on integral Ricci curvature.Wang and Wang [24] concentrated on investigating a class of nonlinear elliptic equations on a collapsed complete Riemannian manifold and its parabolic counterpart under integral curvature conditions.
In this context, for the scenario where the boundary of M is not convex, Olivé [25] formulated a Li-Yau type gradient estimate for positive solutions to the heat equation with Neumann boundary conditions under the integral Ricci curvature assumption.To be more precise, the principal theorem proposed by Olivé was: , and R > 0 small enough, there exists K(n, p) > 0 such that if M n is a compact Riemannian submanifold with boundary of a Riemannian manifold N n with the properties: (2) II ≥ −H, where II is the second fundamental form of ∂M; (3) M satisfies the interior rolling R-ball condition, then any positive solution u(x, t) to: satisfies the Li-Yau type gradient estimate: where, given any 0 < ξ < 1, we can choose any 0 < α < 1−ξ (1+H) 2 and any 0 < β ≤ ξ 2 (1−ξ) 2ξ 2 +n 2 (1+H) 2 , and where .
Our objective is to derive a gradient estimate for the positive solutions to the nonlinear heat equation with Neumann conditions, under the integral Ricci curvature assumption.The methodology utilized to establish the primary results fundamentally draws upon Olivé [25].However, in contrast to the previously mentioned points, our preliminary task entails identifying a methodology that can proficiently manage the nonlinear terms.Moreover, the function Q(t) can be inferred through analysis, leveraging the Sobolev inequality [20] and volume doubling results [17].
The remainder of this paper is structured as follows.In Section 2, we revisit some definitions and prior works that will serve a fundamental role.Our principal results are presented in Section 3.

Preliminaries
In this section, we will effectively derive and present the definitions and lemmas that were instrumental in shaping the arguments and conclusions of this paper.Initially, let τ(s) denote a nonnegative function on [0, ∞) and the properties of τ(s) as For convenience, τ( D(s) R ) is written as σ(s), where D(s) denotes the distance from a point s ∈ M to the boundary of M.Then, we introduce two functions, namely ω(s) and ω(s), as proposed in references [13,14].The specific expressions of two functions are defined as where b satisfies the conditions outlined in the following lemma.Lemma 2.1.( [13,14]) The function ω satisfies (2.1) Let ∂ ∂ν be the outward pointing unit normal vector to ∂M.To estimate a unique smooth solution Q(x, t), we need the following lemma provided in [25] and the key tools are referenced in [26][27][28].Lemma 2.2.( [25]) There exists a unique smooth solution Q(x, t) to the problem where c > 1 is a constant, and satisfies and r) be a geodesic ball centered at x with raduis r.For any y ∈ ∂N, there exists x ∈ N such that O(x, r) ⊆ N, and O(x, r) ∩ ∂N = y.Then, this condition is commonly referred to as the interior rolling r-ball condition on a Riemannian manifold N with boundary.
Based on the above important definitions and lemmas, for q > n 2 and n > 0, we need to point out that this paper investigates manifold M n in light of the scaling invariant curvature condition where K = K(n, q) > 0 small enough satisfies the volume doubling proprty [17] and diam(M) ≤ D.

Results
The objective of this section is to establish estimates on the positive solution of the equation and we will give the proof of the main theorem.In order to use the maximum principle, by substituting u = ln f into (3.1),we can transform equation (3.1).
Theorem 3.1.Let M n be a compact Riemannian submanifold with boundary of a Riemannian manifold N n and K, H, and r small enough be nonnegative constants.Given H > 0, n > 0, p > n 2 , there exists K(n, p) > 0. The second fundamental form denoted as II and the boundary of the manifold ∂M is bounded from below by −H.Suppose that M satisfies the interior rolling r-ball conditon.If u(x, t) is a positive solution of equation (3.1) on M × (0, ∞), then the Li-Yau type gradient estimate is given by where Proof.Before we give our proof of main results in this section, for simplicity of notation, we define where u = ln f , ε > 0, and Q satisfies the lemma corresponding to c = (3 For T > 0, we suppose that P takes the maximum at point (x, t 0 ) in M × [0, T ], thus ∆P − ∂ t P ≤ 0 at point (x, t 0 ).The difference from [25] is that we need to deal with the terms involving the coefficient of au p .Now we deal with the case that, for x ∈ • M = M\∂M, P takes the maximum at point (x, t 0 ).Without loss of generality, we can assume that P(x, t 0 ) > 0, then where h = |∇u| 2 + ε.
Obviously, the first step of this section should calculate the evolution equation.Computing the evolution equation of P and using Bochner's formula, the following result is straightforward: and it is obvious that (3.9) Towards the same direction, for the term of −2|∇u||∇Q| on the right-hand side, we need to find the upper bound.Taking −2|∇u||∇Q| ≥ − |∇Q| 2 dQ − dQ|∇u| 2 and the Cauchy-Schwarz inequality into consideration, the derivation of (3.9) is as follows: By employing a similar calculation as previously described, we can determine the upper bound of 2∇ ω∇Q, where 2∇ ω∇Q ≥ − 1 bd ω|∇Q| 2  Q − dQ|∇ ω| 2 .We arrive at Subsequently, from (3.11) we obtain where Q satisfies Lemma 2.2 corresponding to c = (3 + 1 b ) 1 d .Since G achieves its maximum value at the point (x, t 0 ) in M n × (0, T ], then The following identities are widely recognized and will be utilized in (3.13): and Therefore, this can be demonstrated by rearranging the inequality as follows: Here, we expand the terms containing h = |∇u| 2 + ε in (3.14) and note the merging of the terms involving the coefficient of a f p .Since we have the inequality 2∇ ω∇u ≥ −(|∇ ω| 2 + |∇u| 2 ), we check where θ(ε) denotes a function that goes to zero as ε goes to zero.
To apply the maximum principle, we further simplify the right-hand side of the inequality.It is necessary to compute the relationship between the term (|∇u| 2 − ∂ t u) 2 and P2 .That is to say that By substituting ∂ t u = ωQ(|∇u| + ε) + au p − P, we derive the transformation of (3.16) to be The negativity of 2(1 − ωQ)|∇u| 2 can lead to the inequality Combining the above two identities, (3.17) and (3.18), we deduce, at point (x, t 0 ), that In the following step, we appropriately select values for b, d and estimate each of the coefficient terms.We choose b ≤ 1−ξ (1−H) 2 for some 0 < ξ < 1 and choose d > 0 so that 2(1−d) ω−b n 2 > 0. The coefficient of |∇u| 4 in the inequality derived from (3.19) can be further refined.For some G > 0, we have Given the evolution of (3.21), the bound of d is given by (3.22) In (3.22), we assume that d is a positive and consider choosing G = bξ 2 n 2 .Then, (3.23) As Lemma 2.1 holds for ω, it allows that the coefficients of |∇u| 2 satisfies where
Since the calculations are similar to that of Wang [13] and Olivé [25], we will only sketch the key steps.Calling x = |∇u| 2 , the coefficient of t 0 J is transformed such that (3.30) By substituting the result of (3.30) into (3.29),we derive that (3.29) becomes Multiplying both sides of (3.31) by t 0 gives rise to Obviously, the right side of (3.32) is quadratic in P. So, if it is nonpositive, we need, at point (p, t 0 ), where P takes the maximum at point (p, t For any point x ∈ M, under the conditions satisfied by Q outlined in Lemma 2.2, (3.33) can be transformed into Expanding the definition of P, we observe where the inequality is being evaluated at (x, T ).At this point, the inequality does not depend on the (x, t 0 ), so we have Then, at any point (x, t) ∈ M × (0, ∞), we derive Hence, dividing by t on both sides of (3.36) leads to Finally, we aim to bound the case of x ∈ ∂M.Following [13], here we employ an idea that has been developed in [19] and [18] to deal with manifolds.We will only sketch the key steps.Choosing an orthonormal frame {e 1 , e 2 ......e n } at x so that e n = ν and ∂ ν P(x, t 0 ) = 0, we obtain Using the fact that ∂ ν u = 0 and ∂ ν J = 0 on ∂M × (0, ∞), and dividing by t 0 ω(|∇u| + ε), at point (x, t 0 ), one can show that This leads to a contradiction.Then, this case can not occur, which completes the proof.□ Remark 3.1.In Table 1, there have been observed inconsistencies in the assignment of constants among the theorems.However, in the case where a = 0, equation (3.1) can indeed be simplified to equation (1) using appropriate transformations.The constant C 1 in Theorem A is equivalent to the constant C 2 in Theorem 3.1, and the remaining constants in the theorems exhibit a similar relationship.The equivalence of these constants underscores the close relationship between the two theorems, further reinforcing the interconnectedness and consistency of the results.

Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.