曲線の変分論

同義語:calculus of variations for curves

概要

曲線の変分論(calculus of variations for curves)とは、リーマン多様体上の曲線を滑らかな変分の族の中で動かし、エネルギー $E$ や長さ $L$ がどう変化するかを調べる理論である。第1変分公式により、両端点固定でエネルギーが停留する曲線はちょうど測地線であり、正則曲線の長さが停留するのは弧長で再パラメータ化すると測地線になるときであって、最短曲線は測地線になる。測地線に沿うエネルギーの第2変分(指数形式)には曲率が現れ、最短な測地線では非負になるので、単位球面上の長さが $\pi$ を超える大円弧は測地線でも最短でないことが分かる。区分的に滑らかな曲線では折点項が残るため、停留曲線に角は現れない。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 多様体, リーマン幾何学, 内積, Riemann積分

曲線を動かす

曲線の変分論(calculus of variations for curves)は、一本の曲線を滑らかな族の中で動かし、長さやエネルギーが1次にどう変化するかを調べる理論である。有限次元の関数 $f(x)$ では微分 $df_x$ が停留点を検出する。曲線の変分論では、変数 $x$ の代わりに曲線 $\gamma$ 自身を動かし、汎関数の微分に相当する第1変分を考える。
本記事では、リーマン多様体上のエネルギーと長さの第1変分を導き、その停留曲線が測地線であることを示す。後半では、測地線に沿うエネルギーの第2変分公式を導き、第1変分だけでは分からない「停留曲線が最短かどうか」を曲率で判定する入口を示す。定式化は Mil63 の第 III 部(§12–13)と Nic26 の §5.2 に沿う。

変分と変分ベクトル場

曲線の変分

$(M,g)$ をリーマン多様体、$a< b$ とし、$\gamma:[a,b]\to M$ を滑らかな曲線とする。$\varepsilon>0$ に対する滑らかな写像
$$ F:(-\varepsilon,\varepsilon)\times[a,b]\to M $$
が $F(0,t)=\gamma(t)$ を満たすとき、$F$ を $\gamma$ の変分という。各 $s$ に対する曲線を $\gamma_s(t):=F(s,t)$ と書く。
変分が両端点を固定するとは、すべての $s$ に対して
$$ F(s,a)=\gamma(a),\qquad F(s,b)=\gamma(b) $$
が成り立つことである。
$F$ から得られる $\gamma$ に沿うベクトル場
$$ V(t):=\left.\frac{\partial F}{\partial s}\right|_{(0,t)} $$
を変分ベクトル場という。両端点を固定する変分では $V(a)=V(b)=0$ である。

写像 $F$ は、曲線 $\gamma_s$ を $s$ の方向に並べた小さな曲面とみなせる。その曲面上に
$$ T:=\frac{\partial F}{\partial t},\qquad S:=\frac{\partial F}{\partial s} $$
を定める。$\nabla$ を リーマン接続 とすると、接続が捩れをもたず $[\partial_s,\partial_t]=0$ であることから
$$ \nabla_S T=\nabla_T S $$
が成り立つ。この微分の交換が第1変分公式の要点である。

沿うベクトル場の変分としての実現

$\gamma:[a,b]\to M$ に沿う滑らかなベクトル場 $V$ は、ある変分の変分ベクトル場として実現できる。$V(a)=V(b)=0$ ならば、その変分を両端点固定に選べる。

$\exp$ をリーマン指数写像とする。零断面の近傍で $\exp$ は定義され、$[a,b]$ はコンパクトであるから、十分小さい一定の $\varepsilon>0$ に対して
$$ F(s,t):=\exp_{\gamma(t)}(sV(t)) $$
が $|s|<\varepsilon$ とすべての $t\in[a,b]$ で定義される。$F(0,t)=\gamma(t)$ なので $F$ は変分である。また指数写像の定義から
$$ \left.\frac{\partial}{\partial s}\right|_{s=0} \exp_{\gamma(t)}(sV(t))=V(t) $$
である。$V(a)=V(b)=0$ なら
$F(s,a)=\exp_{\gamma(a)}(0)=\gamma(a)$ および $F(s,b)=\gamma(b)$ となり、両端点も固定される。$\blacksquare$

