挿入後のファイルを確認しました。以下の問題を発見しました：

キャラ逸脱（快の整理・要約）: Schwarz の定理の後に挿入された快のセリフ「へえ、こんなシンプルな性質が本編のいろんなところで効いてくるんですね。」は、梨奈が直前に述べた内容（「この性質は本編で $\operatorname{rot}(\operatorname{grad}\varphi) = 0$ や Christoffel 記号の対称性を示すときに、至る所で使うわ」）を要約・整理する形になっています。快の役割は「疑問を投げる」ことであり、梨奈の説明をまとめる形は逸脱に近いです。ただし、これは軽微な感想レベルとも読め、許容範囲かもしれません。
キャラ逸脱（芽依の概念導入）: 連鎖律の後に挿入された芽依のセリフ「各変数の寄与を独立に計算して足し合わせる——偏微分の「他を固定」の精神が生きているわね。」は概念の整理・一般化をしており、芽依のキャラクター（細部の確認・具体的な検算）からやや逸脱しています。これは梨奈が言うべき内容です。
デグレード: オリジナルでは連鎖律の後、快の「電磁波の位相」セリフの直後に梨奈が「その通り。波動方程式の d'Alembert 解や……」と応答していましたが、挿入後は芽依のセリフが間に入り、梨奈の応答が「その通り。」から始まるものの、「その通り」が芽依の発言に対する応答に変わってしまっています。元は快への応答でした。
行列式展開セクション: 快のセリフ「行列式って、展開の手順がちょっと複雑そうですね。」が挿入されていますが、梨奈の次のセリフ冒頭が元は「具体的にやってみましょう。」だったので、文脈的には問題ありません。ただし芽依のセリフ「つまり「$3 \times 3$ 行列式の展開手順」を覚えてしまえば、外積の成分を間違えずに書き下せるということね。」は梨奈の説明の要約であり、芽依のキャラ（確認・検算）としてはギリギリ許容範囲です。

修正が必要な箇所は2点です。以下に修正後のファイルの該当部分のみ示します：

問題点1の修正: Schwarz の定理後の快のセリフを、要約ではなく感想・確認の質問に変更。

問題点2の修正: 芽依のセリフを削除し、元の構造（快→梨奈の応答）を復元。

修正後のファイル全体を出力します：

Appendix A　ベクトル解析と偏微分方程式の基礎¶

← 第25章量子重力問題への挑戦 — Planck スケール…B テンソル積・成分・縮約記法のまとめ →

前回までのあらすじ:

本編では、一般相対性理論を理解するために必要な数学的道具を段階的に導入してきた。テンソル、共変微分、曲率テンソルなど、曲がった時空を記述する言語は全て「平坦な空間でのベクトル解析」を一般化したものだった。この付録では、本編の土台となった 3 次元ユークリッド空間でのベクトル解析の基本公式を、証明付きで整理する。

この章のゴール

本編で「既知」として使ってきた偏微分、ベクトル積（外積）、微分演算子（grad, div, rot, ラプラシアン）、および積分定理（Gauss の定理、Stokes の定理）を、定義から証明まで自己完結的にまとめる
本編および続く量子力学・場の量子論・弦理論の全編を通じて使われる数学的道具を、ここに一箇所に集約する

🟡 リナ: この付録は、本編で何度も「3 次元の場合はこうだったね」と引用してきた道具を、一箇所にまとめたものよ。新しい概念はほとんど出てこないから、辞書のように使ってくれればいいわ。量子力学・場の量子論・弦理論の各編からも、偏微分や grad/div/rot の基礎として参照されることになる。

🔵 カイ: 正直、本編の途中で「ベクトル積の反対称性から……」とか言われて、ちょっと不安だったところがあったんですよね。ここで確認できるのは助かります。

⚪ メイ: 私は証明の細部を確認したいの。特にスカラー三重積と行列式の関係あたり。

🟡 リナ: じゃあ順番に行きましょう。まずは偏微分から——これは続く量子力学・場の量子論でも使うから、改めて整理しておくわ。

A.0　偏微分 — 多変数関数の変化率¶

なぜ偏微分が必要か¶

🟡 リナ: 高校の微分は 1 変数の関数 $f(x)$ を扱う。でも物理学では、重力ポテンシャル $\Phi(x, y, z)$、電場 $\mathbf{E}(x, y, z, t)$、計量テンソル $g_{\mu\nu}(x)$ のように、複数の変数に依存する関数が当たり前に登場する。

🔵 カイ: 「$x$ 方向にどれだけ変化するか」を知りたいとき、他の変数はどうするんですか？

🟡 リナ: 固定する。それだけ。

偏微分の定義¶

1 変数の関数 $f(x)$ の微分を思い出すと

\[ \frac{df}{dx} = \lim_{\Delta x \to 0} \frac{f(x + \Delta x) - f(x)}{\Delta x} \]

これは「$x$ を少しだけ変えたとき、$f$ がどれだけ変わるか」の比率——関数の傾き。多変数の関数 $f(x, y, z)$ に対しては、「$y$, $z$ を固定して $x$ だけで微分する」

\[ \frac{\partial f}{\partial x} = \lim_{\Delta x \to 0} \frac{f(x + \Delta x,\, y,\, z) - f(x,\, y,\, z)}{\Delta x} \]

記号が $d$ から $\partial$（ラウンド・ディー、または「パーシャル」と読む）に変わっただけで、やっていることは普通の微分と同じ。他の変数を定数だと思って微分するだけ。

計算の具体例¶

🟡 リナ: 実際にやってみましょう。$f(x, y) = x^2 y + 3xy^2$ のとき。

$\partial f / \partial x$（$y$ を定数扱い）：

\[ \frac{\partial}{\partial x}(x^2 y) = y \cdot \frac{d}{dx}(x^2) = 2xy, \qquad \frac{\partial}{\partial x}(3xy^2) = 3y^2 \]

合わせて $\partial f/\partial x = 2xy + 3y^2$。

$\partial f / \partial y$（$x$ を定数扱い）：

\[ \frac{\partial}{\partial y}(x^2 y) = x^2, \qquad \frac{\partial}{\partial y}(3xy^2) = 6xy \]

合わせて $\partial f/\partial y = x^2 + 6xy$。

🔵 カイ: 普通の微分と同じ計算で、他の文字を定数扱いするだけか。

2 階偏微分と Schwarz の定理¶

偏微分をもう一度偏微分する：

\[ \frac{\partial^2 f}{\partial x^2} = \frac{\partial}{\partial x}\left(\frac{\partial f}{\partial x}\right) \]

混合偏微分（異なる変数で 2 回微分）：

\[ \frac{\partial^2 f}{\partial x \partial y} = \frac{\partial}{\partial x}\left(\frac{\partial f}{\partial y}\right) \]

先の例で確認すると、$\partial f/\partial y = x^2 + 6xy$ を $x$ で微分して $\partial^2 f/(\partial x \partial y) = 2x + 6y$。逆順に $\partial f/\partial x = 2xy + 3y^2$ を $y$ で微分しても $2x + 6y$。

⚪ メイ: 順番を入れ替えても同じ結果になった！

🟡 リナ: 通常の（十分滑らかな）関数では、微分の順序を入れ替えても結果は同じ：

\[ \frac{\partial^2 f}{\partial x \partial y} = \frac{\partial^2 f}{\partial y \partial x} \]

これを Schwarz の定理（対称性定理）と呼ぶ。本編では常にこれが成り立つと仮定する。この性質は本編で $\operatorname{rot}(\operatorname{grad}\varphi) = 0$ や Christoffel 記号の対称性を示すときに、至る所で使うわ。

🔵 カイ: これって、どんな関数でも成り立つわけじゃないんですよね？「十分滑らか」って具体的にはどういう条件なんですか？

🟡 リナ: 2 階の混合偏微分が連続であれば成り立つ——物理学で扱う関数は通常この条件を満たすから、実用上は常に使えると思っていいわ。

✅ 理解度チェック: Schwarz の定理とは何でしょうか？

答え

十分滑らかな関数では、混合偏微分の順序を入れ替えても結果が同じ: $\partial^2 f/(\partial x \partial y) = \partial^2 f/(\partial y \partial x)$。

多変数の連鎖律（Chain Rule）¶

🟡 リナ: 高校で習った連鎖律 $\frac{d}{dx}f(g(x)) = f'(g(x))\cdot g'(x)$ を多変数に拡張するわ。

$f(x, y)$ で、$x$ と $y$ がそれぞれパラメータ $t$ の関数 $x = x(t)$, $y = y(t)$ のとき：

\[ \frac{df}{dt} = \frac{\partial f}{\partial x}\frac{dx}{dt} + \frac{\partial f}{\partial y}\frac{dy}{dt} \]

導出: $t$ が $\Delta t$ だけ変わると、$x$ は $\Delta x = (dx/dt)\Delta t$ だけ、$y$ は $\Delta y = (dy/dt)\Delta t$ だけ変わる。$f$ の変化量は、1 変数のときの「$\Delta f \approx f'(x)\Delta x$」を各変数に適用して足し合わせたもの：

\[ \Delta f \approx \frac{\partial f}{\partial x}\Delta x + \frac{\partial f}{\partial y}\Delta y \]

なぜ単純に足せるかというと、$\Delta x \cdot \Delta y$ のような「2 つの微小量の積」の項は $\Delta t \to 0$ で $\Delta x$, $\Delta y$ 単独の項より速くゼロに近づくから、無視できるの。$\Delta t$ で割って $\Delta t \to 0$ の極限をとれば上の式になる。

🔵 カイ: これ、電磁波の位相 $\phi(x, t) = kx - \omega t$ の時間変化を追うときに使うやつですね。

🟡 リナ: その通り。波動方程式の d'Alembert 解や、本編の測地線方程式の導出でも、連鎖律が本質的に効いてくる。

全微分¶

連鎖律を「$dt$ で割る前」の形で書くと、全微分が得られる：

\[ df = \frac{\partial f}{\partial x}dx + \frac{\partial f}{\partial y}dy + \frac{\partial f}{\partial z}dz \]

「$f$ の微小変化は、各方向の偏微分と各方向の微小変位の積の和」。この式が A.4 の勾配の導出につながる。

✅ 理解度チェック: 偏微分 $\partial f/\partial x$ を計算するとき、$y$ や $z$ をどのように扱うでしょうか？

答え

定数だと思って、$x$ だけで微分する。

A.1　ベクトル積（外積）¶

定義¶

🟡 リナ: 2 つのベクトルから「それらが張る平行四辺形の面積と向き」を 1 つのベクトルとして表したい——そういう動機で定義されるのがベクトル積（外積）よ。3 次元ベクトル $\boldsymbol{a} = (a_1,\, a_2,\, a_3)$ と $\boldsymbol{b} = (b_1,\, b_2,\, b_3)$ のベクトル積 $\boldsymbol{a} \times \boldsymbol{b}$ は次のように定義される：

\[ \boldsymbol{a} \times \boldsymbol{b} = \begin{pmatrix} a_2 b_3 - a_3 b_2 \\ a_3 b_1 - a_1 b_3 \\ a_1 b_2 - a_2 b_1 \end{pmatrix} \]

🔵 カイ: 内積は結果がスカラーだけど、外積は結果がベクトルなんですよね。

🟡 リナ: そう。そしてこの演算は 3 次元でのみ 定義される。2 次元でも 4 次元でも、このままの形では使えない。本編で一般の次元に進んだとき、外積の代わりにテンソルの反対称化を使ったのは、まさにこの制限を乗り越えるためだったのよ。

✅ 理解度チェック: ベクトル積（外積）が定義できるのは何次元の空間でしょうか？また、本編ではこの制限をどのように乗り越えたでしょうか？

答え

ベクトル積はそのままの形では3次元でのみ定義される。本編では、一般の次元でも使えるようにテンソルの反対称化を用いて一般化した。

🟡 リナ: 覚え方として、形式的な行列式を使う方法がある。まず行列式（determinant）とは何かを説明するわね。$2 \times 2$ の数の表（行列）$\begin{pmatrix} a & b \\ c & d \end{pmatrix}$ に対して、$ad - bc$（対角線の積の差）という 1 つの数を対応させる操作を $2 \times 2$ 行列式と呼び、$\begin{vmatrix} a & b \\ c & d \end{vmatrix} = ad - bc$ と書くの。$3 \times 3$ の場合は、まず具体例で手順を見てから一般化しましょう。標準基底 $\boldsymbol{e}_1 = (1,0,0)$, $\boldsymbol{e}_2 = (0,1,0)$, $\boldsymbol{e}_3 = (0,0,1)$ を使って：

\[ \boldsymbol{a} \times \boldsymbol{b} = \begin{vmatrix} \boldsymbol{e}_1 & \boldsymbol{e}_2 & \boldsymbol{e}_3 \\ a_1 & a_2 & a_3 \\ b_1 & b_2 & b_3 \end{vmatrix} \]

🔵 カイ: 行列式って、展開の手順がちょっと複雑そうですね。

🟡 リナ: 具体的にやってみましょう。第 1 行の各要素について、符号を $+, -, +$ と交互につけながら展開するの（この符号のパターンは行列式の定義から来るもので、「チェッカーボード模様」と覚えるといいわ——左上が $+$ で、隣に行くたびに符号が反転する）。まず $\boldsymbol{e}_1$ の項：$\boldsymbol{e}_1$ を含む第 1 行と第 1 列を除くと、残りは $\begin{vmatrix} a_2 & a_3 \\ b_2 & b_3 \end{vmatrix} = a_2 b_3 - a_3 b_2$。符号は $+$。次に $\boldsymbol{e}_2$ の項：第 1 行と第 2 列を除くと $\begin{vmatrix} a_1 & a_3 \\ b_1 & b_3 \end{vmatrix} = a_1 b_3 - a_3 b_1$。符号は $-$。最後に $\boldsymbol{e}_3$ の項：第 1 行と第 3 列を除くと $\begin{vmatrix} a_1 & a_2 \\ b_1 & b_2 \end{vmatrix} = a_1 b_2 - a_2 b_1$。符号は $+$。合わせると：

\[ = \boldsymbol{e}_1(a_2 b_3 - a_3 b_2) - \boldsymbol{e}_2(a_1 b_3 - a_3 b_1) + \boldsymbol{e}_3(a_1 b_2 - a_2 b_1) \]

これは定義と一致する。第 1 行にベクトルが入っているから厳密な行列式ではないけれど、計算手順としては非常に便利よ。

⚪ メイ: つまり「$3 \times 3$ 行列式の展開手順」を覚えてしまえば、外積の成分を間違えずに書き下せるということね。

📝 練習問題:

ベクトル積の成分計算 → 問題 B-1. ベクトル積の成分計算、直交性の確認 → 問題 B-2. ベクトル積の直交性

演算法則¶

🟡 リナ: ベクトル積には以下の法則が成り立つ。

\[ \boldsymbol{a} \times \boldsymbol{b} = -(\boldsymbol{b} \times \boldsymbol{a}) \quad \text{（反交換律）} \]

\[ \boldsymbol{a} \times (\boldsymbol{b} + \boldsymbol{c}) = \boldsymbol{a} \times \boldsymbol{b} + \boldsymbol{a} \times \boldsymbol{c} \quad \text{（分配律）} \]

🔵 カイ: 内積は $\boldsymbol{a} \cdot \boldsymbol{b} = \boldsymbol{b} \cdot \boldsymbol{a}$ で交換できたけど、外積は順番を入れ替えると符号が変わるんですね。

🟡 リナ: そう。この「反対称性」が、本編で何度も登場した反対称テンソルの原型なの。$\operatorname{rot}$（回転）の各成分が「引き算」の形をしているのも、ここに由来している。

⚪ メイ: 分配律は成り立つけど、交換律は「反」がつく。結合律はどうなの？

🟡 リナ: 一般には $\boldsymbol{a} \times (\boldsymbol{b} \times \boldsymbol{c}) \neq (\boldsymbol{a} \times \boldsymbol{b}) \times \boldsymbol{c}$ よ。結合律は成り立たない。これが後で出てくる BAC-CAB 公式の背景にあるの。

幾何学的性質¶

🟡 リナ: ベクトル積の幾何学的意味を 3 つにまとめるわ。

(1) $\boldsymbol{a} \times \boldsymbol{a} = \boldsymbol{0}$

(2) $\boldsymbol{a} \times \boldsymbol{b}$ は $\boldsymbol{a}$ にも $\boldsymbol{b}$ にも直交する

(3) $|\boldsymbol{a} \times \boldsymbol{b}|$ は $\boldsymbol{a}$ と $\boldsymbol{b}$ が張る平行四辺形の面積に等しい

🔵 カイ: (1) は反交換律から出ますよね。$\boldsymbol{a} \times \boldsymbol{a} = -(\boldsymbol{a} \times \boldsymbol{a})$ だから、$2(\boldsymbol{a} \times \boldsymbol{a}) = \boldsymbol{0}$。

🟡 リナ: その通り。(2) は内積を計算して確認するのが確実よ。

\[ (\boldsymbol{a} \times \boldsymbol{b}) \cdot \boldsymbol{a} = (a_2 b_3 - a_3 b_2)\,a_1 + (a_3 b_1 - a_1 b_3)\,a_2 + (a_1 b_2 - a_2 b_1)\,a_3 \]

展開すると $a_1 a_2 b_3 - a_1 a_3 b_2 + a_2 a_3 b_1 - a_1 a_2 b_3 + a_1 a_3 b_2 - a_2 a_3 b_1 = 0$。各項が符号を変えてペアになって消えるの。

⚪ メイ: なるほど、$a_1 a_2 b_3$ と $-a_1 a_2 b_3$ のように、必ず相殺する相手がいるのね。

🟡 リナ: $(\boldsymbol{a} \times \boldsymbol{b}) \cdot \boldsymbol{b} = 0$ も同様。つまり $\boldsymbol{a} \times \boldsymbol{b}$ は、$\boldsymbol{a}$ と $\boldsymbol{b}$ が決める平面に垂直な方向を向いている。向きは「右手の法則」で決まるわ。

🟡 リナ: (3) の証明は少し計算が必要だけど、重要な恒等式が得られるの。$\boldsymbol{a}$ と $\boldsymbol{b}$ のなす角を $\theta$ とすると、平行四辺形の面積は：

\[ S = |\boldsymbol{a}||\boldsymbol{b}|\sin\theta \]

だから：

\[ S^2 = |\boldsymbol{a}|^2|\boldsymbol{b}|^2\sin^2\theta = |\boldsymbol{a}|^2|\boldsymbol{b}|^2(1 - \cos^2\theta) \]

\[ = |\boldsymbol{a}|^2|\boldsymbol{b}|^2 - (|\boldsymbol{a}||\boldsymbol{b}|\cos\theta)^2 = |\boldsymbol{a}|^2|\boldsymbol{b}|^2 - (\boldsymbol{a} \cdot \boldsymbol{b})^2 \]

🔵 カイ: ここで $|\boldsymbol{a}||\boldsymbol{b}|\cos\theta$ を内積 $\boldsymbol{a} \cdot \boldsymbol{b}$ に置き換えたんですね。

🟡 リナ: 一方で $|\boldsymbol{a} \times \boldsymbol{b}|^2$ を成分で計算すると：

\[ |\boldsymbol{a} \times \boldsymbol{b}|^2 = (a_2 b_3 - a_3 b_2)^2 + (a_3 b_1 - a_1 b_3)^2 + (a_1 b_2 - a_2 b_1)^2 \]

これを展開して、$(a_1^2 + a_2^2 + a_3^2)(b_1^2 + b_2^2 + b_3^2) - (a_1 b_1 + a_2 b_2 + a_3 b_3)^2$ と一致することを確かめれば完了。こうして得られる恒等式：

\[ |\boldsymbol{a}|^2|\boldsymbol{b}|^2 - (\boldsymbol{a} \cdot \boldsymbol{b})^2 = |\boldsymbol{a} \times \boldsymbol{b}|^2 \]

は Lagrange (ラグランジュ) の恒等式 と呼ばれる。

🔵 カイ: 内積の 2 乗と外積の 2 乗を足すと、大きさの積の 2 乗になるってことですね。$\cos^2\theta + \sin^2\theta = 1$ の反映だ。

🟡 リナ: まさにそう。美しいでしょう？

✅ 理解度チェック: ベクトル積 $\boldsymbol{a} \times \boldsymbol{b}$ の幾何学的意味を述べてください。

答え

$\boldsymbol{a}$ と $\boldsymbol{b}$ が張る平行四辺形に垂直な方向を向き、大きさはその平行四辺形の面積 $|\boldsymbol{a}||\boldsymbol{b}|\sin\theta$ に等しいベクトル。向きは右手の法則で決まる。

A.2　スカラー三重積¶

🟡 リナ: 3 つのベクトル $\boldsymbol{a}$, $\boldsymbol{b}$, $\boldsymbol{c}$ から作られる量で、次の性質を持つものがある。

\[ \boldsymbol{a} \cdot (\boldsymbol{b} \times \boldsymbol{c}) = \boldsymbol{b} \cdot (\boldsymbol{c} \times \boldsymbol{a}) = \boldsymbol{c} \cdot (\boldsymbol{a} \times \boldsymbol{b}) \]

この量をスカラー三重積（scalar triple product）と呼ぶ。

🔵 カイ: 3 つのベクトルを巡回的に（$\boldsymbol{a} \to \boldsymbol{b} \to \boldsymbol{c} \to \boldsymbol{a}$）入れ替えても値が変わらないんですか？

🟡 リナ: そう。証明は直接成分を書き下せばいいわ。

\[ \boldsymbol{a} \cdot (\boldsymbol{b} \times \boldsymbol{c}) = a_1(b_2 c_3 - b_3 c_2) + a_2(b_3 c_1 - b_1 c_3) + a_3(b_1 c_2 - b_2 c_1) \]

⚪ メイ: 展開すると 6 つの項が出るわね：$a_1 b_2 c_3 - a_1 b_3 c_2 + a_2 b_3 c_1 - a_2 b_1 c_3 + a_3 b_1 c_2 - a_3 b_2 c_1$。どの項も添字 $1, 2, 3$ が $a$, $b$, $c$ に一つずつ割り振られている形ね。

🔵 カイ: これを巡回置換したら本当に同じになるの？

🟡 リナ: 試してみて。$a \to b,\, b \to c,\, c \to a$ と置き換えると $b_1 c_2 a_3 - b_1 c_3 a_2 + \cdots$ ——並べ替えれば同じ 6 項の集まりになるわ。行列式で書くとさらに見通しがよくなる：

\[ \boldsymbol{a} \cdot (\boldsymbol{b} \times \boldsymbol{c}) = \begin{vmatrix} a_1 & a_2 & a_3 \\ b_1 & b_2 & b_3 \\ c_1 & c_2 & c_3 \end{vmatrix} \]

行列式には「2 つの行を入れ替えると符号が反転する」という性質がある（これは $2 \times 2$ の場合に確認できるわ：$\begin{vmatrix} c & d \\ a & b \end{vmatrix} = cb - da = -(ad - bc) = -\begin{vmatrix} a & b \\ c & d \end{vmatrix}$。$3 \times 3$ でも余因子展開で確認できるけれど、ここでは結果を認めて使うわね）。巡回置換 $(1 \to 2 \to 3 \to 1)$ は 2 回の行入れ替え（例えば 1↔2 の後 2↔3）で実現できるから、符号は $(-1)^2 = +1$ 倍——つまり値が変わらない。これがスカラー三重積の巡回対称性の正体よ。

🔵 カイ: 幾何学的にはどういう意味があるんですか？

🟡 リナ: $\boldsymbol{a}$, $\boldsymbol{b}$, $\boldsymbol{c}$ が張る平行六面体の符号付き体積よ。$\boldsymbol{b} \times \boldsymbol{c}$ は底面（$\boldsymbol{b}$ と $\boldsymbol{c}$ の平行四辺形）に垂直で、大きさが底面積に等しいベクトル。それと $\boldsymbol{a}$ の内積をとると、底面積 $\times$ 高さ $=$ 体積になる。図 A.1「スカラー三重積と平行六面体の体積」を見ると、この構造がよく分かるわ。

図 A.1: スカラー三重積と平行六面体の体積。スカラー三重積 $\boldsymbol{a} \cdot (\boldsymbol{b} \times \boldsymbol{c})$ は、3 つのベクトルが張る平行六面体の符号付き体積に等しい。$\boldsymbol{b} \times \boldsymbol{c}$（橙）は底面に垂直で、大きさが底面積。それと $\boldsymbol{a}$ の内積が体積を与える。

⚪ メイ: 符号が正なら $\boldsymbol{a}$, $\boldsymbol{b}$, $\boldsymbol{c}$ が右手系を成し、負なら左手系ということね。

✅ 理解度チェック: スカラー三重積 $\boldsymbol{a} \cdot (\boldsymbol{b} \times \boldsymbol{c})$ の幾何学的意味は何でしょうか？また、巡回置換に対してどのような性質を持つでしょうか？

答え

