Let τ(n) be the Dirichlet divisor function and k⩾2 be a fixed integer. We give an asymptotic formula of the mean square of
Δk(x)=∑n1,⋯,nk⩽xτ(n1⋯nk)−xkPk(logx).
Citation: Zhen Guo. On mean square of the error term of a multivariable divisor function[J]. AIMS Mathematics, 2024, 9(10): 29197-29219. doi: 10.3934/math.20241415
[1] | Marco Cirant, Kevin R. Payne . Comparison principles for viscosity solutions of elliptic branches of fully nonlinear equations independent of the gradient. Mathematics in Engineering, 2021, 3(4): 1-45. doi: 10.3934/mine.2021030 |
[2] | Isabeau Birindelli, Kevin R. Payne . Principal eigenvalues for k-Hessian operators by maximum principle methods. Mathematics in Engineering, 2021, 3(3): 1-37. doi: 10.3934/mine.2021021 |
[3] | Edgard A. Pimentel, Miguel Walker . Potential estimates for fully nonlinear elliptic equations with bounded ingredients. Mathematics in Engineering, 2023, 5(3): 1-16. doi: 10.3934/mine.2023063 |
[4] | Aleksandr Dzhugan, Fausto Ferrari . Domain variation solutions for degenerate two phase free boundary problems. Mathematics in Engineering, 2021, 3(6): 1-29. doi: 10.3934/mine.2021043 |
[5] | Antonio Vitolo . Singular elliptic equations with directional diffusion. Mathematics in Engineering, 2021, 3(3): 1-16. doi: 10.3934/mine.2021027 |
[6] | Isabeau Birindelli, Giulio Galise . Allen-Cahn equation for the truncated Laplacian: Unusual phenomena. Mathematics in Engineering, 2020, 2(4): 722-733. doi: 10.3934/mine.2020034 |
[7] | Youchan Kim, Seungjin Ryu, Pilsoo Shin . Approximation of elliptic and parabolic equations with Dirichlet boundary conditions. Mathematics in Engineering, 2023, 5(4): 1-43. doi: 10.3934/mine.2023079 |
[8] | Menita Carozza, Luca Esposito, Raffaella Giova, Francesco Leonetti . Polyconvex functionals and maximum principle. Mathematics in Engineering, 2023, 5(4): 1-10. doi: 10.3934/mine.2023077 |
[9] | Juan-Carlos Felipe-Navarro, Tomás Sanz-Perela . Semilinear integro-differential equations, Ⅱ: one-dimensional and saddle-shaped solutions to the Allen-Cahn equation. Mathematics in Engineering, 2021, 3(5): 1-36. doi: 10.3934/mine.2021037 |
[10] | Lucio Boccardo . A "nonlinear duality" approach to W1,10 solutions in elliptic systems related to the Keller-Segel model. Mathematics in Engineering, 2023, 5(5): 1-11. doi: 10.3934/mine.2023085 |
Let τ(n) be the Dirichlet divisor function and k⩾2 be a fixed integer. We give an asymptotic formula of the mean square of
Δk(x)=∑n1,⋯,nk⩽xτ(n1⋯nk)−xkPk(logx).
Tanti auguri caro amico Sandro! Con molta stima, Italo.
The purpose of this paper is to provide a review of a few classical and some quite recent results about the validity of the weak Maximum Principle in a rather general setting involving non smooth functions satisfying fully nonlinear degenerate elliptic partial differential inequalities in the viscosity sense in an unbounded domain Ω. Since no regularity assumption will be made on ∂Ω the classical methods based on the construction of suitable distance-like barrier functions are not viable in our framework.
Consider the set USC(¯Ω) of upper semicontinuous functions on ¯Ω and a mapping
F:Ω×R×Rn×Sn→R |
where Sn is the space of n×n symmetric matrices and F is degenerate elliptic, that is
F(x,s,p,Y)≥F(x,s,p,X) if Y−X is non-negative definite | (1.1) |
In this context, the weak Maximum Principle for (F,Ω) is the validity of the following sign propagation property:
u∈USC(Ω):F(x,u,Du,D2u)≥0in Ω,u≤0on ∂Ω |
implies
u≤0inΩ |
The partial differential inequality above involving in a nonlinear way the gradient and the Hessian matrix of u will be always understood in the Crandall-Ishii-Lions viscosity sense, see [15]. The motivation for considering non-smooth functions is of course a strong one when dealing with fully nonlinear partial differential inequalities or equations but is also a relevant one in the linear case, as pointed out by E. Calabi in his 1958 paper [10] on entire solutions of some elliptic partial differential inequalities. In that paper, brought to my attention by N. Garofalo, some version of the Hopf Maximum Principle is indeed proved to hold for upper semicontinuous functions u satisfying the linear partial differential inequality
Tr(A(x)D2u)+b(x)⋅Du+c(x)u≥0, |
in an appropriately defined weak sense.
It is worth to note from an historical point of view that the weak notion introduced by Calabi is in fact similar although somehow stronger than the Crandall-Lions viscosity notion, see [17] for comments in this respect.
In Section 2 I recall a few results concerning Alexandrov-Bakelman-Pucci estimate, in short (ABP), and as a by-product the weak Maximum Principle in the uniformly elliptic case for possibly unbounded domains satisfying a measure-type condition. Section 3 reviews more recent results concerning a class of directionally degenerate elliptic operators in domains satisfying geometric conditions related to the directions of ellipticity. Section 4 is the announcement of a recently obtained result in collaboration with A. Vitolo concerning the weak Maximum Principle in the case of cooperative systems.
Due to the expository character of the paper the proofs of the reviewed results are just sketched; for full details we refer the interested reader to the original papers [11,12,13,14].
It is well-known that the weak Maximum Principle may not hold in unbounded domains: just observe that
u(x)=1−1|x|n−2 |
with n≥3 satisfies Δu=0 in the exterior domain Ω=Rn∖¯B1(0), u≡0 on ∂Ω but u>0 in Ω.
Some remarkable results concerning the validity of (ABP) estimates and the weak Maximum Principle for linear uniformly elliptic operators in unbounded domains are due to X. Cabré [7]. He considered domains satisfying the following measure-geometric condition:
(G) for fixed numbers σ,τ∈(0,1), there exists a positive real number R(Ω) such that for any y∈Ω there exists an n-dimensional ball BRy of radius Ry≤R(Ω) satisfying
y∈BRy,|BRy∖Ωy,τ|≥σ|BRy| |
where Ωy,τ is the connected component of Ω∩BRy/τ containing y.
The above measure-type condition introduced in [3] requires, roughly speaking, that there is enough boundary near every point in Ω allowing, so to speak, to carry over the information on the sign of u from the boundary to the interior of the domain.
Note that (G) holds in particular in the following cases:
● Ω is bounded: R(Ω)=C(n)diam(Ω)
● Ω is unbounded with finite Lebesgue measure: R(Ω)=C(n)|Ω|1n
● Ω is a cylinder with diam(Ω)=|Ω|=+∞ such as {x=(x1,x′):a<x1<b} : R(Ω)=b−a
Since (G) implies the metric condition supΩdist(y,∂Ω)<+∞, (G) does not hold on cones. Note also that the periodically perforated plane
R2per=R2∖⋃(i,j)∈Z2¯Br(i,j) |
as well as its n−dimensional analogues or variants are (G) domains while of course this is not the case for exterior domains such as Rn∖¯Br(0).
For domains satisfying condition (G), X.Cabré [7] proved an Alexandrov-Bakelman-Pucci estimate for strong solutions:
Theorem 2.1. Let Ω satisfy (G), f∈Ln(Ω) and c≤0. If u∈W2,n(Ω) satisfies almost everywhere the uniformly elliptic partial differential inequality
Tr(A(x)D2u)+b(x)⋅Du+c(x)u≥f(x)inΩ |
with
A(x)ξ⋅ξ≥λ|ξ|2,λ>0, |
then
(ABP)supΩu≤sup∂Ωu++CR(Ω)||f||Ln(Ω) | (2.1) |
If f≡0 and u≤0 on ∂Ω, the validity of weak Maximum Principle follows from (2.1) in the case of linear uniformly elliptic operators. Observe also that (2.1) extends to the case of (G) domains the classical (ABP) estimate for bounded domains, see [16].
Some of the results of [7] have been later generalised in two directions in the paper [11]:
● viscosity solutions of fully nonlinear uniformly elliptic inequalities
● more general domains satisfying a weaker form of (G)
The above mentioned condition, which differs from (G) since does not require a priori the boundedness of the radii Ry reads as
(wG) there exist constants σ,τ∈(0,1) such that for all y∈Ω there is a ball BRy of radius Ry containing y such that
|BRy∖Ωy,τ|≥σ|BRy| |
where Ωy,τ is the connected component of Ω∩BRy/τ containing y.
If supy∈ΩRy<+∞, then (wG) boils down to condition (G). Relevant examples of unbounded domains satisfying condition (wG) but not (G) are cones of Rn and their unbounded subsets. Indeed, condition (wG) is satisfied in this case of cones with Ry=O(|y|) as |y|→∞.
A less standard example which shows that there is no a priori restriction on the growth at infinity of Ry is the plane domain described in polar coordinates as
Ω=R2∖{ϱ=eθ, θ≥0} |
Here (wG) holds with Ry=O(e|y|) as |y|→∞.
In order to state the main result in [11] concerning a nonlinear version of the Alexandrov-Bakelman-Pucci we assume that the structure conditions (SC) hold:
● F is continuous with respect to all variables x,t,p,X
● uniform ellipticity: λTr(Q)≤F(x,t,p,X+Q)−F(x,t,p,X)≤ΛTr(Q)
for all Q positive semidefinite and some 0<λ≤Λ
● properness: t↦F(x,t,p,X) is nonincreasing
● linear growth with respect to the gradient F(x,0,p,O)≤β(x)|p| for some bounded β
Under the assumptions above we proved the following form of the (ABP) estimate for bounded above viscosity subsolutions:
Theorem 2.2. Let u∈USC(¯Ω) with supΩu<+∞ satisfy in the viscosity sense
F(x,u,Du,D2u)≥f(x),x∈Ω |
where f∈C(Ω)∩L∞(Ω). If Ω satisfies (wG) for some σ,τ∈(0,1), F satisfies (SC) and, moreover,
(C)supy∈ΩRy‖ |
then
\sup\limits_{\Omega}u \leq \sup\limits_{ \partial \Omega}u^+ +C\, \sup\limits_{y\in \Omega} R_y\, \Vert f^- \Vert_{L^n(\Omega_{y, \tau})}\, |
for some positive constant C depending on n, \ \lambda, \ \Lambda, \sigma, \tau and \sup_{y\in \Omega} R_y\, \Vert \beta \Vert _{L^\infty (\Omega_{y, \tau})}\, .
For u\leq 0 on \partial\Omega and f\geq 0 the above yields the weak Maximum Principle.
To obtain the {\bf(ABP)} estimate in this more general case we will assume, besides condition {\bf(wG)} on the domain, the extra condition (C) which can be seen as a control of the interplay between the geometrical size of the domain (the opening, in the case of cones) and the growth of the first order transport term.
This condition is trivially satisfied if \sup_{y\in\Omega} R_y\leq R_0 < +\infty in {\bf(wG)}, i.e., if \Omega satisfies {\bf(G)} or, for example, when b\equiv 0 , namely when F does not depend on the first order transport term. Let us observe also that condition (C) is unavoidable in general, as shown by by the next
Example 2.3. The function u(x) = u(x_1, x_2) = \left(1-{\rm e}^{1-x_1^\alpha}\right)\, \left(1-{\rm e}^{1-x_2^\alpha}\right)\, , with 0 < \alpha < 1 , is bounded and strictly positive in the cone \Omega = \left\{x = (x_1, x_2)\in {\mathbb R}^2\, :\ x_1 > 1, \ x_2 > 1\right\} and satisfies
u\equiv 0\quad\ on\;{\partial\Omega }, \quad \Delta u+b(x)\cdot Du = 0\quad in \;\Omega |
where the vectorfield b is given by b(x, y) = \left(\frac{\alpha}{x_1^{1-\alpha}}+\frac{1-\alpha}{x_1}\; , \; \frac{\alpha}{x_2^{1-\alpha}}+\frac{1-\alpha}{x_2}\right)\, .
As observed above, \Omega satisfies {\bf(wG)} with R_z = O(|z|) as |z|\to \infty . Since ||b||_{L^\infty (\Omega_{z, \tau})} = 1 the interplay condition (C) fails in this example.
Some non trivial cases in which condition (C) is fulfilled are:
(a) the cylinder \Omega = \left\{x = (x', x_n)\in {\mathbb R}^{n-1}\times {\mathbb R}\; :\; |x'| < 1\, , \ x_n > 0\right\} is a (G) domain and therefore condition (C) is obviously satisfied,
(b) \Omega = \left\{ (x', x_n)\in {\mathbb R}^{n-1}\times {\mathbb R}\; :\; x_n > |x'|^q\right\}\, with q > 1 , This non convex conical set is a genuine (wG) domain with radii R_y = O\left(|y|^{1/q}\right) as |y|\to \infty .
In this case, (C) imposes to the function \beta a rate of decay \beta(y) = O\left(1/|y|^{1/q}\right) as |y|\to \infty .
(c) \Omega is the strictly convex cone \left\{ x\in {\mathbb R}^n\setminus \{0\}\; :\; x/|x|\in \Gamma\right\} where \Gamma is a proper subset of the unit half-sphere S^{n-1}_+ = \left\{ x = (x', x_N)\in {\mathbb R}^{n-1}\times {\mathbb R} \; :\; |x| = 1\, , \ x_n > 0\right\} .
In this case (wG) is satisfied with R_y = O(|y|) for |y|\to \infty and condition (C) requires the rate of decay \beta(y) = O\left(1/|y|\right) as |y|\to \infty .
Note that cases (a) and (c) can be seen as limiting cases of situation (b) when, respectively, q\to +\infty and q = 1 .
The proof of Theorem 2.2 is based on several viscosity calculus tools but not on the Comparison Principle for viscosity sub and supersolutions. It relies indeed on rescaled versions of the weak Harnack and of the boundary weak Harnack inequality for operators including first order terms, see [11] and [8] for the case of Pucci operators.
Let us briefly sketch the main steps of the proof with the aid of the following three lemmas. The minimal Pucci operator with parameters 0 < \lambda\leq \Lambda is defined by
P^-_{\lambda, \Lambda}(X) = \lambda\, {\rm Tr}(X^+) -\Lambda\, {\rm Tr}(X^-) |
where X^+ and X^- are nonnegative definite matrices such that X = X^+-X^- and X^+X^- = O . Under our assumption the following holds:
Lemma 2.4. If w\in USC(\Omega) satisfies
F(x, w(x), Dw(x), D^2w(x)) \geq 0 \;, \; x\in\Omega |
in the viscosity sense, then for each M\in {\mathbb R} the function u = M-w^+ \in LSC(\Omega) satisfies
P^-_{\lambda, \Lambda}(D^2u) -\beta(x)|Du|\leq f^-(x) \;, \; x\in\Omega |
The proof is a direct application of the viscosity notions.
Lemma 2.5 (the rescaled weak Harnack inequality). Let R > 0 , \tau\in(0, 1) and f\in C(B_{R/\tau})\cap L^\infty(B_{R/\tau}) . If u\in LSC(\overline B_{R/\tau}), \, u\geq 0 is a viscosity solution of
P^-_{\lambda, \Lambda}(D^2u) -\beta(x)|Du|\leq f^-(x) \;, \; x\in B_{ R/\tau} |
then
\begin{equation} \left(\frac{1}{|B_R|}\int_{B_R} u^p\right)^{1/p} \leq C^* \, \left( \inf\limits_{B_{R}} u + R\, \|f^-\|_{L^N(B_{R/\tau})}\;\right)\, \end{equation} | (2.2) |
where p and C^* are positive constants depending on \lambda, \Lambda, n, \tau and on the product R\, \|\beta\|_{L^\infty(B_{R/\tau})} .
Inequality (2.2) is a consequence of the weak Harnack inequality for nonnegative viscosity solutions in a cube Q_1 with side-length 1
\Vert u\Vert_{L^{\hat p}(Q_{1/4})}\leq \hat{C}\left( \inf\limits_{Q_{1/2}}\, u+\Vert f^-\Vert_{L^N(Q_1)}\right) \, , |
where \hat{p} and \hat{C} depend on \lambda, \ \Lambda, \ N as well as on \|\beta\|_{L^\infty(Q_1)} , see [8]. Inequality (2.2) is indeed obtained from the above by rescaling and using the positive homogeneity of the operator {P}^-_{\lambda, \Lambda} and a standard covering argument.
Let A be a bounded domain in {\mathbb R}^n and B_R\, , B_{R/\tau} be concentric balls such that
A \cap B_R \not = \emptyset\;\;, \;\;B_{R/\tau} \backslash A \not = \emptyset\, \;. |
For u \in LSC(\bar A) , u\geq 0 , consider the following lower semicontinuous extension u^-_m of function u
u^-_m (x) = \left\{ \begin{array}{ll} \min(u(x); m) & {\rm if }\; x \in A \\ m & {\rm if }\; x \not \in A \end{array}\right. \, |
where m = \inf_{x \in \partial A \cap B_{R/\tau}} u(x)\, .
Lemma 2.6. With the above notations, if g \in C(A)\cap L^\infty (A) and u \in LSC(\bar A) satisfy
u \geq 0\;\;, \;\; {P}^-_{\lambda, \Lambda}(D^2 u)-b(x)|Du|\leq g(x)\;\;\mathit{\mbox{in}}\;\;A |
in the viscosity sense, then
\Big ({\frac{1}{ |B_{R}|}} \int_{B_{R}} (u^-_m) ^{p} \Big )^{1/p} \leq C^*\left(\inf\limits_{A \cap B_{R}} u + R\, \|g^+\|_{L^n (A \cap B_{R/\tau})}\right)\, |
where p and C^* are positive constants depending on \lambda, \Lambda, n, \tau and on the product R\, \|\beta\|_{L^\infty(B_{R/\tau})} .
It is not hard to check, see [8] for the case b\equiv 0 , that u^-_m is a viscosity solution of
u^-_m \geq 0\;\;, \;\; {P}^-(D^2 u^-_m)-\chi_{A}(x)\, b(x)\, |Du^-_m|\leq \chi_{A}(x)\, g^+(x)\;\;\mbox{in}\;\; B_{R/\tau}\, |
where \chi_{A} is the characteristic function of A . The statement follows then by applying the weak Harnack inequality (2.2) to u^-_m\, , observing that \inf _{B_{R}} u^-_m\leq \inf _{A \cap B_{R}} u and that \chi_A\, g^+\equiv 0 outside A .
The role of Lemma 2.6 in the proof of Theorem 2.2 is a crucial one, in connection with the {\bf(wG)} condition, in establishing a localized form of the {\bf(ABP)} estimate. Indeed, if w\in USC(\Omega) satisfies F(x, w, Dw, D^2w)\geq f then, as we have seen above, for u(x) = \sup_{\Omega}w^+-w(x)\geq 0 we have
P^-_{\lambda, \Lambda}(D^2u) -\beta(x)|Du|\leq f^-(x) \;, \; x\in\Omega |
Apply now Lemma 2.6 with A = \Omega_{y, \tau} and g = f^- to obtain
\left( \frac{1}{|B_{R_y}|}\int_{B_{R_y}}(u^-_m)^p\right)^{1/p}\leq C^*_y\, \left( \inf\limits_{\Omega_{y, \tau}\cap B_{R_y}}u +R_y\, \Vert f^-\Vert_{L^N(\Omega_{y, \tau})}\right)\, |
for positive constants p and C^*_y depending on N, \ \lambda, \ \Lambda, \ \tau and R_y\, \Vert b\Vert _{L^\infty (\Omega_{y, \tau})} .
The left-hand side of the above inequality can be of course estimated from below as follows
\left(\frac{1}{|B_{R_y}|}\int_{B_{R_y}}(u^-_m)^p\right)^{1/p}\geq \left( \frac{1}{|B_{R_y}|}\int_{B_{R_y}\setminus \Omega_{y, \tau}}(u^-_m)^p\right)^{1/p} = m\, \left( \frac{|B_{R_y}\setminus \Omega_{y, \tau}|}{|B_{R_y}|}\right)^{1/p}\, |
where m = \inf_{x \in \partial \Omega_{y, \tau} \cap B_{R/\tau}} u(x)\, . Hence, using condition {\bf(wG)},
\left(\frac{1}{|B_{R_y}|}\int_{B_{R_y}}(u^-_m)^p\right)^{1/p}\geq m\, \sigma^{1/p}\, |
At this point a few technicalities lead to estimate the value of w^+ at any y\in \Omega as follows:
\begin{equation} w^+(y) \leq (1-\theta_y)\, \sup\limits_{\Omega}w^+ + \theta_y\, \sup\limits_{ \partial \Omega}w^+ + R_y\, \Vert f^- \Vert_{L^N(\Omega_{y, \tau})}\, \end{equation} | (2.3) |
for suitable choice of \theta_y \in(0, 1) depending on n, \ \lambda, \ \Lambda, \ \sigma, \ \tau and on y through the quantity R_y\, \Vert b\Vert _{L^\infty (\Omega_{y, \tau})} which appears in the statement of Theorem 2.2. The proof of the theorem is easily deduced from (2.3) through the use of condition (C).
This section is about the validity of the weak Maximum Principle for degenerate elliptic operators F which are strictly elliptic on unbounded domains \Omega of \mathbb R^n whose geometry is related to the direction of ellipticity. Some results of that kind for one-directional elliptic operators in bounded domains have been previously established, among other qualitative properties, by Caffarelli-Li-Nirenberg [9].
We assume:
\begin{equation} F \;\; {\rm \, is \, continuous \, with\, respect\, to\, } (x, t, p, X) \end{equation} | (3.1) |
the following standard monotonicity conditions on F
\begin{equation} F(x, s, p, Y) \ge F(x, s, p, X) \quad \hbox{ if} \ Y\ge X \end{equation} | (3.2) |
\begin{equation} F(x, s, p, X)\le F(x, r, p, X) \quad \hbox{ if} \ s \gt r \end{equation} | (3.3) |
and, just for simplicity,
\begin{equation} F(x, 0, 0, O) = 0 \quad \forall\, x \in \Omega\, \end{equation} | (3.4) |
where O is the null matrix. Our last requirement is on a Lipschitz behavior of F with respect to the gradient variable,
\begin{equation} \mbox{there exists}\, \gamma \gt 0: |F(x, 0, p, X)-F(x, 0, 0, X)| \le \gamma \, |p|\, \mbox{for all}\, p\in \mathbb R^n \end{equation} | (3.5) |
The class of domains that will be considered in this section is conveniently described by decomposing {\mathbb R}^n as \mathbb R^n = U\bigoplus U^\bot , where U is a k -dimensional subspace and U^\bot is its orthogonal complement and denoting by P and P^\bot the projection matrices on U and U^\bot , respectively.
We will consider the open domains \Omega which are contained in unbounded slabs whose k -dimensional sections are cube of edge d , namely
\begin{equation} \Omega \subset \{x \in \mathbb R^n : a \le x\cdot \nu^h \le a+d_h, \; h = 1, \dots, k\}: = C\quad \mbox{for some} \; a \in \mathbb R, d_h \gt 0\, , \end{equation} | (3.6) |
where \{\nu^1, \dots, \nu^k\} is an orthonormal system for the subspace U . We will refer to vectors in U^\bot as unbounded directions.
Such domains may of course be unbounded and of infinite Lebesgue measure but they do satisfy the measure-geometric (wG) condition of the previous section. Note once more that no regularity requirement is made on \partial \Omega .
The next assumptions are crucial ones in establishing our results on the weak maximum and on Phragmèn-Lindelöf principles
\begin{equation} \mbox{ there exist} \, \nu \in U \, \mbox{ and} \, \lambda \gt 0: \, F(x, 0, p, X+t\nu \otimes \nu)-F(x, 0, p, X) \ge \lambda\, t \quad \end{equation} | (3.7) |
for all t > 0 .
Observe that the matrix \nu \otimes \nu is the orthogonal projection over the one dimensional subspace generated by \nu : the strict directional ellipticity condition (3.7) related to the geometry of \Omega will play a crucial role in our results.
\begin{equation} \mbox{there exists}\, \Lambda \gt 0\, : F(x, 0, 0, X+tP^\bot)-F(x, 0, 0, X) \le \Lambda t\, |x| \quad\\ \end{equation} | (3.8) |
for all t > 0 , as |x|\to \infty .
The above one-sided Lipschitz condition is satisfied in the linear case if the coefficients corresponding to second derivatives in the unbounded directions (i.e., belonging to U^\bot ) have at most linear growth with respect to x . It is worth to observe that \nu \otimes \nu and P^\bot are symmetric and positive semidefinite matrices and that (3.7), (3.8) form a much weaker condition than uniform ellipticity. Indeed, they comprise a control from below only with respect to a single direction \nu\in U together with control from above in the directions of U^\bot , a much weaker condition on F than uniform ellipticity which would indeed require a uniform control of the difference quotients both from below and from above with respect to all possible increments with positive semidefinite matrices.
We will refer collectively to conditions (3.1) to (3.8) above as (\bf {SC}) .
The following very basic two dimensional example is useful to clarify our structural assumptions. Let U = \{x = (x_1, x_2): x_2 = 0\} and consider the linear operator
F(x, u, Du, D^2u) = \frac{\partial ^2u}{\partial x_1^2 }+H(x, Du)+c(x)u |
where H , c are continuous and bounded with respect to x . Obviously, F satisfies (3.1), (3.2), (3.4) and (3.3) if c(x)\leq 0 and H(x, 0) = 0 . The strict ellipticity condition (3.7) is trivially satisfied by \nu = (1, 0) and \lambda = 1 . Condition (3.7) is fulfilled with \Lambda = 1 since
F(x, 0, 0, X+tP^\bot)-F(x, 0, 0, X) = t\, Tr (P^\bot) |
and P^\bot is the matrix
\left(\begin{array}{cc}0 & 0 \\0 & 1\end{array}\right) |
Condition (3.5) is trivially satisfied if H is Lipschitz continuous.
Operators of this kind arise in the dynamic programming approach to the optimal control of a deterministic system in {\mathbb R}^n which is perturbed by a lower dimensional Brownian motion. A more general example is
F(x, u, Du, D^2u) = {\rm Tr}(AD^2u)+H(x, Du)+c(x)u |
with
A = \left(\begin{array}{cc} \mathbb I_k & \bf 0 \hskip0.3cm\\ \bf 0 & \varphi(x) \, \mathbb I_{n-k} \end{array}\right) |
The structure condition (3.7) requires in this case a sublinear growth of the coefficient \varphi in the unbounded directions of the domain.
The Bellman-Isaacs operators are an important class of fully nonlinear operators arising in the theory of differential games. They have the following form
F(x, u, Du, D^2u) = \sup\limits_{\alpha}\inf\limits_{\beta}\, \big[{\rm Tr}(A^{\alpha\beta}D^2u)+b^{\alpha\beta}(x)Du)+ c^{\alpha\beta}u(x)\big] |
with constant coefficients depending \alpha and \beta running in some sets of indexes \mathcal{A} , \mathcal {B} .
It is not hard to check that condition (\bf{SC}) is satisfied in any domain \Omega contained in a (n-k) -infinite slab as in (3.6) if all matrices A^{\alpha\beta} , whose entries are denoted by a_{ij}^{\alpha\beta}(x) , are positive semidefinite and
\begin{align*} &A^{\alpha\beta} \nu^h \cdot \nu^h\geq \lambda \gt 0 \;\;\; h = 1, \dots, k, \\ &\sum\limits_{i, j = 1}^n a_{ij}^{\alpha\beta}\nu^h_i\nu^h_j a_{ij}^{\alpha\beta}\nu^h_i\nu^h_j \le \Lambda|x|, \;\;\; h = k+1, \dots, n \\ &|b^{\alpha\beta}(x) |\leq \gamma, \quad c^{\alpha\beta}(x)\leq 0, \end{align*} |
where \lambda, \Lambda, \gamma are independent on \alpha, \beta and \{\nu^1, \dots, \nu^k\} is an orthonormal basis of U .
The main results concerning the validity of (wMP) are stated in the following theorems:
Theorem 3.1. Let \Omega be a domain of \mathbb R^n satisfying condition (3.6):
\Omega \subset \{x \in \mathbb R^n : a \le x\cdot \nu^h \le a+d_h, \; h = 1, \dots, k\} \quad \mathit{\mbox{for some}} \; a \in \mathbb R, d_h \gt 0\, , |
and assume that F satisfies the structure condition (\bf SC) .
Then (wMP) holds for any u \in USC(\overline\Omega) such that u^+(x) = o(|x|) as |x| \to \infty .
Note that some restriction on the behaviour of u at infinity is unavoidable. Observe indeed that u(x_1, x_2, x_3) = e^{x_1}\sin x_2\sin x_3 solves the degenerate Dirichlet problem
\frac{\partial^2 u}{\partial x_1^2} +\frac{\partial^2 u}{\partial x_2^2} = 0 \ \mbox{in}\ \ \Omega, \ \ u(x_1, x_2, x_3) = 0\, \, \mbox{on}\, \, \partial\Omega |
in the 1 -infinite cylinder \Omega = \mathbb R\times (0, \pi)^2\subset \mathbb R^3 and u(x_1, x_2, x_3) > 0 in \Omega so (wMP) fails in this case.
As a consequence of the previous Theorem a quantitative estimate can be obtained:
Theorem 3.2. Let \Omega be a domain of \mathbb R^n satisfying condition (3.6) and assume that F satisfies the structure condition (\bf SC) .
If
F(x, u, Du, D^2u) \ge f(x)\quad \ in \; \Omega\, |
where f is continuous and bounded from below and u^+(x) = o(|x|) as |x| \to \infty , then
\sup\limits_{\overline\Omega} u \le \sup\limits_{\partial\Omega}u^+ + e \frac{\sup f^-(x)}{\lambda}\, d^2 |
where d = \min _h d_h and f^-(x) = -\min(f(x), 0) .
The method of proof of the above theorems is considerably different from that of the analogous results in Section 2 since tools as Harnack inequalities are not available in the present degenerate elliptic setting.
In order to simplify the notations a sketch of the proof of Theorem 3.1 is given below for the case F = F(x, \nabla^2u)+ c(x)u where c\leq 0 and \Omega = \{x\in {\mathbb R}^2: -d/2 < x_1 < d/2\} . So we take \nu = (1, 0) and the projection matrix on U is
\nu\otimes\nu = P = \left(\begin{array}{cc}1& 0 \\0 & 0\end{array}\right) |
Arguing by contradiction we suppose that u(x^0) = k > 0 at some x^0 \in \Omega .
Let us consider then, for \varepsilon > 0 , the function u^\varepsilon(x) = u(x)-\varepsilon \phi(x) where \phi(x) = \sqrt {x_2^2+1} is smooth function penalising the unbounded variable x_2 . The prescribed behavior of u at infinity implies then that there exists a bounded domain \Omega^\varepsilon \subset \Omega such that u^\varepsilon \leq 0 outside \Omega^\varepsilon . Hence, for \varepsilon < \frac{k}{2\phi(x^0)} ,
\sup\limits_{\Omega} u^\varepsilon = \max\limits_{\overline\Omega^\varepsilon} u^\varepsilon = u^\varepsilon (x^\varepsilon)\geq k/2 |
Since \phi is smooth it is easy to check that
\begin{equation} F(x, \nabla^2u^\varepsilon(x)+ \varepsilon \nabla^2 \phi(x) )+ c(x)(u^\varepsilon (x)+\varepsilon \phi(x)) \geq 0 \end{equation} | (3.9) |
in the viscosity sense. The next step is to construct a smooth function h^\varepsilon , depending on a positive parameter \alpha to be chosen later, to be used as a test function replacing u^\varepsilon in 3.9. At this purpose, set
h^\varepsilon(x) = k^\varepsilon (1+e^{-\alpha d})-k^\varepsilon e^{\alpha(x_1-x^\varepsilon_1)-\alpha d} |
where k^\varepsilon = u^\varepsilon (x^\varepsilon) . We have
h^\varepsilon(x^\varepsilon) = u^\varepsilon(x^\varepsilon) \;\mbox{and} \; h^\varepsilon \geq \frac{ k}{2} e^{-2\alpha d} \; \;\mbox{on}\; \overline \Omega |
Modulo a translation by a nonnegative constant \delta^\varepsilon we can find some point \overline x^\varepsilon such that h^\varepsilon+\delta^\varepsilon +\varepsilon \phi touches u^\varepsilon from above at \overline x^\varepsilon . So, by definition of viscosity subsolution and using the assumption that the zero-order coefficient c is \leq 0 and that h ^\varepsilon , \delta^\varepsilon and \varepsilon \phi are positive,
\begin{equation} F(\overline x^\varepsilon, \nabla^2 h ^\varepsilon+\varepsilon \nabla^2 \phi)\geq 0 \; \;\mbox{at}\;\overline x^\varepsilon. \end{equation} | (3.10) |
Observe now that \nabla^2 h^\varepsilon(\overline x^\varepsilon) = -\alpha^2 k^\varepsilon e^{\alpha(\overline x^\varepsilon_1-x^\varepsilon_1)-\alpha d}\, \nu\otimes\nu and apply the directional ellipticity assumption (3.7) with X = \nabla^2h^\varepsilon(\overline x^\varepsilon)+\varepsilon\nabla^2\phi(\overline x^\varepsilon) and t = \frac{ k}{2} e^{-2\alpha d} to obtain, taking (3.10) into account,
\begin{equation} 0\leq F(\overline x^\varepsilon, \nabla^2 h ^\varepsilon (\overline x^\varepsilon)+\varepsilon \nabla^2 \phi(\overline x^\varepsilon))\leq F(\overline x^\varepsilon, \varepsilon\nabla^2 \phi(\overline x^\varepsilon) )-\lambda\alpha^2\frac{ k}{2} e^{-2\alpha d} \end{equation} | (3.11) |
In order to let \varepsilon \to 0 in the above let us consider first the case where \overline x^\varepsilon is bounded and therefore, at least for a subsequence \overline x^\varepsilon \to \overline x as \varepsilon \to 0^+ .
By continuity of F and \nabla^2 \phi we obtain thanks to (3.4) the contradiction
0 \lt \lambda\alpha^2\frac{ k}{2} e^{-2\alpha d}\leq F( \overline x, O) = 0 |
and the Theorem is proved in this case.
Suppose now that \overline x^\varepsilon is unbounded and observe that \frac{1}{\phi(x) }P^\bot \geq \nabla^2 \phi(x) in the matrix ordering as in (3.2) where
P^\bot = \left(\begin{array}{cc}0 & 0 \\0 & 1\end{array}\right) |
Hence, by degenerate ellipticity and (3.11),
\lambda\alpha^2\frac{ k}{2} e^{-2\alpha d}\leq F(\overline x^\varepsilon, \varepsilon\nabla^2 \phi(\overline x^\varepsilon) )\leq F(\overline x^\varepsilon, \varepsilon \frac{1}{\phi(\overline x^\varepsilon) }P^\bot ) |
Using condition (3.8) with t = \frac{\varepsilon}{\phi(\overline x^\varepsilon)} and X = O we have
F(\overline x^\varepsilon, \varepsilon \frac{1}{\phi(\overline x^\varepsilon) }P^\bot ) = F(\overline x^\varepsilon, \varepsilon \frac{1}{\phi(\overline x^\varepsilon) }P^\bot) -F(\overline x^\varepsilon, O) \leq \Lambda \frac{\varepsilon}{\phi(\overline x^\varepsilon)}|\overline x^\varepsilon| |
Since \frac{|\overline x^\varepsilon|} {\phi(\overline x^\varepsilon)} remains bounded as \varepsilon \to 0 , using the above estimate we derive the contradiction
0 \lt \lambda\alpha^2\frac{ k}{2} e^{-2\alpha d}\leq\liminf\limits_{\varepsilon \to 0}F(\overline x^\varepsilon, \varepsilon \frac{1}{\phi(\overline x^\varepsilon) }P^\bot) = 0 \, . |
Let us now sketch the proof of Theorem 3.2 in the setting adopted for the proof of Theorem 3.1 and assume as a further simplification that \Omega is the strip (0, 1)\times {\mathbb R} . The first step is to show by simple viscosity calculus that if u satisfies F(x, D^2u)+c(x)u \ge f(x) in \Omega with c\leq 0 , then u^+ = {\rm max}(u, 0) satisfies
\begin{equation} F(x, D^2u^+) \ge -f^-(x)\, . \end{equation} | (3.12) |
Consider the auxiliary function
w(x) = u^+(x) - \sup\limits_{\partial\Omega}u^+ +C(e^{ x_1} -e)\, |
where C is a positive constant to be chosen in the sequel.
The structure condition (\bf SC) yields
F(x, \nabla^2 w) \ge F(x, \nabla^2u^+)+ \lambda C e^{ x_1} \ge -f^-(x) + \lambda C \;\hbox{ in} \, \Omega\, |
Choosing C = \sup_{\Omega}\frac{f^-}{\lambda} we obtain
F(x, \nabla^2w) \ge 0 \ \, \hbox{in}\, \ \Omega, \ \ w \le 0 \, \hbox{on}\, \partial\Omega\, |
By Theorem 3.1 we conclude that w \le 0 in \Omega , yielding
u(x) \le u^+(x) \le \sup\limits_{\partial\Omega}u^+ +e \, \sup\limits_{\Omega}\frac{f^-}{\lambda}\, \, |
which proves the statement in this slightly simplified case. The case d\neq 1 can be dealt by working with the rescaled variable y = \frac{x}{d} since the function v(x) = u(dy) satisfies the inequality
G(y, \nabla^2v(y)\geq d^2 f(dy) |
for y\in d^{-1} \Omega \subset (0, 1)\times {\mathbb R} where G(y, X) = d^2 F(dx, d^{-2} X) .
It is well-known that the weak Maximum Principle for linear uniformly elliptic operators
{\rm Tr}(A(x) D^{2}u)+b(x)\cdot Du+c(x)u\;\; \hbox{in}\; \Omega\, |
in bounded domains may hold by replacing the standard assumption c\leq 0 by suitable alternative conditions, see for example [3,18].
The next result which easily follow from the previous Theorem go in this direction.
Theorem 3.3. Let \Omega satisfy (3.6) and assume that F satisfies (\bf SC) with (3.3) replaced by the weaker condition
\begin{equation} F(x, s, p, M)-F(x, r, p, M) \le c(x)\, (s-r) \end{equation} | (3.13) |
for some continuous bounded function c(x) > 0 and for all s > r .
Then ({\bf wMP}) holds for u \in USC(\overline \Omega) , u bounded above, provided \frac{1}{\lambda} \min d_h^2 \sup_{\Omega} c(x) is small enough.
As already observed F(x, u, Du, D^2u) \ge 0 implies F(x, u^+, Du^+, D^2u^+) \ge 0 and so, by (3.13) with r = 0, s = u^+ ,
F(x, 0, Du^+, D^2u^+)\ge F(x, u^+, Du^+, D^2u^+)-cu^+\ge -cu^+. |
Since u^+ = 0 on \partial \Omega , Theorem (3.2) applies with f = -c\, u^+ leading to
\sup\limits_{\overline\Omega} u^+ \le \frac{e d^2 C\ }{\lambda}\sup\limits_{\overline\Omega} u^+ |
where C\geq c(x) . From this estimate the statement follows if K = \frac{e d^2 C\ }{\lambda} < 1 .
For fixed C > 0 this results applies to narrow domains, that is those whose thickness d is sufficiently small. Conversely, for fixed d > 0 (wMP) holds provided C is a sufficiently small positive number. The estimate above shows also that (wMP) holds if the directional ellipticity constant \lambda is large.
In a work in progress with A.Vitolo [14] we consider systems of elliptic partial differential inequalities of the form
\begin{equation} F[u]+C(x)u \geq 0 \end{equation} | (4.1) |
in a bounded domain \Omega\subset {\mathbb R}^n . Here u = (u_1, ..., u_N) is a vector function u: {\mathbb R}^n\to {\mathbb R}^N , C(x) = (c_{ij}(x)) is a N\times N matrix and F = (F_1, ..., F_N) are second order operators acting on u of the form
F_i[u] = F_i(x, \nabla u_i, \nabla^2 u_i)\;\;(i = 1, ..., N) |
It is well-known that the validity of the weak Maximum Principle, namely the sign propagation property
u_i\leq 0 \;\mbox{in} \;\partial \Omega \quad \mbox{implies}\; u_i\leq 0 \; \mbox{in}\; \Omega \;\;(i = 1, ..., N) |
may fail for general matrices C . Several papers, see for example [1,5,19] for linear F_i and [6] for fully nonlinear differential operators, considered such systems in the case C(x) is a cooperative matrix, often called essentially positive, see [4], that is
\begin{equation} c_{ij}(x)\geq 0\; \mbox {for all}\; i \neq j\;, \; \sum\limits_{j = 1} ^N c_{ij}(x)\leq 0\, \end{equation} | (4.2) |
in order to establish the validity of the weak Maximum Principle. The results obtained in those papers require uniform ellipticity of the operators F_i .
In our work in progress we select a sufficient condition guaranteing the validity of a vector form of the weak Maximum Principle for upper semicontinuous viscosity solutions of systems (4.1), where the interaction matrix satisfies (4.2), in the more general context of operators F_i which are just degenerate elliptic, namely
0\leq F_i(x, p, X+Q)-F_i(x, p, X) \;\;(i = 1, ..., N) \; |
for all matrix Q\geq O . A basic example of such operators are Bellman operators
F_i(x, \nabla^2u_i) = \sup\limits_{\gamma \in \Gamma}{\rm Tr}\big(A_i^{\gamma}(x) \nabla^2 u_i\big)+b_i^\gamma(x)\cdot \nabla u_i+c_i^\gamma(x) u_i |
where A_i^{\gamma}(x) are symmetric positive semidefinite and \gamma is a parameter running in an arbitrary set \Gamma .
For the proof of that result, which seems to be new even in the linear case
F_i(x, \nabla u_i, \nabla^2u_i) = {\rm Tr} \big(A_i(x) \nabla^2 u_i\big)+b_i(x)\cdot \nabla u_i+c_i(x) u_i |
where A^i(x) are symmetric positive semidefinite matrices, the reader is referred to [14].
It relies on a suitable combination of results in [2] concerning the validity of the weak Maximum Principle for scalar degenerate elliptic operators with ideas in [6] and some viscosity calculus to exploit conveniently the cooperativity assumption.
The author declares no conflict of interest.
[1] |
G. F. Voronoï, Sur une fonction transcendante et ses applications à la sommation de quelques séries, Ann. Sci. École Norm. Sup., 21 (1904), 207–267. https://doi.org/10.24033/asens.539 doi: 10.24033/asens.539
![]() |
[2] | J. G. van der Corput, Zum teilerproblem, Math. Ann., 98 (1928), 697–716. https://doi.org/10.1007/BF01451619 |
[3] | G. Kolesnik, On the order of \zeta ({1\over 2}+it) and \Delta (R), Pacific J. Math., 98 (1982), 107–122. |
[4] |
M. N. Huxley, Exponential sums and lattice points III, Proc. London Math. Soc., 87 (2003), 591–609. https://doi.org/10.1112/S0024611503014485 doi: 10.1112/S0024611503014485
![]() |
[5] |
H. Cramér, Über zwei Sätze des Herrn G. H. Hardy, Math. Z., 15 (1922), 201–210. https://doi.org/10.1007/BF01494394 doi: 10.1007/BF01494394
![]() |
[6] | K. C. Tong, On division problems I, Acta Math. Sinica, 5 (1955), 313–324. |
[7] | E. Preissmann, Sur la moyenne quadratique du terme de reste du problème du cercle, C. R. Acad. Sci. Paris Sér. I Math., 306 (1988), 151–154. |
[8] |
Y. K. Lau, K. M. Tsang, On the mean square formula of the error term in the Dirichlet divisor problem, Math. Proc. Cambridge Philos. Soc., 146 (2009), 277–287. https://doi.org/10.1017/S0305004108001874 doi: 10.1017/S0305004108001874
![]() |
[9] |
L. Tóth, W. Zhai, On multivariable averages of divisor functions, J. Number Theory, 192 (2018), 251–269. https://doi.org/10.1016/j.jnt.2018.04.015 doi: 10.1016/j.jnt.2018.04.015
![]() |
[10] | A. P. Ivić, The Riemann zeta-function, Wiley-Interscience Publication, 1985. |
[11] |
K. M. Tsang, Higher-power moments of \Delta(x), \; E(t) and P(x), Proc. London Math. Soc., 65 (1992), 65–84. https://doi.org/10.1112/plms/s3-65.1.65 doi: 10.1112/plms/s3-65.1.65
![]() |
[12] |
W. Zhai, On higher-power moments of \Delta(x) II, Acta Arith., 114 (2004), 35–54. https://doi.org/10.4064/aa114-1-3 doi: 10.4064/aa114-1-3
![]() |
[13] | L. Tóth, Short proofs, generalizations, and applications of certain identities concerning multiple Dirichlet series, J. Integer Seq., 26 (2023), 15. |
[14] |
R. Heyman, L. Tóth, Hyperbolic summation for functions of the GCD and LCM of several integers, Ramanujan J., 62 (2023), 273–290. https://doi.org/10.1007/s11139-022-00681-2 doi: 10.1007/s11139-022-00681-2
![]() |