Loading [MathJax]/jax/output/SVG/jax.js
Research article

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

    tf=Δf+af(lnf)p

    with Neumann boundary conditions on a compact Riemannian manifold satisfying the integral Ricci curvature assumption, where p0 is a real constant. This contrasts Olivé's gradient estimate, which works mainly for the heat equation rather than nonlinear parabolic equations and the result can be regarded as a generalization of the Li-Yau [P. Li, S. T. Yau, On the parabolic kernel of the Schrödinger operator, Acta Math., 156 (1986), 153–201] and Olivé [X. R. Olivé, Neumann Li-Yau gradient estimate under integral Ricci curvature bounds, Proc. Amer. Math. Soc., 147 (2019), 411–426] gradient estimates.

    Citation: Hao-Yue Liu, Wei Zhang. Neumann gradient estimate for nonlinear heat equation under integral Ricci curvature bounds[J]. AIMS Mathematics, 2024, 9(2): 3881-3894. doi: 10.3934/math.2024191

    Related Papers:

    [1] Yuwen He, Jun Li, Lingsheng Meng . Three effective preconditioners for double saddle point problem. AIMS Mathematics, 2021, 6(7): 6933-6947. doi: 10.3934/math.2021406
    [2] Snezhana Hristova, Kremena Stefanova . p-moment exponential stability of second order differential equations with exponentially distributed moments of impulses. AIMS Mathematics, 2021, 6(3): 2886-2899. doi: 10.3934/math.2021174
    [3] Waleed Mohamed Abd-Elhameed, Anna Napoli . New formulas of convolved Pell polynomials. AIMS Mathematics, 2024, 9(1): 565-593. doi: 10.3934/math.2024030
    [4] Akbar Ali, Ivan Gutman, Boris Furtula, Abeer M. Albalahi, Amjad E. Hamza . On chemical and mathematical characteristics of generalized degree–based molecular descriptors. AIMS Mathematics, 2025, 10(3): 6788-6804. doi: 10.3934/math.2025311
    [5] Ze-Miao Dai, Jia-Bao Liu, Kang Wang . Analyzing the normalized Laplacian spectrum and spanning tree of the cross of the derivative of linear networks. AIMS Mathematics, 2024, 9(6): 14594-14617. doi: 10.3934/math.2024710
    [6] Jia-Bao Liu, Kang Wang . The multiplicative degree-Kirchhoff index and complexity of a class of linear networks. AIMS Mathematics, 2024, 9(3): 7111-7130. doi: 10.3934/math.2024347
    [7] Mohamed Abdalla . On Hankel transforms of generalized Bessel matrix polynomials. AIMS Mathematics, 2021, 6(6): 6122-6139. doi: 10.3934/math.2021359
    [8] Emrah Polatlı . On some properties of a generalized min matrix. AIMS Mathematics, 2023, 8(11): 26199-26212. doi: 10.3934/math.20231336
    [9] Abdulrahman B. Albidah, Ibraheem M. Alsulami, Essam R. El-Zahar, Abdelhalim Ebaid . Advances in mathematical analysis for solving inhomogeneous scalar differential equation. AIMS Mathematics, 2024, 9(9): 23331-23343. doi: 10.3934/math.20241134
    [10] Yang Zhang, Jizhu Nan . A note on the degree bounds of the invariant ring. AIMS Mathematics, 2024, 9(5): 10869-10881. doi: 10.3934/math.2024530
  • In this paper, we consider a Li-Yau gradient estimate on the positive solution to the following nonlinear parabolic equation

    tf=Δf+af(lnf)p

    with Neumann boundary conditions on a compact Riemannian manifold satisfying the integral Ricci curvature assumption, where p0 is a real constant. This contrasts Olivé's gradient estimate, which works mainly for the heat equation rather than nonlinear parabolic equations and the result can be regarded as a generalization of the Li-Yau [P. Li, S. T. Yau, On the parabolic kernel of the Schrödinger operator, Acta Math., 156 (1986), 153–201] and Olivé [X. R. Olivé, Neumann Li-Yau gradient estimate under integral Ricci curvature bounds, Proc. Amer. Math. Soc., 147 (2019), 411–426] gradient estimates.



    Throughout this paper, we write Cp×q (Rp×q, resp.) for the set of all p×q complex (real, resp.) matrices. We use the symbol C+ to stand for the open upper half complex plane. For a Hermitian matrix A=ACp×p, the number of negative eigenvalues (counting algebraic multiplicities) of A is denoted by ν(A). For convenience, the symbol Rn[z] (R0n[z], resp.) represents the set of all real coefficient polynomials of degree at most n (real coefficient polynomials of degree n, resp.) of the variable z. For a rational function r(z) of the form r(z)=p(z)/q(z), in which p(z),q(z) are nonzero complex polynomials such that gcd(p(z),q(z))=1, the McMillan degree of r(z) is defined by degr(z)=max{degp(z),degq(z)}.

    Let f(z) be a function meromorphic on CR. The domain of f(z) is denoted by D(f). For a nonnegative integer κ, f(z) is called a generalized Nevanlinna function with negative index κ if the following statements hold: (i) f(z) satisfies the symmetry condition f(¯z)=¯f(z); (ii) for each choice of a positive integer m and m distinct points z1,,zmC+D(f), we have

    ν(Pf(z1,,zm))κ, (1.1)

    and, for some particular choice, the equality in (1.1) holds, where

    Pf(z1,,zm)=(f(zi)¯f(zj)zi¯zj)mi,j=1

    is a Hermitian Loewner matrix. We denote by Nκ the class of all generalized Nevanlinna functions with negative index κ.

    The indefinite Hamburger moment problem in the class Nκ (short for the HM(Nκ) problem) can be formulated in the following manner: Given a sequence of real numbers s0,,s2n2, it is required to find all functions f(z)Nκ such that the following asymptotic expansion at infinity

    f(z)=s0zs1z2s2n2z2n1+o(z2n+1) (1.2)

    holds when z tends to in the sector πϵ(0)={zCϵargzπϵ} (0<ϵ<π/2).

    The classical Hamburger moment problem (e.g., [1,15]) is identical to the HM(Nκ) problem with κ=0. In comparison with the classical case, the HM(Nκ) problem with κ>0 seems to be much more complicated. In 2003, Derevyagin [7] applied the step-by-step Schur algorithm to give a description of the solutions of the HM(Nκ) problem when the Hankel matrix

    H=(si+j)n1i,j=0 (1.3)

    determined by the asymptotic expansion (1.2) is nonsingular. In 2012, Derkach et al. [8] gave the solvability criterion of the HM(Nκ) problem and a complete parametrization description of the solutions by using the same algorithm. In this paper, we derive the solvability criterion for the HM(Nκ) problem by a new approach, which is more algebraic and different from the existing methods. As a by-product of this approach, we obtain a concrete rational solution of the HM(Nκ) problem with the least McMillan degree when the solvability conditions are met.

    We remark that, starting from the Hankel matrix H given by (1.3), we can derive the solvability criterion of the HM(Nκ) problem and the concrete formula of the solutions when the solvability conditions are met. For this reason, we say H is the Hankel matrix of the HM(Nκ) problem. Moreover, the HM(Nκ) problem is said to be non-degenerate (degenerate, resp.) if its Hankel matrix is nonsingular (singular, resp.). In this paper, we divide the HM(Nκ) problem into the non-degenerate case and the degenerate case to derive the solvability criterion by using the structural characteristics of the Hankel matrix H, such as the characteristic degrees and characteristic polynomial quadruple (see, e.g., [4,5]), the quasidirect decomposition (see, e.g., [9]), and the relation among the Hankel, Loewner, Bezout and some other structured matrices (see, e.g., [4,6,11]).

    A brief synopsis of this paper is as follows. In Section 2, we introduce the characteristic degrees, characteristic polynomial quadruple and quasidirect decomposition of the Hankel matrix of the HM(Nκ) problem, and we list some basic results about these structural characteristics without proofs. In Section 3, we first recall some known properties of the generalized Nevanlinna functions in [17] given by the first two authors of this paper and their collaborators, and then, we prove several new properties by using the structural characteristics of the Hankel matrix and the relation among the Hankel, Loewner, Bezout and some other structured matrices. The last section is devoted to the solvability criterion of the HM(Nκ) problem and a concrete rational solution with the least McMillan degree when the solvability conditions are met.

    Let H=(si+j)n1i,j=0Rn×n be the Hankel matrix of the HM(Nκ) problem. The first and the second characteristic degrees of H are defined by n1=rankH, n2=2nn1, respectively. Clearly, n1nn2. For a pair of positive integers k,l such that k+l=2n, we write Hkl=(si+j)k1,l1i,j=0. In particular, H=Hnn. It follows from [13] that rank Hkl=min{k,l,n1}. We use the symbol Al to stand for the subspace of Rl1[z]:

    Al=(1,z,,zl1)KerHkl,k+l=2n.

    In the case n1=n2=n, H is nonsingular, and there exists a unique monic polynomial p(z)Rn1[z] and a polynomial q(z)R0n[z], forming a basis of the space An+1. In the case n1<n<n2, H is singular and there exists uniquely a monic polynomial p(z)Rn1[z], forming a basis of An1+1 and, moreover, a polynomial q(z)Rn2[z] such that p(z),zp(z),,zn2n1p(z),q(z) forms a basis of the space An2+1. For convenience, in this paper we always assume that either degp(z)=n1 and degq(z)<n2 or degp(z)<n1 and degq(z)=n2. Such a pair of polynomials p(z) and q(z) are referred to as the first and the second characteristic polynomials of H, respectively.

    Let p(z)=pn1zn1++p0 and q(z)=qn2zn2++q0. We define two real coefficient polynomials γ(z) and δ(z) by

    γ(z)=(1,z,,zn11)(p1pn1...pn1)(s0sn11)Rn11[z],δ(z)=(1,z,,zn21)(q1qn2...qn2)(s0sn21)Rn21[z].

    Hereafter, [p(z),q(z),γ(z),δ(z)] is called the characteristic polynomial quadruple of H. Such a quadruple, together with the first and second characteristic degrees, plays an important role in our discussion. By using the definitions, we check easily that the following asymptotic expansions at infinity hold:

    p(z)(s0z+s1z2++s2n2z2n1+o(z2n+1))=γ(z)+o(zn2+1)(z),q(z)(s0z+s1z2++s2n2z2n1+o(z2n+1))=δ(z)+o(zn1+1)(z).

    The following lemma comes as a direct consequence of the definitions and the last two equations (e.g., [4]).

    Lemma 2.1. Let n1 and [p(z),q(z),γ(z),δ(z)] be the first characteristic degree and the characteristic polynomial quadruple of the Hankel matrix H=(si+j)n1i,j=0 given by (1.3), respectively. Then,

    δ(z)p(z)γ(z)q(z)=σ, (2.1)

    where σ is a nonzero constant. Moreover, if degp(z)=n1, then fγ,p(z)=γ(z)/p(z) admits the following asymptotic expansion at infinity:

    fγ,p(z)=s0zs1z2s2n2z2n1+o(z2n+1)(z); (2.2)

    if degp(z)<n1 then fδ,q(z)=δ(z)/q(z) admits the following asymptotic expansion at infinity:

    fδ,q(z)=s0zs1z2s2n2z2n1+o(z2n+1)(z). (2.3)

    Remark 2.2. By (2.1) and the fact that degp(z)+degq(z)2n1, we can show that the asymptotic expansions (2.2) and (2.3) cannot hold simultaneously. This means that if the asymptotic expansion (2.2) ((2.3), resp.) holds, then degp(z)=n1 (degp(z)<n1, resp.).

    Now, we introduce the quasidirect decomposition of a singular Hankel matrix. Let H=(si+j)n1i,j=0 given by (1.3) be singular. We say H is a proper Hankel matrix if Δn10 and Δk=0 (k=n1+1,,n), in which Δi stands for the ith leading principle minor of H. Moreover, we say that H is a degenerate Hankel matrix if s0=s1==sn1=0. By definitions, a n×n zero matrix is both a proper Hankel matrix and a degenerate Hankel matrix. In [9], Fielder showed that under certain conditions, each singular Hankel matrix can be uniquely decomposed into the sum of a proper Hankel matrix and a degenerate Hankel matrix.

    Lemma 2.3. Let H=(si+j)n1i,j=0 given by (1.3) be singular. Then, H has a unique decomposition of the form:

    H=Hp+Hd,rank(H)=rank(Hp)+rank(Hd), (2.4)

    in which Hp is a proper Hankel matrix and Hd is a degenerate Hankel matrix.

    The formula (2.4) is called the quasidirect decomposition of the Hankel matrix H. By definition, together with Lemmas 2.1, 2.3 and Remark 2.2, we can give a characterization of the proper Hankel matrix by using the structural characteristics.

    Lemma 2.4. Let H=(si+j)n1i,j=0 given by (1.3) be singular, and let n1 and [p(z),q(z),γ(z),δ(z)] be the first characteristic degree and the characteristic polynomial quadruple of H, respectively. Then, H is a proper Hankel matrix if, and only if, one of the following statements holds:

    (1) degp(z)=n1;

    (2) fγ,p(z)=γ(z)/p(z) admits the asymptotic expansion (2.2) at infinity.

    The following lemma shows that if the Hankel matrix of the HM(Nκ) problem is singular, then it is equivalent and congruent to a block diagonal matrix, where the upper left corner and the lower right corner blocks are Hankel matrices whose orders coincide with the ranks of its proper part and degenerate part, respectively, and the other blocks are zero matrices. Such a structural characteristic of the singular Hankel matrix plays an important role in deducing the solvability criterion of the HM(Nκ) problem for the degenerate case.

    Lemma 2.5. Let H=(si+j)n1i,j=0 be a singular Hankel matrix with the quasidirect decomposition (2.4), let n1 and p(z)=zr+pr1zr1++p1z+p0 be the first characteristic degree and the first characteristic polynomial of H, respectively and let

    Q=(Irp0pr11p0pr11)Rn×n. (2.5)

    Then, Q is a nonsingular matrix satisfying

    QHpQ=diag(ˆHp,0nr),QHdQ=diag(0nn1+r,ˆHd),QHQ=diag(ˆHp,0nn1,ˆHd), (2.6)

    in which

    ˆHp=(si+j)r1i,j=0,ˆHd=(0ˆs1...ˆs1ˆsn1r) (2.7)

    are nonsingular Hankel matrices of small sizes.

    Hereafter, we always assume that the Hankel matrix ˆHp (ˆHd, resp.) does not appear in the case r=0 (r=n1, resp.).

    In [17], Song et al. presented an equivalent definition of the generalized Nevanlinna function with negative index κ in terms of the generalized Loewner matrix. On the basis of this definition, they have shown several properties of such a kind of function. Here, we list two interesting properties of them. The first property is stated as follows, which gives a necessary condition for the HM(Nκ) problem to have a solution.

    Lemma 3.1. [17] Let H=(si+j)n1i,j=0 be the Hankel matrix of the HM(Nκ) problem. If f(z)Nκ admits the asymptotic expansion (1.2) at infinity, then κν(H).

    The second property is actually a generalization of the first one. Starting from this property, we can derive the solvability criterion of the HM(Nκ) problem in the degenerate case.

    Lemma 3.2. [17] If f(z)Nκ admits the asymptotic expansion (1.2) at infinity, then for each positive integer m and m distinct points z1,,zmC+D(f),

    ν(Lf(z1,,zm))κ,

    in which

    Lf(z1,,zm)=(HC(z1)C(zm)C(z1)Pf(z1,,zm)C(zm)), (3.1)

    H=(si+j)n1i,j=0 is the Hankel matrix of the HM(Nκ) problem and

    C(z)=(f(z),z(f(z)+s0z),,zn1(f(z)+s0z+s1z2++sn2zn1)) (3.2)

    for arbitrary zC+D(f).

    We observe that H is a principle submatrix of Lf(z1,,zm). By the interlacing relation between the eigenvalues of a Hermitian matrix and its principle submatrices (e.g., [14, Theorem 4.3.28]), we have ν(H)ν(Lf(z1,,zm)), then Lemma 3.1 is a direct consequence of Lemma 3.2. If the Hankel matrix H of the HM(Nκ) problem is singular, by Lemma 2.3 it has a unique quasidirect decomposition. In that case, applying Lemmas 2.5 and 3.2, we can prove the following property of the functions in the class Nκ.

    Theorem 3.3. Let H=(si+j)n1i,j=0 given by (1.3) be singular, n1 and [p(z),q(z),γ(z),δ(z)] be the first characteristic degree and the characteristic polynomial quadruple of H, respectively and let f(z)Nκ admit the asymptotic expansion (1.2) at infinity. Then, for each positive integer m and m distinct points z1,,zmC+D(f), the structured matrix Lf(z1,,zm) defined by (3.1) and (3.2) is equivalent and congruent to

    ˜Lf(z1,,zm)=(ˆHpE0ˆHdFEFPf(z1,,zm)), (3.3)

    where ˆHp, ˆHd are the same as in (2.7), and

    E=(f(z1)z1(f(z1)+s0z1)zr11(f(z1)+s0z1++sr2zr11)f(zm)zm(f(zm)+s0zm)zs1m(f(zm)+s0zm++sr2zr1m))Cm×r,F=(p(z1)f(z1)+γ(z1)znr11(p(z1)f(z1)+γ(z1))p(zm)f(zm)+γ(zm)znr1m(p(zm)f(zm)+γ(zm)))Cm×(nr). (3.4)

    Proof. Here, we give a proof only for the case m=1. The proof of the case m>1 is completely analogous and, thus, omitted. Let p(z)=zr+pr1zr1++p1z+p0 and ˜Q=diag(Q,1), in which Q is given by (2.5). Clearly, ˜Q is also a nonsingular matrix. It follows from Lemma 2.5 and (2.6) that ˜Lf(z1,,zm)=˜QLf(z1,,zm)˜Q is of the form (3.3), in which

    (E,F)=(f(z1),z1(f(z1)+s0z1),,zn11(f(z1)+s0z1++sn2zn11))Q.

    By a direct calculation, we have

    E=(f(z1),z1(f(z1)+s0z1),,zr11(f(z1)+s0z1++sr2zr11)),

    and the (k+1)-th element of F is

     p0zk1(f(z1)+s0z1++sk1zk1)++pr1zr+k11(f(z1)+s0z1++sr+k2zr+k11)+zr+k1(f(z1)+s0z1++sr+k1zr+k1)= zk1(p(z1)f(z1)+γ(z1)),k=0,1,,nr1,

    then (3.4) holds. Therefore, we complete the proof of Theorem 3.3.

    To derive the solvability criterion of the HM(Nκ) problem, we need some properties of the rational generalized Nevanlinna function. To introduce these properties, we recall the concept of the Bezout matrix (see, e.g., [10,16]). For a pair of complex polynomials a(z),b(z) with the maximal degree n, the Bezout matrix B(a,b) is defined by the bilinear form

    a(z)b(w)a(w)b(z)zw=(1,z,,zn1)B(a,b)(1,w,,wn1)T.

    It is well known that the Bezout matrix has many applications in the theory of system and control (e.g., [2,3,12]). The following result shows that each real rational function is a generalized Nevanlinna function, whose negative index coincides with the number of negative eigenvalues of the Bezout matrix of its denominator and numerator polynomials.

    Theorem 3.4. Let a(z)R0n[z], b(z)Rn[z]. Then, the real rational function fb,a(z)=b(z)/a(z)Nκ, in which κ=ν(B(a,b)).

    Proof. Since a(z),b(z) are real coefficient polynomials, fb,a(z) is meromorphic in CR such that for each zD(fb,a), ¯zD(fb,a) and fb,a(¯z)=¯fb,a(z). Moreover, for each choice of a positive integer m and m distinct points z1,,zmC+D(fb,a), we have

    Pfb,a(z1,,zm)=(fb,a(zi)fb,a(¯zj)zi¯zj)mi,j=1=Λ(a(zi)b(¯zj)b(zi)a(¯zj)zi¯zj)mi,j=1Λ=ΛVB(a,b)VΛ,

    where Λ=diag(a(z1)1,,a(zm)1) is a nonsingular diagonal matrix and V=(zj1i)m,ni,j=1 is a Vandermonde matrix. This implies that ν(Pfb,a(z1,,zm))ν(B(a,b)). Particularly, in the case of m=n, V is also a nonsingular matrix and, thus, ν(Pfb,a(z1,,zn))=ν(B(a,b)). By the definition of generalized Nevanlinna functions in the class Nκ, we have fb,a(z)=b(z)/a(z)Nκ, in which κ=ν(B(a,b)), then the proof of Theorem 3.4 is completed.

    We remark that there are many interesting connections between Bezout and Hankel matrices. The following lemma shows that the Bezout matrix B(a,b) in Theorem 3.4 is equivalent and congruent to a real Hankel matrix generated by the rational function b(z)/a(z) (see, e.g., [6,11] for the general case).

    Lemma 3.5. Let a(z)=ni=0aizi (an0) and b(z)=n1i=0bizi be real coefficient polynomials, and let the asymptotic expansion of b(z)/a(z) at infinity be of the form

    b(z)a(z)=h0z+h1z2++h2n2z2n1+o(z2n+1)(z).

    Then, B(a,b)=S(a)H(a,b)S(a), in which H(a,b)=(hi+j)n1i,j=0 and

    S(a)=(a1an...an)

    is a nonsingular real symmetric matrix.

    By Lemma 3.5, the negative index κ in Theorem 3.4 can be formulated in terms of the number of negative eigenvalues of the Hankel matrix H(a,b)=(hi+j)n1i,j=0.

    Corollary 3.6. Let a(z),b(z) and H(a,b) be the same as in Lemma 3.5. Then, the real rational function fb,a(z)=b(z)/a(z)Nκ, in which κ=ν(H(a,b)).

    In this section, we apply the structural characteristics of the Hankel matrix of the HM(Nκ) problem and the properties of the generalized Nevanlinna functions to deduce the solvability criterion of the HM(Nκ) problem and a concrete rational solution with the least McMillan degree for both the non-degenerate and degenerate cases. We first derive the solvability criterion for the non-degenerate HM(Nκ) problem.

    Theorem 4.1. Let H=(si+j)n1i,j=0 given by (1.3) be nonsingular. Then, the HM(Nκ) problem is solvable if, and only if, κν(H).

    Proof. The "only if" part is a direct consequence of Lemma 3.1. For the proof of the "if" part, we consider two cases.

    Case I: κ=ν(H). Let n1,n2 and [p(z),q(z),γ(z),δ(z)] be the first characteristic degree, the second characteristic degree and the characteristic polynomial quadruple of H, respectively. Since H is nonsingular, we have n1=n2=n, degp(z)<n1=n and degq(z)=n. By Lemma 2.1, the rational function fδ,q(z)=δ(z)/q(z) admits the asymptotic expansion (1.2) at infinity, and moreover, H=H(q,δ). On the other hand, by Corollary 3.6, we obtain that fδ,q(z)Nκ, in which κ=ν(H(q,δ))=ν(H)=κ. Then, fδ,q(z) is a solution of the non-degenerate HM(Nκ) problem.

    Case II: κ>ν(H). For convenience, we denote by Hn=H and m=κν(H)>0. Define a sequence of Hankel matrices recursively by

    Hn+k=(si+j)n+k1i,j=0, (4.1)

    in which

    s2n+2k3=0,  s2n+2k2=(sn+k,,s2n+2k3)H1n+k1(sn+k,,s2n+2k3)T1,k=1,2,,m. (4.2)

    Then, Hn+m=(si+j)n+m1i,j=0 is equivalent and congruent to diag(H,Im), which implies that Hn+m is a nonsingular Hankel matrix and ν(Hn+m)=ν(H)+m=κ. Let [u(z),v(z),α(z),β(z)] be the characteristic polynomial quadruple of Hn+m. Then degu(z)<n+m and degv(z)=n+m. According to the analysis in Case I, we have that fβ,v(z)Nκ and the asymptotic expansion of fβ,v(z) at infinity is of the form:

    fβ,v(z)=s0zs1z2s2n+2m2z2n+2m1+o(z2n2m+1)(z),

    then fβ,v(z) is a solution of the HM(Nκ) problem. Summarizing the analysis above, we complete the proof of the "if" part.

    From the proof of Theorem 4.1, we obtain immediately a concrete rational solution of the non-degenerate HM(Nκ) problem with the least McMillan degree when the solvability conditions are met.

    Theorem 4.2. Let H=(si+j)n1i,j=0 given by (1.3) be nonsingular. If κν(H), then r(z)=β(z)/v(z) is a rational solution with the least McMillan degree among all rational solutions of the HM(Nκ) problem, in which [u(z),v(z),α(z),β(z)] is the characteristic quadruple of the Hankel matrix Hn+κν(H) defined recursively by (4.1)(4.2).

    Proof. Let m=κν(H)0. We can easily see from the proof of Theorem 4.1 that r(z) presented in Theorem 4.2 is a rational solution of the HM(Nκ) problem with McMillan degree n+m. It remains to prove that the HM(Nκ) problem has not any rational solution whose McMillan degree is less than n+m. If there exists a rational solution s(z) of the non-degenerate HM(Nκ) problem such that degs(z)<n+m, then s(z)=b(z)/a(z), in which a(z),b(z) are two co-prime real polynomials satisfying degb(z)<dega(z)=t<n+m. Let a(z)=a0+a1z++atzt and

    b(z)a(z)=s0zs1z2s2n+2m2z2n+2m1+o(z2n2m+1)(z),

    in which si=si, i=0,1,,2n2. If t<n, then H(a0,,at,0,,0)T=0. It contradicts to the nonsingularity of H. If nt<n+m, we denote by Ht=(si+j)t1i,j=0. In this case, H is a nonsingular principle submatrix of Ht, and then ν(Ht)ν(H)+tn<ν(H)+m=κ. On the other hand, by Corollary 3.6, we have κ=ν(Ht). It is a contradiction, so r(z) given in Theorem 4.2 is a rational solution of the HM(Nκ) problem with the least McMillan degree.

    To derive the solvability criterion for the HM(Nκ) problem for the degenerate case, we need the following result, which can be verified by a direct computation.

    Lemma 4.3. Let ACm×m be nonsingular, B=BCm×m and

    C=(0AAB)C2m×2m.

    Then, ν(C)=m.

    Now, we apply the results above to deduce the solvability criterion of the HM(Nκ) problem for the degenerate case.

    Theorem 4.4. Let H=(si+j)n1i,j=0 given by (1.3) be singular and n1 be the first characteristic degree of H. Then, the HM(Nκ) problem is solvable if, and only if, one of the following statements holds:

    (1) κ=ν(H) and H is a proper Hankel matrix;

    (2) κν(H)+nn1.

    Proof. Let n1, n2 and [p(z),q(z),γ(z),δ(z)] be the first characteristic degree, the second characteristic degree and the characteristic polynomial quadruple of H, respectively, and let degp(z)=rn1. We first prove the "only if" part. Suppose that f(z)Nκ is a solution of the degenerate HM(Nκ) problem. By Lemma 3.1, we have κν(H). In the case of κ=ν(H), it follows from Lemma 3.2 that for each zC+D(f), we have ν(Lf(z))κ, in which

    Lf(z)=(HC(z)C(z)Pf(z)),

    and C(z) is the same as in (3.2). Since H is a principle submatrix of Lf(z), we have ν(Lf(z))ν(H)=κ and thus ν(Lf(z))=κ for all zC+D(f). By Theorem 3.3, Lf(z) is equivalent and congruent to

    ˜Lf(z)=(ˆHpE0nn1ˆHdFEFPf(z)),

    in which ˆHp, ˆHd are the same as in (2.7), and

    E=(f(z),z(f(z)+s0z),,zr1(f(z)+s0z++sr2zr1)),F=(p(z)f(z)+γ(z),,znr1(p(z)f(z)+γ(z))).

    We check easily that ˜Lf(z) is furtherly equivalent and congruent to the following block diagonal matrix:

    ˆLf(z)=diag(ˆHp,ˆHd,0nn11,(0¯d(z)d(z)e(z))),

    in which d(z)=p(z)f(z)+γ(z) and e(z)=¯e(z). Note that ν(Lf(z))=ν(˜Lf(z))=ν(ˆLf(z))=κ=ν(H)=ν(ˆHp)+ν(ˆHd), then by Lemma 4.3, we have d(z)=p(z)f(z)+γ(z)=0, zC+D(f). This implies that f(z)=fγ,p(z)=γ(z)/p(z) for a sufficiently large |z| and zC+D(f). Since f(z) admits the asymptotic expansion (1.2) at infinity, we have

    fγ,p(z)=γ(z)p(z)=s0zs1z2s2n2z2n1+o(z2n+1)(z). (4.3)

    By Lemma 2.4, we have degp(z)=n1 and, thus, H=Hp is a proper Hankel matrix.

    When κ>ν(H), f(z)fγ,p(z)=γ(z)/p(z). Otherwise, fγ,p(z)Nκ and the asymptotic expansion (4.3) holds. In this case, by Lemma 2.4, H=Hp is a proper Hankel matrix. Moreover, by Corollary 3.6, we have κ=ν(H(p,γ))=ν(ˆHp)=ν(Hp)=ν(H), which contradicts the assumption κ>ν(H). We write m=nn1 and g(z)=p(z)f(z)+γ(z) for short, then there exist m distinct points z1,,zmC+D(f) such that g(zk)0, k=1,,m. By Lemma 3.2 and Theorem 3.3, we have ν(˜Lf(z1,,zm))κ, in which

    ˜Lf(z1,,zm)=(ˆHp0000m0A00ˆHdAPf(z1,,zm))C(n+m)×(n+m),
    A=(g(z1)zm11g(z1)g(zm)zm1mg(zm))Cm×m.

    We check easily that ˜Lf(z1,,zm) is equivalent and congruent to the block diagonal matrix diag(ˆHp,ˆHd,D), where

    D=(OAAB)

    and B=BCm×m. Since detA=g(z1)g(zm)1i<jm(zjzi)0, A is a nonsingular m×m matrix. Applying Lemma 4.3, we obtain that ν(Lf(z1,,zm))=ν(ˆHp)+ν(ˆHd)+m=ν(H)+nn1. Therefore, κν(H)+nn1, as needed.

    Let us turn to prove the "if" part. First we suppose that κ=ν(H) and H is a proper Hankel matrix. Then, degp(z)=n1<n, and by Lemma 2.1, fγ,p(z)=γ(z)/p(z) admits the asymptotic expansion (1.2) at infinity. By Corollary 3.6, we have fγ,p(z)Nκ, in which κ=ν(H(p,γ))=ν(ˆHp)=ν(Hp)=ν(H)=κ, then fγ,p(z) is a solution of the degenerate HM(Nκ) problem.

    Now, we suppose that κν(H)+nn1. The proof is divided into two cases.

    Case I: degp(z)=n1. In this case, by Lemma 2.4, H=Hp is a proper Hankel matrix, and the asymptotic expansion (2.2) holds. Assume that

    fγ,p(z)=s0zs1z2s2n2z2n1s2n22z2n21+o(z2n2+1)(z). (4.4)

    We define an n2×n2 Hankel matrix by

    Hn2=(˜si+j)n21i,j=0,˜sk={sk,k2n1;sk+1,k=2n1. (4.5)

    Analogous to the proof of Lemma 2.5, there exists a nonsingular matrix ˜Q of order n2 such that

    ˜QHn2˜Q=(Hn10000A0AB),

    in which

    Hn1=(si+j)n11i,j=0,A=(01...1)C(nn1)×(nn1).

    Clearly, the Hankel matrix Hn2 defined by (4.4) and (4.5) is nonsingular. Moreover, by Lemma 4.3, we have ν(Hn2)=ν(Hn1)+nn1=ν(H)+nn1κ, then by Theorem 4.1, there exists a function f(z)Nκ such that the following asymptotic expansion at infinity

    f(z)= ˜s0z˜s1z2˜s2n2z2n1˜s2n22z2n21+o(z2n2+1)

    holds when z tends to in the sector πϵ(0). This means that f(z) is a solution of the degenerate HM(Nκ) problem.

    Case II: degp(z)<n1. By Lemma 2.1, fd,q(z)=δ(z)/q(z) admits the asymptotic expansion (2.3) at infinity. We assume that

    fd,q(z)=s0zs1z2s2n2z2n1s2n22z2n21+o(z2n2+1)(z). (4.6)

    We define an n2×n2 Hankel matrix by

    Hn2=(si+j)n21i,j=0. (4.7)

    Analogous to the proof of Lemma 2.5, there exists a nonsingular matrix ˜Q of order n2 such that

    ˜QHn2˜Q=(ˆHp00000nn10A00ˆHdB0ABC),

    in which ˆHp and ˆHd are the same as in (2.7), and

    A=(Oˆs1...ˆs1ˆsnn1) (ˆs10),B=(ˆsi+j)n2n,ti,j=1,C=(ˆsi+j+t)n2ni,j=1

    are real Hankel matrices. We check easily that the Hankel matrix Hn2 defined by (4.6) and (4.7) is nonsingular. Moreover, it is furtherly equivalent and congruent to the following block diagonal matrix:

    Λ=diag(ˆHp,ˆHd,(0nn1AA˜C)),

    in which ˜C is a real symmetric matrix of order n2n. Applying Lemma 4.3, we obtain that ν(Hn2)=ν(Λ)=ν(ˆHp)+ν(ˆHd)+nn1=ν(H)+nn1κ, then by Theorem 4.1, there exists a function f(z)Nκ such that the following asymptotic expansion at infinity

    f(z)= s0zs1z2s2n2z2n1s2n22z2n21+o(z2n2+1)

    holds when z tends to in the sector πϵ(0). This implies that f(z) is a solution of the degenerate HM(Nκ) problem, and the proof of Theorem 4.4 is completed.

    From the proofs of Theorems 4.1 and 4.4, we can obtain a concrete rational solution with the least McMillan degree among all rational solutions of the degenerate HM(Nκ) problem, when the solvability conditions are met.

    Theorem 4.5. Let H=(si+j)n1i,j=0 given by (1.3) be singular, and let n1,n2 and [p(z),q(z),γ(z),δ(z)] be the first characteristic degree, the second characteristic degree and the characteristic quadruple of H, respectively.

    (1) If κ=ν(H) and H is a proper Hankel matrix, then r(z)=γ(z)/p(z) is a rational solution with the least McMillan degree among all rational solutions of the HM(Nκ) problem;

    (2) If κν(H)+nn1, then r(z)=β(z)/v(z) is a rational solution with the least McMillan degree among all rational solutions of the HM(Nκ) problem, where [u(z),v(z),α(z),β(z)] is the characteristic quadruple of the Hankel matrix Hκ+nν(H)=(si+j)κ+nν(H)1i,j=0 defined by (4.4) if κ+nν(H)=n2 and degp(z)=n1, by (4.6) if κ+nν(H)=n2 and degp(z)<n1, and defined recursively by

    s2n2+2k3=0,s2n2+2k2=(sn2+k,,s2n2+2k3)H1n2+k1(sn2+k,,s2n2+2k3)T1,  k=1,2,

    if κ+nν(H)>n2.

    Proof. We see from the proofs of Theorems 4.1 and 4.4 that the rational function r(z) presented in Theorem 4.5 for each case is a solution of the HM(Nκ) problem. Now it remains to prove that such a rational solution has the least McMillan degree among all rational solutions of the HM(Nκ) problem, when the solvability conditions are met. We divide the proof into two cases.

    Case I: κ=ν(H) and H is a proper Hankel matrix. In this case, degr(z)=degp(z)=n1 and ˆHp=(si+j)n11i,j=0 is nonsingular. Assume that the HM(Nκ) problem has a rational solution f(z) such that degf(z)<n1, then there exist two co-prime real polynomials a(z),b(z) such that f(z)=b(z)/a(z), dega(z)n11. Let a(z)=a0+a1z++an11zn11. In view of the fact that

    f(z)=b(z)a(z)=s0zs1z2s2n2z2n1+o(z2n+1),

    we have ˆHp(a0,a1,,an11)T=0. It contradicts to the nonsingularity of ˆHp, then r(z)=γ(z)/p(z) is a rational solution of the HM(Nκ) problem with the least McMillan degree n1.

    Case II: κν(H)+nn1. In this case, degr(z)=degv(z)=κ+nν(H). If the HM(Nκ) problem has a rational solution f(z) such that degf(z)=t<κ+nν(H), then there exist two coprime real polynomials a(z),b(z) such that f(z)=b(z)/a(z) and dega(z)=t. Let

    f(z)=b(z)a(z)=s0zs1z2s2k2z2k1+(z)

    in which si=si, i=1,,2n2, and let Hk=(si+j)k1i,j=0, k=1,2,. By Corollary 3.6, we have ν(Ht)=κν(H)+nn1. If t<n2, then Hn2 is singular and, thus, ν(Ht)ν(Hn2)<ν(H)+(n2n1)/2=ν(H)+nn1. It is a contradiction. Therefore, tn2 and Hn2 is a principle submatrix of Ht. By the interlacing relation between the eigenvalues of a Hermitian matrix and its principle submatrices (see, e.g., [14, Theorem 4.3.28]),

    ν(Ht)ν(Hn2)+tn2ν(H)+n2n12+tn2=ν(H)+tn<κ.

    It is also a contradiction, and r(z)=β(z)/v(z) is a rational solution of the HM(Nκ) problem with the least McMillan degree κ+nν(H).

    In this paper, we introduced some basic structural characteristics of the Hankel matrix, such as the first and second characteristic degrees, the characteristic polynomial quadruple and the quasidirect decomposition for the singular case, and then, we applied these structural characteristics and the relation among the Hankel, Loewner, Bezout and some other structured matrices to deduce several new properties of the functions in the class Nκ and the solvability criterion of the HM(Nκ) problem for both the non-degenerate and degenerate cases. As a by-product, we simultaneously obtained a rational solution of the HM(Nκ) problem with the least McMillan degree when the solvability conditions were met.

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

    The corresponding author is supported by the Scientific Research Fund from Beijing Normal University at Zhuhai (Grant No. 111032119).

    Yongjian Hu is the Guest Editor of special issue "Matrix theory and its applications" for AIMS Mathematics. Yongjian Hu was not involved in the editorial review and the decision to publish this article. All authors declare no conflicts of interest in this paper.



    [1] P. Li, S. T. Yau, On the parabolic kernel of the Schrödinger operator, Acta Math., 156 (1986), 153–201. https://doi.org/10.1007/BF02399203 doi: 10.1007/BF02399203
    [2] R. S. Hamilton, A matrix Harnack estimate for the heat equation, Comm. Anal. Geom., 1 (1993), 113–126. https://doi.org/10.1090/S0002-9939-99-04967-9 doi: 10.1090/S0002-9939-99-04967-9
    [3] G. Perelman, The entropy formula for the Ricci flow and its geometric applications, arXiv preprint math/0211159, 2002.
    [4] J. F. Li, X. J. Xu, Differential Harnack inequalities on Riemannian manifolds I: linear heat equation, Adv. Math., 226 (2011), 4456–4491. https://doi.org/10.1016/j.aim.2010.12.009 doi: 10.1016/j.aim.2010.12.009
    [5] B. Peng, Y. D. Wang, G. D. Wei, Yau type gradient estimates for Δu+au(logu)p+bu=0 on Riemannian manifolds, J. Math. Anal. Appl., 498 (2021), 124963. https://doi.org/10.1016/j.jmaa.2021.124963 doi: 10.1016/j.jmaa.2021.124963
    [6] J. Y. Wu, Elliptic gradient estimates for a nonlinear heat equation and applications, Nonlinear Anal., 151 (2017), 1–17. https://doi.org/10.1016/j.na.2016.11.014 doi: 10.1016/j.na.2016.11.014
    [7] L. V. Dai, N. T. Dung, N. D. Tuyen, L. Zhao, Gradient estimates for weighted p-Laplacian equations on Riemannian manifolds with a Sobolev inequality and integral Ricci bounds, Kodai Math. J., 45 (2022), 19–37. https://doi.org/10.2996/kmj/kmj45102 doi: 10.2996/kmj/kmj45102
    [8] X. J. Xu, Gradient estimates for ut=ΔF(u) on manifolds and some Liouvilletype theorems, J. Differential Equations, 252 (2012), 1403–1420. https://doi.org/10.1016/j.jde.2011.08.004 doi: 10.1016/j.jde.2011.08.004
    [9] X. B. Zhu, Gradient estimates and Liouville theorems for nonlinear parabolic equations on noncompact Riemannian manifolds, Nonlinear Anal., 74 (2011), 5141–5146. https://doi.org/10.1016/j.na.2011.05.008 doi: 10.1016/j.na.2011.05.008
    [10] Y. Li, X. R. Zhu, Harnack estimates for a heat-type equation under the Ricci flow, J. Differential Equations, 260 (2016), 3270–3301. https://doi.org/10.1016/j.jde.2015.10.024 doi: 10.1016/j.jde.2015.10.024
    [11] I. P. de Jesus, A. M. Oliveira, M. R. Clark, P. P. A. Oliveira, Controllability for semilinear heat equation with globally Lipschitz nonlinearities and memory term, Nonlinear Anal. Real World Appl., 60 (2021), 103277. https://doi.org/10.1016/j.nonrwa.2020.103277 doi: 10.1016/j.nonrwa.2020.103277
    [12] I. P. de Jesus, S. B. de Menezes, On the approximate controllability of Stackelberg-Nash strategies for linear heat equations in RN with potentials, Appl. Anal., 94 (2015), 780–799. https://doi.org/10.1080/00036811.2014.901503 doi: 10.1080/00036811.2014.901503
    [13] J. P. Wang, Global heat kernel estimates, Pacific J. Math., 178 (1997), 377–398. https://doi.org/10.2140/pjm.1997.178.377 doi: 10.2140/pjm.1997.178.377
    [14] R. Chen, Neumann eigenvalue estimate on a compact Riemannian manifold, Proc. Amer. Math. Soc., 108 (1990), 961–970. https://doi.org/10.2307/2047954 doi: 10.2307/2047954
    [15] P. Petersen, G. F. Wei, Relative volume comparison with integral curvature bounds, Geom. Funct. Anal., 7 (1997), 1031–1045. https://doi.org/10.1007/s000390050036 doi: 10.1007/s000390050036
    [16] D. Yang, Convergence of Riemannian manifolds with integral bounds on curvature. I, Ann. Sci. École Norm. Sup., 25 (1992), 77–105. https://doi.org/10.24033/asens.1644
    [17] P. Petersen, G. F. Wei, Analysis and geometry on manifolds with integral Ricci curvature bounds. II, Trans. Amer. Math. Soc., 353 (2001), 457–478. https://doi.org/10.1090/S0002-9947-00-02621-0 doi: 10.1090/S0002-9947-00-02621-0
    [18] Q. S. Zhang, M. Zhu, Li-Yau gradient bounds on compact manifolds under nearly optimal curvature conditions, J. Funct. Anal., 275 (2018), 478–515. https://doi.org/10.1016/J.JFA.2018.02.001 doi: 10.1016/J.JFA.2018.02.001
    [19] Q. S. Zhang, M. Zhu, Li-Yau gradient bound for collapsing manifolds under integral curvature condition, Proc. Amer. Math. Soc., 145 (2017), 3117–3126. https://doi.org/10.1090/proc/13418 doi: 10.1090/proc/13418
    [20] X. Z. Dai, G. F. Wei, Z. L. Zhang, Local Sobolev constant estimate for integral Ricci curvature bounds, Adv. Math., 325 (2018), 1–33. https://doi.org/10.1016/j.aim.2017.11.024 doi: 10.1016/j.aim.2017.11.024
    [21] W. Wang, Harnack inequality, heat kernel bounds and eigenvalue estimates under integral Ricci curvature bounds, J. Differential Equations, 269 (2020), 1243–1277. https://doi.org/10.1016/j.jde.2020.01.003 doi: 10.1016/j.jde.2020.01.003
    [22] W. Wang, Elliptic-type gradient estimates under integral Ricci curvature bounds, Proc. Amer. Math. Soc., 150 (2022), 4965–4979. https://doi.org/10.1090/proc/14774 doi: 10.1090/proc/14774
    [23] J. H. Wang, Elliptic gradient estimate for heat equation with integral Ricci curvature condition, Acta Math. Sinica (Chinese Ser.), 65 (2022), 763–774. https://doi.org/10.12386/A20210040 doi: 10.12386/A20210040
    [24] J. Wang, Y. D. Wang, Gradient estimates for Δu+a(x)ulogu+b(x)u=0 and its parabolic countepart under integral Ricci curvature bounds, arXiv preprint arXiv: 2109.05235, 2022.
    [25] X. R. Olivé, Neumann Li-Yau gradient estimate under integral Ricci curvature bounds, Proc. Amer. Math. Soc., 147 (2019), 411–426. https://doi.org/10.1090/proc/14213 doi: 10.1090/proc/14213
    [26] C. Rose, Heat kernel upper bound on Riemannian manifolds with locally uniform Ricci curvature integral bounds, J. Geom. Anal., 27 (2017), 1737–1750. https://doi.org/10.1007/S12220-016-9738-3 doi: 10.1007/S12220-016-9738-3
    [27] M. Choulli, L. Kayser, E. M. Ouhabaz, Observations on Gaussian upper bounds for Neumann heat kernels, Bull. Aust. Math. Soc., 92 (2015), 429–439. https://doi.org/10.1017/S0004972715000611 doi: 10.1017/S0004972715000611
    [28] S. Gallot, Isoperimetric inequalities based on integral norms of Ricci curvature, Astérisque, 157 (1988), 191–216.
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1036) PDF downloads(57) Cited by(0)

Figures and Tables

Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog