On the heat content functional and its critical domains

We study and classify smooth bounded domains in an analytic Riemannian manifold which are critical for the heat content at all times t>0. We do that by first computing the first variation of the heat content, and then showing that a domain is critical if and only if it has the so-called constant flow property, so that we can use a previous classification result established by the author. The outcome is that a domain is critical for the heat content at all times if and only if it admits an isoparametric foliation, that is, a foliation whose leaves are all parallel to the boundary and have constant mean curvature. Then, we consider the sequence of functionals given by the exit-time moments, which generalize the torsional rigidity. We prove that a domain is critical for all exit time moments if and only it is critical for the heat content at all times, and then we get a classification as well. The main purpose of the paper is to understand the variational properties of general isoparametric foliations and their role in PDE's theory; in some respect isoparametric foliations generalize the properties of the foliation of Euclidean space by round spheres.


Introduction and main results
In this paper we consider smooth bounded domains Ω in a complete Riemannian manifold M , and we study those domains which are critical for the heat content H Ω (t) of the domain at any time fixed t, under deformations which preserve the inner volume (see its definition in Section 1.1). The heat content has been amply studied in the literature and there are many important contributions. Of special relevance to this paper are the following aspects: its asymptotic behaviour for small times (see [2], [3], [29], [30], [14], [15]), comparison theorems of isoperimetric nature ( [5]), its relation with the Dirichlet spectrum and the Brownian motion ( [20], [19]), and the geometric rigidity of an associated overdetermined problem ( [31], [32]). The main scope of this paper is to understand the geometry of critical domains, and possibly give a classification. We observe that that it is easy to prove a classification when Ω is a domain in Euclidean space R n , hyperbolic space H n and the hemisphere S n + . In all these cases the only critical domains are geodesic balls, which in fact are absolute maxima for the heat content, by known comparison results ( [5]). So, our main focus will be on arbitrary ambient Riemannian manifolds M . After computing the formula for the first variation of the heat content, we see that critical domains are characterized by the so-called constant flow property (defined in section 1.2), which basically says that the normal derivative of the temperature function (the heat flow) is constant on the boundary of the domain. The overdetermined heat equation which follows has been studied in the papers [31], [32] and there, under the assumption that the ambient manifold is analytic, it is shown that a domain supports a solution to such overdetermined heat equation if and only if it is an isoparametric tube, that is, it admits a foliation by smooth, parallel, constant mean curvature hypersurfaces, degenerating to a unique (possibly singular) leaf which is a closed, minimal submanifold of M (see Section 1.3). The fact that the only (compact) isoparametric tubes of R n , H n and S n + are the geodesic balls (in which case the leaves are geodesic spheres, and the singular leaf is simply a point, the center of the ball) is a consequence of a classical result due to Cartan, or else, it follows from the classical Alexandrov theorem and from its generalization to H n and S n + (see [21]), to the effect that the only embedded mean curvature hypersurfaces of those spaces are geodesic spheres. That is why the full geometric meaning of criticality is perhaps not seen on those spaces.
It is already seen, however, on the whole sphere, which hosts a reasonably large, and quite interesting, variety of isoparametric tubes which are not geodesic balls: the simplest such family is given by domains bounded by Clifford tori (products of spheres of appropriate radii). More generally (at least when the boundary is connected) it follows from the classical theory that isoparametric tubes in the sphere are precisely those bounded by an isoparametric hypersurface, that is, a hypersurface having constant principal curvatures. For general facts about isoparametric hypersurfaces, see for example [38]. The classification of isoparametric hypersurfaces in the sphere is a classical problem in differential geometry which started from Cartan in the 30's ( [7]) and which, after several important intermediate results, has been completed only recently ( [8]).
Thus, as a result of our findings and thanks to the collective efforts leading to the final classification [8], the family of critical domains in the whole sphere can be completely classified, and is in fact much larger than the family of geodesic balls. Criticality for the heat content, in S n , implies in particular that any component of the boundary is a smooth algebraic variety, more precisely, it is the zero set of the restriction to S n of a harmonic polynomial in R n+1 satisfying precise algebraic conditions (Cartan-Müntzner polynomials). We will briefly summarize the main facts about isoparametric foliations on the sphere in section 1.4. At the moment, there is no classification (in the strict sense of the word) of isoparametric foliations in general, although there is some progress. It is clear however that this condition is quite strict, and that a random manifold will not support any isoparametric foliation.
The probabilistic aspect of the heat content, its relation with the Dirichlet spectrum and with the so-called exit-time moments {T 1 (Ω), T 2 (Ω), . . . } have been studied in several papers ( [9], [10], [20], [19], [16], [17]). We point out that T 1 (Ω) has also an intepretation in mechanics, being the torsional rigidity of the domain Ω; isoperimetric inequalities for the functional T 1 (Ω) are classical, and were recently extended to the higher exit time moments. For example, it is shown in [19] that, in the spaces R n , H n and S n + the only domains which are critical for T k (Ω), for at least one k ∈ N, are geodesic balls. We point out that the same rigidity definitely does not extend to spherical domains not contained in a hemisphere; in fact all isoparametric tubes (and not only geodesic balls) are critical for T k (Ω) for all k ∈ N. So, the problem we want to bring to attention is the following: • classify the domains which are critical for all exit-time moments T k (Ω) when Ω is a domain in an arbitrary Riemannian manifold. The outcome is that this family coincides with the family of domains having the constant flow property, hence also with the family of domains which are critical for the heat content functional at all times: these domains are precisely the isoparametric tubes. At this point, we summarize the main results and state the main theorem of this paper. In what follows, "critical" means "critical under volume preserving deformations".
Let Ω be a smooth bounded domain in a complete Riemannian manifold (M, g). The following are equivalent: a) Ω is critical for the heat content H Ω (t) at every fixed time t > 0. b) Ω is critical for the k-th exit time moment T k (Ω), for all k ≥ 1. c) Ω has the constant flow property. If the ambient Riemannian manifold M is analytic, a), b) and c) are in turn all equivalent to: d) Ω is an isoparametric tube over a closed, minimal submanifold of M . This paper is a natural continuation of [31] and [32], and clarifies the variational meaning of the constant flow property. Its main theme is to stress the relation between variational geometry, overdetermined PDE's and the isoparametric property, which is not fully evident when studying problems in the spaces R n , H n and S n + due to the Alexandrov rigidity theorem. That the isoparametric condition is sufficient to guarantee the existence of solutions to many overdetermined problems was perhaps first discussed in detail by Shklover in [34]. That this property is actually also necessary for the constant flow property was the main outcome of [31] and [32]. The main theorem shows that certain variational problems, and the overdetermined PDE's they generate, naturally lead to the consideration of isoparametric foliations. For previous facts on PDE's and the isoparametric theory, we mention the works [35] and [36], where the spectrum of the Laplace-Beltrami operator was computed for cubic isoparametric minimal hypersurfaces of the Euclidean sphere, and the paper [37] where the first eigenvalue of any closed minimal isoparametric hypersurface of S n+1 was shown to be equal to n, thus confirming Yau's conjecture in these cases. Another important functional is the Dirichlet heat trace, which is known to be a spectral invariant of a domain Ω. For results about its critical domains, we refer to [11]: there it is shown, using the asymptotic expansion of the heat trace and the Alexandrov theorem, that the only critical domains in R n , H n and S n + are the geodesic balls. At the moment, there seems to be no classification in other manifolds, in particular, in the sphere S n of arbitrary dimension. In the rest of this introduction, we give the precise definitions and discuss the problem in greater detail. Proofs will be given starting from Section 2, where we prove the first variation formula of the heat content, Theorem 10, which is the main step.