エネルギーと長さ

エネルギー汎関数と長さ汎関数

滑らかな曲線 $\gamma:[a,b]\to M$ のエネルギーと長さをそれぞれ
$$ E(\gamma):=\frac12\int_a^b\|\dot\gamma(t)\|^2\,dt, \qquad L(\gamma):=\int_a^b\|\dot\gamma(t)\|\,dt $$
と定める。長さの変分を微分するときは $\dot\gamma(t)\ne0$ をすべての $t$ で仮定し、このような曲線を正則曲線という。

長さは向きを保つ再パラメータ化で変わらないが、エネルギーは速度の選び方に依存する。この違いは停留曲線の特徴に現れる。エネルギーはアフィン・パラメータ付きの測地線を、長さは再パラメータ化を除いた測地線の軌跡を検出する。

エネルギーと長さの不等式

任意の滑らかな曲線 $\gamma:[a,b]\to M$ に対して
$$ E(\gamma)\ge \frac{L(\gamma)^2}{2(b-a)} $$
が成り立つ。等号は速度 $\|\dot\gamma(t)\|$ が定数であるとき、またそのときに限り成り立つ。

$[a,b]$ 上の関数 $1$ と $\|\dot\gamma\|$ にCauchy--Schwarzの不等式を適用すると
$$ L(\gamma)^2 =\left(\int_a^b 1\cdot\|\dot\gamma(t)\|\,dt\right)^2 \le (b-a)\int_a^b\|\dot\gamma(t)\|^2\,dt =2(b-a)E(\gamma) $$
を得る。これを変形すれば所望の不等式になる。Cauchy--Schwarzの等号条件は $\|\dot\gamma(t)\|$ が関数 $1$ の定数倍であること、すなわち速度が定数であることである。$\blacksquare$

第1変分公式

エネルギーの第1変分

$F$ を $\gamma$ の変分、$V$ をその変分ベクトル場とする。$D_tV:=\nabla_{\dot\gamma}V$ と書くと
$$ \left.\frac{d}{ds}\right|_{s=0}E(\gamma_s) =\bigl[g(V,\dot\gamma)\bigr]_a^b -\int_a^b g(V,D_t\dot\gamma)\,dt $$
が成り立つ。特に両端点固定の変分では
$$ \left.\frac{d}{ds}\right|_{s=0}E(\gamma_s) =-\int_a^b g(V,D_t\dot\gamma)\,dt $$
である。

リーマン接続 の計量両立性と $\nabla_S T=\nabla_T S$ を用いると
$$ \begin{aligned} \left.\frac{d}{ds}\right|_{s=0}E(\gamma_s) &=\frac12\int_a^b \left.\frac{\partial}{\partial s}\right|_{s=0}g(T,T)\,dt\\ &=\int_a^b g(\nabla_S T,T)|_{s=0}\,dt\\ &=\int_a^b g(D_tV,\dot\gamma)\,dt \end{aligned} $$
となる。一方、再び計量両立性から
$$ \frac{d}{dt}g(V,\dot\gamma) =g(D_tV,\dot\gamma)+g(V,D_t\dot\gamma) $$
である。これを積分して項を移せば最初の公式を得る。両端点固定なら $V(a)=V(b)=0$ なので境界項が消える。$\blacksquare$

エネルギーの停留曲線

両端点を固定したすべての変分に対してエネルギーの第1変分が $0$ であることと、
$$ D_t\dot\gamma=0 $$
が成り立つことは同値である。すなわちエネルギーの停留曲線は、アフィン・パラメータをもつ測地線に他ならない。

$D_t\dot\gamma=0$ なら、thm-first-variation-energy の両端点固定公式からすべての第1変分が $0$ になる。
逆にすべての両端点固定変分に対して第1変分が $0$ と仮定する。prop-variation-field-realization により、端点で消える任意の滑らかな $V$ を試験できるから
$$ \int_a^b g(V,D_t\dot\gamma)\,dt=0 $$
がすべてのそのような $V$ に対して成り立つ。もし $D_t\dot\gamma(t_0)\ne0$ となる内点 $t_0$ があれば、$t_0$ の近傍ではそのノルムは正である。その近傍に台をもつ非負の滑らかな関数 $\chi$ で $\chi(t_0)>0$ を満たすものを取り、$V=\chi D_t\dot\gamma$ と置けば
$$ \int_a^b g(V,D_t\dot\gamma)\,dt =\int_a^b\chi\|D_t\dot\gamma\|^2\,dt>0 $$
となって矛盾する。よって $D_t\dot\gamma=0$ である。$\blacksquare$

