Here the geometry of a smooth mapping between Riemannian manifolds is considered via local coordinates and, equivalently, via linear algebra. Such a mapping defines both a first fundamental form and second fundamental form. The Laplacian (also called tension field) is defined via the second fundamental form, and its vanishing is the condition for the map to be harmonic. The definitions extend without modification to the setting of pseudo-Riemannian manifolds.

Local coordinates edit

Let U be an open subset of ℝ^m and let V be an open subset of ℝⁿ. For each i and j between 1 and n, let g_ij be a smooth real-valued function on U, such that for each p in U, one has that the m × m matrix [g_ij (p)] is symmetric and positive-definite. For each α and β between 1 and m, let h_αβ be a smooth real-valued function on V, such that for each q in V, one has that the n × n matrix [h_αβ (q)] is symmetric and positive-definite. Denote the inverse matrices by [g^ij (p)] and [h^αβ (q)].

For each i, j, k between 1 and n and each α, β, γ between 1 and m define the Christoffel symbols Γ(g)^k_ij : U → ℝ and Γ(h)^γ_αβ : V → ℝ by^[4]

${\displaystyle {\begin{aligned}\Gamma (g)_{ij}^{k}&={\frac {1}{2}}\sum _{\ell =1}^{m}g^{k\ell }{\Big (}{\frac {\partial g_{j\ell }}{\partial x^{i}}}+{\frac {\partial g_{i\ell }}{\partial x^{j}}}-{\frac {\partial g_{ij}}{\partial x^{\ell }}}{\Big )}\\\Gamma (h)_{\alpha \beta }^{\gamma }&={\frac {1}{2}}\sum _{\delta =1}^{n}h^{\gamma \delta }{\Big (}{\frac {\partial h_{\beta \delta }}{\partial y^{\alpha }}}+{\frac {\partial h_{\alpha \delta }}{\partial y^{\beta }}}-{\frac {\partial h_{\alpha \beta }}{\partial y^{\delta }}}{\Big )}\end{aligned}}}$ 

Given a smooth map f from U to V, its second fundamental form defines for each i and j between 1 and m and for each α between 1 and n the real-valued function ∇(df)^α_ij on U by^[5]

${\displaystyle \nabla (df)_{ij}^{\alpha }={\frac {\partial ^{2}f^{\alpha }}{\partial x^{i}\partial x^{j}}}-\sum _{k=1}^{m}\Gamma (g)_{ij}^{k}{\frac {\partial f^{\alpha }}{\partial x^{k}}}+\sum _{\beta =1}^{n}\sum _{\gamma =1}^{n}{\frac {\partial f^{\beta }}{\partial x^{i}}}{\frac {\partial f^{\gamma }}{\partial x^{j}}}\Gamma (h)_{\beta \gamma }^{\alpha }\circ f.}$ 

Its laplacian defines for each α between 1 and n the real-valued function (∆f)^α on U by^[6]

${\displaystyle (\Delta f)^{\alpha }=\sum _{i=1}^{m}\sum _{j=1}^{m}g^{ij}\nabla (df)_{ij}^{\alpha }.}$ 

Bundle formalism edit

Let (M, g) and (N, h) be Riemannian manifolds. Given a smooth map f from M to N, one can consider its differential df as a section of the vector bundle T ^*M ⊗ f ^*TN over M; this is to say that for each p in M, one has a linear map df_p between tangent spaces T_pM → T_f(p)N.^[7] The vector bundle T ^*M ⊗ f ^*TN has a connection induced from the Levi-Civita connections on M and N.^[8] So one may take the covariant derivative ∇(df), which is a section of the vector bundle T ^*M ⊗ T ^*M ⊗ f ^*TN over M; this is to say that for each p in M, one has a bilinear map (∇(df))_p of tangent spaces T_pM × T_pM → T_f(p)N.^[9] This section is known as the hessian of f.

Using g, one may trace the hessian of f to arrive at the laplacian of f, which is a section of the bundle f ^*TN over M; this says that the laplacian of f assigns to each p in M an element of the tangent space T_f(p)N.^[10] By the definition of the trace operator, the laplacian may be written as

${\displaystyle (\Delta f)_{p}=\sum _{i=1}^{m}{\big (}\nabla (df){\big )}_{p}(e_{i},e_{i})}$ 