The heat content
Let Ω be a smooth, bounded domain in a complete Riemannian manifold M . The heat content of Ω is the function H Ω (t) of time t > 0 which measures the total heat inside Ω at time t, assuming that initially the temperature is uniformly distributed and equal to 1, and that the boundary is subject to absolute refrigeration at all times (Dirichlet boundary conditions). The heat content is expressed as where u(t, x) is the temperature function, solution of the relevant heat equation: It is natural to expect that the extrinsic geometry of ∂Ω plays a major role in the asymptotic behavior of H Ω (t) for small times t; in fact, for small times, only the points near the boundary will feel the sudden drop of the temperature due to the boundary refrigeration (in fact, points far away from the boundary obey the so-called principle of not feeling the boundary). This behavior is reflected in an asymptotic series in powers of √ t: H as t → 0 (2) whose coefficients are, for what we have just said, supported on the boundary and depend on the extrinsic curvature of ∂Ω, precisely, the second fundamental form and its covariant derivatives, which in turn involve the ambient curvature tensor together with its covariant derivatives. After some works in Euclidean space, the existence of the asymptotic series in the general Riemannian case, together with the calculation of the coefficients β k (Ω) up to k = 2 was carried out by van den Berg and Gilkey [2]. For more general operators of Laplace-type, with various boundary conditions see the review papers [14] and [15]. An algorithm for the computation of the whole asymptotic series (2) has been developped by the author in [29] and [30]; this approach has been recently adapted in [28] to treat the sub-Riemannian case. On the other hand, for large times the heat content is governed by the Dirichlet spectrum of the domain; in particular, H Ω (t) decays exponentially to zero as t → ∞, with speed proportional to the lowest Dirichlet eigenvalue λ 1 (Ω): with c Ω being the integral over Ω of a unit L 2 -norm first Dirichlet eigenfunction. This can be easily seen by writing the Fourier series of the heat content.

