We establish a new on-diagonal lower estimate of continuous-time heat kernels for large time on graphs. To achieve the goal, we first introduce an upper estimate of heat kernels in natural graph metric, then use the upper estimate and the volume growth condition to show the validity of the on-diagonal lower estimate.

Let G = (V, E) be a locally finite graph, whose measure μ(x) has positive lower bound, and Δ be the usual graph Laplacian. Applying the mountain-pass theorem due to Ambrosetti and Rabinowitz (1973), we establish existence results for some nonlinear equations, namely Δu + hu = f(x, u), x ∈ V. In particular, we prove that if h and f satisfy certain assumptions, then the above-mentioned equation has strictly positive solutions. Also, we consider existence of positive solutions of the perturbed equation Δu + hu = f(x, u) + ϵg. Similar problems have been extensively studied on the Euclidean space as well as on Riemannian manifolds.

We derive a gradient estimate for positive functions, in particular for positive solutions to the heat equation, on finite or locally finite graphs. Unlike the well known Li-Yau estimate, which is based on the maximum principle, our estimate follows from the graph structure of the gradient form and the Laplacian operator. Though our assumption on graphs is slightly stronger than that of Bauer et al. (J Differ Geom 99:359–405, 2015), our estimate can be easily applied to nonlinear differential equations, as well as differential inequalities. As applications, we estimate the greatest lower bound of Cheng’s eigenvalue and an upper bound of the minimal heat kernel, which is recently studied by Bauer et al. (Preprint, 2015) by the Li-Yau estimate. Moreover, generalizing an earlier result of Lin and Yau (Math Res Lett 17:343–356, 2010), we derive a lower bound of nonzero eigenvalues by our gradient estimate.

Let G=(V,E) be a connected finite graph and Δ be the usual graph Laplacian. Using the calculus of variations and a method of upper and lower solutions, we give various conditions such that the Kazdan–Warner equation Δu=c−he^u has a solution on V, where c is a constant, and h:V→R is a function. We also consider similar equations involving higher order derivatives on graph. Our results can be compared with the original manifold case of Kazdan and Warner (Ann. Math. 99(1):14–47, 1974).

Continuing our previous work (arXiv:1509.07981v1), we derive another global gradient estimate for positive functions, particularly for positive solutions to the heat equation on finite or locally finite graphs. In general, the gradient estimate in the present paper is independent of our previous one. As applications, it can be used to get an upper bound and a lower bound of the heat kernel on locally finite graphs. These global gradient estimates can be compared with the Li–Yau inequality on graphs contributed by Bauer et al. [J. Differential Geom., 99, 359–409 (2015)]. In many topics, such as eigenvalue estimate and heat kernel estimate (not including the Liouville type theorems), replacing the Li–Yau inequality by the global gradient estimate, we can get similar results.