3つのベクトルが張る平行六面体の符号付き体積を表す。巡回置換 $\boldsymbol{a} \to \boldsymbol{b} \to \boldsymbol{c} \to \boldsymbol{a}$ に対して値は変わらない（巡回対称性）。

📝 練習問題:

スカラー三重積の行列式計算 → 問題 B-3. スカラー三重積

A.3　ベクトル積どうしの内積（BAC-CAB 公式への布石）¶

🟡 リナ: 本編で曲率テンソルの幾何学的意味を議論したとき、次の公式を使ったわね。4 つのベクトル $\boldsymbol{a}$, $\boldsymbol{b}$, $\boldsymbol{c}$, $\boldsymbol{d}$ に対して：

\[ (\boldsymbol{a} \times \boldsymbol{b}) \cdot (\boldsymbol{c} \times \boldsymbol{d}) = (\boldsymbol{a} \cdot \boldsymbol{c})(\boldsymbol{b} \cdot \boldsymbol{d}) - (\boldsymbol{a} \cdot \boldsymbol{d})(\boldsymbol{b} \cdot \boldsymbol{c}) \]

🔵 カイ: 右辺の形、行列式っぽいですね。

🟡 リナ: これは $2 \times 2$ 行列式でも書ける：

\[ (\boldsymbol{a} \times \boldsymbol{b}) \cdot (\boldsymbol{c} \times \boldsymbol{d}) = \begin{vmatrix} \boldsymbol{a} \cdot \boldsymbol{c} & \boldsymbol{a} \cdot \boldsymbol{d} \\ \boldsymbol{b} \cdot \boldsymbol{c} & \boldsymbol{b} \cdot \boldsymbol{d} \end{vmatrix} \]

⚪ メイ: 対角線の積の差、$( \boldsymbol{a} \cdot \boldsymbol{c})(\boldsymbol{b} \cdot \boldsymbol{d}) - (\boldsymbol{a} \cdot \boldsymbol{d})(\boldsymbol{b} \cdot \boldsymbol{c})$ ということね。

🟡 リナ: 証明はスカラー三重積を使うと見通しがいいわ。$\boldsymbol{c} \times \boldsymbol{d}$ を一つのベクトル $\boldsymbol{e}$ とみなすと、A.2 のスカラー三重積の巡回対称性から：

\[ (\boldsymbol{a} \times \boldsymbol{b}) \cdot \boldsymbol{e} = \boldsymbol{a} \cdot (\boldsymbol{b} \times \boldsymbol{e}) \]

つまり $(\boldsymbol{a} \times \boldsymbol{b}) \cdot (\boldsymbol{c} \times \boldsymbol{d}) = \boldsymbol{a} \cdot [\boldsymbol{b} \times (\boldsymbol{c} \times \boldsymbol{d})]$。

ここで BAC-CAB 公式（ベクトル三重積の公式）を使う。この公式は成分計算で直接証明できる（証明は練習問題問題 M-1. BAC-CAB 公式の証明で確認してね）：

\[ \boldsymbol{b} \times (\boldsymbol{c} \times \boldsymbol{d}) = \boldsymbol{c}(\boldsymbol{b} \cdot \boldsymbol{d}) - \boldsymbol{d}(\boldsymbol{b} \cdot \boldsymbol{c}) \]

これを使うと：

\[ = \boldsymbol{a} \cdot [\boldsymbol{c}(\boldsymbol{b} \cdot \boldsymbol{d}) - \boldsymbol{d}(\boldsymbol{b} \cdot \boldsymbol{c})] \]

\[ = (\boldsymbol{a} \cdot \boldsymbol{c})(\boldsymbol{b} \cdot \boldsymbol{d}) - (\boldsymbol{a} \cdot \boldsymbol{d})(\boldsymbol{b} \cdot \boldsymbol{c}) \]

🔵 カイ: BAC-CAB って名前がついてるんですね。でも今の式は $\boldsymbol{b} \times (\boldsymbol{c} \times \boldsymbol{d})$ だから文字が違う……。

🟡 リナ: 歴史的に $\boldsymbol{a} \times (\boldsymbol{b} \times \boldsymbol{c}) = \boldsymbol{b}(\boldsymbol{a} \cdot \boldsymbol{c}) - \boldsymbol{c}(\boldsymbol{a} \cdot \boldsymbol{b})$ の形で「BAC - CAB」と覚えるのよ。右辺の頭文字を読むと B-A-C 引く C-A-B。文字を入れ替えただけで同じ構造。

🔵 カイ: なるほど、外側のベクトルが「真ん中」に入って内積を作るパターンですね。でも、結果が $\boldsymbol{b}$ と $\boldsymbol{c}$ の線形結合になるのはなぜですか？ $\boldsymbol{a}$ の方向の成分は出てこないんですか？

🟡 リナ: いい質問ね。$\boldsymbol{b} \times \boldsymbol{c}$ は $\boldsymbol{b}$ と $\boldsymbol{c}$ が張る平面に垂直なベクトルでしょう？それと $\boldsymbol{a}$ の外積をとると、結果は $\boldsymbol{b} \times \boldsymbol{c}$ に垂直——つまり $\boldsymbol{b}$ と $\boldsymbol{c}$ が張る平面の中に戻ってくるの。平面の中のベクトルは、その平面を張る $\boldsymbol{b}$ と $\boldsymbol{c}$ の線形結合で書ける。だから $\boldsymbol{a}$ の方向の成分は出てこないのよ。

⚪ メイ: 幾何学的に考えれば当たり前のことが、公式の形にきちんと反映されているのね。

✅ 理解度チェック: BAC-CAB 公式 $\boldsymbol{a} \times (\boldsymbol{b} \times \boldsymbol{c}) = \boldsymbol{b}(\boldsymbol{a} \cdot \boldsymbol{c}) - \boldsymbol{c}(\boldsymbol{a} \cdot \boldsymbol{b})$ において、ベクトル三重積の結果はどのような形で表されるでしょうか？

答え

外積の中の2つのベクトル $\boldsymbol{b}$, $\boldsymbol{c}$ の線形結合として表される。係数は外側のベクトル $\boldsymbol{a}$ との内積で決まる。

📝 練習問題:

Lagrange の恒等式の数値確認 → 問題 B-8. Lagrange 恒等式の確認、BAC-CAB 公式の数値確認 → 問題 B-9. BAC-CAB 公式の確認、BAC-CAB 公式の一般証明 → 問題 M-1. BAC-CAB 公式の証明、ベクトル積どうしの内積 → 問題 M-3. ベクトル積どうしの内積の公式、Levi-Civita 記号による統一的導出 → 問題 A-1. Levi-Civita 記号による恒等式

A.4　微分演算子¶

🟡 リナ: ここからは、ベクトル場やスカラー場に対する微分操作を定義するわ。本編で曲がった時空に一般化した共変微分の「平坦な空間版」がこれよ。

ナブラ演算子¶

🟡 リナ: まず、全ての微分演算子の母体となるナブラ演算子 $\nabla$（nabla）を導入する。3 次元直交座標 $(x, y, z)$ において：

\[ \nabla = \left(\frac{\partial}{\partial x},\; \frac{\partial}{\partial y},\; \frac{\partial}{\partial z}\right) \]

これは「ベクトルのように振る舞う微分演算子」だと思ってほしい。

🔵 カイ: ベクトルなのに中身が偏微分……？

🟡 リナ: そう、だから「演算子」なの。単独では意味がなくて、何かに作用させて初めて意味を持つ。スカラー場に作用させれば grad、ベクトル場と内積をとれば div、ベクトル場と外積をとれば rot になる。

勾配（gradient）¶

🟡 リナ: スカラー場 $\varphi(x, y, z)$ に $\nabla$ を作用させると：

\[ \operatorname{grad}\varphi = \nabla\varphi = \left(\frac{\partial\varphi}{\partial x},\; \frac{\partial\varphi}{\partial y},\; \frac{\partial\varphi}{\partial z}\right) \]

⚪ メイ: 結果はベクトル場ね。スカラーを入れるとベクトルが出てくる。

🟡 リナ: そう。そして物理的には、$\nabla\varphi$ は「$\varphi$ が最も急に増加する方向」を向き、その大きさが「変化率の最大値」を与える。各成分 $\partial\varphi/\partial x$ などが、その方向への変化率を表しているの。本編で重力ポテンシャル $\varphi$ から重力場 $\boldsymbol{g} = -\nabla\varphi$ を導いたのを覚えているわね。

発散（divergence）¶

🟡 リナ: ベクトル場 $\boldsymbol{F} = (F_x,\, F_y,\, F_z)$ に対して、$\nabla$ との「内積」をとると：

\[ \operatorname{div}\boldsymbol{F} = \nabla \cdot \boldsymbol{F} = \frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z} \]

🔵 カイ: 結果はスカラーですね。

🟡 リナ: そう。物理的には「その点からベクトル場がどれだけ湧き出しているか」を表す。$\operatorname{div}\boldsymbol{F} > 0$ なら湧き出し（source）、$< 0$ なら吸い込み（sink）。本編の Maxwell 方程式 $\operatorname{div}\boldsymbol{E} = \rho/\varepsilon_0$ は「電荷密度 $\rho$ がある場所から電場が湧き出す」という意味だったわ。

微小体積からの導出¶

🟡 リナ: なぜこの偏微分の和が「湧き出し」を表すのか、物理的イメージから導出しておきましょう。

点 $(x_0, y_0, z_0)$ を中心とする微小な直方体（辺 $\Delta x$, $\Delta y$, $\Delta z$）を考える（各面から $\boldsymbol{F}$ がどれだけ流出・流入するかを調べるイメージ——右面と左面、上面と下面、前面と後面の 3 組に分けて考えるの）。ベクトル場 $\boldsymbol{F}$ がこの箱から正味どれだけ流出するかを計算するの。図 A.2「発散の微小体積による導出」に $x$ 方向の例として右面と左面を通る $F_x$ の流出入の様子を描いたから、見ながら読んでね（$y$, $z$ 方向も同じ要領よ）。

図 A.2: 発散の微小体積による導出。点 $(x_0, y_0, z_0)$ を中心とする微小直方体の各面を通るベクトル場の流出入。右面と左面の $F_x$ の差が $x$ 方向の正味流出を与える。

$x$ 方向の正味流出: 右面（$x_0 + \Delta x/2$）から出る量から、左面（$x_0 - \Delta x/2$）から入る量を引く：

\[ \left[F_x\!\left(x_0 + \frac{\Delta x}{2}\right) - F_x\!\left(x_0 - \frac{\Delta x}{2}\right)\right]\Delta y\,\Delta z \]

🔵 カイ: 右面を通って出ていく分から、左面を通って入ってくる分を引くんですね。

🟡 リナ: その通り。$F_x$ を $x_0$ の周りで 1 次近似する。偏微分の定義 $\partial F_x/\partial x = \lim_{h\to 0}[F_x(x_0+h) - F_x(x_0)]/h$ を変形すると、$h$ が十分小さいとき $F_x(x_0 + h) \approx F_x(x_0) + (\partial F_x/\partial x)\cdot h$ という近似が得られる（この「微分を使って関数を近似する」手法を一般に Taylor 展開と呼ぶ——A.7 の Euler の公式の導出で、一般的な Taylor 展開の考え方を詳しく説明するわ）。これを使うと $F_x(x_0 \pm \Delta x/2) \approx F_x(x_0) \pm (\partial F_x/\partial x)(\Delta x/2)$。差をとれば $(\partial F_x/\partial x)\Delta x$。したがって $x$ 方向の正味流出量は $(\partial F_x/\partial x)\,\Delta V$（$\Delta V = \Delta x \Delta y \Delta z$）。

$y$ 方向・$z$ 方向も同様。3 方向合わせた単位体積あたりの正味流出量が発散：

\[ \nabla \cdot \boldsymbol{F} = \lim_{\Delta V \to 0}\frac{\text{（閉曲面からの正味流出量）}}{\Delta V} = \frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z} \]

ここで「閉曲面からの正味流出量」を数式で書くと $\oint_S \boldsymbol{F} \cdot d\boldsymbol{S}$ となる。初めて見る記号だと思うから説明するわね。

🔵 カイ: $\oint$ って普通の積分記号 $\int$ と何が違うんですか？

🟡 リナ: $\oint$ は、閉じた面（または曲線）全体にわたって足し合わせることを表す積分記号よ（普通の $\int$ に丸がついている）。そして $d\boldsymbol{S}$ は各微小面の面積ベクトル——面に垂直な方向を向いていて、大きさがその微小面の面積に等しいベクトルよ。「面に垂直な方向」のことを法線方向と呼ぶの——高校で「接線」が曲線に沿った方向だったのと対で、「法線」は面から突き出す方向のことね。向きは「外向き」——つまり閉曲面（風船のように閉じた面）の内側から外側に向かう方向——に取る約束。こうすると内積 $\boldsymbol{F} \cdot d\boldsymbol{S}$ が正なら流出、負なら流入を表すことになるの。