Constant flow property
The function ∂u ∂N (t, ·) : ∂Ω → R is called the (pointwise) heat flow at time t; when integrated on the boundary, it measures the speed at which the domain is losing heat due to boundary refrigeration: where dσ is the Riemannian surface measure. Now, for a general domain the heat flow is not uniform (i.e. constant) across the boundary. If the heat flow is constant on ∂Ω, at every fixed value of time, then we say that Ω has the constant flow property. In other words: Definition. The domain Ω has the constant flow property if there exists a smooth function for all y ∈ ∂Ω.
Adding the condition (4) to the heat equation (1) defining u(t, x), we obtain an overdetermined problem, and this overdetermination imposes strict conditions on the geometry of the domain. In the papers [31] and [32] we actually give a geometric characterization of the domains with the constant flow property when the ambient manifold M is analytic: they are isoparametric tubes (see the definition below) that is, they admit an isoparametric foliation or, in other words, a foliation by parallel hypersurfaces having constant mean curvature. We will discuss this in more details in the next section. Before doing that, we remark that for domains in R n another interesting overdetermined condition on the function u(t, x) was studied in [18] and the works that follow: namely, the existence of a stationary isothermic hypersurface in Ω. Under mild assumptions, it is proved there that existence of such S forces Ω to be a ball.

Isoparametric tubes
Here is the definition.
Definition 2. Let P be a smooth, closed submanifold of the domain Ω. We say that Ω is a smooth tube around P if : a) Ω is the set of points at distance at most R from P , b) For each s ∈ (0, R], the equidistant is a smooth hypersurface of Ω.
We say that the smooth tube Ω is an isoparametric tube if every equidistant Σ s as above has constant mean curvature.
The submanifold P is called the soul of Ω, and can have dimension dim P = 0, . . . , n − 1. The soul is then an embedded submanifold, which is always minimal (see [13]). Any isoparametric tube has at most two boundary components (for the easy proof see [32]). For example, a solid revolution torus in R 3 with radii a > b > 0 is a smooth tube (the soul P is a circle), but is not an isoparametric tube because equidistants have variable mean curvature. In fact: Proposition 3. The only (compact) isoparametric tubes in the spaces R n , H n and S n + are the geodesic balls (in which case the soul is a point).
For the proof, just observe that the Alexandrov theorem holds true in any of these spaces, so that any compact embedded hypersurface must be a geodesic sphere. If the boundary has two components, then by definition they must be spheres with the same (constant) mean curvature (hence the same radius) and must be parallel, which is impossible in the cases at hand. Thus, the boundary consists of one piece, which is then a sphere.
• Observe that a domain in a whole sphere S n , bounded by two geodesic spheres, is an isoparametric tube if and only if the two boundary spheres are isometric and have equal (or antipodal) centers : in that case, the soul is an equatorial (i.e. totally geodesic) hypersurface. In fact, if the centers are neither equal nor antipodal then the region is not even a smooth tube.
• More generally, any geodesic ball in a locally harmonic manifold is (more or less by definition) an isoparametric tube around its center.
• Now assume that the metric of Ω is smooth and that there is a distinguished point p ∈ Ω such that (Ω\{p}, g) is isometric to (0, R]×S n−1 endowed with the metric g = dr 2 +θ 2 (r)g S n−1 (r being the radial parameter). Then we say that Ω is a revolution manifold. Clearly any revolution manifold is an isoparametric tube around its soul, the point p. Note that its boundary ∂Ω has only one component, namely {R} × S n−1 .
In the next subsection we will discuss the main class of examples of isoparametric tubes, namely, spherical domains bounded by isoparametric hypersurfaces.