長さの第1変分

$\gamma$ を正則曲線、$F$ を $\gamma$ の変分、$V$ を変分ベクトル場とし、
$$ U:=\frac{\dot\gamma}{\|\dot\gamma\|} $$
を単位接ベクトル場とする。このとき
$$ \left.\frac{d}{ds}\right|_{s=0}L(\gamma_s) =\bigl[g(V,U)\bigr]_a^b-\int_a^b g(V,D_tU)\,dt $$
が成り立つ。従って両端点固定の変分で長さが停留することと $D_tU=0$ は同値である。これは、$\gamma$ が弧長で再パラメータ化すると測地線になることを意味する。

$T=\partial_tF$ とする。正則性により $s=0$ で $\|T\|>0$ であるから、積分の微分を計算すると
$$ \begin{aligned} \left.\frac{d}{ds}\right|_{s=0}L(\gamma_s) &=\int_a^b \left.\frac{1}{2\|T\|}\frac{\partial}{\partial s}g(T,T)\right|_{s=0}dt\\ &=\int_a^b g(D_tV,U)\,dt. \end{aligned} $$
ここで
$$ \frac{d}{dt}g(V,U)=g(D_tV,U)+g(V,D_tU) $$
を積分すれば主張の公式を得る。両端点固定なら境界項は消える。さらに prop-variation-field-realization と cor-energy-critical-geodesic と同じ台関数の論法を用いると、すべての両端点固定変分に対する消失は $D_tU=0$ と同値である。
弧長パラメータ $r$ において $d\gamma/dr=U$ であるから、$D_tU=0$ と $D_r(d\gamma/dr)=0$ は同値である。後者が測地線方程式である。$\blacksquare$

第 1 変分公式の最も直接の応用は、最短の曲線が測地線であることである。

最短曲線は測地線

$\gamma\colon[a,b]\to M$ を滑らかな曲線とする。

  1. $\gamma$ が、同じ端点を結ぶ $[a,b]$ 上の滑らかな曲線の中でエネルギーを最小にするなら、$\gamma$ は測地線である(とくに速さは一定)。
  2. $\gamma$ が正則で、同じ端点を結ぶ滑らかな曲線の中で長さを最小にするなら、$\gamma$ を弧長で再パラメータ化したものは測地線である。

1:両端点固定の任意の変分 $F$ について、$s\mapsto E(\gamma_s)$ は $s=0$ で最小値をとる滑らかな関数なので、その微分は $s=0$ で $0$ である。すべての両端点固定の変分で第 1 変分が $0$ なので、cor-energy-critical-geodesic により $D_t\dot\gamma=0$ である。
2:$\gamma$ は正則なので、$s$ が $0$ に十分近ければ $\gamma_s$ も正則であり($[a,b]$ がコンパクトで $\partial_tF$ が連続なので、$\|\partial_tF\|$ は $s=0$ の近くで正のまま)、$s\mapsto L(\gamma_s)$ は $s=0$ の近くで滑らかで、$s=0$ で最小になる。よって長さの第 1 変分が $0$ であり、thm-first-variation-length により $D_tU=0$、すなわち弧長で再パラメータ化すると測地線である。$\blacksquare$

逆は成り立たない。測地線は停留曲線であるが最短とは限らない(後の ex-sphere-long-great-circle)。

例と反例

Euclid空間の線分

$M=\mathbb R^n$ に標準計量を入れ、$p,q\in\mathbb R^n$ に対して
$$ \gamma(t)=(1-t)p+tq\qquad(0\le t\le1) $$
と置く。$\dot\gamma=q-p$ は定ベクトルなので $D_t\dot\gamma=0$ である。従ってこの線分は両端点固定のエネルギーと長さの双方に対する停留曲線である。さらに
$$ E(\gamma)=\frac12\|q-p\|^2, \qquad L(\gamma)=\|q-p\| $$
であり、prop-energy-length-inequality の等号が成り立つ。

