\(% 自動抽出されたマクロ定義
% 元ファイル: Geometry3.tex
% 生成日時: 2026-09-20 23:02:45
% MathJax用の標準コマンド定義
\def\slash{/\mkern-5mu}
\def\v#1{{\bf #1}}
\def\d{\mathrm{d}}
\def\pd{\partial}
\def\vf#1{\bar{\bf #1}}
\def\t#1{\tilde{#1}}
\def\ttheta{\t{\theta}}
\def\tomega{\t{\omega}}
\def\w{\wedge}
\def\tt{\ttheta}
\def\to{\tomega}
\def\dq#1{\d q^{#1}}
\def\pq#1{\frac{\pd}{\pd q^{#1}}}
\def\tv#1{\tilde{{\bf #1}}}
\def\vov#1{\stackrel{#1}{\vee}}
\def\nabsla{\nabla\!\!\!\!\slash}
\def\delsla{\pd\!\!\!\slash}
\def\sla#1{#1\!\!\!\slash}
\def\diffop#1{\frac{\d}{\d #1}}
\def\diffopn#1#2{\frac{\d^#1}{\d #2^#1}}
\def\diff#1#2{\frac{\d #1}{\d #2}}
\def\diffn#1#2#3{\frac{\d^#1 #2}{\d #3^#1}}
\def\pdiffop#1{\frac{\pd}{\pd #1}}
\def\pdiffopn#1#2{\frac{\pd^#1}{\pd #2^#1}}
\def\pdiff#1#2{\frac{\pd #1}{\pd #2}}
\def\pdiffn#1#2#3{\frac{\pd^#1 #2}{\pd #3^#1}}
\def\cdiffop#1{\frac{d}{d #1}}
\def\cdiffopn#1#2{\frac{d^#1}{d #2^#1}}
\def\cdiff#1#2{\frac{d #1}{d #2}}
\def\cdiffn#1#2#3{\frac{d^#1 #2}{d #3^#1}}
\def\tr{\mathrm{tr}}
\def\Tr{\mathrm{Tr}}
\def\dbra#1{\langle\!\bra{#1}}
\def\dket#1{\ket{#1}\!\rangle}
\def\dbraket#1{\langle\!\braket{#1}\!\rangle}
\def\Dbra#1{\left.\Left\langle #1 \Right.\right|}
\def\Dket#1{\left|\Left. #1 \Right\rangle\right.}
\def\Dbraket#1{\Left\langle #1 \Right\rangle}
\def\dlangle{\langle\!\langle}
\def\drangle{\rangle\!\rangle}
\def\Dlangle{\Big\langle\!\!\Big\langle}
\def\Drangle{\Big\rangle\!\!\Big\rangle}
\def\sdual#1{#1^*}
\def\dsdual#1{#1^{**}}
\def\dual{{}^\star}
\def\exdual{{}^\dagger}
\def\combi#1#2{{}_{#1}C_{#2}}
\def\grad{\mathrm{grad}}
\def\rot{\mathrm{rot}}
\def\divergent{\mathrm{div}}
\def\inner#1#2{\Braket{{}#1{},{}#2{}}}
\def\dinner#1#2{\Left\langle {}#1{},{}#2{}\Right\rangle}
\def\tensorUD#1#2#3{#1^{#2}_{\,#3}}
\def\tensorDU#1#2#3{#1_{#2}^{\,#3}}
\def\gv#1{\boldsymbol{#1}}
\def\Lie#1{\mathcal{L}_{#1}}
\def\iTM#1{i_{#1}}
\def\iE#1{\iota_{#1}}
\def\a{a_{\mathrm{n}}}
\def\qmark{\mathrm{n}}
\def\qket#1{\ket{#1}_{\!\qmark}}
\def\qbra#1{{}_{\qmark\!\!}\bra{#1}}
\def\qbraket#1{\braket{#1}_{\!\qmark}}
\def\Left#1#2
\def\ts@r{\nulldelimiterspace=0pt \mathsurround=0pt}
\def\sht@im{#2}
\def\@t{{\mathchoice{\def\@fen{\displaystyle}\k@fel}
{\def\@fen{\textstyle}\k@fel}
{\def\@fen{\scriptstyle}\k@fel}
{\def\@fen{\scriptscriptstyle}\k@fel}}}
\def\g@rin{\ts@r\left\@hat\vphantom{\sht@im}\right.}
\def\k@fel{\setbox0=\hbox{$\@fen\g@rin$}\hbox{
$\@fen \kern.3875\wd0 \copy0 \kern-.3875\wd0
\llap{\copy0}\kern.3875\wd0$}}
\def\pt@h{\mathopen\@t}
\def\Right#1{\let\@hat=#1
\def\st@m{\mathclose\@t}
\st@m\endgroup}
\)
INTRODUCTION
このnoteでは"幾何のお話2"ノート[1]の続編となる、一般的な次元をもつ空間の幾何学を取り上げます。 このnoteでは物理学における一般相対性理論で扱う幾何学的な概念や数学の理解、 それと曲がった空間の幾何学のイメージを掴むこと、またそこで現れる各数学的対象の直観的なイメージを掴むこと。 これらを主な目標とします。 ものの本ではあまり詳しく踏み込んで説明されることが少ないように思われる 正規座標系や曲率などの自然な導入を行って、それらの幾何学的なイメージが自然に持てるような 議論の展開を目指しました。 形式には拘らずに、また数学的な完全な記述もあまり目指さずに書きました。 数学的に完全なあるいは形式的な記述が好みの読者は必要に応じて数学の本を参照することをおすすめします。
ユークリッド空間の中の超曲面を扱う
最初は考えている空間\(M\)がより大きな次元を持つ\(N\)次元ユークリッド空間\(\mathbb {R}^N\)の部分空間の中の ある\(n\)次元の超曲面であるとして論じる。 一般的に\(n\)次元の多様体と呼ばれる空間はそのような大きな次元をもつユークリッド空間の部分空間として扱うことができることが知られている。 ここではそのことを証明しないが、その事実だけを流用する。 最後の方でそのような大きな空間の部分空間としてではなく、それそのものとして多様体として扱うことができることに少しだけ 触れる予定ではあるが、そこから出発するとこのnoteの主目的である曲がった空間の幾何学的なイメージや 各数学的対象の直観的なイメージを掴むことから少し遠ざかってしまう。 曲がった空間に対する素直なイメージを持つためには空間を大きなユークリッド空間に埋め込まれた 超曲面として議論をすることが一番の近道のように思う。
定義 .1.
\(N\)次元ユークリッド空間\(\mathbb {R}^N\)とは\(N\)個の実数の組\(\v {x}=(x_1,x_2,\cdots ,x_N)\)全体からなる集合であり、 それらに通常の実ベクトル空間の構造を入れたものである。
定義 .2.
\(\mathbb {R}^N\)の部分集合\(M\)が\(n\)次元の超曲面であるとは\(M\)の任意の点\(p\)に対して、\(M\)の中での\(p\)の近傍\(U\ (\subset M)\)と \(\mathbb {R}^n\)のある開集合\(V\)との間の滑らかな全単射写像をとることができるときをいう。 この時\(\v {u}=(u_1,\cdots ,u_n)\in V\)を\(U\)の局所座標と呼ぶ。 こう書いた時は\(p\)は\(V\)の中のある点\(\v {u}_p\)に対応していて、\(\v {u}\)は\(p\)のまわりの局所的な座標として捉える。
\(\v {u}=(u_1,\cdots ,u_n)\)を\(p\in M\)のまわりの局所座標とする。 一方、\(M\subset \mathbb {R}^N\)なので\(\v {u}\)の対応する点\(\v {x}=(x_1,x_2,\cdots ,x_N)\in \mathbb {R}^N\)を 取ることができる。この対応を\(\v {x}=\v {x}(\v {u})\)のようにして\(\v {x}\)を\(\v {u}\)の関数として表す。 \(M\)の\(\v {x}\)を始点とする接ベクトルは
\begin{align} \pdiff {\v {x}}{u_1},\ \ \pdiff {\v {x}}{u_2},\ \ \cdots ,\ \ \pdiff {\v {x}}{u_n} \end{align}
の\(n\)個ベクトルの線形結合によって表すことができる。 これらひとつひとつがベクトルであり、
\begin{align} \pdiff {\v {x}}{u_i} = \left ( \pdiff {x_1}{u_i},\pdiff {x_2}{u_i},\cdots ,\pdiff {x_N}{u_i} \right ) \ \ \ (i=1,2,\dots ,n) \end{align}
である。ここでベクトルは縦ベクトルのように書いたり横ベクトルのように書いたりするが適宜読み直してほしい。 必ずしも接ベクトルの基底をこれらに取る必要もないので一般的に接ベクトルの基底を \(\v {e}_1,\v {e}_2,\cdots ,\v {e}_n\)と表すことにする。ただしこれらはどれも直交しているものとする。またその大きさも\(1\)であるように取ることにする。 こうすると\(\v {x}\)を始点とするベクトルで\(M\)に接するものはすべて表すことができるが、それらと直交するベクトルはこれらでは表すことができない。\(\mathbb {R}^N\)の中のベクトルなので、残り\(N-n\)個の基底ベクトルをとることができる。これらを\(\v {n}_1,\v {n}_2,\cdots ,\v {n}_{N-n}\)で表すことにする。これらもお互いに直交していて大きさが\(1\)であるとする。 こうして\(\v {x}\)を始点とする\(\mathbb {R}^N\)の中のベクトルの正規直交基底\(\v {e}_1,\cdots ,\v {e}_n,\v {n}_1,\cdots ,\v {n}_{N-n}\)が出そろった。 一般的に\(\v {x}\)を始点とするどの接ベクトルとも直交するベクトルを\(\v {x}\)を始点とする\(M\)の法線ベクトルと呼ぶ。
超曲面\(M\)の中の曲線\(C\)上の位置は\(\mathbb {R}^N\)の中での座標を用いて\(\v {x}=\v {x}_t\)のように表される。速度ベクトルは\(\v {v}=\diff {\v {x}}{t}\)であり、その曲線の長さは
\begin{align} s=\int ^t_0|\v {v}|\d t=\int ^t_0\sqrt {\dot {x}_1^2+\dot {x}_1^2+\cdots +\dot {x}_N^2}\d t\ \ \ \ (\d s=|\v {v}|\d t) \end{align}
である。速度ベクトル\(\v {v}\)に対して、パラメータを\(t\)の代わりに\(s\)で取り直したものを\(\v {e}\)とする。即ち
\begin{align} \v {e}=\diff {\v {x}}{s}=\frac {\v {v}}{|\v {v}|} \end{align}
これは曲線\(C\)に接した大きさが\(1\)のベクトルである。\(|\v {e}|^2=\inner {\v {e}}{\v {e}}=1\)を\(s\)で微分すると \(\inner {\v {e}}{\diff {\v {e}}{s}}=0\)が得られるので\(\v {e}\)とその\(s\)での微分\(\diff {\v {e}}{s}\)とは直交している。 従って
\begin{align} \diff {\v {e}}{s} = \kappa _g\v {e}_g+\kappa _n\v {n}_n \end{align}
のように表すことができる。ここで\(\v {e}_g\)は\(\v {e}\)と直交する\(M\)の単位接ベクトルであり、\(\v {n}_n\)は単位法線ベクトルである。 \(2\)次元の曲面の中の曲線と同じで、\(\kappa _g=0\)となる曲線を測地線という。 \(\kappa _g\v {e}_g\)が加速度ベクトルであり、\(\kappa _n\v {n}_n\)が拘束力と解釈できる。 \(\kappa _g\)が加速度の大きさであり、\(\kappa _n\)が拘束力の大きさを表す。
曲線の線素\((\d s)^2\)は
\begin{align} (\d s)^2 &= |\v {v}|^2(\d t)^2 \nonumber \\ &= \inner {\diff {\v {x}}{t}}{\diff {\v {x}}{t}}(\d t)^2 \nonumber \\ &= \sum _{i,j=1,\cdots ,n}\inner {\pdiff {\v {x}}{u^i}}{\pdiff {\v {x}}{u^j}} \diff {u^i}{t}\diff {u^j}{t}(\d t)^2 \end{align}
であることから形式的な線素の表式は
\begin{align} (\d s)^2 = \sum _{i,j=1,\cdots ,n}\inner {\pdiff {\v {x}}{u^i}}{\pdiff {\v {x}}{u^j}} \d u^i \d u^j \end{align}
と表される。 この時
\begin{align} g_{ij}=\inner {\pdiff {\v {x}}{u^i}}{\pdiff {\v {x}}{u^j}} \end{align}
を計量と呼ぶ。ここで計量は\(\v {x}\)と\(\v {u}\)との間の関係のみによって決まるということに注意されたい。つまり\(s\)や\(t\)といった他の要素には依存しない。
さて\(\v {e}\)は\(\v {x}\)を始点とした接ベクトルなので \((\v {u}_1,\cdots ,\v {u}_n):=\left (\pdiff {\v {x}}{u^1},\cdots ,\pdiff {\v {x}}{u^n}\right )\)の線形結合で表すことができる。即ち
\begin{align} \v {e} = v^1\v {u}_1+v^2\v {u}_2+\cdots +v^n\v {u}_n, \ \ \ \ \ \ v^i=\diff {u^i}{s} \end{align}
である。ここで\(\v {e}\)の\(s\)での微分は
\begin{align} \diff {\v {e}}{s} = \sum _{i=1,\cdots ,n}\left (\diff {v^i}{s}\v {u}_i+v^i\diff {\v {u}_i}{s}\right ) \label {eq:de/ds} \end{align}
である。また\(\diff {\v {u}_i}{s}\)は一般的に\(M\)の接ベクトルとも法線ベクトルとも限らないので、
\begin{align} \diff {\v {u}_i}{s} = \sum _{j=1,\cdots ,n}\tensorDU {w}{i}{j}\v {u}_j +\sum _{\alpha =1,\cdots ,N-n}\tensorDU {\sigma }{i}{\alpha }\v {n}_\alpha \label {eq:du/ds} \end{align}
と書ける。従って\(\diff {\v {e}}{s}\)は
\begin{align} \diff {\v {e}}{s} = \sum _{i,j=1,\cdots ,n}\left (\diff {v^i}{s}+v^j\tensorDU {w}{j}{i}\right )\v {u}_i +\sum _{\substack {i=1,\cdots ,n\\ \alpha =1,\cdots ,N-n}}v^i\tensorDU {\sigma }{i}{\alpha }\v {n}_\alpha \label {eq:de/ds2} \end{align}
この右辺の\(\v {u}_i\)成分を取り出して書くと
\begin{align} \diff {v^i}{s}+\sum _{j=1,\cdots ,n}v^j\tensorDU {w}{j}{i}\ \ \ \ (i=1,\cdots ,n) \label {eq:kasokudoseibun} \end{align}
一方で\(\v {u}_i\)の\(s\)での微分を具体的に書くと
\begin{align} \diff {\v {u}_i}{s} &= \diffop {s}\left (\pdiff {\v {x}}{u^i}\right ) \nonumber \\ &= \sum _{j=1,\cdots ,n}\frac {\pd ^2 \v {x}}{\pd u^i\pd u^j}\diff {u^j}{s} \nonumber \\ &= \sum _{j=1,\cdots ,n}\frac {\pd ^2\v {x}}{\pd u^i\pd u^j}v^j \label {eq:du/ds2} \end{align}
ここで(15)式と(11)式と見比べるとそれぞれの右辺が等しいことが分かる。\(\diff {\v {u}_i}{s}\)と\(\v {u}_1,\cdots ,\v {u}_n\)との内積をとると
\begin{align} \begin {pmatrix} \inner {\v {u}_1}{\v {u}_1} & \cdots & \inner {\v {u}_1}{\v {u}_n} \\ \vdots & \ddots & \vdots \\ \inner {\v {u}_n}{\v {u}_1} & \cdots & \inner {\v {u}_n}{\v {u}_n} \end {pmatrix} \begin {pmatrix} \tensorDU {w}{i}{1} \\ \vdots \\ \tensorDU {w}{i}{n} \end {pmatrix} = \begin {pmatrix} \inner {\v {u}_1}{\frac {\pd ^2\v {x}}{\pd u^i\pd u^1}} & \cdots & \inner {\v {u}_1}{\frac {\pd ^2\v {x}}{\pd u^i \pd u^n}} \\ \vdots & \ddots & \vdots \\ \inner {\v {u}_n}{\frac {\pd ^2\v {x}}{\pd u^i\pd u^1}} & \cdots & \inner {\v {u}_n}{\frac {\pd ^2\v {x}}{\pd u^i \pd u^n}} \end {pmatrix} \begin {pmatrix} v^1 \\ \vdots \\ v^n \end {pmatrix} \end{align}
となる。ここで左辺の最初の行列は計量\(g_{ij}=\inner {\v {u}_i}{\v {u}_j}=\inner {\pdiff {\v {x}}{u^i}}{\pdiff {\v {x}}{u^j}}\)になっている。 計量の逆行列を\(g^{ij}\)で表すと、 \(\pdiffop {u^i}\inner {\v {u}_r}{\v {u}_k}=\inner {\v {u}_r}{\frac {\pd ^2\v {x}}{\pd u^i\pd u^k}}+\inner {\v {u}_k}{\frac {\pd ^2\v {x}}{\pd u^i\pd u^r}}\)に注意すると
\begin{align} \tensorDU {w}{i}{j} &= \sum _{r,k=1,\cdots ,n} g^{jr}\inner {\v {u}_r}{\frac {\pd ^2\v {x}}{\pd u^i\pd u^k}}v^k \nonumber \\ &= \frac {1}{2}\sum _{r,k=1,\cdots ,n} g^{jr}\left ( \pdiff {g_{rk}}{u^i}+\pdiff {g_{ri}}{u^k}-\pdiff {g_{ik}}{u^r} \right )v^k \end{align}
であることが分かる。 ここで
\begin{align} \tensorDU {w}{i}{j} = \sum _{k=1,\cdots ,n}\tensorUD {\Gamma }{j}{ik}v^k \label {ref:setsuzoku} \end{align}
のように書くようにすれば
\begin{align} \tensorUD {\Gamma }{i}{jk} = \frac {1}{2}\sum _{r=1,\cdots ,n} g^{ir}\left ( \pdiff {g_{rk}}{u^j}+\pdiff {g_{rj}}{u^k}-\pdiff {g_{jk}}{u^r} \right ) \end{align}
と表すことができる。\(\tensorUD {\Gamma }{i}{jk}\)はChristoffelの記号と呼ばれる。 また(18)の左辺の\(\tensorDU {w}{i}{j}\)はLevi-Civita接続と呼ばれる。特に誤解のないときには接頭語を省略して単に接続と呼ぶこともある。 このChristoffelの記号は計量\(g_{ij}\)とその逆行列\(g^{ij}\)のみによって決まることに注意されたい。 計量と従ってその逆行列は\(\v {x}\)と\(\v {u}\)との間の関係のみによって決まるのであった。 従ってChristoffelの記号もまた\(\v {x}\)と\(\v {u}\)との間の関係のみによって決まる。 このことはLevi-Civita接続\(\tensorDU {w}{i}{j}\)が単にベクトル\(\v {e}\)の成分\(v^1,\cdots ,v^n\)の 一次結合によって表すことができて、\(v^1,\cdots ,v^n\)への依存がそれのみであることを意味する。 Christoffelの記号は下付きの添え字に関して対称である。即ち
\begin{align} \tensorUD {\Gamma }{i}{jk}=\tensorUD {\Gamma }{i}{kj} \end{align}
である。
\(\v {n}_\alpha \)と\(\diff {\v {e}}{s}\)との内積をとると(10)式と(15)とを使うと
\begin{align} \inner {\v {n}_\alpha }{\diff {\v {e}}{s}} =\sum _{j,k=1,\cdots ,n} \inner {\v {n}_\alpha }{\frac {\pd ^2\v {x}}{\pd u^j\pd u^k}}v^jv^k \ \ \ \ (\alpha =1,\cdots ,N-n) \label {ref:nofde/ds} \end{align}
を得る。従って\(\diff {\v {e}}{s}\)の法線成分\(\v {n}\)は
\begin{align} \v {n} = \sum _{\alpha =1,\cdots ,N-n}\inner {\v {n}_\alpha }{\diff {\v {e}}{s}}\v {n}_\alpha =\sum _{\substack {j,k=1,\cdots ,n\\ \alpha =1,\cdots ,N-n}} \inner {\v {n}_\alpha }{\frac {\pd ^2\v {x}}{\pd u^j\pd u^k}}v^jv^k\v {n}_\alpha \end{align}
となる。 これは\(2\)次元の曲面上の曲線の法曲率に対応するもので、第二基本形式と呼ばれる。 これは上記で述べたように曲線上を運動する質点に働く拘束力を表している。 (13)式と(21)式とを見比べると
\begin{align} \tensorDU {\sigma }{i}{\alpha } = \sum _{j=1,\cdots ,n} \tensorUD {\sigma }{\alpha }{ij}v^j, \ \ \ \ \ \ \
\tensorUD {\sigma }{\alpha }{ij} = \inner {\v {n}_\alpha }{\frac {\pd ^2\v {x}}{\pd u^i\pd u^j}} \label {eq:housenvector} \end{align}
の関係があることが分かる。
加速度が\(0\)の場合は\(\diff {\v {e}}{s}\)の接線成分が\(0\)である場合なので、 (14)式より
\begin{align} \diff {v^i}{s}+\sum _{j=1,\cdots ,n}v^j\tensorDU {w}{j}{i}=0\ \ \ \ (i=1,\cdots ,n) \end{align}
が得られる。これは測地線の方程式と呼ばれる。 あるいは(18)を使って書き下すと
\begin{align} \diff {v^i}{s}+\sum _{j,k=1,\cdots ,n}\tensorUD {\Gamma }{i}{jk}v^jv^k=0\ \ \ \ (i=1,\cdots ,n) \label {eq:sokuchisen} \end{align}
と書ける。 この条件を満たす場合はその曲線は測地線である。 言い換えると測地線に沿って運動する質点には加速度が働かず、拘束力のみが働いている。 測地線の方程式は超曲面\(M\)の座標\(u^i\)に関しては\(2\)階の微分方程式なので(\(v^i=\diff {u^i}{s}\)に注意)初期条件として\(s=0\)の時の初期位置\(p\)と、初速度ベクトル\(\v {v}_0\)が与えられれば方程式を一意に解くことができる。測地線の方程式の解はパラメータ\(s\)に関してどこまでも延長できるとは限らないが、\(s=0\)の近傍には延長できるので初期値\(\v {v}_0\)の大きさが十分に小さければ\(s=1\)まで延長できる。このことは測地線の方程式で\(v^i\)を\(1/C\)倍にすればパラメータ \(s\)は\(C\)倍に取り直すことで同じ測地線の方程式となることからすぐに理解できる。そうすると\(p\)での超曲面\(M\)の接ベクトル\(\v {v}_0\)を、初期値が\(\v {v}_0\)の測地線の方程式の解の\(s=1\)での位置に対応させる写像\(\exp _p:\v {v}_0\rightarrow \v {u}\)を考えることができる。 ここまでの話からもちろん点\(p\)の全ての接ベクトルに対して定義できるとは限らないので、一般には点\(p\)の近傍にある接ベクトルに対して定義される。 この対応関係を正規座標あるいは測地座標という。
次に超曲面\(M\)に接している一般的なベクトル\(\v {a}\)について同じように微分をすることを考える。 ベクトル\(\v {a}\)は\(M\)に接しているので接ベクトル\(\v {u}_1,\cdots ,\v {u}_n\)の線形結合で書ける。
\begin{align} \v {a}=\sum _{i=1,\cdots ,n}a^i\v {u}_i \end{align}
今度はこのベクトルは曲線\(C\)に沿っているとは限らない。 このベクトルを曲線\(C\)に沿った\(s\)方向への微分をとると
\begin{align} \diff {\v {a}}{s} &= \sum _{i,j=1,\cdots ,n} \left ( \diff {a^i}{s}+a^j\tensorDU {w}{j}{i} \right )\v {u}_i + \sum _{\substack {i=1,\cdots ,n\\ \alpha =1,\cdots ,N-n}} a^i\tensorDU {\sigma }{i}{\alpha }\v {n}_\alpha \end{align}
と書ける。ここで\(\tensorDU {w}{j}{i}\)は(18)と同一のものである。 この時測地線と同じように接線成分が\(0\)であるとき、即ち
\begin{align} \diff {a^i}{s}+\sum _{j,k=1,\cdots ,n}\tensorUD {\Gamma }{i}{jk}a^jv^k=0\ \ \ \ (i=1,\cdots ,n) \label {eq:heikouidou} \end{align}
の時、ベクトル\(\v {a}\)は曲線\(C\)に沿って平行であるという。 \(\v {a}\)が曲線\(C\)に沿って平行であるとき、\(\diffop {s}\inner {\v {a}}{\v {a}}=0\)となるので、 曲線\(C\)に沿ってベクトル\(\v {a}\)を動かしていってもその大きさが変わらない。 逆にある点\(\v {p}\in M\)において超曲面\(M\)に接する接ベクトル\(\v {a}_0\)が与えられた時に、 平行移動の微分方程式(30)式と初期条件\(s=0\)で\(\v {a}=\v {a}_0\)で\(s\)の一階微分方程式を解くことができる。こうして与えられるベクトル\(\v {a}\)は曲線\(C\)に沿って平行なベクトルとなっている。 一般的に任意に与えられた超曲面\(M\)に接しているベクトル場\(\v {a}\)に対しては、\(s\)方向の微分は平行移動になっているとも限らず、また測地線の接ベクトルとも限らないので\(\v {u}_i\)の成分は\(0\)とはならない。 平行移動という意味を考えてみる。超曲面\(M\)に接しているベクトル\(\v {a}\)を超曲面内の曲線\(C\)に沿って移動させることを考える時、その向きや大きさが変えずに移動するというのは\(\diff {\v {a}}{s}\)の\(M\)に接する成分がすべて\(0\)であることを意味する。そのような移動は平面\(\mathbb {R}^2\)内であれば平行移動と呼ぶので、そのような平行移動の概念の曲がった超曲面\(M\)への拡張とも言える。平面内であればベクトルを平行移動させてまた同じ位置に戻ってくれば、その大きさや向きが変化していないので最初と同じベクトルと同じ向きと大きさになっている。しかし一般的に曲がった超曲面\(M\)であれば一般的にそうとは限らない。移動の際にその向きや大きさを変えないように移動させていたとしてもである。その辺のところを以下で見ていくことにする。
その前にまず以下の計算で必要なものを揃えておく。接ベクトルの基底\(\v {u}_i\)の\(u^a\ (a=1,\cdots ,n)\)方向微分は
\begin{align} \pdiff {\v {u}_i}{u^a} = \sum _{j=1,\cdots ,n}\tensorUD {\Gamma }{j}{ia}\v {u}_j + \sum _{\alpha =1,\cdots ,N-n}\tensorUD {\sigma }{\alpha }{ia}\v {n}_\alpha \ \ \ \ (i,a=1,\cdots ,n) \end{align}
のようになるのであった。 ここで、今度は\(\v {n}_\alpha \)の\(u^a\)方向の微分を考える。\(u^a\)で微分すると
\begin{align} \pdiff {\v {n}_\alpha }{u^a} = \sum _{j=1,\cdots ,n}\tensorUD {\rho }{j}{\alpha a}\v {u}_j + \sum _{\beta =1,\cdots ,N-n}\tensorUD {\kappa }{\beta }{\alpha a}\v {n}_\beta \ \ \ \ (a=1,\cdots ,n,\ \alpha =1,\cdots ,N-n) \end{align}
のようになる。ここで\(\inner {\v {u}_i}{\v {n}_\alpha }=0\)より、これを\(u^a\)で微分すると\(\inner {\v {u}_i}{\v {u}_j}=g_{ij}\)に注意すれば
\begin{align} \tensorUD {\rho }{i}{\alpha a} &= -\sum _{j=1,\cdots ,n}g^{ij}\tensorUD {\sigma }{\alpha }{ja} \label {ref:rho-sigma} \end{align}
であることが分かる。 また\(\inner {\v {n}_\alpha }{\v {n}_\beta }=\delta _{\alpha \beta }\)を\(u^a\)で微分すると
\begin{align} \tensorUD {\kappa }{\alpha }{\beta a}=-\tensorUD {\kappa }{\beta }{\alpha a}\ \ \ \ (a=1,\cdots ,n,\ \alpha ,\beta =1,\cdots ,N-n) \end{align}
であることが分かる。 ここで計量の\(g_{ij}\)や計量の逆\(g^{ij}\)で添え字の上げ下げを行うという規約を導入する。即ち\(v_i=\sum _{j}g_{ij}v^j\)や \(v^i=\sum _{j}g^{ij}v_j\)のようなものを表すと決める。こう決めると (33)の右辺は\(-\tensorUD {\sigma }{\alpha i}{a}\)と簡単に表すことができる。
さて超曲面\(M\)に接するベクトルを平行移動して元の位置に戻った時のずれについて考察をしてみる。そのためには接ベクトルの\(2\)回微分を見てみると良い。\(\v {v}=v^1\v {u}_1+\cdots +v^n\v {u}_n\)に対して\(u^a\)の方向微分を行った後に\(u^b\)方向微分をとると
\begin{align} \pdiffop {u^b}\left (\pdiff {\v {v}}{u^a}\right ) &=\!\!\!\! \sum _{i=1,\cdots ,n} \pdiffop {u^b} \left [ \left ( \pdiff {v^i}{u^a} + \sum _{j=1,\cdots ,n}\tensorUD {\Gamma }{i}{ja}v^j \right ) \v {u}_i \right ] +\!\!\!\! \sum _{\substack {i=1,\cdots ,n\\ \alpha =1,\cdots ,N-n}} \pdiffop {u^b} \left ( \tensorUD {\sigma }{\alpha }{ia}v^i\v {n}_\alpha \right ) \nonumber \\ &=\!\!\!\! \sum _{i=1,\cdots ,n}\! \left ( \frac {\pd ^2 v^i}{\pd u^b\pd u^a} +\!\!\!\!\!\! \sum _{j=1,\cdots ,n}\!\!\!\! \pdiff {v^j}{u^a} \tensorUD {\Gamma }{i}{jb} +\!\!\!\!\!\! \sum _{j=1,\cdots ,n}\!\!\!\! \pdiff {v^j}{u^b} \tensorUD {\Gamma }{i}{ja} +\!\!\!\!\!\! \sum _{j=1,\cdots ,n}\!\!\!\! \pdiff {\tensorUD {\Gamma }{i}{ja}}{u^b}v^j +\!\!\!\!\!\!\!\! \sum _{j,k=1,\cdots ,n}\!\!\!\!\!\! \tensorUD {\Gamma }{k}{ja}\tensorUD {\Gamma }{i}{kb}v^j -\!\!\!\!\!\!\!\!\!\!\! \sum _{\substack {j=1,\cdots ,n\\ \alpha =1,\cdots ,N-n}}\!\!\!\!\!\!\!\! \tensorUD {\sigma }{\alpha }{ja}\tensorUD {\sigma }{\alpha i}{b} v^j \right ) \v {u}_i \nonumber \\ & \hspace {2cm} +\!\!\!\!\!\! \sum _{\alpha =1,\cdots ,N-n} \left [ \sum _{i=1,\cdots ,n}\!\! \left ( \pdiff {v^i}{u^a} +\!\!\!\!\!\! \sum _{j=1,\cdots ,n}\!\!\!\! \tensorUD {\Gamma }{i}{ja}v^j \right ) \tensorUD {\sigma }{\alpha }{ib} +\!\!\!\!\!\! \sum _{i=1,\cdots ,n}\!\!\!\! \pdiffop {u^b} \left ( \tensorUD {\sigma }{\alpha }{ia}v^i \right ) +\!\!\!\!\!\!\!\! \sum _{\substack {i=1,\cdots ,n\\ \beta =1,\cdots ,N-n}}\!\!\!\!\!\! \tensorUD {\sigma }{\beta }{ia}\tensorUD {\kappa }{\alpha }{\beta b}v^i \right ] \v {n}_\alpha \label {eq:d^2v} \end{align}
のようになる。ここで\(\v {v}\)の\(u^a\)での微分と\(u^b\)での微分はその微分をする順番を変えても変わらないので\(a\)と\(b\)を入れ替えたものとの差をとると\(0\)となる。一方で接線成分、即ち(38)式の右辺の\(2\)段目の\(\v {u}_i\)の係数のカッコの中の最初の\(3\)項は\(a\)と\(b\)を入れ替えたものとの差をとると消えてしまう。従って
\begin{align} \pdiff {\tensorUD {\Gamma }{i}{lj}}{u^k} - \pdiff {\tensorUD {\Gamma }{i}{kj}}{u^l} + \sum _{m=1,\cdots ,n} \left ( \tensorUD {\Gamma }{m}{jl}\tensorUD {\Gamma }{i}{km} - \tensorUD {\Gamma }{m}{jk}\tensorUD {\Gamma }{i}{lm} \right ) = \sum _{\alpha =1,\cdots ,N-n} \left ( \tensorUD {\sigma }{\alpha }{jl}\tensorUD {\sigma }{\alpha i}{k} - \tensorUD {\sigma }{\alpha }{jk}\tensorUD {\sigma }{\alpha i}{l} \right ) \end{align}
が分かる。この式の左辺を(従って右辺を)\(\tensorUD {R}{i}{jkl}\)と置いてRiemann曲率テンソルと呼ぶ。即ち
\begin{align} \tensorUD {R}{i}{jkl} &= \pdiff {\tensorUD {\Gamma }{i}{lj}}{u^k} - \pdiff {\tensorUD {\Gamma }{i}{kj}}{u^l} + \sum _{m=1,\cdots ,n} \left ( \tensorUD {\Gamma }{m}{jl}\tensorUD {\Gamma }{i}{km} - \tensorUD {\Gamma }{m}{jk}\tensorUD {\Gamma }{i}{lm} \right ) \nonumber \\ &= \sum _{\alpha =1,\cdots ,N-n} \left ( \tensorUD {\sigma }{\alpha }{jl}\tensorUD {\sigma }{\alpha i}{k} - \tensorUD {\sigma }{\alpha }{jk}\tensorUD {\sigma }{\alpha i}{l} \right ) \label {eq:kyokuritsu} \end{align}
である。従って添え字の\(k\)と\(l\)の入れ替えで反対称である。即ち
\begin{align} \tensorUD {R}{i}{jkl} = -\tensorUD {R}{i}{jlk} \end{align}
である。また(40)の右辺の二段目の式を見ると、上付きの添え字を計量で下げたもの、即ち\(R_{ijkl}=\sum _{m}g_{im}\tensorUD {R}{m}{jkl}\)に対してもうひとつの反対称性と対称性
\begin{align} R_{ijkl}=-R_{jikl}\ , \ \ \ \ \ \
R_{ijkl}=R_{klij} \end{align}
の関係があることが分かる。 (40)の右辺の\(2\)段目の等式はGaussの方程式と呼ばれる。 ここで式の見通しをよくするために\(\nabla _a\)を
\begin{align} \nabla _a\v {v} = \sum _{i=1,\cdots ,n} \left ( \pdiff {v^i}{u^a} + \sum _{j=1,\cdots ,n} \tensorUD {\Gamma }{i}{ja}v^j \right )\v {u}_i \end{align}
として定義する。\(M\)に沿ったベクトル\(\v {a}\)の積分曲線のパラメータを\(s\)とすれば
\begin{align} \nabla _{\v {a}}\v {v} = \sum _{i=1,\cdots ,n} \left ( \diff {v^i}{s} + \sum _{j,k=1,\cdots ,n} \tensorUD {\Gamma }{i}{jk}v^ja^k \right ) \v {u}_i \end{align}
のように表すこともある。この記法を使うとRiemannの曲率テンソルは
\begin{align} \sum _{i,j=1\cdots ,n} \tensorUD {R}{i}{jkl}v^j\v {u}_i = \left ( \nabla _k\nabla _l-\nabla _l\nabla _k \right ) \v {v} \end{align}
と書くこともできる。 また(38)の右辺の\(2\)段目の法線成分に関して同様に\(a\)と\(b\)を入れ替えたものの差をとることで
\begin{align} \pdiff {\tensorUD {\sigma }{\alpha }{ia}}{u^b} - \pdiff {\tensorUD {\sigma }{\alpha }{ib}}{u^a} + \sum _{j=1,\cdots ,n} \left ( \tensorUD {\Gamma }{j}{ia}\tensorUD {\sigma }{\alpha }{jb} - \tensorUD {\Gamma }{j}{ib}\tensorUD {\sigma }{\alpha }{ja} \right ) = \sum _{\beta =1,\cdots ,N-n} \left ( \tensorUD {\sigma }{\beta }{ib}\tensorUD {\kappa }{\alpha }{\beta a} - \tensorUD {\sigma }{\beta }{ia}\tensorUD {\kappa }{\alpha }{\beta b} \right ) \end{align}
を得る。この関係式をCodazzi–Mainardiの方程式という。
さて、超曲面\(M\)上の点\(p\)での接ベクトル\(\v {v}_p\)と、\(p\)の近くにある別の\(M\)上の点\(q\)での接ベクトル\(\v {v}_q\)とを比較する時に\(M\)上の座標系\(\v {u}=(u^1,\cdots ,u^n)\)を基準に比較したときにはベクトルの成分を比較することになる。また比較する際にはベクトルを一方から他方へ平行移動する。従って点\(p\)から点\(q\)へ平行移動したときのベクトルの成分の差分は
\begin{align} v^i(q)-v^i(p) &= \sum _{j=1,\cdots ,n}\pdiff {v^i}{u^j}(p)\Delta u^j \nonumber \\ &= -\sum _{j,k=1,\cdots ,n}\tensorUD {\Gamma }{i}{jk}(p)v^j(p)\Delta u^k \end{align}
となる。ここで\(\Delta \v {u}=(\Delta u^1,\cdots ,\Delta u^n)\)は点\(p\)と点\(q\)との\(M\)上の座標の差分である。 次にさらに\(q\)の近くの\(p\)とも\(q\)とも異なる点\(r\)を取る。 点\(r\)と点\(q\)との\(M\)上の座標の差分を\(\delta \v {u}=(\delta u^1,\cdots ,\delta u^n)\)と表すと、 点\(q\)から点\(r\)へベクトルを平行移動したときのベクトルの成分は
\begin{align} v^i(r) &= v^i(q) - \sum _{j,k=1,\cdots ,n} \tensorUD {\Gamma }{i}{jk}(q)v^j(q)\delta u^k \nonumber \\ &= v^i(p) - \sum _{j,k=1,\cdots ,n} \left ( \tensorUD {\Gamma }{i}{jk}(q)v^j(q)\delta u^k + \tensorUD {\Gamma }{i}{jk}(p)v^j(p)\Delta u^k \right ) \nonumber \\ &= v^i(p) - \sum _{j,k=1,\cdots ,n} \left ( \tensorUD {\Gamma }{i}{jk}(q)v^j(p)\delta u^k -\!\!\!\! \sum _{k,l=1,\cdots ,n} \tensorUD {\Gamma }{i}{jk}(p)\tensorUD {\Gamma }{j}{lm}(p)v^l(p) \delta u^k\Delta u^m + \tensorUD {\Gamma }{i}{jk}(p)v^j(p)\Delta u^k \right ) \nonumber \\ &= v^i(p) -\!\!\!\!\!\!\!\!\!\! \sum _{j,k=1,\cdots ,n}\!\! \left ( \tensorUD {\Gamma }{i}{jk}(p)v^j(p)\delta u^k +\!\!\!\!\!\! \sum _{l=1,\cdots ,n}\!\!\!\! \pdiff {\tensorUD {\Gamma }{i}{jk}}{u^l}(p)v^j(p)\delta u^k\Delta u^l -\!\!\!\!\!\!\!\! \sum _{l,m=1,\cdots ,n}\!\!\!\!\!\! \tensorUD {\Gamma }{i}{jk}(p)\tensorUD {\Gamma }{j}{lm}(p) v^l(p) \delta u^k\Delta u^m + \tensorUD {\Gamma }{i}{jk}(p)v^j(p)\Delta u^k \right ) \end{align}
となる。これに対して点\(p\)から\(\delta \v {u}\)の分だけ異なる点を新たに点\(s\)とする。点\(s\)と点\(r\)とは\(\Delta \v {u}\)だけ座標が異なる点となる。 この時ベクトルを\(p\)から\(s\)へ平行移動してそれから\(r\)へ平行移動したベクトル\(\v {v}^\prime \)の成分も同様になるが\(\Delta u^i\)と\(\delta u^i\)とが対称的に現れる。 これらは一致せずにその差は、最初の点\(q\)を経由して平行移動したベクトルを\(v_1^i\)とし、点\(s\)を経由して平行移動したベクトルを\(v_2^i\)とすると
\begin{align} v_1^i(r)-v_2^i(r) = \sum _{j,k=1,\cdots ,n} \left ( \pdiff {\tensorUD {\Gamma }{i}{jk}}{u^l}(p) - \pdiff {\tensorUD {\Gamma }{i}{jl}}{u^k}(p) + \sum _{l,m=1,\cdots ,n} \left ( \tensorUD {\Gamma }{i}{ml}(p)\tensorUD {\Gamma }{m}{jk}(p) - \tensorUD {\Gamma }{i}{mk}(p)\tensorUD {\Gamma }{m}{jl}(p) \right ) \right ) v^j(p)\delta u^k\Delta u^l \label {eq:vector-zure} \end{align}
となる。右辺のカッコの中に現れたのはまさにRiemannの曲率テンソルである。 これはベクトルを超曲面の閉路に沿って平行移動して元の位置まで戻した時に元のベクトルからベクトルがずれてしまうことを意味する。またそのずれがRiemannの曲率テンソルで計ることができることを意味する。 [2]
ここまでの結果を見返してみる。\(3\)次元ユークリッド空間の中の\(2\)次元の曲面の時は、これはそのまま幾何のお話2[1]ノートの内容になっている。その場合(40)式の右辺の二段目はガウス曲率(正確にはそれを計量の行列式で割ったもの)に等しい。 こうしてみると(40)式の右辺の一段目の式はChristoffelの記号で、従って計量のみを用いて表すことができる。即ち曲面のガウス曲率は曲面上の計量が与えられると完全に決まる。 これはガウスが "Theorema Egregium"(驚異の定理)と呼んだものである。 1
このことは一般的な\(N\)次元ユークリッド空間の中の\(2\)次元の曲面であっても同じようなことが言える。以下でそれを詳しく見ていく。 \(2\)次元の曲面の時は曲率の添え字は\(1\)と\(2\)のみであることと、曲率\(R_{ijkl}\)の添え字に関する反対称性を考慮すると実質的に独立な成分が\(1\)つしかないことが分かる。 即ちそのことを表すと \begin{align} R_{ijkl} = \sigma \cdot \varepsilon _{ij}\varepsilon _{kl} \ \ \ \ (i,j,k,l=1,2) \label {eq:2d-kyokuritsu} \end{align}
のように表される。ここで\(\sigma \)は\(\tensorUD {\sigma }{\alpha }{ij}\)からなる行列の\(\alpha \)に関する和、即ち
\begin{align} \sigma = \sum _{\alpha =1,\cdots ,N-2}\det \left (\tensorUD {\sigma }{\alpha }{ij}\right )_{ij} \ \ \ \ (i,j=1,2) \end{align}
である。また\(\varepsilon _{ij}\ (i,j=1,2)\)は
\begin{align} \varepsilon _{ij} = \left \{ \begin {array}{ll} 1 & :i=1,j=2\\ -1 &:i=2,j=1\\ 0 & :i=j \end {array} \right . \end{align}
である。いったんまた話を\(n\)次元の超曲面\(M\)に戻して、 ここまでは基底ベクトルとして\(\v {u}_i=\pdiff {\v {x}}{u^i}\)を選んでいたが、代わりに局所正規直交基底を選ぶ。 即ち点\(p\in M\)のでの接ベクトルの基底として\(\v {e}_1,\cdots ,\v {e}_n\)をそれらが互いに直交していて、かつ大きさが\(1\)であるベクトルを選ぶ。\(\v {u}_i\)と\(\v {e}_i\)との間の関係はベクトルの基底の取り換えなので次のような線形な関係にある。
\begin{align} (\v {u}_1\cdots \v {u}_n) = (\v {e}_1\cdots \v {e}_n) \mathcal {A} \end{align}
ここで\((\v {e}_1\cdots \v {e}_n)\)や\((\v {u}_1\cdots \v {u}_n)\)は各ベクトル\(\v {e}_i\)、\(\v {u}_i\)を縦ベクトルとして扱って並べてできる\(n\times n\)行列を表す。また\(\mathcal {A}\)は基底の変換行列であり、一般に\(\v {u}=(u^1,\cdots ,u^n)\)に依存する\(n\times n\)行列である。 \(\v {E}=(\v {e}_1\cdots \v {e}_n)\)と置いて
\begin{align} \v {E}^t\v {E} &= \begin {pmatrix} \v {e}_1^t\\ \vdots \\ \v {e}^t_n \end {pmatrix} \begin {pmatrix} \v {e}_1\cdots \v {e}_n \end {pmatrix} \nonumber \\ &= \begin {pmatrix} \inner {\v {e}_1}{\v {e}_1}&\cdots &\inner {\v {e}_1}{\v {e}_n}\\ \vdots &\ddots &\vdots \\ \inner {\v {e}_n}{\v {e}_1}&\cdots &\inner {\v {e}_n}{\v {e}_n} \end {pmatrix} \nonumber \\ &= 1 \end{align}
なので(ここでの\(1\)は\(n\times n\)の単位行列を意味する)、
\begin{align} \v {g} = \mathcal {A}^t\mathcal {A} \label {eq:det-g-AA} \end{align}
となることが分かる。ここで\(\v {g}\)は計量\(g_{ij}\)を行列として扱ったものを表す。ベクトル\(\v {v}=v^1\v {u}_1+\cdots +v^n\v {u}_n\)の基底を\(\v {e}_1,\cdots ,\v {e}_n\)で置き換えると
\begin{align} \v {v} &= \sum _{i,j=1,\cdots ,n}\tensorUD {A}{j}{i}v^i\v {e}_j \nonumber \\ &= \sum _{i=1,\cdots ,n}v^{\prime i}\v {e}_i \ \ \ \ \ \ \left (v^{\prime i}=\sum _{j=1,\cdots ,n}\tensorUD {A}{i}{j}v^j\right ) \end{align}
と書ける。従って第二基本形式は
\begin{align} \tensorUD {\sigma }{\alpha }{ij} = \sum _{k,l=1,\cdots ,n} \tensorUD {\sigma }{\prime \alpha }{kl} \tensorUD {A}{k}{i}\tensorUD {A}{l}{j} \end{align}
の関係が分かる。 するとRiemannの曲率テンソルもまた
\begin{align} R_{ijkl} &= \sum _{a,b,c,d=1,\cdots ,n} R^\prime _{abcd} \tensorUD {A}{a}{i} \tensorUD {A}{b}{j} \tensorUD {A}{c}{k} \tensorUD {A}{d}{l} \end{align}
となることが分かる。 ここでまた\(n=2\)の\(2\)次元曲面の場合を考える。その場合にはRiemannの曲率テンソルは(50)の関係式から
\begin{align} R_{ijkl} &= \sum _{a,b,c,d=1,2} \sigma ^\prime \cdot \varepsilon _{ab} \varepsilon _{cd} \tensorUD {A}{a}{i} \tensorUD {A}{b}{j} \tensorUD {A}{c}{k} \tensorUD {A}{d}{l} \nonumber \\ &= \sigma ^\prime \left (\det A\right )^2 \varepsilon _{ij}\varepsilon _{kl} \label {eq:2-riemann-kyokuritsu} \end{align}
さらに(50)と(55)式から
\begin{align} \sigma ^\prime = \frac {\sigma }{\det \v {g}} \end{align}
となる。これは\(N=3\)の時はちょうど幾何のお話2[1]で得られたガウス曲率に対する一般式と一致している。即ち一般の\(N\)次元ユークリッド空間の中の\(2\)次元曲面への拡張となっている。 ここでまた一般的な\(n\)次元の超曲面\(M\)の場合を考える。超曲面\(M\)の点\(p\)の\(2\)つの異なる接ベクトルを\(\v {v}_1=v_1^1\v {u}_1+\cdots +v_1^n\v {u}_n\)と\(\v {v}_2=v_2^1\v {u}_1+\cdots +v_2^n\v {u}_n\)として、これらの張る\(p\)での\(M\)の\(2\)次元接平面を\(\Pi \)とする。これまでと同じ議論を行っう。 \(\v {v}_a\ (a=1,2)\)の微分は
\begin{align} \diff {\v {v}_a}{s} &= \sum _{b=1,2} \tensorDU {w}{a}{b}\v {v}_a + \sum _{i=1,\cdots ,n-2} \tensorDU {w}{a}{i}\v {u}_i + \sum _{\alpha =1,\cdots ,N-n} \tensorDU {\sigma }{a}{\alpha } \v {n}_\alpha \nonumber \\ &= \sum _{b=1,2} \tensorDU {\tilde {w}}{a}{b}\v {v}_b + \sum _{\alpha =1,\cdots ,N-2} \tensorDU {\tilde {\sigma }}{a}{\alpha }\v {n}_\alpha \end{align}
のようになる。ここで右辺の\(1\)段目は\(M\)の接ベクトルの基底を \(\v {v}_1,\v {v}_2,\v {u}_1,\cdots ,\v {u}_{n-2}\)と選んで表示したものである。 \(2\)段目は\(\Pi \)の接ベクトルと法線ベクトルとで分けて表したものである。 こうして並べてかき分けてみると、\(1\)段目の第一項の\(\tensorDU {w}{a}{b}\) と\(2\)段目の第一項の\(\tensorDU {\tilde {w}}{a}{b}\)とは等しいことが分かる。 また\(1\)段目の第三項の\(\tensorDU {\sigma }{a}{\alpha }\)は \(2\)段目の第二項に現れる\(\alpha =1,\cdots ,N-n\)までの \(\tensorDU {\tilde {\sigma }}{a}{\alpha }\)と等しいことが分かる。 従って\(1\)段目の\(\tensorDU {w}{a}{i}\)は \(2\)段目の\(\alpha =N-n+1,\cdots ,N-2\)までの \(\tensorDU {\tilde {\sigma }}{a}{\alpha }\)が等しいことが分かる。 ここで\(\Pi \)を測地座標とする\(2\)次元曲面\(S\)を考えると、 \(S\)のRiemannの曲率テンソル\(\tilde {R}_{abcd}\)を考えることができる。 それは\(\tensorDU {\tilde {\sigma }}{a}{\alpha }\)を用いて (40)のように表すことができた。 一方で超曲面\(M\)のRiemannの曲率テンソル\(R_{ijkl}\)は \(\tensorDU {\sigma }{a}{\alpha }\)を用いて同様に(40) のように表される。ベクトル\(\v {v}_1\)と\(\v {v}_2\)が測地線の方程式を満たすとき、 \(\tensorDU {w}{a}{b}\)と\(\tensorDU {w}{a}{i}\)は全て\(0\)となる。 従って\(\tensorDU {\tilde {\sigma }}{a}{\alpha }\)の\(\alpha =N-n+1,\cdots ,N-2\) は全て\(0\)となる。従ってその場合にはそれぞれのRiemannの曲率テンソルは等しい。 即ち
\begin{align} \tilde {R}_{abcd} = \sum _{i,j,k,l=1,\cdots ,n} R_{ijkl} v^i_a v^j_b v^k_c v^l_d \end{align}
この左辺は(59)式の左辺と等しいものである。 この場合\(S\)の計量は\(g_{ab}=\inner {\v {v}_a}{\v {v}_b}\)なので 従って\(S\)のガウス曲率は
\begin{align} K = \sum _{i,j,k.l=1,\cdots ,n} \frac { R_{ijkl} v^i_1 v^j_2 v^k_1 v^l_2 }{ \inner {\v {v}_1}{\v {v}_1} \inner {\v {v}_2}{\v {v}_2} -\inner {\v {v}_1}{\v {v}_2}^2 } \end{align}
のように表すことができる。この時\(K\)を\(M\)の\(\Pi \)に対する 断面曲率とも呼ぶ。[3]
ここまで\(N\)次元ユークリッド空間\(\mathbb {R}^N\)の中の\(n\)次元の超曲面\(M\)について 考えてきた。ここでは超曲面\(M\)が\(\mathbb {R}^N\)の座標\(\v {x}\)の関数 \(f(\v {x})\)を用いて\(f(\v {x})=0\)で定められている場合を考える。 \(M\)の座標\(\v {u}\)の関数として見た時は\(f(\v {x}(\v {u}))=0\)という式は 恒等式となる。これを\(u^i\)で微分すると
\begin{align} \inner {\pdiff {f}{\v {x}}}{\pdiff {\v {x}}{u^i}}=0 \end{align}
を得る。 この式から\(\pdiff {f}{\v {x}}\)と\(\v {u}_i=\pdiff {\v {x}}{u^i}\)とは直交している。 即ち\(\pdiff {f}{\v {x}}\)は法線ベクトルとなっていることが分かる。 点\(\v {x}\)から出る一次独立なベクトルの組は\(N\)個あるので、 そのうちの\(1\)つを\(\pdiff {f}{\v {x}}\)と平行にとれば、他の\(N-1\)個は \(\pdiff {f}{\v {x}}\)に直交するように取れる。 即ち\(f=0\)を満たす自由度は\(N-1\)個ある。従って\(M\)の座標としては \(\v {u}=(u^1,\cdots ,u^{N-1})\)のように取ることができて、\(M\)は\(N-1\)次元空間をなす。 さらに上記の関係式も\(\v {u}\)に関しては恒等式なのでさらに\(u^j\)で微分しても\(0\)である。 即ち
\begin{align} \sum _{a,b=1,\cdots ,N} \frac {\pd ^2 f}{\pd x^a\pd x^b}\pdiff {x^a}{u^i}\pdiff {x^b}{u^j} + \inner {\pdiff {f}{\v {x}}}{\frac {\pd ^2 \v {x}}{\pd u^i \pd u^j}} = 0 \end{align}
を得る。(24)式から、法線成分は\(1\)つしかないので 単位法線ベクトルを\(\v {n}\)と書くと
\begin{align} \inner {\pdiff {f}{\v {x}}}{\v {n}} \sigma _{ij} = - \inner {\v {u}_i}{\frac {\pd ^2 f}{\pd \v {x}\pd \v {x}}\cdot \v {u}_j} \end{align}
の関係が得られる。従って\(\pdiff {f}{\v {x}}=\left |\pdiff {f}{\v {x}}\right |\v {n}\) に注意すれば、\(\sigma _{ij}\)は\(f\)を用いて
\begin{align} \sigma _{ij} = - \frac {1}{\left |\pdiff {f}{\v {x}}\right |} \inner {\v {u}_i}{\frac {\pd ^2 f}{\pd \v {x}\pd \v {x}}\cdot \v {u}_j} \end{align}
と表すことができる。\(\sigma _{ij}\)を用いてRiemannの曲率テンソルを表すことが できるので\(f\)を用いてRiemannの曲率テンソルを書き表す関係式が得られたことになる。 具体的に書き表すと
\begin{align} R_{ijkl} = \frac {1}{\left |\pdiff {f}{\v {x}}\right |^2} \left ( \inner {\v {u}_i}{\frac {\pd ^2 f}{\pd \v {x}\pd \v {x}}\cdot \v {u}_k} \inner {\v {u}_j}{\frac {\pd ^2 f}{\pd \v {x}\pd \v {x}}\cdot \v {u}_l} - \inner {\v {u}_i}{\frac {\pd ^2 f}{\pd \v {x}\pd \v {x}}\cdot \v {u}_l} \inner {\v {u}_j}{\frac {\pd ^2 f}{\pd \v {x}\pd \v {x}}\cdot \v {u}_k} \right ) \end{align}
となる。これらの関係式を見ると\(\sigma _{ij}\)は\(f\)のヘッシアンに関係しているし、 曲率はその部分\(2\times 2\)行列の行列式に関係している。 例として簡単な具体例で計算してみる。\(n\)次元の半径\(r\)の球面を\(M\)とする。 \(M\)を規定する関数は\(f=\frac {1}{2}(\v {x}^2-r^2)\)である。微分すると \(\d f=\inner {\v {x}}{\d \v {x}}\)なので\(\v {x}_0=(r,0,0,\cdots ,0)\) では\(\pdiff {f}{\v {x}}(\v {x}_0)=(r,0,0,\cdots ,0)\)となる。 従って\(\left |\pdiff {f}{\v {x}}(\v {x}_0)\right |=r\) ヘッシアン\(\frac {\pd ^2 f}{\pd \v {x}\pd \v {x}}\)は単位行列\(1\)に等しい。 従って
\begin{align} R_{ijkl} = \frac {1}{r^2} \left ( \delta _{ik}\delta _{jl} - \delta _{il}\delta _{jk} \right ) \end{align}
となり、よく知られた\(n\)次元球面のRiemannの曲率テンソルが得られた。
\(M\)を定める関数が複数あり、\(N-n\)個の 条件\(f^1(\v {x})=0,\ f^2(\v {x})=0,\cdots ,\ f^{N-n}(\v {x})=0\) によって定められ、 点\(\v {x}\)の近傍で各\(\pdiff {f^\alpha }{\v {x}}\)が全て一次独立である時は、 点\(\v {x}\)での\(M\)の接ベクトルで一次独立なものが 各\(\pdiff {f^\alpha }{\v {x}}\)と直交する\(n\)個が取れるので\(M\)は\(n\) 次元の空間となる。 この場合にも各\(f^\alpha \)に対して上記の計算が従い、
\begin{align} \sum _{\beta =1\cdots ,N-n} \inner {\pdiff {f^\alpha }{\v {x}}}{\v {n}^\beta } \tensorUD {\sigma }{\beta }{ij} = - \inner {\v {u}_i}{\frac {\pd ^2 f^\alpha }{\pd \v {x}\pd \v {x}}\cdot \v {u}_j} \end{align}
となる。従って
\begin{align} \tensorUD {\sigma }{\alpha }{ij} = - \sum _{\beta =1,\cdots ,N-n} \inner {\pdiff {f^\alpha }{\v {x}}}{\v {n}_\beta }^{-1} \inner {\v {u}_i}{\frac {\pd ^2 f^\beta }{\pd \v {x}\pd \v {x}}\cdot \v {u}_j} \end{align}
として得られる。ここで\(\inner {\pdiff {f^\alpha }{\v {x}}}{\v {n}_\beta }^{-1}\) は\(\inner {\pdiff {f^\alpha }{\v {x}}}{\v {n}_\beta }\)の逆行列を表す。 従ってこれからRiemannの曲率テンソルを\(f^\alpha \)を用いて表すことができる。 多くの場合で、超曲面\(M\)をこのように関数を用いて規定することがよくある。 その場合には関数を使って直接Riemannの曲率を計算することができる。
ここで正規座標について少し考察する。 正規座標系は点\(p\)の接ベクトル\(\v {u}\)を初期ベクトルとした測地線方程式(26) 式の解は\(\v {x}_s\ (\diff {\v {x}}{s}=\v {v}_s)\)の\(s=1\)での位置座標を\(\v {u}\)とするものであった。 従って\(\v {x}_1=\v {x}(\v {u})\)である。初期ベクトルを\(\v {u}\)から \(C\v {u}\)に置き直した測地線の方程式の解\(\v {x}^\prime _s\)はもとの測地線の方程式 の解\(\v {x}_s\)の\(s\)を\(Cs\)に置き換えたものに等しい。従って\(\v {x}_s=\v {x}(s\v {u})\) となることが分かる。これから
\begin{align} \diff {\v {x}_s}{s} &= \v {v}_s \nonumber \\ &= \sum _{i=1,\cdots ,n} u^i\v {u}_i \ \ \ \ \left (\v {u}_i=\pdiff {\v {x}}{u^i}(s\v {u})\right ) \end{align}
これをさらに\(s\)で微分すると\(\v {v}_s\)が測地線の方程式を満たすことから \(\diff {\v {x}_s}{s}\)の接線成分が\(0\)になることに注意すれば
\begin{align} \sum _{i,j=1,\cdots .n} \frac {\pd ^2\v {x}}{\pd u^i\pd u^j} u^iu^j &= \sum _{\substack {i,j=1,\cdots ,n\\ \alpha =1,\cdots .N-n}} \tensorUD {\sigma }{\alpha }{ij}u^iu^j\v {n}_\alpha \end{align}
を得る。右辺は\(\v {v}_s\)を\(s\)で微分したときの法線成分になる。 従って左辺は接線成分を含まないことが分かる。 従って正規座標系では
\begin{align} \sum _{j,k=1,\cdots ,n}\tensorUD {\Gamma }{i}{jk}u^ju^k = \sum _{j,k=1,\cdots ,n} \frac {1}{2} \left ( \tensorUD {\Gamma }{i}{jk}+\tensorUD {\Gamma }{i}{kj} \right ) u^ju^k = 0 \end{align}
である。 これは恒等式なのでさらに\(s\)で微分すると
\begin{align} \sum _{j,k,l=1,\cdots ,n} \frac {1}{3} \left ( \pdiff {\tensorUD {\Gamma }{i}{jk}}{u^l} + \pdiff {\tensorUD {\Gamma }{i}{kl}}{u^j} + \pdiff {\tensorUD {\Gamma }{i}{lj}}{u^k} \right ) u^ju^ku^l =0 \end{align}
となる。ここでこれらの\(\tensorUD {\Gamma }{i}{jk}\)は\(s\v {u}\)に依存している ためそのままこれらの和の係数が\(0\)であるとは言えない。一点、\(s=0\)の時のみ\(\v {u}\)に依存しなくなるため 点\(p\)においてのみ成分が\(0\)であるということができる。 従って正規座標系では\(\tensorUD {\Gamma }{i}{jk}\)の下付きの添え字の対称性から
\begin{align} \tensorUD {\Gamma }{i}{jk}(0)=0 \end{align}
および
\begin{align} \pdiff {\tensorUD {\Gamma }{i}{jk}}{u^l}(0) + \pdiff {\tensorUD {\Gamma }{i}{kl}}{u^j}(0) + \pdiff {\tensorUD {\Gamma }{i}{lj}}{u^k}(0) = 0 \label {eq:cyclic-Gamma} \end{align}
となることを意味する。 さてRiemannの曲率テンソルは
\begin{align} \tensorUD {R}{i}{jkl}(0) = \pdiff {\tensorUD {\Gamma }{i}{jl}}{u^k}(0) - \pdiff {\tensorUD {\Gamma }{i}{jk}}{u^l}(0) \end{align}
一方、添え字を入れ替えて
\begin{align} \tensorUD {R}{i}{lkj}(0) = \pdiff {\tensorUD {\Gamma }{i}{lj}}{u^k}(0) - \pdiff {\tensorUD {\Gamma }{i}{lk}}{u^j}(0) \end{align}
これらを足して、さらに(78)式を使うと
\begin{align} \pdiff {\tensorUD {\Gamma }{i}{jl}}{u^k}(0) = \frac {1}{3} \left ( \tensorUD {R}{i}{jkl}(0) + \tensorUD {R}{i}{lkj}(0) \right ) \label {eq:Gamma-R} \end{align}
を得る。 この式を使うと正規座標系での計量をRiemannの曲率テンソルを用いた表示を 得ることができる。計量\(g_{ij}(\v {u})\)を点\(p\)、即ち\(\v {u}=0\)まわりで展開すると
\begin{align} g_{ij}(\v {u}) &= g_{ij}(0) + \sum _{k=1,\cdots ,n} \pdiff {g_{ij}}{u^k}(0)u^k + \frac {1}{2!} \sum _{k,l=1,\cdots ,n} \frac {\pd ^2 g_{ij}}{\pd u^k\pd u^l}(0) u^ku^l +\cdots \end{align}
ここで
\begin{align} \pdiff {g_{ij}}{u^k} &= \pdiffop {u^k}\inner {\v {u}_i}{\v {u}_j} \nonumber \\ &= \inner {\v {u}_i}{\frac {\pd ^2\v {x}}{\pd u^k\pd u^j}} + \inner {\frac {\pd ^2\v {x}}{\pd u^k\pd u^i}}{\v {u}_j} \nonumber \\ &= \Gamma _{ijk}+\Gamma _{jik} \end{align}
さらに
\begin{align} \frac {\pd ^2g_{ij}}{\pd u^k\pd u^l} &= \pdiff {\Gamma _{ijk}}{u^l} + \pdiff {\Gamma _{jik}}{u^l} \end{align}
となる。従って点\(p\)では\(\v {u}=0\)なので(81)式を使うと
\begin{align} \frac {\pd ^2g_{ij}}{\pd u^k\pd u^l}(0) &= \pdiff {\Gamma _{ijk}}{u^l}(0) + \pdiff {\Gamma _{jik}}{u^l}(0) \nonumber \\ &= \frac {1}{3} \left ( R_{iklj}(0) + R_{jlki}(0) \right ) \end{align}
従って
\begin{align} g_{ij}(\v {u}) = g_{ij}(0) - \frac {1}{3} \sum _{k,l=1,\cdots ,n} R_{ikjl}(0)u^ku^l + \cdots \end{align}
を得る。ここでさらに点\(p\)での計量を計算する。
\begin{align} \inner {\v {v}_s}{\v {v}_s} &= \inner {\diff {\v {x}_s}{s}}{\diff {\v {x}_s}{s}} \nonumber \\ &= \sum _{i,j=1,\cdots ,n} \inner {\pdiff {\v {x}_s}{u^i}}{\pdiff {\v {x}_s}{u^j}}u^iu^j \nonumber \\ &= \sum _{i,j=1,\cdots ,n} g_{ij}(s\v {u})u^iu^j \end{align}
となるが、\(\lim _{s=0}\inner {\v {v}_s}{\v {v}_s}=\v {u}^2\)であるので \(g_{ij}(0)=\delta _{ij}\)となる。 ここで\(\delta _{ij}\)は\(i=j\)のとき\(1\)でそれ以外で\(0\)となる。 従ってまとめると 正規座標系では、計量は
\begin{align} g_{ij}(\v {u}) = \delta _{ij} - \frac {1}{3} \sum _{k,l=1,\cdots ,n} R_{ikjl}(0)u^ku^l + \cdots \label {eq:keiryou-seikizahyou} \end{align}
と書けることが分かった。 従って正規座標系で表した時には計量が、従って曲線の長さが ユークリッド空間の時からどれだけずれたかをRiemannの曲率テンソル によって測ることができる。 さらに、体積要素は幾何のお話ノート[4]によれば 計量の行列式のルート\(\sqrt {g}\)で表された。 従って(88)の行列式を取れば
\begin{align} g = \det g_{ij} = 1 - \frac {1}{3} \sum _{i,k,l=1,\cdots ,n} \tensorUD {R}{i}{kil}(0)u^ku^l + \cdots \end{align}
となるため
\begin{align} \sqrt {g} = 1 - \frac {1}{6} \sum _{i,j=1,\cdots ,n} R_{ij}(0)u^iu^j + \cdots \label {eq:taisekiyouso-seikizahyou} \end{align}
となる。ここで
\begin{align} R_{ij} = \sum _{k=1,\cdots ,n} \tensorUD {R}{k}{ikj} \end{align}
をRicciテンソルと呼ぶ。 添え字の上げ下げは計量\(g_{ij}\)やその逆行列\(g^{ij}\)で行う。 ここでは高次のオーダーは無視しているので実質的に\(\delta _{ij}\) で上げ下げを行ってよい。 即ち(90)式から 正規座標系で表した時に、体積がユークリッド空間の時の体積から どれくらいずれているかをRicciテンソルで計ることができる。
この正規座標系での\(|\v {u}|<r\)の領域\(D\)の体積を計算してみる。
\begin{align} V &= \int _D\d ^n\v {x} \nonumber \\ &= \int _D\sqrt {g}\cdot \d ^n\v {u} \nonumber \\ &= \int _D\d ^n\v {u} - \frac {1}{6}R_{ij}(0) \int _Du^iu^j\d ^n\v {u} + \cdots \nonumber \\ &= \frac {2\pi ^{\frac {n}{2}}}{n\Gamma (\frac {n}{2})}r^n - \frac {1}{6}R_{ij}(0) \int _Du^iu^j\d ^n\v {u} + \cdots \end{align}
ここで最後の段の第二項の計算は\(i\neq j\)の時には被積分関数が奇関数になることと、 積分領域が\(\v {u}=0\)を中心に対称であることから\(0\)となる。 残るのは\(i=j\)の積分であり、対称性から添え字の\(i\)に依存しない。 従って
\begin{align} \int _D(u^i)^2\d ^n\v {u} &= \frac {1}{n}\int _D\v {u}^2\d ^n\v {u} \nonumber \\ &= \frac {2\pi ^{\frac {n}{2}}}{n\Gamma (\frac {n}{2})} \int _0^rr^{n+1}\d r \nonumber \\ &= \frac {2\pi ^{\frac {n}{2}}}{n(n+2)\Gamma (\frac {n}{2})} r^{n+2} \end{align}
となる。従って
\begin{align} V = \frac {2\pi ^{\frac {n}{2}}}{n\Gamma (\frac {n}{2})}r^n \left ( 1 - \frac {R(0)}{6(n+2)} r^2 + \cdots \right ) \end{align}
となる。ここで
\begin{align} R=\sum _{i=1,\cdots ,n}\tensorUD {R}{i}{i} \end{align}
をスカラー曲率という。 即ち正規座標系で点\(p\)からの距離が一定の球体の体積が ユークリッド空間の場合からどれくらいずれるかをスカラー曲率で計ることが できる。
曲がった空間を考察する
ここまではユークリッド空間の中にある超曲面を考えてきたので、自然に超曲面の 曲がりを計る量として接続や曲率といったものが自然に誘導された。 ここでは空間が曲がっていることを空間の外のユークリッド空間を持ち出すことなく 計ることができるか?という問いに関して考察したい。 最初に考察したいのは空間\(M\)に\(2\)つの点が与えられた時に、それらの間の最短ルートと そのルートの長さ(=距離)である。 これを数学的に扱う手法は確立されている。 即ち点\(p\)から点\(q\)までを結ぶ可能な曲線\(C\)の長さ
\begin{align} L = \int _0^t\sqrt { \left (\diff {x_1}{t}\right )^2 + \left (\diff {x_2}{t}\right )^2 + \cdots + \left (\diff {x_N}{t}\right )^2 }\cdot \d t \end{align}
を最小にする曲線\(C\)を求める問題である。 しかしこれは\(\mathbb {R}^N\)に埋め込まれた超曲面であればこれで良いが、今回は そのような描像を採用しない。従って空間\(M\)の外の情報を使えない。 ここまでの議論から空間\(M\)に計量\(g_{ij}\)が与えられていれば曲線\(C\)の長さは
\begin{align} L = \int _0^t \sqrt { \sum _{i,j=1,\cdots ,n}g_{ij}(\v {u}_t) \diff {u^i}{t}\diff {u^j}{t} }\cdot \d t \end{align}
で与えられることが理解できる。 端点の\(p\)と\(q\)を固定してその間を結ぶ曲線を動かした時にその曲線の長さが最小になる問題は 曲線を動かした時に\(L\)がどう変わるかを計ればよい。 曲線は\(\v {u}_t\)のように\(t\)を与えればそれに対して\(M\)の中の位置\(\v {u}\)が決まるという 形で与えられるので、曲線を動かすというのはそれに対して \(\v {u}_t\rightarrow \v {u}_t+\delta \v {u}_t\)のような変動 が加わることを意味する。従って\(u^i\)に関して\(L\)の被積分関数の微分を行えば、 それが\(L\)の変分となる。曲線の長さが最小になるというのはその変分に対して\(L\)の変化が 停留するというよく使われる幾何学の手法をとることで計算することができる。 このことは
\begin{align} \mathcal {L} = \sqrt { \sum _{i,j=1,\cdots ,n}g_{ij}(\v {u}_t) \diff {u^i}{t}\diff {u^j}{t} } \end{align}
に対して \(\delta u^i\)の変分\(\delta \mathcal {L}\) を計算すればよい。実際に計算すると
\begin{align} \delta \mathcal {L} &= \sum _{i,j=1,\cdots ,n} \frac {1}{2\mathcal {L}} \left ( 2g_{ij}\diff {\delta u^i}{t}\diff {u^j}{t} + \sum _{k=1,\cdots ,n} \pdiff {g_{jk}}{u^i}\delta u^i \diff {u^j}{t}\diff {u^k}{t} \right ) \nonumber \\ &= - \frac {1}{\mathcal {L}} \sum _{i,j=1,\cdots ,n} \left [ \left ( g_{ij}\diffn {2}{u^j}{t} + \sum _{k=1,\cdots ,n} \left ( \pdiff {g_{ij}}{u^k}\diff {u^j}{t}\diff {u^k}{t} - \frac {1}{2} \pdiff {g_{jk}}{u^i}\diff {u^j}{t}\diff {u^k}{t} \right ) \right )\delta u^i - \diffop {t} \left ( g_{ij}\diff {u^j}{t}\delta u^i \right ) \right ] \nonumber \\ &= - \frac {1}{\mathcal {L}} \sum _{i,j=1,\cdots ,n} g_{ij} \left [ \left ( \diffn {2}{u^j}{t} + \sum _{k,l,m=1,\cdots ,n} \frac {1}{2}g^{jk} \left ( \pdiff {g_{kl}}{u^m} + \pdiff {g_{km}}{u^l} - \frac {1}{2} \pdiff {g_{lm}}{u^k} \right ) \diff {u^l}{t}\diff {u^m}{t} \right )\delta u^i - \diffop {t} \left ( g_{ij}\diff {u^j}{t}\delta u^i \right ) \right ] \nonumber \\ &= - \frac {1}{\mathcal {L}} \sum _{i,j=1,\cdots ,n} g_{ij} \left ( \diffn {2}{u^j}{t} + \sum _{l,m=1,\cdots ,n} \tensorUD {\Gamma }{j}{lm} \diff {u^l}{t}\diff {u^m}{t} \right ) \delta u^i + \sum _{i,j=1,\cdots ,n} \left ( \diffop {t} \left ( \frac {g_{ij}\diff {u^j}{t}\delta u^i}{\mathcal {L}} \right ) - \frac {g_{ij}\diff {u^j}{t}\delta u^i}{\mathcal {L}^2} \diff {\mathcal {L}}{t} \right ) \end{align}
ここで最後の等式の最初の項は測地線の方程式の形となっている。 第二項は\(t\)に関する全微分の形である。 第三項は\(\mathcal {L}\)が曲線の長さの\(t\)での微分\(\diff {s}{t}\)なので パラメータ\(t\)を曲線の長さそのもの\(s\)にとっておけば \(\diff {s}{t}=\mathcal {L}=1\) となり定数となるためその場合には \(t\)を\(s\)で取り直して\(\diff {\mathcal {L}}{s}=0\)となる。 従って\(s\)で積分すれば
\begin{align} \delta L &= -\int _0^L \left [ \inner {\nabla _\v {v} \v {v}}{\delta \v {u}_s} - \sum _{i,j=1,\cdots ,n} \diffop {s} \left ( g_{ij}\diff {u^j}{s}\delta u^i \right ) \right ]\d s \end{align}
となる。ここで変分は点\(p\)と点\(q\)を固定して行うため、 \(s=0\)と\(s=L\)での\(\v {u}_s\)の変分\(\delta \v {u}_s\)は\(0\)である。 従ってこの積分の第二項は\(0\)となる。結果として第一項のみが 残り、この変分が停留するという条件は\(\nabla _{\v {v}}\v {v}=0\) という測地線の方程式となる。 従って空間\(M\)の\(2\)点間の最短ルートは測地線である ことが分かる。
空間の中の点\(p\)を基準にしたとき、\(p\)のまわりの座標の取り方 は任意であるが、通常のユークリッド空間と同じような直交座標系 をとろうと思うのは自然なことである。 点\(p\)を始点とした方向をひとつ決めると、その方向へまっすぐに直線を 伸ばしたものを\(x\)軸とみなすことになる。 この時、この空間\(M\)の中での"まっすぐ"というのは測地線をを指すというのは とても自然である。このようにして点\(p\)を始点とする、互いに直交する \(n\)個の方向ベクトルを選び、それらを初期値とした測地線を引くことで 座標系を作ることができる。 空間\(M\)の中での直線とは測地線を意味することになる。 このようにしてできる座標系は正規座標系である。 正規座標系では点\(p\)のまわりではChristoffelの記号は\(0\)となるため、 点\(p\)の近傍では測地線の方程式は近似的に単なる直線の方程式である。 その直線からのずれは曲率で与えられる。
ユークリッド空間での幾何学の基本的な性質で、\(2\)本の平行な直線は どこまでいっても互いに交わらないというのがある。 このことは次のように考察できる。 点\(p\)を始点とするベクトル\(\v {v}_1\)を点\(q\)まで平行移動させたものを \(\v {v}_2\)とする。それら\(2\)つのベクトルを初期値とする測地線の方程式 から\(2\)つの測地線を得る。ユークリッド空間ではこれらは交わることが なかったが、今考えているような空間の幾何学ではこれは一般に成り立たない。 その平行からのずれもまた曲率テンソルで与えられる。 始点\(p\)のベクトル\(\v {v}\)を終点\(r\)まで測地線に沿って平行移動して さらに点\(s\)まで平行移動する。 一方、点\(p\)から点\(q\)にベクトル\(\v {v}\)を平行移動して、 \(q\)から終点\(s\)まで測地線に沿って平行移動したベクトル とのずれは、(49)式より
\begin{align} \Delta v^i = -\tensorUD {R}{i}{jkl}v^j\delta u^k\Delta u^l \end{align}
がそのずれである。 従ってこのずれがあるために平行な\(2\)直線の間はずっと平行とならずに いずれ交わったり逆に離れていったりする。 Christoffelの記号やRiemannの曲率は計量が与えられれば 決定することができる。この意味で空間の中での直線(=測地線)や 平行性、長さ、体積といったものがユークリッド空間からどれくらい ずれているかが決定できる。 ここまでの議論を見返してみると次のことが言える。 空間における計量のユークリッド計量からのずれは空間の曲がりを表す。 ある場所にある物差しの長さが別の場所に持ってきた時に同じ物差しなので 同じ目盛りを刻んでいるが、実際の長さが場所によって変わっていたり あるいは実は曲がっていたりといったことを検知する方法があるということを 示唆している。 アインシュタインの一般相対性理論はこれらの事情をちょうど含んだ理論となっている。 相対性理論はこの空間の曲がりともうひとつ特殊相対性理論からくる 時間と空間との関係まで含めた物理の幾何学理論である。 現実の物理世界では理想的な直線を用意することが難しいが、 相対性理論では光の進む道が直線すなわち測地線の役割を果たしている。 このノートは幾何学の数学的な理解や直観を得ることを目的としているため 相対性理論の話題にこれ以上踏み込むことはしないでおきたい。 興味のある方には[5]などをおすすめしたい。
参考文献