Isoparametric tubes in the standard sphere
By definition, the boundary of an isoparametric tube, and all of its regular equidistants, are isoparametric hypersurfaces of M , which allows for a classification when M is the round sphere. Let us recall that a hypersurface Σ in a space form R n , H n or S n is isoparametric if it has constant principal curvatures, that is, if the characteristic polynomial of the second fundamental form is the same at all points of Σ. While on R n , H n and the hemisphere S n + the only compact isoparametric hypersurfaces are the geodesic spheres, in S n there is a much larger variety of them. The classification of such hypersurfaces started from Cartan and was a major problem in Differential Geometry, which came to a complete solution only very recently ( [8]). Let us review the main properties of an isoparametric hypersurface Σ of S n .
• The number g of distinct principal curvatures of Σ can be only 1, 2, 3, 4 and 6 ( [22]). The case g = 1 corresponds to the family of geodesic spheres, and g = 2 corresponds to Clifford tori; these are tubes around a totally geodesic submanifold and are hypersurfaces of type: which admit a natural embedding into S n with constant principal curvatures λ = b a (counted p times) and µ = − a b (counted q times). It is a remarkable and perhaps surprising fact that when g = 4 there exists non-homogenous isoparametric hypersurfaces.
• For each Σ there exist two regular, connected submanifolds Σ + , Σ − of S n such that Σ is the surface of the tube with radius r + (resp. r − ) around Σ + (resp. Σ − ). These submanifolds are called the focal submanifolds of Σ, and are minimal in S n .
• Every isoparametric hypersurface Σ belongs to a one-parameter family, which gives rise to what is known to be an isoparametric foliation of S n . This foliation has precisely two singular leaves (the focal submanifolds Σ + and Σ − ) and contains exactly one minimal isoparametric hypersurface: when g = 1 it is the unique equator of the family (which is totally geodesic), and when g = 2, for fixed p and q, it is the minimal Clifford torus defined by the identity pb 2 = qa 2 .
• The work of Münzner [22], [23] established a connection between isoparametric hypersurfaces in S n and homogeneous polynomials with n + 1 variables. Let Σ be isoparametric in S n . Then Σ is a regular level set of the restriction f to S n of a homogeneous polynomial F : R n+1 → R of degree g, where g is the number of distinct principal curvatures of Σ. One has f (S n ) = [−1, 1] and Σ = f −1 (t) for some t ∈ (−1, 1) (the focal submanifolds correspond to the critical values −1 and 1). The polynomial F satisfies the conditions for a suitable rational number c, which is determined by the multiplicity of the principal curvatures of Σ and∆,∇ are the Laplacian and the gradient in R n+1 . The polynomial F is an example of Cartan-Münzner polynomial. Conversely, any Cartan-Münzner polynomial of degree g with c = (n−1)g gives rise to an isoparametric foliation of S n with g distinct principal curvatures. Having that, the geometric classification of isoparametric foliations reduces to the (difficult) algebraic problem of classifying all Cartan-Münzner polynomials. Now we can easily classify the isoparametric tubes in the round sphere. Recall that an isoparametric tube can have at most two boundary components.

Theorem 4.
Let Ω be an isoparametric tube in S n . Then either Ω is the domain bounded by a (connected) isoparametric hypersurface of S n , or Ω is a tube around a minimal isoparametric hypersurface Σ of S n . In the first case the soul is a focal submanifold of ∂Ω, in the second it is the hypersurface Σ (of codimension 1) and the boundary has two components. In low dimensions: Corollary 5. a) An isoparametric tube in S 2 is either a geodesic disk or a tube around an equator. b) An isoparametric tube in S 3 is congruent to one of the following: a geodesic disk, a tube around the equator, a domain bounded by a Clifford torus or a tube around the minimal Clifford torus