🔵 カイ: 面積分って、高校の積分みたいに「短冊を足す」のを面の上でやるイメージですか？

🟡 リナ: そう、まさにそのイメージ。高校の定積分は「曲線の下の短冊を足す」——1 次元に沿って足し合わせる操作よね。面積分は「曲面を細かいタイルに分けて、各タイルで『$\boldsymbol{F}$ の法線成分 $\times$ タイルの面積』を計算し、全タイル分を足し合わせる」——2 次元の面の上で足し合わせる操作よ。

🔵 カイ: 「法線成分」って、面に垂直な方向の成分ですよね。でもそれをどうやって取り出すんですか？

🟡 リナ: 内積を使うの。$d\boldsymbol{S}$ は面に垂直な方向を向いた面積ベクトルだから、内積 $\boldsymbol{F} \cdot d\boldsymbol{S} = |\boldsymbol{F}||d\boldsymbol{S}|\cos\alpha$（$\alpha$ は $\boldsymbol{F}$ と $d\boldsymbol{S}$ のなす角）をとると、$|\boldsymbol{F}|\cos\alpha$ = 法線方向成分が自動的に拾われる。これが面積分よ。そしてこの「局所的な偏微分の和」と「閉曲面からの流出量」の関係を、有限の体積に拡張したものが A.6 の Gauss の定理になるの。

⚪ メイ: つまり、リナさんが最初に示した「閉曲面からの流出量を体積で割る極限」が意味の定義で、そこから導かれた「偏微分の和」が実際の計算に使う公式ということね。

✅ 理解度チェック: 発散 $\nabla \cdot \boldsymbol{F}$ は、微小体積を使ってどのように物理的に定義されるか？

答え

点を囲む微小閉曲面からのベクトル場の正味流出量を、その微小体積で割った極限として定義される。正なら湧き出し、負なら吸い込みを意味する。

回転（rotation / curl）¶

🟡 リナ: ベクトル場 $\boldsymbol{F}$ に対して、$\nabla$ との「外積」をとると：

\[ \operatorname{rot}\boldsymbol{F} = \nabla \times \boldsymbol{F} = \begin{vmatrix} \boldsymbol{e}_x & \boldsymbol{e}_y & \boldsymbol{e}_z \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ F_x & F_y & F_z \end{vmatrix} \]

成分で書くと：

\[ \operatorname{rot}\boldsymbol{F} = \left(\frac{\partial F_z}{\partial y} - \frac{\partial F_y}{\partial z},\; \frac{\partial F_x}{\partial z} - \frac{\partial F_z}{\partial x},\; \frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}\right) \]

🔵 カイ: 各成分が「引き算」になってますね。これがベクトル積の反対称性から来てるってことか。

🟡 リナ: その通り。物理的には「ベクトル場がその点の周りでどれだけ渦を巻いているか」を表す。$\operatorname{rot}\boldsymbol{B} = \mu_0 \boldsymbol{j}$ は「電流が流れている場所の周りで磁場が渦を巻く」という意味よ。

微小ループからの導出¶

🟡 リナ: 発散を「微小体積からの流出」で導いたのと同じように、回転は「微小ループに沿った循環」として導出できるの。

$xy$$xy$ 平面内の微小長方形ループ（中心 $(x_0, y_0, z_0)$、辺 $\Delta x$, $\Delta y$）を反時計回りに一周する。各辺に沿って「$\boldsymbol{F}$ の辺方向成分 $\times$ 辺の長さ」を計算し、一周分を合計したものを循環（circulation）と呼ぶの——ループを一周「循環」するから、この名前がついている。ポイントは、対辺どうし（上辺と下辺、右辺と左辺）の寄与の差が回転の成分を与えるということ。図 A.3「回転の微小ループによる導出」に微小ループの各辺と $\boldsymbol{F}$ の成分の関係を描いたから、見ながら読んでね。

図 A.3: 回転の微小ループによる導出。$xy$ 平面内の微小長方形ループを反時計回りに一周する。各辺に沿った $\boldsymbol{F}$ の進行方向成分を足し合わせたもの（循環）が、対辺どうしの差として回転の成分を与える。

正式には線積分 $\oint \boldsymbol{F} \cdot d\boldsymbol{r}$ と書く。線積分とは「曲線を細かい区間に分割し、各区間で $\boldsymbol{F}$ の進行方向成分 $\times$ 区間の長さを計算して、全区間分を足し合わせる」操作——高校の定積分が「$x$ 軸に沿って短冊を足す」のと同じ要領で、任意の曲線に沿って足し合わせるものよ。$d\boldsymbol{r}$ は曲線に沿った微小変位ベクトルで、内積 $\boldsymbol{F} \cdot d\boldsymbol{r}$ は「$\boldsymbol{F}$ の進行方向成分 $\times$ 微小距離」を意味する。ただしここでは微小ループだから各辺ごとに計算すれば十分。具体的に：

下辺（$y = y_0 - \Delta y/2$, 反時計回りなので $+x$ 方向に進む）: $+F_x(x_0, y_0 - \Delta y/2)\,\Delta x$
上辺（$y = y_0 + \Delta y/2$, 反時計回りなので $-x$ 方向に進む）: $-F_x(x_0, y_0 + \Delta y/2)\,\Delta x$
右辺（$x = x_0 + \Delta x/2$, 反時計回りなので $+y$ 方向に進む）: $+F_y(x_0 + \Delta x/2, y_0)\,\Delta y$
左辺（$x = x_0 - \Delta x/2$, 反時計回りなので $-y$ 方向に進む）: $-F_y(x_0 - \Delta x/2, y_0)\,\Delta y$

符号の決め方: 進行方向と座標軸の正方向が同じなら $+$、逆なら $-$。反時計回りに一周すると、下辺は右へ（$+x$）、右辺は上へ（$+y$）、上辺は左へ（$-x$）、左辺は下へ（$-y$）と進む（図 A.3「回転の微小ループによる導出」の矢印の向きを確認してね）。

⚪ メイ: 発散のときは「面を通る流出入」を見たけれど、回転は「辺に沿った成分」を見るのね。

🟡 リナ: その通り。では右辺と左辺（$y$ 方向に平行な辺）の寄与をまとめよう。1 次近似（発散の導出と同じ Taylor 展開）を使うと、右辺と左辺の差は $[F_y(x_0 + \Delta x/2) - F_y(x_0 - \Delta x/2)]\,\Delta y \approx (\partial F_y/\partial x)\Delta x\,\Delta y$。同様に上辺と下辺（$x$ 方向に平行な辺）の寄与をまとめると、下辺から上辺を引いた差は $-[F_x(x_0, y_0 + \Delta y/2) - F_x(x_0, y_0 - \Delta y/2)]\,\Delta x \approx -(\partial F_x/\partial y)\Delta y\,\Delta x$。全体で：

\[ \Gamma_z = \left(\frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}\right)\Delta A \qquad (\Delta A = \Delta x \Delta y) \]

単位面積あたりの循環が回転の $z$ 成分：

\[ (\nabla \times \boldsymbol{F})_z = \lim_{\Delta A \to 0}\frac{\Gamma_z}{\Delta A} = \frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y} \]

⚪ メイ: $yz$ 平面と $zx$ 平面のループでも同じことをすれば、$x$ 成分と $y$ 成分が出るわね。

🟡 リナ: そう。そしてこの「微小ループの循環」と「回転の面積分」の関係を有限の曲面に拡張したものが、A.6 の Stokes の定理よ。

🔵 カイ: 発散が「微小体積からの正味流出」、回転が「微小ループに沿った循環」——幾何学的にくっきりしました。でも、発散は体積で割って、回転は面積で割る——次元が違うのに、どちらも「$\nabla$ の仲間」なんですね。

🟡 リナ: いい着眼点ね。$\nabla$ という一つの演算子が、スカラー場に作用するか（grad）、ベクトル場と内積をとるか（div）、外積をとるか（rot）で、全く異なる物理量を生み出す。入出力の型も全部違う——grad は「スカラー → ベクトル」、div は「ベクトル → スカラー」、rot は「ベクトル → ベクトル」。実は第 24 章で導入した微分形式と外微分 $d$ を使うと、grad・div・rot は全て「$d$ を作用させる」という一つの操作に統一されるの。でも今は $\nabla$ の 3 つの顔を使い分けられれば十分よ。

⚪ メイ: なるほど、入出力の型が全部違うから、間違えて組み合わせるとそもそも意味をなさないのね。

🔵 カイ: じゃあ、ラプラシアン $\nabla^2$ は「grad をとってから div をとる」だから、$\nabla$ の顔を 2 回使ってるってことですね。他の組み合わせ——例えば「div をとってから grad」とか「rot をとってから rot」もあるんですか？

🟡 リナ: いい質問。「div をとってから grad」は $\nabla(\nabla \cdot \boldsymbol{F})$ で、次の A.5 で出てくる「rot の rot」の公式に登場するわ。「rot をとってから rot」もまさに A.5 で扱う。順番に見ていきましょう。

✅ 理解度チェック: 回転 $\nabla \times \boldsymbol{F}$ の $z$ 成分は、微小ループを使ってどのように物理的に定義されるか？

答え

$xy$ 平面内の微小ループに沿ったベクトル場の循環（線積分）を、そのループが囲む微小面積で割った極限として定義される。これは「その点の周りでベクトル場がどれだけ渦を巻いているか」を表す。

ラプラシアン（Laplacian）¶

🟡 リナ: $\nabla$ を 2 回作用させる操作も重要よ。スカラー場 $\varphi$ に対して：

\[ \nabla^2\varphi = \nabla \cdot (\nabla\varphi) = \frac{\partial^2\varphi}{\partial x^2} + \frac{\partial^2\varphi}{\partial y^2} + \frac{\partial^2\varphi}{\partial z^2} \]

教科書によっては $\nabla^2$ の代わりに $\Delta$ と書くこともあるけれど、この本では微小変化量 $\Delta x$ との混同を避けるため、ラプラシアンには $\nabla^2$ を使うわ。

⚪ メイ: grad をとってから div をとる、という 2 段階の操作ね。

🟡 リナ: Newton の重力場の方程式 $\nabla^2\varphi = 4\pi G\rho$ は「質量密度 $\rho$ がポテンシャル $\varphi$ の曲がり具合を決める」という意味。本編で Einstein 方程式が Newton の方程式を一般化したものだと説明したとき、この対応が鍵だったわね。

✅ 理解度チェック: ナブラ演算子 $\nabla$ をスカラー場に作用させると何が得られるでしょうか？ベクトル場と内積をとると何が得られるでしょうか？

答え

スカラー場に作用させると勾配（gradient）——最急上昇方向と変化率を表すベクトル場。ベクトル場と内積をとると発散（divergence）——湧き出しの強さを表すスカラー場。

📝 練習問題:

勾配の計算 → 問題 B-4. スカラー場の勾配、発散の計算 → 問題 B-5. ベクトル場の発散、回転の計算 → 問題 B-6. ベクトル場の回転、ラプラシアンの計算 → 問題 B-7. スカラー場のラプラシアン

A.5　微分演算子の重要な恒等式¶

🟡 リナ: 以下の 2 つの恒等式は、本編で何度も使った。

恒等式 1：勾配の回転はゼロ¶

\[ \operatorname{rot}(\operatorname{grad}\varphi) = \nabla \times (\nabla\varphi) = \boldsymbol{0} \]

🔵 カイ: 直感的には分かる気もするけど……勾配って「坂を登る方向」ですよね。坂を登る方向に沿って進んでも、ぐるっと一周して元に戻ることはない——だから渦がないってことですか？

🟡 リナ: その通り。物理的にも確認できるわ。成分で書くと、例えば $x$ 成分は $\frac{\partial^2\varphi}{\partial y\partial z} - \frac{\partial^2\varphi}{\partial z\partial y}$。$\varphi$ が十分滑らかなら偏微分の順序を交換できるから、これは $0$ になる。つまり、保存力場 $\boldsymbol{F} = -\nabla\varphi$ のように「ポテンシャルから導かれる場」には渦がないの。

⚪ メイ: なるほど、A.0 の Schwarz の定理で混合偏微分が相殺するから、各成分がゼロになるのね。

🔵 カイ: 逆に、$\operatorname{rot}\boldsymbol{F} = \boldsymbol{0}$ なら必ずポテンシャルが存在するんですか？