半円は平面では停留しない

$\gamma(t)=(\cos t,\sin t)$($0\le t\le\pi$)をEuclid平面上の半円とする。両端点を固定する変分
$$ F(s,t)=(1+s\sin t)\gamma(t) $$
の変分ベクトル場は $V(t)=\sin t\,\gamma(t)$ である。$D_t\dot\gamma=-\gamma$ なので
$$ \left.\frac{d}{ds}\right|_{s=0}E(\gamma_s) =-\int_0^\pi\langle\sin t\,\gamma(t),-\gamma(t)\rangle\,dt =\int_0^\pi\sin t\,dt=2. $$
よって半円は平面のエネルギーの停留曲線ではない。円弧が「まっすぐ」ではないことが、第1変分の非消失として定量的に現れている。

反例:長さの停留性からアフィン・パラメータは従わない

$\gamma:[1,2]\to\mathbb R^2$ を $\gamma(t)=(t^2,0)$ とする。この曲線の単位接ベクトルは $U=(1,0)$ であり $D_tU=0$ なので、両端点固定の長さの第1変分はすべて消える。しかし
$$ D_t\dot\gamma=(2,0)\ne0 $$
なので、エネルギーの停留曲線ではない。この反例は、長さの停留性が測地線の軌跡を特徴づけるだけで、与えられたパラメータがアフィンであることまでは保証しないことを示す。

区分的に滑らかな曲線

曲線が折点
$$ a=t_0< t_1<\cdots< t_m=b $$
をもち、各開区間で滑らかな場合は、thm-first-variation-energy を各区間に適用して足し合わせればよい。このとき内部の境界項が打ち消し合わず、
$$ \sum_{i=1}^{m-1} g\!\left(V(t_i),\dot\gamma(t_i^-)-\dot\gamma(t_i^+)\right) $$
という折点項が残る。したがって、折点を動かせるすべての変分に対してエネルギーが停留するなら、各区間で測地線方程式が成り立つだけでなく、折点で左右の速度が一致しなければならない。真の停留曲線に角は残らない。
正則な区分的に滑らかな曲線に長さの第1変分公式を各区間で適用すると、内部の折点から
$$ \sum_{i=1}^{m-1} g\!\left(V(t_i),U(t_i^-)-U(t_i^+)\right) $$
が生じる。したがって、折点を動かせるすべての変分に対して長さが停留するなら、左右の単位接ベクトルが一致しなければならない。これにより、2本の直線を折れ曲がって繋いだ道は、各部分では直線であっても全体として長さの停留曲線ではない。

第2変分公式

第 1 変分は停留曲線を検出するが、その曲線が長さやエネルギーを最小にするかどうかはそれだけでは判定できない。通常の 1 変数関数で $f'(x)=0$ が極小を保証しないのと同じである。極小かどうかを見るには 2 階の微分が要り、測地線に沿うエネルギーの 2 階微分には曲率が現れる。
以下、曲率テンソルを $R(X,Y)Z:=\nabla_X\nabla_YZ-\nabla_Y\nabla_XZ-\nabla_{[X,Y]}Z$ とする(リーマン曲率テンソル の記事の定義)。変分 $F$ に沿うベクトル場 $W(s,t)$ に対し、$D_s$、$D_t$ を $s$ 方向・$t$ 方向の共変微分とする。

共変微分の交換と曲率

変分 $F$ に沿う滑らかなベクトル場 $W$ について
$$ D_sD_tW-D_tD_sW=R(S,T)W\qquad\bigl(S=\partial_sF,\ T=\partial_tF\bigr) $$
が成り立つ。

交換の補題の出典

この補題はこの記事では証明しない。局所座標で両辺を書いて比べる計算であり、Mil63 の Lemma 9.2(p. 52)にある(そこでは曲率テンソルの符号が上と逆の規約なので、左辺の順序が入れ替わった形で書かれている)。

エネルギーの第2変分