Exit time moments
The function E 1 : Ω → R, unique solution of the Dirichlet problem is known in the literature as the torsion function of the domain Ω, and its integral defines the so-called torsional rigidity of Ω. Isoperimetric inequalities for the torsion function and the torsional rigidity are by now classical (see for example [1], [26]). It is also wellknown (but see below) that a domain is critical for torsional rigidity (under volume preserving deformations) if and only if its torsion function has constant normal derivative, that is, if and only if Ω supports a solution to the overdetermined problem known in the literature as Serrin problem. As a consequence of Serrin's rigidity theorem (see [33]) we know that the only domains in R n which are critical for T 1 under volume preserving deformations are balls, and these are all maxima by the classical isoperimetric result by Polya ([25], [26]). The same rigidity holds in the spaces H n and S n + : the only critical domains are geodesic balls (also in this case, these are absolute maxima, see for example [5] and [6]). On the sphere there are many critical domains which are not geodesic balls, for example, domains bounded by any connected isoparametric hypersurface, as proved in [34] (see also [31]). Yet more generally, isoparametric tubes in general Riemannian manifolds are critical for T 1 ( [32]). However, the condition of being critical for T 1 seems to be too weak to guarantee a reasonable classification: see a recent example in [12] of a spherical domain which is critical for T 1 and is not even an isoparametric tube (actually, the boundary has variable mean curvature). We will prove in this paper a classification result under criticality for the whole family of exit time moments, which we are going to define. Now, the function E 1 has also a probabilistic intepretation, being the mean exit time associated to the Brownian motion of Ω. As such, it is part of a hierarchy of exit time moments.Precisely, define the k-th exit time function E k on Ω inductively as follows. We set E 0 = 1 and, for k ≥ 1, we let E k be the unique solution of Note that E 1 is just the torsion function. The k-th exit time moment of the bounded domain Ω is now defined as The sequence is known as the exit time moment spectrum of Ω. These invariants have been studied in the papers [9], [10], [17], [16], [19], [20].
The following expressions hold, for k ≥ 1 (see for example [9]): where H Ω (t) is the heat content of Ω. These identities show the strict relation between the exit time moments and the heat content. The relation between the moment spectrum and the Dirichlet spectrum has been studied in [10], [9] and [16].
In this paper we characterize the geometry of Riemannian domains which are critical for all exit time moments T k (Ω), see Theorem 6 below.

Domain deformations and critical domains
Let V be a smooth vector field defined in a neighborhood U of the domain Ω in M . Define a one-parameter domain deformation f ǫ : Ω → M by: For ǫ small enough, f ǫ restricts to a diffeomorphism: We call Ω ǫ the one-parameter deformation of Ω associated to the vector field V .
Given a geometric functional F = F (Ω) depending smoothly on the domain Ω we define its first variation F ′ (Ω) as follows: Note that F ′ (Ω) depends on the vector field V which defines the deformation, hence it would be more correct to write F ′ (Ω) = DF (Ω, V ), interpreting such expression as the directional derivative of the functional F at Ω in the direction V . We will say that Ω is critical for the functional F if DF (Ω, V ) = 0 for all deformations of Ω hence, for all vector fields V . However, to have a meaningful geometric problem one should impose (and we will) that the deformation f ǫ is volume preserving: for ǫ small enough. To preserve volume, the vector field V must satisfy the condition Ω divV = 0 and then, by Green formula: Under condition (11) we can study the first variation of the heat content when the deformed domains have the same volume.

The main result
Here is the main result of this paper, as stated at the beginning of the introduction. The word critical means really critical under volume preserving deformations.

Remarks
We first recall a comparison result proved by Burchard and Shmuckenschlager [5]: let Ω be a domain in a constant curvature space form (hence, up to homotheties, R n , H n and S n ) and let Ω ⋆ be the geodesic ball with the same volume: |Ω ⋆ | = |Ω|. Then, at all times t > 0 one has Therefore geodesic balls in constant curvature space forms are absolute maxima for the heat content functional, if one restricts to deformations which preserve the inner volume.
We have seen in Theorem 3 that the only isoparametric tubes in R n , H n and S n + are geodesic balls: hence, by our main theorem, we have: The only bounded domains in R n , H n and S n + which are critical for the heat content are geodesic balls, and these are absolute maxima.
Of course the same conclusion holds for the exit time moments; this was first proved in [19] by the Alexandrov-Serrin moving plane argument. We remark that the moving plane method cannot be applied in our general case, and in fact the situation changes drastically already in the whole sphere, as there are many critical domains which are not geodesic balls. This follows from part d) of the main theorem and the classification of the isoparametric tubes given in Theorem 4 of the previous section: Corollary 8. A domain Ω in S n is critical for the heat content at all times t and for the k-th exit time moment T k , for all k, if and only if: Ω is bounded by a (connected) isoparametric hypersurface or is a tube around a minimal isoparametric hypersurface.
For example, in S 3 , the critical domains are: geodesic balls, tubes around an equator, domains bounded by a Clifford torus and tubes around a minimal Clifford torus. The classification exists in higher dimension, but it gets more complicated due to the large variety of isoparametric foliations in higher dimensions.
• As a final remark, we ask the following question: is it possible to weaken the assumption "Ω is critical for all T k " in the statement of the main theorem ? For example, if we assume criticality for only one k can we get the required rigidity? Well, the answer is negative, at least if k = 1: Remark 9. There exist (analytic) Riemannian domains which are critical for torsional rigidity T 1 but are not isoparametric tubes.
A first such example was constructed in [32], and consists of any minimal free boundary immersion in the 3-dimensional unit Euclidean ball B 3 having more than two boundary components: the normal derivative of the torsion function E 1 is constant on the boundary, but Ω cannot be a smooth tube. Other "exotic" examples exist even in the round sphere; in [12] one can find domains which are critical for torsional rigidity but have boundary with non-constant mean curvature. These domains are in fact suitable perturbations of a tubular neighborhood of the equator.