where e₁, ..., e_m is any g_p-orthonormal basis of T_pM.

From the perspective of local coordinates, as given above, the energy density of a mapping f is the real-valued function on U given by^[11]

${\displaystyle {\frac {1}{2}}\sum _{i=1}^{m}\sum _{j=1}^{m}\sum _{\alpha =1}^{n}\sum _{\beta =1}^{n}g^{ij}{\frac {\partial f^{\alpha }}{\partial x^{i}}}{\frac {\partial f^{\beta }}{\partial x^{j}}}(h_{\alpha \beta }\circ f).}$ 

Alternatively, in the bundle formalism, the Riemannian metrics on M and N induce a bundle metric on T ^*M ⊗ f ^*TN, and so one may define the energy density as the smooth function 1/2 | df |² on M.^[12] It is also possible to consider the energy density as being given by (half of) the g-trace of the first fundamental form.^[13] Regardless of the perspective taken, the energy density e(f) is a function on M which is smooth and nonnegative. If M is oriented and M is compact, the Dirichlet energy of f is defined as

${\displaystyle E(f)=\int _{M}e(f)\,d\mu _{g}}$ 

where dμ_g is the volume form on M induced by g.^[14] Since any nonnegative measurable function has a well-defined Lebesgue integral, it is not necessary to place the restriction that M is compact; however, then the Dirichlet energy could be infinite.

The variation formulas for the Dirichlet energy compute the derivatives of the Dirichlet energy E(f) as the mapping f is deformed. To this end, consider a one-parameter family of maps f_s : M → N with f₀ = f for which there exists a precompact open set K of M such that f_s|_{M − K} = f|_{M − K} for all s; one supposes that the parametrized family is smooth in the sense that the associated map (−ε, ε) × M → N given by (s, p) ↦ f_s(p) is smooth.

• The first variation formula says that^[15]

${\displaystyle \int _{M}{\frac {\partial }{\partial s}}{\Big |}_{s=0}e(f_{s})\,d\mu _{g}=-\int _{M}h\left({\frac {\partial }{\partial s}}{\Big |}_{s=0}f_{s},\Delta f\right)\,d\mu _{g}}$ 

There is also a version for manifolds with boundary.^[16]

• There is also a second variation formula.^[17]

Due to the first variation formula, the Laplacian of f can be thought of as the gradient of the Dirichlet energy; correspondingly, a harmonic map is a critical point of the Dirichlet energy.^[18] This can be done formally in the language of global analysis and Banach manifolds.

Let (M, g) and (N, h) be smooth Riemannian manifolds. The notation g_stan is used to refer to the standard Riemannian metric on Euclidean space.

• Every totally geodesic map (M, g) → (N, h) is harmonic; this follows directly from the above definitions. As special cases:
  • For any q in N, the constant map (M, g) → (N, h) valued at q is harmonic.
  • The identity map (M, g) → (M, g) is harmonic.
• If f : M → N is an immersion, then f : (M, f ^*h) → (N, h) is harmonic if and only if f is minimal relative to h. As a special case:
  • If f : ℝ → (N, h) is a constant-speed immersion, then f : (ℝ, g_stan) → (N, h) is harmonic if and only if f solves the geodesic differential equation.

Recall that if M is one-dimensional, then minimality of f is equivalent to f being geodesic, although this does not imply that it is a constant-speed parametrization, and hence does not imply that f solves the geodesic differential equation.

• A smooth map f : (M, g) → (ℝⁿ, g_stan) is harmonic if and only if each of its n component functions are harmonic as maps (M, g) → (ℝ, g_stan). This coincides with the notion of harmonicity provided by the Laplace-Beltrami operator.
• Every holomorphic map between Kähler manifolds is harmonic.
• Every harmonic morphism between Riemannian manifolds is harmonic.

Well-posedness edit

Let (M, g) and (N, h) be smooth Riemannian manifolds. A harmonic map heat flow on an interval (a, b) assigns to each t in (a, b) a twice-differentiable map f_t : M → N in such a way that, for each p in M, the map (a, b) → N given by t ↦ f_t (p) is differentiable, and its derivative at a given value of t is, as a vector in T_{ft (p)}N, equal to (∆ f_t )_p. This is usually abbreviated as:

${\displaystyle {\frac {\partial f}{\partial t}}=\Delta f.}$ 

Eells and Sampson introduced the harmonic map heat flow and proved the following fundamental properties:

• Regularity. Any harmonic map heat flow is smooth as a map (a, b) × M → N given by (t, p) ↦ f_t (p).

Now suppose that M is a closed manifold and (N, h) is geodesically complete.

• Existence. Given a continuously differentiable map f from M to N, there exists a positive number T and a harmonic map heat flow f_t on the interval (0, T) such that f_t converges to f in the C¹ topology as t decreases to 0.^[19]
• Uniqueness. If { f_t : 0 < t < T } and { f _t : 0 < t < T } are two harmonic map heat flows as in the existence theorem, then f_t = f _t whenever 0 < t < min(T, T).

As a consequence of the uniqueness theorem, there exists a maximal harmonic map heat flow with initial data f, meaning that one has a harmonic map heat flow { f_t : 0 < t < T } as in the statement of the existence theorem, and it is uniquely defined under the extra criterion that T takes on its maximal possible value, which could be infinite.

Eells and Sampson's theorem edit

The primary result of Eells and Sampson's 1964 paper is the following:^[1]

Let (M, g) and (N, h) be smooth and closed Riemannian manifolds, and suppose that the sectional curvature of (N, h) is nonpositive. Then for any continuously differentiable map f from M to N, the maximal harmonic map heat flow { f_t : 0 < t < T } with initial data f has T = ∞, and as t increases to ∞, the maps f_t subsequentially converge in the C^∞ topology to a harmonic map.

In particular, this shows that, under the assumptions on (M, g) and (N, h), every continuous map is homotopic to a harmonic map.^[1] The very existence of a harmonic map in each homotopy class, which is implicitly being asserted, is part of the result. Shortly after Eells and Sampson's work, Philip Hartman extended their methods to study uniqueness of harmonic maps within homotopy classes, additionally showing that the convergence in the Eells−Sampson theorem is strong, without the need to select a subsequence.^[20] Eells and Sampson's result was adapted by Richard Hamilton to the setting of the Dirichlet boundary value problem, when M is instead compact with nonempty boundary.^[21]

Singularities and weak solutions edit

For many years after Eells and Sampson's work, it was unclear to what extent the sectional curvature assumption on (N, h) was necessary. Following the work of Kung-Ching Chang, Wei-Yue Ding, and Rugang Ye in 1992, it is widely accepted that the maximal time of existence of a harmonic map heat flow cannot "usually" be expected to be infinite.^[22] Their results strongly suggest that there are harmonic map heat flows with "finite-time blowup" even when both (M, g) and (N, h) are taken to be the two-dimensional sphere with its standard metric. Since elliptic and parabolic partial differential equations are particularly smooth when the domain is two dimensions, the Chang−Ding−Ye result is considered to be indicative of the general character of the flow.

Modeled upon the fundamental works of Sacks and Uhlenbeck, Michael Struwe considered the case where no geometric assumption on (N, h) is made. In the case that M is two-dimensional, he established the unconditional existence and uniqueness for weak solutions of the harmonic map heat flow.^[23] Moreover, he found that his weak solutions are smooth away from finitely many spacetime points at which the energy density concentrates. On microscopic levels, the flow near these points is modeled by a bubble, i.e. a smooth harmonic map from the round 2-sphere into the target. Weiyue Ding and Gang Tian were able to prove the energy quantization at singular times, meaning that the Dirichlet energy of Struwe's weak solution, at a singular time, drops by exactly the sum of the total Dirichlet energies of the bubbles corresponding to singularities at that time.^[24]

Struwe was later able to adapt his methods to higher dimensions, in the case that the domain manifold is Euclidean space;^[25] he and Yun Mei Chen also considered higher-dimensional closed manifolds.^[26] Their results achieved less than in low dimensions, only being able to prove existence of weak solutions which are smooth on open dense subsets.

The main computational point in the proof of Eells and Sampson's theorem is an adaptation of the Bochner formula to the setting of a harmonic map heat flow { f_t : 0 < t < T }. This formula says^[27]

${\displaystyle {\Big (}{\frac {\partial }{\partial t}}-\Delta ^{g}{\Big )}e(f)=-{\big |}\nabla (df){\big |}^{2}-{\big \langle }\operatorname {Ric} ^{g},f^{\ast }h{\big \rangle }_{g}+\operatorname {scal} ^{g}{\big (}f^{\ast }\operatorname {Rm} ^{h}{\big )}.}$ 