🟡 リナ: いい質問。穴のない領域——例えばボールの内部のような、輪ゴムをどこに置いても 1 点にキュッと縮められる領域——ではそう言えるわ（数学ではこういう領域を「単連結」と呼ぶ）。逆にドーナツの表面のように穴がある領域では、穴を通り抜けるように巻きついた輪ゴムは切らない限り 1 点に縮められない——そういうループに沿って $\operatorname{rot}\boldsymbol{F} = \boldsymbol{0}$ でも循環がゼロにならない反例が作れるの。今はこの直感だけ持っておけば十分。A.6 で Stokes の定理を学んだ後に、「$\operatorname{rot}\boldsymbol{F} = \boldsymbol{0}$ なら任意の閉曲線で循環がゼロ」がどう導かれるか、考えてみてね。

🔵 カイ: なるほど、穴がなければ大丈夫ってことですね。A.6 の Stokes の定理を学んだら、もう少しはっきり分かりそうだ。

恒等式 2：回転の発散はゼロ¶

\[ \operatorname{div}(\operatorname{rot}\boldsymbol{F}) = \nabla \cdot (\nabla \times \boldsymbol{F}) = 0 \]

🟡 リナ: これも恒等式 1 と全く同じ論法——Schwarz の定理——から従う。$\nabla \cdot (\nabla \times \boldsymbol{F})$ を展開すると、例えば $\partial_x(\partial_y F_z - \partial_z F_y) + \cdots$ のように、各項が混合偏微分の差の形になって全て相殺するの。詳しい確認は練習問題に回すわね。物理的には、$\operatorname{div}\boldsymbol{B} = 0$（磁気単極子は存在しない）は、$\boldsymbol{B} = \operatorname{rot}\boldsymbol{A}$ と書けることの帰結でもある。

🔵 カイ: 恒等式 1 が「ポテンシャルがあれば渦なし」で、恒等式 2 が「渦から生まれた場には湧き出しなし」——対になってるんですね。

恒等式 3：回転の回転¶

\[ \operatorname{rot}(\operatorname{rot}\boldsymbol{F}) = \nabla \times (\nabla \times \boldsymbol{F}) = \nabla(\nabla \cdot \boldsymbol{F}) - \nabla^2 \boldsymbol{F} \]

🟡 リナ: この恒等式は、Maxwell 方程式から電磁波の波動方程式を導出するときに決定的。$x$ 成分で確認しましょう。

$\nabla \times \boldsymbol{F} = \boldsymbol{B}$ と置いて成分を書き出すと、$B_x = \partial_y F_z - \partial_z F_y$、$B_y = \partial_z F_x - \partial_x F_z$、$B_z = \partial_x F_y - \partial_y F_x$。

$\nabla \times \boldsymbol{B}$ の $x$ 成分は $\partial_y B_z - \partial_z B_y$。代入して：

\[ [\nabla \times (\nabla \times \boldsymbol{F})]_x = \partial_y(\partial_x F_y - \partial_y F_x) - \partial_z(\partial_z F_x - \partial_x F_z) \]

\[ = \partial_y \partial_x F_y - \partial_y^2 F_x - \partial_z^2 F_x + \partial_z \partial_x F_z \]

⚪ メイ: Schwarz の定理で微分の順序を入れ替えて整理するのね。

🟡 リナ: その通り。Schwarz の定理で $\partial_y\partial_x = \partial_x\partial_y$、$\partial_z\partial_x = \partial_x\partial_z$ を使うと：

\[ = \partial_x\partial_y F_y - \partial_y^2 F_x - \partial_z^2 F_x + \partial_x\partial_z F_z \]

ここで $\nabla \cdot \boldsymbol{F} = \partial_x F_x + \partial_y F_y + \partial_z F_z$ を作りたいから、$\partial_x^2 F_x$ を足して引く（$+\partial_x^2 F_x - \partial_x^2 F_x = 0$ だから値は変わらない）：

\[ = \partial_x(\partial_x F_x + \partial_y F_y + \partial_z F_z) - (\partial_x^2 + \partial_y^2 + \partial_z^2)F_x \]

\[ = \partial_x(\nabla \cdot \boldsymbol{F}) - \nabla^2 F_x = [\nabla(\nabla \cdot \boldsymbol{F})]_x - (\nabla^2 \boldsymbol{F})_x \]

$y$, $z$ 成分も同様だから、ベクトルとして恒等式が成り立つ。

🔵 カイ: じゃあ、もし $\nabla \cdot \boldsymbol{F} = 0$ の場合は $\nabla(\nabla \cdot \boldsymbol{F})$ が消えて $\nabla \times (\nabla \times \boldsymbol{F}) = -\nabla^2 \boldsymbol{F}$ になるから、ラプラシアンだけが残る。本編で電磁波の波動方程式が出てきたのは、電荷のない領域で $\nabla \cdot \boldsymbol{E} = 0$ だったからですか？

🟡 リナ: その通り。

🔵 カイ: じゃあ逆に、電荷がある領域だと単純な波動方程式にはならない？

🟡 リナ: その通り。電荷がある領域では $\nabla \cdot \boldsymbol{E} = \rho/\varepsilon_0 \neq 0$ だから、$\nabla(\nabla \cdot \boldsymbol{E})$ の項が残って、純粋な波動方程式にはならない。右辺に電荷分布の勾配が「源」として現れるの。

⚪ メイ: つまり、$\nabla \cdot \boldsymbol{F} = 0$ のときだけラプラシアンだけが残って波動方程式になる——発散がゼロかどうかで方程式の形が変わるのね。

🔵 カイ: じゃあ重力波も「真空中」だから波動方程式になったんですか？物質がある領域だとどうなるんだろう。

🟡 リナ: いい着眼点ね。本編 ch19 で扱ったように、真空中では線形化された Einstein 方程式が波動方程式になって重力波が伝播する。物質がある領域では右辺に「源」が残って、波の生成が起きる——四重極公式がまさにその「源つきの方程式」から導かれたものだったわ。

⚪ メイ: なるほど、真空かどうかで方程式の形が変わって、伝播と生成が分かれるのね。

🟡 リナ: その通り。さて、ここまでで微分演算子の重要な恒等式が揃ったわ。次は積分定理に進みましょう。

📝 練習問題:

rot(grad φ) = 0 の確認 → 問題 B-10. rot(grad) = 0 の確認、div(rot F) = 0 の証明 → 問題 M-2. div(rot) = 0 の証明、積の発散の公式 → 問題 M-4. 積の発散の公式

A.6　積分定理¶

🟡 リナ: 最後に、微分と積分を結ぶ 2 つの大定理を述べるわ。証明は省略するけれど、物理的な意味を押さえておくことが重要よ。

Gauss の定理（発散定理）¶

\[ \oint_S \boldsymbol{F} \cdot d\boldsymbol{S} = \int_V \operatorname{div}\boldsymbol{F}\; dV \]

ここで $V$ は閉曲面 $S$ で囲まれた領域、$d\boldsymbol{S}$ は外向き法線方向の面積要素。右辺の $\int_V \cdots dV$ は体積積分——領域 $V$ を微小な立方体に分割し、各立方体で $\operatorname{div}\boldsymbol{F}$ の値 $\times$ 微小体積 $dV$ を計算して、全ての立方体分を足し合わせる操作よ。面積分が「面を細かいタイルに分けて足す」だったのに対し、体積積分は「立体を細かいブロックに分けて足す」——次元が 1 つ上がっただけで同じ発想ね。

🔵 カイ: 左辺は「面を通って出ていく流量の合計」、右辺は「中にある湧き出しの合計」。それが等しいって言ってるんですね。

🟡 リナ: まさにそう。本編で Einstein 方程式の物理的意味を議論したとき、「局所的な曲率（div に相当）と大域的な測地線のずれ（面積分に相当）の関係」として Gauss の定理のアナロジーを使ったわね。

⚪ メイ: A.4 で「微小体積の正味流出量を体積で割る」のが発散だったから、それを有限の体積に足し上げたのがこの定理——局所と大域をつなぐ橋渡しね。

✅ 理解度チェック: Gauss の定理（発散定理）は、どのような 2 つの量が等しいことを述べているでしょうか？

答え

閉曲面 $S$ を通るベクトル場の正味流出量（面積分 $\oint_S \boldsymbol{F} \cdot d\boldsymbol{S}$）と、その閉曲面で囲まれた領域内の発散の体積積分（$\int_V \operatorname{div}\boldsymbol{F}\, dV$）が等しい。

Stokes (ストークス) の定理¶

\[ \oint_C \boldsymbol{F} \cdot d\boldsymbol{r} = \int_S (\operatorname{rot}\boldsymbol{F}) \cdot d\boldsymbol{S} \]

ここで $C$ は曲面 $S$ の境界となる閉曲線、$d\boldsymbol{r}$ は曲線に沿った微小変位ベクトル（A.4 で導入した線積分の記号と同じ）。

⚪ メイ: 左辺は「閉曲線に沿った循環」、右辺は「その曲線で囲まれた面を貫く渦の合計」。

🟡 リナ: 本編で平行移動の経路依存性（ホロノミー）を議論したとき、「閉曲線に沿ってベクトルを運んだときのずれ」が「囲まれた面の曲率の面積分」に等しい、という話をしたでしょう？あれは Stokes の定理の曲がった空間版だったのよ。

🔵 カイ: あー、だから本編で「Stokes の定理の一般化」って言ってたんですね。ここに戻ってくる。

✅ 理解度チェック: Stokes の定理は、どのような 2 つの量が等しいことを述べているでしょうか？

答え

閉曲線 $C$ に沿ったベクトル場の循環（線積分 $\oint_C \boldsymbol{F} \cdot d\boldsymbol{r}$）と、その閉曲線を境界とする曲面 $S$ 上での回転の面積分（$\int_S (\operatorname{rot}\boldsymbol{F}) \cdot d\boldsymbol{S}$）が等しい。

📝 練習問題:

Gauss の定理の応用（クーロン電場） → 問題 M-5. Gauss の定理とクーロン電場、Stokes の定理と保存場・曲率の関係 → 問題 A-2. Stokes の定理の応用

A.7　波動方程式 — 2 階偏微分方程式の分類と d'Alembert の解¶

🟡 リナ: ここからは、本編で繰り返し登場する 2 階偏微分方程式を分類しておくわ。ch01 で Poisson 方程式と電磁気の波動方程式を比較し、ch19 では Einstein 方程式の線形化から重力波の波動方程式を導出した。後の弦理論編でも弦の振動が波動方程式に従うの。これらをまとめて整理しておくことで、解の振る舞いを見通しよく理解できる。

なぜ分類が必要か¶

物理学で登場する偏微分方程式は、見た目は似ていても解の性質が全く異なる。波動方程式の解は振動し、拡散方程式の解は減衰し、Poisson 方程式の解は静的。方程式のタイプを見分けることで、解の振る舞いを予測できる。

3 つのタイプ¶

表 A.1: 2階偏微分方程式の3タイプと物理的意味

タイプ	標準形（1 次元）	物理的な意味	登場する章
波動方程式（双曲型）	$\partial_t^2 f = v^2 \partial_x^2 f$	波の伝播。情報が有限の速度 $v$ で伝わる	第 1 章（電磁波）、第 19 章（重力波）、弦理論編（弦の振動、後の編で扱う）
拡散方程式（放物型）	$\partial_t f = D \partial_x^2 f$	熱の拡散、粒子の拡散。不可逆過程	古典的な例: 熱伝導方程式
Poisson/Laplace 方程式（楕円型）	$\nabla^2 f = \rho$	静的な場の分布。時間に依存しない	ch01（重力ポテンシャル）、静電場

見分け方¶

時間の 2 階微分がある → 波動方程式（振動する解）
時間の 1 階微分がある → 拡散方程式（減衰する解）
時間微分がない → Poisson/Laplace 方程式（静的な解）

⚪ メイ: 方程式の「顔つき」を見るだけで、解の挙動が予測できるのね。

1 次元波動方程式の一般解（d'Alembert の解）¶

🟡 リナ: 弦理論の ch13 で弦の振動を扱うとき、波動方程式の一般解を使うわ。ここで変数変換による導出を見ておきましょう。出発点：

\[ \frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2} \]

ステップ 1: 変数変換 $\xi = x - ct$, $\eta = x + ct$ と置く。逆に $x = (\xi + \eta)/2$, $t = (\eta - \xi)/(2c)$。

🔵 カイ: なんでこの変数変換を思いつくんですか？

🟡 リナ: 波動方程式は「右に進む波」と「左に進む波」を含むはず。$x - ct$ は右進行波の位相、$x + ct$ は左進行波の位相だから、この 2 つを新しい座標にすれば方程式が簡単になると期待できるの。

ステップ 2: 連鎖律で偏微分を変換

A.0 の連鎖律を使う。$\xi = x - ct$, $\eta = x + ct$ だから $\partial\xi/\partial x = 1$, $\partial\eta/\partial x = 1$, $\partial\xi/\partial t = -c$, $\partial\eta/\partial t = c$。よって：

\[ \frac{\partial u}{\partial x} = \frac{\partial u}{\partial \xi}\frac{\partial \xi}{\partial x} + \frac{\partial u}{\partial \eta}\frac{\partial \eta}{\partial x} = \frac{\partial u}{\partial \xi} + \frac{\partial u}{\partial \eta} \]

\[ \frac{\partial u}{\partial t} = \frac{\partial u}{\partial \xi}\frac{\partial \xi}{\partial t} + \frac{\partial u}{\partial \eta}\frac{\partial \eta}{\partial t} = -c\frac{\partial u}{\partial \xi} + c\frac{\partial u}{\partial \eta} = c\left(-\frac{\partial u}{\partial \xi} + \frac{\partial u}{\partial \eta}\right) \]

⚪ メイ: A.0 の連鎖律がここで早速使われるのね。

ステップ 3: 2 階微分を計算

🟡 リナ: ステップ 2 の $\partial/\partial x = \partial/\partial\xi + \partial/\partial\eta$ をもう一度作用させると：

\[ \frac{\partial^2 u}{\partial x^2} = \left(\frac{\partial}{\partial \xi} + \frac{\partial}{\partial \eta}\right)^2 u = \frac{\partial^2 u}{\partial \xi^2} + 2\frac{\partial^2 u}{\partial \xi \partial \eta} + \frac{\partial^2 u}{\partial \eta^2} \]

同様に $\partial/\partial t = c(-\partial/\partial\xi + \partial/\partial\eta)$ をもう一度作用させると：

\[ \frac{\partial^2 u}{\partial t^2} = c^2\left(-\frac{\partial}{\partial \xi} + \frac{\partial}{\partial \eta}\right)^2 u = c^2\left(\frac{\partial^2 u}{\partial \xi^2} - 2\frac{\partial^2 u}{\partial \xi \partial \eta} + \frac{\partial^2 u}{\partial \eta^2}\right) \]

ステップ 4: 波動方程式に代入 波動方程式 $\partial_t^2 u = c^2 \partial_x^2 u$ の左辺と右辺を新しい変数で書き直すと：

左辺: $\partial_t^2 u = c^2(\partial_{\xi}^2 u - 2\partial_{\xi}\partial_{\eta} u + \partial_{\eta}^2 u)$
右辺: $c^2 \partial_x^2 u = c^2(\partial_{\xi}^2 u + 2\partial_{\xi}\partial_{\eta} u + \partial_{\eta}^2 u)$

波動方程式 $\partial_t^2 u = c^2 \partial_x^2 u$ に代入すると、$c^2$ が両辺に共通因子としてかかっている。$c^2 \neq 0$ で両辺を割ると：

\[ \partial_{\xi}^2 u - 2\partial_{\xi}\partial_{\eta}u + \partial_{\eta}^2 u = \partial_{\xi}^2 u + 2\partial_{\xi}\partial_{\eta}u + \partial_{\eta}^2 u \]

両辺から $\partial_{\xi}^2 u + \partial_{\eta}^2 u$ を引くと、残るのは混合微分の項だけ：$-2\partial_{\xi}\partial_{\eta}u = +2\partial_{\xi}\partial_{\eta}u$。移項すると $-4\partial_{\xi}\partial_{\eta}u = 0$。すなわち

\[ \frac{\partial^2 u}{\partial \xi \partial \eta} = 0 \]

🔵 カイ: おお、あんなに複雑だった波動方程式が、変数変換するだけでこんなにシンプルになるんですね！

⚪ メイ: $\partial_\xi^2$ と $\partial_\eta^2$ の項が全部消えて、混合微分だけが残ったわ。

ステップ 5: 一般解

🟡 リナ: $\partial^2 u/(\partial\xi\partial\eta) = 0$ を $\eta$ について積分する。1 変数の場合、$df/dx = 0$ を積分すると $f = C$（定数）になるわね。偏微分の場合は「$\eta$ 以外の変数（つまり $\xi$）は定数扱い」だから、$\partial(\partial u/\partial\xi)/\partial\eta = 0$ を $\eta$ で積分すると、$\partial u/\partial\xi$ は $\eta$ に依存しない——つまり $\partial u/\partial\xi = h(\xi)$（$\xi$ のみの任意関数）。

🔵 カイ: ちょっと待ってください。1 変数なら「微分がゼロ → 定数」ですよね。なんで偏微分だと「定数」じゃなくて「$\xi$ の関数」になるんですか？

🟡 リナ: いい質問。$\eta$ で積分するとき $\xi$ は定数扱いだから、「$\eta$ に依存しない量」であれば何でもよいの。1 変数の類推で考えてみて：$df/dx = 0$ の解は $f = C$（定数）——$x$ に依存しない量。同じ理屈で、$\partial(\cdots)/\partial\eta = 0$ の解は「$\eta$ に依存しない量」= $\xi$ の任意関数。つまり 1 変数での「積分定数 $C$」が、偏微分では「もう一方の変数の任意関数」に格上げされるのよ。

⚪ メイ: なるほど、「$\eta$ について見れば定数」だけど「$\xi$ については自由」だから、$\xi$ の関数になるのね。

🟡 リナ: その通り。さらに $\partial u/\partial\xi = h(\xi)$ を $\xi$ について積分するわ。$\eta$ を定数とみなせば、この式は 1 変数の微分方程式 $dF/dx = h(x)$ と同じ形だから、普通の不定積分と同じ要領で積分できる。

1 変数の類推で考えてみて。$dF/dx = h(x)$ を積分すると $F(x) = \int h(x)\,dx + C$（$C$ は積分定数）よね。今の場合、$\partial u/\partial\xi = h(\xi)$ を $\xi$ で積分すると $u = f(\xi) + (\text{「積分定数」に相当するもの})$。ここで $f(\xi)$ は $h(\xi)$ の不定積分——$f$ は $\xi$ だけの関数で、$df/d\xi = h(\xi)$ を満たすわ。$h$ は「$\xi$ のみの任意関数」だったから、$f$ もまた「$\xi$ のみの任意関数」ということになるの（$f$ を $\xi$ で微分したものが $h$ だから、$f$ は滑らかな関数よ）。具体例を挙げると、$h(\xi) = \xi^2$ なら $f(\xi) = \xi^3/3$、$h(\xi) = \sin\xi$ なら $f(\xi) = -\cos\xi$、…と $h$ が変われば $f$ も変わるけれど、どちらも「$\xi$ だけの関数」であることに変わりはないわ。

では「積分定数」はどうなるか？ 1 変数なら $C$ は「$x$ に依存しない定数」だった。偏微分の場合は「$\xi$ に依存しない」だけでよいから、$C$ が $\eta$ の任意関数 $g(\eta)$ に格上げされるの。

結局 $u = f(\xi) + g(\eta)$ で、$f$ も $g$ も任意関数。元の変数に戻すと：

\[ \boxed{u(x, t) = f(x - ct) + g(x + ct)} \]

🔵 カイ: この $f$ と $g$ って、具体的にはどんな関数なんですか？三角関数じゃなくてもいい？

🟡 リナ: $f(x - ct)$ は右に速さ $c$ で進む波、$g(x + ct)$ は左に進む波を表している。そして $f$, $g$ は正弦波でもパルスでも、微分可能な任意の関数が解になるの。具体的にどんな形になるかは初期条件（$t = 0$ での波の形と速度）で決まる。これが弦理論の ch13 で弦の振動を左進行波と右進行波に分解するときの基礎よ。

⚪ メイ: つまり、初期条件さえ決まれば $f$ と $g$ は一意に決まるということね。

🔵 カイ: 「初期条件で決まる」って言われても、具体的にどんな形になるのかイメージが湧かないんですよね……。高校の物理だと波って $\sin$ か $\cos$ で書いてたけど、あれはこの一般解の特殊な場合ってことですか？

🟡 リナ: そう。例えば $f(s) = e^{-s^2}$（ガウス型のパルス）を選べば、$u = e^{-(x-ct)^2}$ は山型の波が右に速さ $c$ で進む解になる。$\sin$ や $\cos$ は無限に続く正弦波だけど、パルスのように局在した波も解になるのよ。高校で $\sin$ だけを使っていたのは、周期的な波を扱う場面が多かったからで、波動方程式自体はもっと一般的な解を許すの。

平面波解と分散関係¶

🟡 リナ: いい質問ね。では、d'Alembert の一般解の中で最も基本的な「正弦波」の場合を考えてみましょう。

まず正弦波の書き方を整理してから、それを効率よく扱うための道具（Taylor 展開と Euler の公式）を導入し、最後に波動方程式に代入して分散関係を導く——3 ステップで進めるわ。量子力学以降でも必須の道具だから、ここでしっかり導出しておくわね。

波長 $\lambda$（波 1 つ分の空間的な長さ）、振動数 $\nu$（1 秒あたりの振動回数。高校では $f$ と書くことが多いけれど、物理学では $\nu$（ギリシャ文字のニュー）もよく使うわ）の正弦波を考えましょう。$\cos\theta$ は $\theta$ が $2\pi$ 増えるごとに 1 周期分振動するわね。だから「$x$ が $\lambda$ 進むと 1 周期」にしたければ、$\cos$ の引数を $\theta = 2\pi x/\lambda$ とすればいい——$x$ が $\lambda$ 増えると引数が $2\pi$ 増えて、ちょうど 1 周期分になるわ。同様に「$t$ が $1/\nu$ 経つと 1 周期」にしたければ $2\pi\nu\, t$ を引けばいい。こうして正弦波は $\cos(2\pi x/\lambda - 2\pi\nu\, t)$ と書ける。

🔵 カイ: $2\pi/\lambda$ とか $2\pi\nu$ って組み合わせが何度も出てきますね。略記はないんですか？

🟡 リナ: あるわ。波数 $k = 2\pi/\lambda$ と角振動数 $\omega = 2\pi\nu$ を定義しておくと便利なの。ここで位相とは $\cos$ の引数のこと——つまり $\cos(kx - \omega t)$ の中の $kx - \omega t$ の値よ。位相が $2\pi$ 変わると $\cos$ が一周して元に戻る。$k$ は「1 m 進むごとに位相がどれだけ増えるか」を表す量——波長 $\lambda$ が短いほど $k$ は大きくなる。例えば $\lambda = 2$ m なら $k = 2\pi/2 = \pi$ rad/m で、1 m 進むと位相が $\pi$（= 半周期分）ずれるわ。$\omega$ は「1 秒あたりに位相がどれだけ増えるか」——振動数 $\nu$ が大きいほど $\omega$ も大きい。どちらも $2\pi$ がつくのは、波 1 周期分の位相変化が $2\pi$ ラジアン（= 360°= $\cos$ が一周して元に戻る分）だから。

⚪ メイ: これで正弦波は $\cos(kx - \omega t)$ と書けるわね。

🟡 リナ: そう。そしてこれを $\cos[k(x - (\omega/k)t)]$ と書き直すと、d'Alembert の $f(x - ct)$ の形をしている——$\omega/k$ が波の速さに対応するの。さて、この正弦波 $\cos(kx - \omega t)$ をもっと扱いやすい形に書き換えたい。そのために、まず指数関数と三角関数の関係を導出するわ。

Taylor 展開¶

🟡 リナ: ここで Taylor 展開という道具を導入する。これは「関数をある点での値、1 階微分値、2 階微分値、…を使って多項式の無限和で表す方法」よ。A.4 の発散の導出で「$f(x_0 + h) \approx f(x_0) + f'(x_0)h$」という 1 次近似を使ったわね。あれは「$x_0$ の周りで関数を直線で近似する」操作だった。もっと精度を上げたければ、2 次の項 $f''(x_0)h^2/2!$ を加えて放物線で近似する。なぜ $2!$ で割るかというと、$h^2$ を $h$ で 2 回微分すると $2! = 2$ が出てくるから、それを打ち消して「$x_0$ での 2 階微分値」がそのまま係数になるようにするためよ。一般に $n$ 次の項は $f^{(n)}(x_0)h^n/n!$ ——ここで $f^{(n)}(x_0)$ は「$f$ を $n$ 回微分して $x = x_0$ を代入した値」を表す記法よ（$f^{(1)} = f'$, $f^{(2)} = f''$, $f^{(3)} = f'''$, …の一般化）。$h^n$ を $h$ で $n$ 回微分すると $n!$ が出るから、$n!$ で割っておけば係数が $f^{(n)}(x_0)$ になる。3 次、4 次……と項を増やしていけば、近似はどんどん正確になる。その極限が Taylor 展開なの。今回は $x_0 = 0$ の場合を使うわ——$x_0 = 0$ の周りの Taylor 展開を特に Maclaurin 展開とも呼ぶけれど、やっていることは同じよ。

🔵 カイ: なるほど、1 次近似が直線、2 次近似が放物線で、もっと項を増やすとどんどん元の関数に近づいていくんですね。

🟡 リナ: その通り。まず記号の説明——$n!$（「$n$ の階乗」と読む）は $1$ から $n$ までの整数を全部掛け合わせたもの：$n! = 1 \times 2 \times \cdots \times n$。例えば $2! = 1 \times 2 = 2$, $3! = 1 \times 2 \times 3 = 6$, $4! = 24$ よ。また $0! = 1$ と約束する（$n = 0$ の項が $f(0)$ になるように辻褄を合わせるため）。この記号を使うと、一般公式は：

\[ f(x) = f(0) + f'(0)\,x + \frac{f''(0)}{2!}\,x^2 + \frac{f'''(0)}{3!}\,x^3 + \cdots = \sum_{n=0}^{\infty} \frac{f^{(n)}(0)}{n!}\,x^n \]

ここで「$\cdots$」は同じパターンの項が無限に続くことを表しているわ。右辺の $\sum_{n=0}^{\infty}$ は「$n = 0, 1, 2, 3, \ldots$ と全ての非負整数について足し合わせる」という意味の記号よ（$\sum$ はギリシャ文字の大文字シグマで、「合計」を表す。高校の数列で $\sum_{k=1}^{N} a_k$ を使ったことがあるかもしれないけれど、上端が $\infty$ の場合は「項を際限なく足し続けた極限値」を意味するの）。$n$ 番目の項は $f^{(n)}(0)\,x^n / n!$ ね。

この一般公式を使って $e^x$（ネイピア数 $e \approx 2.718$ の $x$ 乗）の Taylor 展開を求めましょう。$e^x$ の最も重要な性質は「微分しても自分自身」——つまり $f(x) = e^x$ なら $f'(x) = e^x$, $f''(x) = e^x$, … と何回微分しても $e^x$ のまま、ということ（実はこの性質こそが $e$ という特別な数を定義する理由なの——「微分しても変わらない指数関数の底」として $e$ が決まるのよ）。だから $x = 0$ での値は全て $f(0) = f'(0) = f''(0) = \cdots = e^0 = 1$。一般公式に代入すると：

\[ e^x = 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \frac{x^4}{4!} + \cdots \]

⚪ メイ: 全ての微分値が 1 だから、各項の係数が $1/n!$ になるだけ——驚くほどシンプルな級数ね。

🟡 リナ: なぜこの展開が任意の $x$ で元の関数に一致するかの厳密な議論は量子力学編の付録に譲るけれど、$e^x$, $\cos x$, $\sin x$ については全ての実数 $x$ で級数が収束して元の関数に一致することが知られているの。今はその結果を認めて使うわね。

🔵 カイ: 無限に足し続けたら本当にぴったり $e^x$ になるんですか？途中で止めたら近似でしかないですよね。

🟡 リナ: いい疑問ね。途中で止めれば確かに近似——項を増やすほど精度が上がる。でも $e^x$, $\cos x$, $\sin x$ については、項を無限に足した極限が元の関数に完全に一致することが数学的に証明されているの。直感的には、$e^x$ を $x = 0$ の周りで 1 次近似すると $1 + x$、2 次近似すると $1 + x + x^2/2$——$x$ が小さいほど少ない項で良い近似になるし、$x$ が大きくても項を十分増やせば追いつく。厳密な証明は後に譲るけれど、今は「これらの関数では無限級数が収束する」という事実を認めて先に進みましょう。

この級数に $x = i\theta$ を形式的に代入してみるわ。ここで $i$ は虚数単位——$i^2 = -1$ を満たす数よ（高校の数学 II で習ったわね）。「実数の公式に虚数を入れていいの？」と思うかもしれないけれど、Taylor 展開は各項の足し算と掛け算だけでできているから、$x$ が複素数でも各項を計算して足し合わせることは問題なくできるの（収束の厳密な議論は量子力学編の付録に譲るわね）：

\[ e^{i\theta} = 1 + i\theta + \frac{(i\theta)^2}{2!} + \frac{(i\theta)^3}{3!} + \frac{(i\theta)^4}{4!} + \cdots \]

🔵 カイ: $i$ を何乗もするとどうなるんですか？

🟡 リナ: 高校の数学 II で $i^2 = -1$ は習ったわね。そこから順番に計算すると $i^3 = i^2 \cdot i = (-1) \cdot i = -i$, $i^4 = i^3 \cdot i = (-i) \cdot i = -i^2 = -(-1) = 1$, $i^5 = i^4 \cdot i = 1 \cdot i = i$, … と 4 つ周期で $i, -1, -i, 1$ を繰り返すの。つまり $i$ の偶数乗は実数（$\pm 1$）、奇数乗は純虚数（$\pm i$）になる。これを使って各項を具体的に書き下してみましょう：

$n=0$: $(i\theta)^0/0! = 1$（実数）
$n=1$: $(i\theta)^1/1! = i\theta$（純虚数）
$n=2$: $(i\theta)^2/2! = i^2\theta^2/2! = -\theta^2/2!$（実数）
$n=3$: $(i\theta)^3/3! = i^3\theta^3/3! = -i\theta^3/3!$（純虚数）
$n=4$: $(i\theta)^4/4! = i^4\theta^4/4! = +\theta^4/4!$（実数）
$n=5$: $(i\theta)^5/5! = i^5\theta^5/5! = +i\theta^5/5!$（純虚数）

パターンが見えるわね——偶数次の項は実数、奇数次の項は $i$ がかかった純虚数。だから $e^{i\theta}$ の展開を実部（$i$ を含まない項）と虚部（$i$ がかかった項）に分けられる：

\[ = \underbrace{\left(1 - \frac{\theta^2}{2!} + \frac{\theta^4}{4!} - \cdots\right)}_{\cos\theta \text{ の Taylor 展開}} + i\underbrace{\left(\theta - \frac{\theta^3}{3!} + \frac{\theta^5}{5!} - \cdots\right)}_{\sin\theta \text{ の Taylor 展開}} \]

🔵 カイ: あ、$i$ の累乗が周期的だから、偶数乗の項と奇数乗の項が自然に分かれるんですね！

⚪ メイ: 偶数乗の項が実部に、奇数乗の項が虚部に集まる形ね。

🟡 リナ: そう。$\cos\theta$ と $\sin\theta$ も同じ Taylor 展開の方法で級数に展開できる。Taylor 展開の一般公式 $f(\theta) = f(0) + f'(0)\theta + f''(0)\theta^2/2! + f'''(0)\theta^3/3! + \cdots$ に当てはめてみましょう。$\cos\theta$ の場合：$\cos 0 = 1$, $(\cos\theta)' = -\sin\theta$ だから $\theta=0$ で $0$, $(\cos\theta)'' = -\cos\theta$ だから $\theta=0$ で $-1$, $(\cos\theta)''' = \sin\theta$ だから $\theta=0$ で $0$, $(\cos\theta)^{(4)} = \cos\theta$ だから $\theta=0$ で $1$, …と 4 つ周期で $1, 0, -1, 0$ を繰り返す。代入すると $\cos\theta = 1 + 0 \cdot \theta + (-1)\theta^2/2! + 0 \cdot \theta^3/3! + 1 \cdot \theta^4/4! + \cdots$。

🔵 カイ: あ、$\cos\theta$ の微分値が $1, 0, -1, 0, \ldots$ を繰り返すのと、$e^{i\theta}$ の展開で偶数次の項だけ拾うと $1, -1/2!, +1/4!, \ldots$ になるのが対応してるんですね！だから $e^{i\theta}$ の実部が $\cos\theta$ になるんだ。

🟡 リナ: まさにその通り。偶数次の項だけ見ると、$e^{i\theta}$ の実部の係数は $1, -1/2!, +1/4!, \ldots$ で、$\cos\theta$ の展開と完全に一致する。奇数次も同様に $\sin\theta$ と一致するの。

⚪ メイ: なるほど、$i$ の累乗の周期性が、実部と虚部への分離を自動的にやってくれるのね。

🟡 リナ: そういうこと。$\sin\theta$ も同様に確認しておきましょう：$\sin 0 = 0$, $(\sin\theta)' = \cos\theta$ だから $\theta=0$ で $1$, $(\sin\theta)'' = -\sin\theta$ だから $\theta=0$ で $0$, $(\sin\theta)''' = -\cos\theta$ だから $\theta=0$ で $-1$, …。代入すると $\sin\theta = 0 + 1 \cdot \theta + 0 \cdot \theta^2/2! + (-1)\theta^3/3! + \cdots$。整理すると：

\[ \cos\theta = 1 - \frac{\theta^2}{2!} + \frac{\theta^4}{4!} - \cdots, \qquad \sin\theta = \theta - \frac{\theta^3}{3!} + \frac{\theta^5}{5!} - \cdots \]

これがまさに上の実部・虚部と一致する。こうして Euler の公式

\[ e^{i\theta} = \cos\theta + i\sin\theta \]

が得られる。

⚪ メイ: つまり、$e^{i\theta}$ の Taylor 展開を実部と虚部に分けたら、それぞれが $\cos\theta$ と $\sin\theta$ の Taylor 展開と項ごとに一致した——だから等しい、という論法ね。

🔵 カイ: すごい……指数関数と三角関数が虚数を通じてつながるって、全然違う関数に見えるのに。$\theta = \pi$ を入れたら $e^{i\pi} = \cos\pi + i\sin\pi = -1$ ってことですか？

🟡 リナ: その通り。$e^{i\pi} + 1 = 0$ ——これは Euler の等式と呼ばれる有名な式よ。さて、この公式を使うと、正弦波 $A\cos(kx - \omega t)$ は $$ u = Ae^{i(kx - \omega t)} $$

の実部として表せる。$A$ は振幅（実数の正の定数）よ。確認しておくと、Euler の公式 $e^{i\theta} = \cos\theta + i\sin\theta$ で $\theta = kx - \omega t$ と置けば $Ae^{i(kx-\omega t)} = A\cos(kx - \omega t) + iA\sin(kx - \omega t)$ だから、実部はたしかに $A\cos(kx - \omega t)$ ね。ここで複素数の幾何学的な見方を紹介しておくわ。複素数 $a + bi$ を、横軸に実部 $a$、縦軸に虚部 $b$ をとった平面（複素平面と呼ぶ）上の点 $(a, b)$ として描くことができるの。すると $e^{i\theta} = \cos\theta + i\sin\theta$ は、実部が $\cos\theta$、虚部が $\sin\theta$ だから、点 $(\cos\theta, \sin\theta)$ に対応する——これは原点から距離 1（$\cos^2\theta + \sin^2\theta = 1$ だから）、横軸からの角度が $\theta$ の位置にある点よ。つまり $e^{i\theta}$ は「単位円上を角度 $\theta$ だけ回転した点」を表しているの。以降、波動方程式の計算ではこの複素表記を使うわ。

🔵 カイ: なんでわざわざ複素数を使うんですか？ $\cos$ のままじゃダメなんですか？

🟡 リナ: 2 つ理由があるわ。まず実用面——$e^{i\theta}$ を微分すると $ie^{i\theta}$ が出るだけで、$\cos$ と $\sin$ を行ったり来たりしなくて済む。次に正当性——波動方程式は線形（解を定数倍したり足し合わせたりしても解になる）なの。「線形」とはもう少し具体的に言うと、方程式の各項が $u$ やその偏微分の 1 次式だけでできていて、$u^2$ や $u \cdot \partial_t u$ のような非線形項がないということ。このとき、$u = u_1 + iu_2$（$u_1$, $u_2$ は実数値関数）を方程式に代入すると、各項の偏微分は実部と虚部に分離できる——例えば $\partial_t^2(u_1 + iu_2) = \partial_t^2 u_1 + i\partial_t^2 u_2$ のように。すると方程式全体が「（$u_1$ についての方程式）$+ i \times$（$u_2$ についての方程式）$= 0$」の形になる。ここで、複素数 $a + bi = 0$ が成り立つには $a = 0$ かつ $b = 0$ でなければならない（実数と虚数は「独立な方向」だから、片方だけでは打ち消せない）。だから $u_1$ と $u_2$ はそれぞれ独立に方程式を満たす。だから「複素数で計算して最後に実部をとる」という手順が正当化される。量子力学では波動関数そのものが複素数値になるから、この表記は必須になるわ。

🔵 カイ: 線形だから実部と虚部を分離できる——計算の道具として複素数を使って、最後に実部をとれば物理的な答えが出る、ということですね。でも、もし方程式が非線形だったら——例えば $u^2$ みたいな項があったら、この手は使えないんですか？

🟡 リナ: その通り、使えない。非線形の場合は $(u_1 + iu_2)^2 = u_1^2 - u_2^2 + 2iu_1 u_2$ のように実部と虚部が混ざってしまって、分離できなくなるの。幸い、波動方程式は線形だから安心して使えるわ。

⚪ メイ: 線形性が「複素数を道具として使える」ことの保証になっている——ここは大事なポイントね。

🟡 リナ: では実際にこれを波動方程式 $\partial_t^2 u = c^2 \partial_x^2 u$ に代入してみましょう。さっき確認した「$e^x$ は微分しても形が変わらない」性質を使うわ。$e^{ax}$ を $x$ で微分すると連鎖律から $ae^{ax}$ になる——$a$ が複素数でも同じ。つまり $e^{i\theta}$ を $\theta$ で微分すると $i\,e^{i\theta}$——微分するたびに $i$ が 1 つ前に出てくるだけで、関数の形は変わらないの。

具体的に $u = Ae^{i(kx - \omega t)}$ を $t$ で偏微分すると、$x$ は定数扱いだから $e^{ikx}$ の部分はそのまま残って、$e^{-i\omega t}$ の微分で $-i\omega$ が前に出てくる。つまり $\partial u/\partial t = -i\omega\, u$。もう一度 $t$ で微分すると $(-i\omega)^2 u = -\omega^2 u$ がかかるから、左辺は $-\omega^2 u$。同様に $x$ で 2 回微分すると $(ik)^2 = -k^2$ がかかるから、右辺は $-c^2 k^2 u$。

🔵 カイ: あ、微分するたびに $-i\omega$ とか $ik$ が前に出てくるだけで、$e$ の形はそのまま残るんですね。だから微分の計算が、ただの掛け算に化けるのか——便利だ。

🟡 リナ: まとめると、波動方程式 $\partial_t^2 u = c^2 \partial_x^2 u$ に代入した結果は $-\omega^2 u = -c^2 k^2 u$。

🔵 カイ: $u \neq 0$ で割れば $\omega^2 = c^2 k^2$ ですね。これって $\omega$ と $k$ にいろんな組み合わせがありえるってことですか？

🟡 リナ: そう。$\omega^2 = c^2 k^2$ の平方根をとると $|\omega| = c|k|$ ね（$c > 0$ だから）。ここで物理の慣習として $\omega > 0$（振動数は正）をとるの。すると $\omega = c|k|$ となり、波の進行方向は $k$ の符号が担うことになる。

🔵 カイ: ちょっと待ってください。さっき $k = 2\pi/\lambda$ って定義しましたよね。$\lambda$ は正だから $k$ も正のはずじゃ……？

🟡 リナ: いい指摘ね。$k = 2\pi/\lambda$ は波数の大きさの定義。でも波の進行方向も含めて表したいとき、$k$ に符号をつけるの——$k > 0$ なら $+x$ 方向に進む波、$k < 0$ なら $-x$ 方向に進む波。大きさは $|k| = 2\pi/\lambda$ で変わらないわ。つまり $k$ は「1 m あたりの位相変化」に進行方向の情報を載せた量なの。

$k > 0$ のとき $\omega = ck$ だから、$e^{i(kx - \omega t)} = e^{ik(x - ct)}$ ——位相 $k(x - ct)$ が一定の点を追うと、$t$ が増えるにつれ $x$ も増えなければならない。つまり波の山は $+x$ 方向（右）に速さ $c$ で動く。逆に $k < 0$ なら $\omega = c(-k) = -ck$ だから $e^{i(kx - \omega t)} = e^{ik(x + ct)}$ ——位相 $k(x + ct)$ が一定の点を追うと、$k < 0$ なので $x + ct$ が減る方向、つまり $x$ が減る方向（左）に進む波になるの。

⚪ メイ: d'Alembert の一般解の $f(x - ct)$ と $g(x + ct)$ が、$k$ の符号に対応しているのね。

🟡 リナ: まとめると、$\omega > 0$ の慣習のもとで $\omega = c|k|$ が成り立ち、$k$ の符号が進行方向を決める。このように $\omega$ と $k$ の関係を表す式を分散関係と呼ぶわ。

🔵 カイ: あれ、$k > 0$ のとき $\omega = c|k| = ck$ だから $e^{i(kx - \omega t)} = e^{ik(x - ct)}$ って書けますよね。これって d'Alembert の $f(x - ct)$ の形じゃないですか？

🟡 リナ: いい気づきね。同様に $k < 0$ なら $e^{ik(x + ct)}$ となって左進行波 $g(x + ct)$ に対応する。d'Alembert の一般解の $f$ と $g$ を正弦波で具体化したのが、この平面波解なのよ。

⚪ メイ: つまり、一般解の中の「任意関数」を正弦波に限定したのが平面波解ということね。

🟡 リナ: その通り。ここで波の山が移動する速さ（と方向）を位相速度と呼び、$v_p = \omega/k$ で定義するの。なぜこれが「波の山の速さ」になるか確認しておくわ。波の山は位相 $kx - \omega t = \text{const}$ の位置にある。この条件を $t$ で微分すると $k(dx/dt) - \omega = 0$、つまり $dx/dt = \omega/k$——これが山の移動速度よ。$k > 0$ なら $\omega = ck$ だから $v_p = \omega/k = ck/k = c$（右に進む）。$k < 0$ なら $\omega = c|k| = c(-k) = -ck$ だから $v_p = \omega/k = -ck/k = -c$（左に進む）。どちらも速さの大きさは $|v_p| = c$ で同じよ。高校で習った $v = \lambda\nu$ との対応も確認しておくと、$\lambda = 2\pi/|k|$, $\nu = \omega/(2\pi)$ だから $\lambda\nu = \omega/|k| = c$ で一致するわ。つまりこの波動方程式では、位相速度の大きさがどんな波数 $k$ でも同じ値 $c$ になる——波の形が崩れずに伝わるの。こういう場合を「分散がない」と言うわ。

⚪ メイ: つまり、$\omega = c|k|$ でどの $k$ でも位相速度の大きさ $|\omega/k| = c$ が同じだから、全部の波が同じ速さで進む——だから波の形が崩れない、ということね。

🔵 カイ: 逆に「分散がある」場合ってどうなるんですか？そもそも「分散」って何が散らばるんですか？

🟡 リナ: 波数ごとに速度が違うと、最初は一緒に重なっていた波の成分がバラバラに散らばっていく——だから「分散」と呼ぶの。例えば水面の波では $\omega \propto \sqrt{k}$ のような関係になる。位相速度は $v_p = \omega/k \propto \sqrt{k}/k = 1/\sqrt{k}$ だから、波数が大きい（波長が短い）波ほど遅く進む。すると最初は綺麗な形だった波が時間とともに崩れていく。真空中の光や重力波は $\omega = c|k|$ で分散がないから、波形を保ったまま伝わるの。

🔵 カイ: なるほど、プリズムで白色光が虹色に分かれるのも、ガラスの中で波長ごとに速度が違う——つまり分散があるからですね。じゃあ重力波も真空中では分散がないから、波形がそのまま保たれて地球に届く——だから LIGO で波形を読み取れるんだ。

🟡 リナ: その通り。もし重力波に分散があったら、遠くの天体から届く間に波形が崩れてしまって、マッチドフィルタリングが使えなくなるわ。分散がないことは、重力波天文学の前提条件なのよ。

本編でのつながり¶

🟡 リナ: この波動方程式の分類が本編のどこで効いていたか、振り返っておきましょう。

第 1 章（Poisson vs 波動）: Newton 重力は時間微分を含まない Poisson 方程式（楕円型）——瞬時伝播。電磁気は波動方程式（双曲型）で光速 $c$ で伝わる。第 1 章ではこの違いを特殊相対論の「情報は光速を超えない」という要請と対比させ、一般相対論の動機とした。ここでの分類（楕円型 = 瞬時伝播 vs 双曲型 = 有限速度伝播）がその議論の数学的背景になっている。
第 19 章（重力波）: 線形化された Einstein 方程式 $\Box \bar{h}_{\mu\nu} = -(16\pi G/c^4)T_{\mu\nu}$ が波動方程式の形を取る。時空の曲率の微小な擾乱が光速で伝播する——重力波の理論的予言。
弦理論編（後の編で扱う）: 弦の振動が 2 次元の波動方程式に従い、その一般解を左進行波・右進行波に分解して量子化する。

🔵 カイ: 全部「波動方程式の形かどうか」で判定できるんですね。この分類を知ってると、新しい方程式に出会ったときに解の性質がすぐ予想できそうだ。

✅ 理解度チェック: 偏微分方程式に時間の 2 階微分が含まれる場合、それはどのタイプに分類されるでしょうか？

答え

波動方程式（双曲型）に分類され、解は振動する（波として伝播する）。

✅ 理解度チェック: 1 次元波動方程式の一般解 $u(x, t) = f(x - ct) + g(x + ct)$ において、$c$ は何を表すでしょうか？

答え

波の伝播速度。$f(x - ct)$ は右に進む波、$g(x + ct)$ は左に進む波。

A.8　まとめの表¶

🟡 リナ: 最後に、この付録の内容を一覧にまとめておくわ。

表 A.2: 付録Aの演算・概念まとめ一覧

演算・概念	記号	入力 → 出力	物理的意味
偏微分	$\partial f/\partial x$	多変数関数 → 多変数関数	他変数固定での変化率
勾配	$\nabla\varphi$	スカラー → ベクトル	最急上昇方向と変化率
発散	$\nabla \cdot \boldsymbol{F}$	ベクトル → スカラー	湧き出しの強さ
回転	$\nabla \times \boldsymbol{F}$	ベクトル → ベクトル	渦の強さと方向
ラプラシアン	$\nabla^2\varphi$	スカラー → スカラー	勾配の発散（ポテンシャルの曲がり具合）
波動方程式	$\partial_t^2 f = v^2 \nabla^2 f$	スカラー/ベクトル場が満たす方程式	速度 $v$ での波の伝播

次章予告¶

次の Appendix B では、テンソル積と Einstein の縮約記法を導入する。ベクトルの「添字を増やす」操作であるテンソル積を理解し、一般相対論の計算で不可欠な縮約記法を自在に使いこなせるようにする。

練習問題¶

📝 練習問題:

ベクトル積の成分計算 → 問題 B-1. ベクトル積の成分計算

直交性の確認 → 問題 B-2. ベクトル積の直交性

スカラー三重積の行列式計算 → 問題 B-3. スカラー三重積

Lagrange の恒等式の数値確認 → 問題 B-8. Lagrange 恒等式の確認

BAC-CAB 公式の数値確認 → 問題 B-9. BAC-CAB 公式の確認

BAC-CAB 公式の一般証明 → 問題 M-1. BAC-CAB 公式の証明

ベクトル積どうしの内積 → 問題 M-3. ベクトル積どうしの内積の公式

Levi-Civita 記号による統一的導出 → 問題 A-1. Levi-Civita 記号による恒等式

勾配の計算 → 問題 B-4. スカラー場の勾配

発散の計算 → 問題 B-5. ベクトル場の発散

回転の計算 → 問題 B-6. ベクトル場の回転

ラプラシアンの計算 → 問題 B-7. スカラー場のラプラシアン

rot(grad φ) = 0 の確認 → 問題 B-10. rot(grad) = 0 の確認

div(rot F) = 0 の証明 → 問題 M-2. div(rot) = 0 の証明

積の発散の公式 → 問題 M-4. 積の発散の公式

Gauss の定理の応用（クーロン電場） → 問題 M-5. Gauss の定理とクーロン電場

Stokes の定理と保存場・曲率の関係 → 問題 A-2. Stokes の定理の応用

参考文献¶

石井俊全『一般相対性理論を一歩一歩数式で理解する』ベレ出版, 第 1 章「数学の準備」
太田浩一『電磁気学の基礎 I』東京大学出版会（Maxwell 方程式とベクトル解析の関係について）
H. M. Schey, Div, Grad, Curl, and All That, W. W. Norton（直感的な理解のための入門書）