First variation of the heat content
The main step in the proof of Theorem 1 is the formula for the first variation of the heat content, which is the following.
Theorem 10. For a fixed value of time t > 0, let F t (Ω) = H Ω (t) be the heat content of Ω at time t. Then its first variation in the direction V is given by: where N is the inner unit normal and u(t, x) is the temperature function defined in (1).
Note that, as t → 0, one has, for all y ∈ ∂Ω (see [30]): which guarantees that the integral on the right hand side is convergent. We remark that Ozawa obtained in [24] the first variation of the Dirichlet heat kernel of Euclidean domains. Our methods are different and are based on the fact that the heat kernel depends smoothly on the deformation parameter (see [27]).

Scheme of proof
The rest of the paper is organized as follows.
In Section 2 we prove the first variation formula, Theorem 10. In Section 3 we prove the equivalence between a) and c) of the main theorem: Theorem 11. Ω is critical for the functional F t (Ω) given by the heat content at time t, for all t > 0, if and only if Ω has the constant flow property.

In Section 3 we show the equivalence between b) and c):
Theorem 12. The domain Ω is critical for the k-th exit time moment, for all k ≥ 1 and for all volume preserving deformations, if and only if Ω has the constant flow property.
Having done that, we finish the proof of the main theorem by recalling that, when (M, g) is analytic, the equivalence between c) and d) has been proved in Theorem 7 of [32]. With this in mind, the proof of Theorem 6 is complete.

Proof of the first variation formula
Let Ω ǫ = f ǫ (Ω) be a smooth deformation of Ω associated to the vector field V , as in (10). We adopt the following point of view: deforming the domain in a manifold with a fixed metric is equivalent to keeping the domain fixed and deforming the metric. In more precise terms, we identify (Ω ǫ , g) with (Ω, f ⋆ ǫ g) and denote g ǫ . = f ⋆ ǫ g. We let ∆ ǫ be the Laplacian associated to the metric g ǫ of Ω, and dv ǫ the corresponding Riemannian measure. We denote: It is a standard fact that: where dv = dv g is the Riemannian measure associated to the original metric g and δ is the adjoint of the gradient operator. We point out that the operator ∆ ′ has been computed by Berger in [4]; the explicit expressions of g ′ and ∆ ′ will be given in Lemma 20 of Appendix 1.
Let us denote u ǫ (t, x) the temperature function in the deformed metric g ǫ . It is the unique where k ǫ (t, x, y) is the Dirichlet heat kernel of (Ω, g ǫ ) (resp. the Dirichlet heat kernel of (Ω, g)).
From Proposition 6.1 in [27], we know that the map ǫ → k ǫ is differentiable 1 . Therefore, from the expression (13) also the map ǫ → u ǫ (t, x) is differentiable and we denote: Accordingly, we denote (we stress that the prime indicates differentation with respect to ǫ along V , and not with respect to time t, which in this discussion is fixed). The proof goes as follows: we first express the function u ′ in terms of the Dirichlet heat kernel of g and ∆ ′ u. Then, using an explicit expression of the operator ∆ ′ , and integration by parts, we obtain the final statement. To avoid discussing, at every step, the convergence of the integrals involved we use a suitable approximation v(t, x) of the function u(t, x) by a small parameter δ and then pass to the limit as δ → 0 to obtain the statement.