This is also of interest in analyzing harmonic maps. Suppose f : M → N is harmonic; any harmonic map can be viewed as a constant-in-t solution of the harmonic map heat flow, and so one gets from the above formula that^[28]

${\displaystyle \Delta ^{g}e(f)={\big |}\nabla (df){\big |}^{2}+{\big \langle }\operatorname {Ric} ^{g},f^{\ast }h{\big \rangle }_{g}-\operatorname {scal} ^{g}{\big (}f^{\ast }\operatorname {Rm} ^{h}{\big )}.}$ 

If the Ricci curvature of g is positive and the sectional curvature of h is nonpositive, then this implies that ∆e(f) is nonnegative. If M is closed, then multiplication by e(f) and a single integration by parts shows that e(f) must be constant, and hence zero; hence f must itself be constant.^[29] Richard Schoen and Shing-Tung Yau noted that this reasoning can be extended to noncompact M by making use of Yau's theorem asserting that nonnegative subharmonic functions which are L²-bounded must be constant.^[30] In summary, according to these results, one has:

Let (M, g) and (N, h) be smooth and complete Riemannian manifolds, and let f be a harmonic map from M to N. Suppose that the Ricci curvature of g is positive and the sectional curvature of h is nonpositive.

• If M and N are both closed then f must be constant.
• If N is closed and f has finite Dirichlet energy, then it must be constant.

In combination with the Eells−Sampson theorem, this shows (for instance) that if (M, g) is a closed Riemannian manifold with positive Ricci curvature and (N, h) is a closed Riemannian manifold with nonpositive sectional curvature, then every continuous map from M to N is homotopic to a constant.

The general idea of deforming a general map to a harmonic map, and then showing that any such harmonic map must automatically be of a highly restricted class, has found many applications. For instance, Yum-Tong Siu found an important complex-analytic version of the Bochner formula, asserting that a harmonic map between Kähler manifolds must be holomorphic, provided that the target manifold has appropriately negative curvature.^[31] As an application, by making use of the Eells−Sampson existence theorem for harmonic maps, he was able to show that if (M, g) and (N, h) are smooth and closed Kähler manifolds, and if the curvature of (N, h) is appropriately negative, then M and N must be biholomorphic or anti-biholomorphic if they are homotopic to each other; the biholomorphism (or anti-biholomorphism) is precisely the harmonic map produced as the limit of the harmonic map heat flow with initial data given by the homotopy. By an alternative formulation of the same approach, Siu was able to prove a variant of the still-unsolved Hodge conjecture, albeit in the restricted context of negative curvature.

Kevin Corlette found a significant extension of Siu's Bochner formula, and used it to prove new rigidity theorems for lattices in certain Lie groups.^[32] Following this, Mikhael Gromov and Richard Schoen extended much of the theory of harmonic maps to allow (N, h) to be replaced by a metric space.^[33] By an extension of the Eells−Sampson theorem together with an extension of the Siu–Corlette Bochner formula, they were able to prove new rigidity theorems for lattices.

• Existence results on harmonic maps between manifolds has consequences for their curvature.
• Once existence is known, how can a harmonic map be constructed explicitly? (One fruitful method uses twistor theory.)
• In theoretical physics, a quantum field theory whose action is given by the Dirichlet energy is known as a sigma model. In such a theory, harmonic maps correspond to instantons.
• One of the original ideas in grid generation methods for computational fluid dynamics and computational physics was to use either conformal or harmonic mapping to generate regular grids.

The energy integral can be formulated in a weaker setting for functions u : M → N between two metric spaces. The energy integrand is instead a function of the form

${\displaystyle e_{\epsilon }(u)(x)={\frac {\int _{M}d^{2}(u(x),u(y))\,d\mu _{x}^{\epsilon }(y)}{\int _{M}d^{2}(x,y)\,d\mu _{x}^{\epsilon }(y)}}}$ 

in which μ^ε
_x is a family of measures attached to each point of M.^[34]

• Geometric flow

Footnotes

Articles

Books and surveys

• MathWorld: Harmonic map
• Harmonic Maps Bibliography
• The Bibliography of Harmonic Morphisms