$\gamma\colon[a,b]\to M$ を測地線、$F$ を $\gamma$ の両端点を固定する変分、$V$ をその変分ベクトル場とする。このとき
$$ \left.\frac{d^2}{ds^2}\right|_{s=0}E(\gamma_s)=\int_a^b\Bigl(\|D_tV\|^2-g\bigl(R(V,\dot\gamma)\dot\gamma,V\bigr)\Bigr)\,dt $$
が成り立つ。右辺を $I(V,V)$ と書き、$\gamma$ の指数形式という。

$S=\partial_sF$、$T=\partial_tF$ とする。prf-first-variation-energy と同じ計算を $s=0$ に限らず行うと、すべての $s$ について
$$ \frac{d}{ds}E(\gamma_s)=\int_a^bg(D_sT,T)\,dt=\int_a^bg(D_tS,T)\,dt $$
である($D_sT=D_tS$)。もう一度 $s$ で微分し、計量両立性と $D_sT=D_tS$ を使うと
$$ \frac{d^2}{ds^2}E(\gamma_s)=\int_a^b\Bigl(g(D_sD_tS,T)+\|D_tS\|^2\Bigr)\,dt $$
となる。lem-second-variation-commutator により $D_sD_tS=D_tD_sS+R(S,T)S$ であり、
$$ \int_a^bg(D_tD_sS,T)\,dt=\bigl[g(D_sS,T)\bigr]_a^b-\int_a^bg(D_sS,D_tT)\,dt $$
である。$s=0$ では $D_tT=D_t\dot\gamma=0$(測地線)なので右辺の積分は消える。端点では $F(s,a)=\gamma(a)$ が $s$ によらないので $S(s,a)=0$ がすべての $s$ で成り立ち、したがって $D_sS(s,a)=0$ である。$t=b$ でも同様なので境界項も消える。最後に $g(R(X,Y)Z,W)=-g(R(X,Y)W,Z)$ という反対称性(リーマン曲率テンソル の記事の命題「曲率テンソルの対称性」の 2。同記事の $\mathrm{Rm}(X,Y,Z,W)=g(R(Z,W)Y,X)$ の前の 2 引数についての反対称性)により $g(R(S,T)S,T)=-g(R(S,T)T,S)$ であるから、$s=0$ で $S=V$、$T=\dot\gamma$ とおいて主張を得る。$\blacksquare$

最短な測地線の指数形式は非負

測地線 $\gamma\colon[a,b]\to M$ が、同じ端点を結ぶ $[a,b]$ 上の滑らかな曲線の中でエネルギーを最小にするなら、$\gamma$ に沿う滑らかなベクトル場 $V$ で $V(a)=V(b)=0$ となるものすべてについて $I(V,V)\geq0$ である。

prop-variation-field-realization により、$V$ を変分ベクトル場とする両端点固定の変分 $F$ がある。$s\mapsto E(\gamma_s)$ は $s=0$ で最小値をとる滑らかな関数なので、その 2 階微分は $s=0$ で $0$ 以上である。thm-second-variation-energy により、それは $I(V,V)$ に等しい。$\blacksquare$

反例:球面の長い大円弧

単位球面 $S^2\subset\mathbb{R}^3$ 上の大円 $\gamma(t)=(\cos t,\sin t,0)$($0\le t\le\ell$)は単位速さの測地線である(測地線 の記事の例「球面の二つの道」)。$W:=(0,0,1)$ は各 $\gamma(t)$ で球面に接し($\gamma(t)$ と直交)、定ベクトルなので $\gamma$ に沿って平行である。$V(t):=\sin(\pi t/\ell)\,W$ とおくと $V(0)=V(\ell)=0$、$D_tV=(\pi/\ell)\cos(\pi t/\ell)\,W$ である。$S^2$ は定曲率 $1$ なので(断面曲率 の記事の例「定曲率空間の代表例」)、同じ記事の系「定曲率の特徴づけ」により $R(V,\dot\gamma)\dot\gamma=g(\dot\gamma,\dot\gamma)V-g(V,\dot\gamma)\dot\gamma=V$ である。したがって
$$ I(V,V)=\int_0^\ell\Bigl(\frac{\pi^2}{\ell^2}\cos^2\frac{\pi t}{\ell}-\sin^2\frac{\pi t}{\ell}\Bigr)dt=\frac{\ell}{2}\Bigl(\frac{\pi^2}{\ell^2}-1\Bigr) $$
である。$\ell>\pi$ ならこれは負なので、cor-second-variation-nonnegative により $\gamma$ は端点を固定したエネルギーの最小曲線でない。
さらに長さについても最短でない。$V$ を変分ベクトル場とする両端点固定の変分 $\gamma_s$ をとると、第 1 変分は $0$、第 2 変分は負なので、$0$ に近い $s\neq0$ で $E(\gamma_s)< E(\gamma)=\ell/2$ である。もし $L(\gamma_s)\geq\ell$ なら prop-energy-length-inequality により $E(\gamma_s)\geq L(\gamma_s)^2/(2\ell)\geq\ell/2$ となり矛盾するので、$L(\gamma_s)< \ell=L(\gamma)$ である。こうして「測地線(停留曲線)ならば最短」という含意が破れる。$\ell=\pi$ では $I(V,V)=0$ となり、ちょうど対蹠点で最短性が失われ始める。

測地線に沿うベクトル場 $J$ で $D_tD_tJ+R(J,\dot\gamma)\dot\gamma=0$ を満たすものを Jacobi 場という。$I(V,V)$ を部分積分すると $-\int_a^bg\bigl(V,D_tD_tV+R(V,\dot\gamma)\dot\gamma\bigr)dt$ となり、Jacobi 場は指数形式の「零方向」を与える。上の例の $V$($\ell=\pi$)は Jacobi 場であり、端点で消える $0$ でない Jacobi 場が存在するとき、両端点は $\gamma$ に沿って共役であるという(Mil63 §14、p. 77)。
リーマン計量の正定値性はここでも重要である。擬リーマン多様体では、時間的・空間的・光的曲線を分ける必要があり、特に光的曲線では $g(\dot\gamma,\dot\gamma)=0$ となる。正定値計量から定まるノルムと長さの公式は、そのままでは使えない。

条件を外したときの反例

外す条件反例成り立たなくなること
両端点の固定$\mathbb{R}^2$ の線分 $\gamma(t)=(t,0)$ と $F(s,t)=((1+s)t,0)$($0\le t\le1$)測地線でも第 1 変分が $0$ になること
曲線の正則性定値曲線 $\gamma\equiv0$ と $F(s,t)=(st,0)$長さ $L(\gamma_s)$ の $s=0$ での微分可能性
パラメータ(長さからエネルギーへ)$\gamma(t)=(t^2,0)$(ex-length-energy-parameter)長さの停留曲線がエネルギーの停留曲線であること
最短性(停留から最短へ)単位球面の長さ $\ell>\pi$ の大円弧(ex-sphere-long-great-circle)測地線が端点間で最短であること
曲線の滑らかさ折れ曲がった 2 本の線分各区間で測地線でも全体が停留すること

1 行目:$\gamma$ は測地線($D_t\dot\gamma=0$)だが、$\gamma_s(t)=((1+s)t,0)$ のエネルギーは $E(\gamma_s)=(1+s)^2/2$ で、$s=0$ での微分は $1\neq0$ である。thm-first-variation-energy の境界項 $[g(V,\dot\gamma)]_0^1=g((1,0),(1,0))=1$ が残るためである。2 行目:$L(\gamma_s)=\int_a^b|s|\,dt=|s|(b-a)$ は $s=0$ で微分できない。5 行目は節「区分的に滑らかな曲線」の折点項で確かめた。

文献

第 1 変分公式と「停留曲線は測地線」は Mil63 の Theorem 12.2・Corollary 12.3(pp. 71–73)と Nic26 の Theorem 5.2.3・Corollary 5.2.5(pp. 213–214)、第 2 変分公式は Mil63 の Theorem 13.1(p. 75)と Nic26 の Theorem 5.2.9(p. 215)、最短な測地線でヘッセ形式が半正定値であることは Mil63 の Lemma 13.6(p. 76)にある。Mil63 はエネルギーを $\int\|\dot\gamma\|^2dt$(係数 $1/2$ なし)と定義しているので、公式の係数が本記事と $2$ 倍異なる。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する