Approximation of u(t, x)
We fix δ > 0 and small and we introduce the function v : (0, ∞) × Ω → R defined by: v(t, x) = u(t + δ, x) (for simplicity of notation, we omit to write explicitly the dependance of v on δ). Then, v satisfies the heat equation on (Ω, g) with initial data v(0, x) = u(δ, x) and Dirichlet boundary conditions. Since the initial data vanishes on the boundary, v(t, x) extends to a smooth function on [0, ∞) ×Ω and we thus can avoid dealing with the distributional behavior of u(t, ·) near the boundary, for small times, in the sense that any derivative of v is uniformly bounded on [0, ∞) ×Ω. Since: we see that u(t, x) is decreasing in t at any point x: this implies that v(t, x) ≤ u(t, x) ≤ 1 on (0, ∞) × Ω, and so, as both functions vanish on the boundary: for all t > 0 and y ∈ ∂Ω. Finally, it is clear that v(t, ·) → u(t, ·) together with all of its derivatives, as δ → 0. Let (Ω, g ǫ ) be a smooth deformation of (Ω, g) and let v ǫ (t, x) be the solution of the Dirichlet heat equation in (Ω, g ǫ ) with initial data u(δ, x) (not depending on ǫ). Hence v ǫ satisfies: for all x ∈ Ω v ǫ (t, y) = 0 for all y ∈ ∂Ω, t > 0 We introduce the following notation for a fixed δ > 0: Lemma 13. For any fixed t > 0 one has: u ′ (t, x) = lim δ→0 v ′ (t, x) and therefore: Proof. First observe that v ǫ (t, x) = Ω k ǫ (t, x, y)u(δ, y) dv ǫ (y) hence, differentiating both sides with respect to ǫ and setting ǫ = 0: v ′ (t, x) = Ω k ′ (t, x, y)u(δ, y) dv(y) + Ω k(t, x, y)u(δ, y)dv ′ (y).
We now let δ → 0 on both sides; as u(δ, x) is uniformly bounded by 1 for all δ and converges to 1 as δ → 0 we get, for all fixed t > 0: which is the first assertion. The second assertion follows from the first by a similar argument, passing to the limit in:

Duhamel principle
We now fix δ and work on a convenient expression of H ′ Ω,δ (t). Derive both sides of the equation ∆ ǫ v ǫ + ∂v ǫ ∂t = 0 and set ǫ = 0 to get: Therefore, since the initial data of v ǫ does not depend on ǫ, the function v ′ (t, x) satisfies the following heat equation in the original metric g: Integrating on Ω with respect to x we see (use Fubini and observe lim δ→0 v(t, x) = u(t, x)): Using (16), knowing that lim δ→0 v(t, x) = u(t, x) and that dv ′ = −δV dv we arrive at the following expression.
Lemma 14. In the above notation: The final step is to deal with the integral involving the operator ∆ ′ . This is done in the following lemma.

Lemma 15.
If the functions f, h vanish on the boundary one has: Proof. See Appendix 1.

Proof of Theorem 10
We use Lemma 14 and apply Lemma 15 taking f = v(τ, ·) and h = v(t − τ , ·). We obtain (integration on Ω is of course taken with respect to the variable y).

Equivalence between parts a) and c) of Theorem 1
In this section we prove the following fact.
Theorem 16. Ω is critical for the functional F t (Ω) given by the heat content at time t, for all t > 0, if and only if Ω has the constant flow property.
Recall the formula for the first variation (Theorem 10).
First assume that Ω has the constant flow property. Then, by assumption, ∂u ∂N (t, y) = c(t) for all y ∈ ∂Ω and for a smooth function c = c(t). Then: because the deformation preserves the volume (see (11)). Hence Ω is critical.
Conversely, assume that DF t (Ω, V ) = 0 for all deformation vector fields V and for all t > 0, so that, interchanging the order of integration, we have by assumption: where we set: As V, N has zero mean on ∂Ω, and is otherwise arbitrary, we see that (24) forces B(t, y) to be constant in y, that is: for a smooth (positive) function γ : (0, ∞) → R. Let the hat denote Laplace tranform with respect to time t ∈ (0, ∞). Applied to u(t, x) we obtain the function We can differentiate in the normal direction under the integral sign in (27), so that, for all y ∈ ∂Ω and s > 0: where on the right we have the Laplace trasform of the heat flow. We see, from (25) for all t > 0 and y 1 , y 2 ∈ ∂Ω. Then, the function ∂u ∂N (t, ·) is constant on the boundary, and Ω has the constant flow property.

Equivalence between parts b) and c) of Theorem 1
In this section we prove: Theorem 17. The domain Ω is critical for the k-th exit time moment, for all k ≥ 1 and for all volume preserving deformations, if and only if Ω has the constant flow property.
The proof uses the calculation of the first variation of the k-th exit time moment T k (Ω) first done [19]. For convenience of the reader we derive the same formula, using our approach, in Appendix 2 in this paper. Lemma 18. The first variation of the k-th exit time moment T k , at Ω, in the direction V , is given by: where and E k (x) is the k-th exit time function defined in (6). In particular, Ω is critical for the k-th exit time moment, for all k ≥ 1 and for all volume preserving deformations, if and only if ∂E k ∂N is constant on the boundary, for all k.
Proof. For the proof of (44) see Proposition 2.1 of [19], or Appendix 2. The last assertion follows from an obvious induction on k, observing that

Proof of Theorem 17
Let For φ ∈ H 0 (Ω), we consider the function: which is easily seen to represent the heat content of Ω with initial temperature distribution given by φ(x) and Dirichlet boundary conditions. It is proved in Theorem 9 of [31] that Ω has the constant flow property if and only if H φ ≡ 0 for all φ ∈ H 0 (Ω). An easy Fourier series argument (see Lemma 8 of [31]) shows that the following upper bound holds for all t ≥ 0: where λ 1 > 0 is the first Dirichlet eigenvalue of Ω. Therefore, the Laplace transform Let us compute the Taylor series ofĤ φ (s) at s = 0. One haŝ by (8); more generally, the k-th derivative ofĤ φ at s = 0 is given by: We will use (30) to get the conclusion. First assume that Ω has the constant flow property. Then H φ ≡ 0 for all φ ∈ H 0 (Ω) (in particular Ω φ = 0) and, from (30): for all n ≥ 0. Integrating by parts we obtain: for all n ≥ 0. As any smooth function φ on ∂Ω admits an harmonic extension to the interior, we see that ∂E n+1 ∂N must be orthogonal to the subspace of zero mean functions on ∂Ω, hence it must be constant on ∂Ω, for all n ≥ 0. By Lemma 23 we conclude that Ω is critical for all T k 's.
Conversely, assume that Ω is critical for all T k 's. Then all normal derivatives ∂E k+1 ∂N are constant on the boundary of Ω, so that Ω φE n dv = 0 for all n (from (31)). From (30), we see that the Taylor series ofĤ φ (s) at s = 0 is zero; asĤ φ is analytic on (−λ 1 , ∞) we must haveĤ φ ≡ 0 hence also H φ ≡ 0. Then Ω has the constant flow property.

Appendix 1: proof of Lemma 15
In this section we consider a deformation f ǫ : Ω → Ω ǫ of the domain Ω along a vector field V and identify (Ω ǫ , g) with (Ω, g ǫ ) where g ǫ = f ⋆ ǫ g. We let ∆ ǫ be the Laplacian associated to the metric g ǫ of Ω, and dv ǫ the corresponding Riemannian measure. Recall the notation: We want to prove the following expression.
Lemma 19. Assume that the functions φ, ψ : Ω → R vanish on the boundary. Then: We first remark that it is enough to prove the lemma when φ = ψ; the general formula will follow by polarization. That is, it is enough to show that, for all functions φ vanishing on the boundary: In what follows, we denote by δ the adjoint of the gradient operator, so that δV = −divV = − n j=1 ∇ ej V , e j where (e 1 , . . . , e n ) is any orthonormal frame. We also recall that if h is a symmetric two-tensor then δh is by definition the one-form: Lemma 20. Let g ′ , ∆ ′ and dv ′ be as above. Let V be the deformation vector field. Then g ′ = L V g and one has: In fact, formulae (33), (34), (35) are standard, while (36) is due to Berger [4]. Lemma (19) will drop from the following two facts.
Lemma 21. One has for every φ vanishing on the boundary: The next one is a Rellich-type identity.
Lemma 22. Let X be any vector field on Ω and assume that φ vanishes on ∂Ω. Then: Lemma 21 and Lemma 22 will be proved below.

Proof of (32)
It follows by taking X = V in Lemma 22 and by inserting the resulting identity in Lemma 21.

Appendix 2: proof of Lemma 23
We want to prove the following fact.
Lemma 23. The first variation of the k-th exit time moment T k , at Ω, in the direction V , is given by: where Φ k = k j=1 c kj ∂E j ∂N ∂E k+1−j ∂N , c kj = k! (k + 1 − j)!j! and E k (x) is the k-th exit time function defined in (6). In particular, Ω is critical for the k-th exit time moment, for all k ≥ 1 and for all volume preserving deformations, if and only if ∂E k ∂N is constant on the boundary, for all k.
It is enough to show (44), as the last assertion will follow from an obvious inductive argument.
It will be simpler to work with the functions u k = 1 k! E k which satisfy u 0 = 1 and ∆u k = u k−1 on Ω u k = 0 on ∂Ω for k ≥ 1. Proving (44) is in turn equivalent to proving the following fact. Let Q k (Ω) be the functional Then: Then (47) follows from (50) and (48).
Step 2. Taking into account (47), to prove the final statement it is enough to show that Once this identity has been established, formula (46) follows because Ω u k dv ′ g = − Ω u k δV dv From Lemma 19 we see that: The sum of the two inner integrals, thanks to the defining relations, and after rearranging terms, equals: Similarly one sees that k j=2 D(j) = 0, the only term surviving in the second sum is when j = 1, and then: