体積要素

数学において、体積要素（たいせきようそ、英: volume element）とは、関数を球面座標系や円柱座標系など様々な座標系において体積について積分する際に現われる概念である。次の式により表現される: d V := ρ ( u 1 , u 2 , u 3 ) d u 1 d u 2 d u 3 . {\displaystyle dV:=\rho (u_{1},u_{2},u_{3})\,du_{1}\,du_{2}\,du_{3}.} ここで、ui は座標であり、任意の集合 B の体積を次のように計算できるものとする: Volume ⁡ ( B ) := ∫ B ρ ( u 1 , u 2 , u 3 ) d u 1 d u 2 d u 3 . {\displaystyle \operatorname {Volume} (B):=\int _{B}\rho (u_{1},u_{2},u_{3})\,du_{1}\,du_{2}\,du_{3}.} たとえば、球面座標系においてはdV = u12 sin u2 du1 du2 du3 であり、従って dV = u12 sin u2 である。 体積要素という概念は三次元に留まるものではない。二次元では面積要素（めんせきようそ、area element）と呼ばれることも多く、面積分を行う際に有用である。座標変換の際、（変数変換公式により）体積要素は座標変換のヤコビ行列の行列式の絶対値だけ変化する。この事実から、体積要素は多様体の一種の測度として定義できることが従う。向き付け可能な可微分多様体においては、典型的には体積要素は体積形式、すなわち最高次の微分形式から導かれる。向き付け不可能な多様体においては、典型的には体積要素は（局所的に定義される）体積要素の絶対値であり、1-密度（英語版）を定義する。 ユークリッド空間においては、体積要素はデカルト座標に沿った微分の積により与えられる。 d V = d x d y d z {\displaystyle \mathrm {d} V=\mathrm {d} x\,\mathrm {d} y\,\mathrm {d} z} 他の座標系においては、x = x(u1, u2, u3), y = y(u1, u2, u3), z = z(u1, u2, u3) とするとヤコビ行列を用いて体積要素を以下のように計算できる。 d V = | ∂ ( x , y , z ) ∂ ( u 1 , u 2 , u 3 ) | d u 1 d u 2 d u 3 {\displaystyle \mathrm {d} V=\left|{\frac {\partial (x,y,z)}{\partial (u_{1},u_{2},u_{3})}}\right|\,\mathrm {d} u_{1}\,\mathrm {d} u_{2}\,\mathrm {d} u_{3}} たとえば、球面座標系では x = ρ cos ⁡ θ sin ⁡ ϕ y = ρ sin ⁡ θ sin ⁡ ϕ z = ρ cos ⁡ ϕ {\displaystyle {\begin{aligned}x&=\rho \cos \theta \sin \phi \\y&=\rho \sin \theta \sin \phi \\z&=\rho \cos \phi \end{aligned}}} であるからヤコビアンは | ∂ ( x , y , z ) ∂ ( ρ , θ , ϕ ) | = ρ 2 sin ⁡ ϕ {\displaystyle \left|{\frac {\partial (x,y,z)}{\partial (\rho ,\theta ,\phi )}}\right|=\rho ^{2}\sin \phi } となり、したがって体積要素は以下のように書ける。 d V = ρ 2 sin ⁡ ϕ d ρ d θ d ϕ . {\displaystyle \mathrm {d} V=\rho ^{2}\sin \phi \,\mathrm {d} \rho \,\mathrm {d} \theta \,\mathrm {d} \phi .} このことは微分形式が引き戻し（英語版） F* により以下のように変換することの例と見ることができる。 F ∗ ( u d y 1 ∧ ⋯ ∧ d y n ) = ( u ∘ F ) det ( ∂ F j ∂ x i ) d x 1 ∧ ⋯ ∧ d x n {\displaystyle F^{*}(u\;\mathrm {d} y^{1}\wedge \cdots \wedge \mathrm {d} y^{n})=(u\circ F)\det \left({\frac {\partial F^{j}}{\partial x^{i}}}\right)\mathrm {d} x^{1}\wedge \cdots \wedge \mathrm {d} x^{n}} n-次元ユークリッド空間 Rn の線形部分空間が次の線形独立なベクトルにより張られるものとする。 X 1 , … , X k {\displaystyle X_{1},\dots ,X_{k}} この部分空間における体積要素を計算する場合、Xi の張る平行多胞体[訳語疑問点]が線形幾何学から Xi のグラム行列の行列式の平方根により与えられることを知っておくと便利である。 det ( X i ⋅ X j ) i , j = 1 … k {\displaystyle {\sqrt {\det(X_{i}\cdot X_{j})_{i,j=1\dots k}}}} この部分空間上の任意の点 p はある座標 (u1, u2, ..., uk) により以下のように表わされる。 p = u 1 X 1 + ⋯ + u k X k {\displaystyle p=u_{1}X_{1}+\cdots +u_{k}X_{k}} 点 p において、辺を dui とする微小平行多胞体を作ると、その体積はグラム行列の行列式の平方根により与えられる。 det ( ( d u i X i ) ⋅ ( d u j X j ) ) i , j = 1 … k = det ( X i ⋅ X j ) i , j = 1 … k d u 1 d u 2 ⋯ d u k {\displaystyle {\sqrt {\det \left((\mathrm {d} u_{i}X_{i})\cdot (\mathrm {d} u_{j}X_{j})\right)_{i,j=1\dots k}}}={\sqrt {\det(X_{i}\cdot X_{j})_{i,j=1\dots k}}}\;\mathrm {d} u_{1}\,\mathrm {d} u_{2}\,\cdots \,\mathrm {d} u_{k}} これにより線形部分空間における体積形式を定義することができる。 向き付け可能な次元 n のリーマン多様体における体積要素は定数関数の f(x) = 1 のホッジ双対に等しい。 ω = ⋆ 1 {\displaystyle \omega =\star 1} これと等価に、体積要素は正確にレヴィ＝チヴィタテンソル ε と正確に一致する。座標を用いて書けば、以下のようになる。 ω = ϵ = | det g | d x 1 ∧ ⋯ ∧ d x n {\displaystyle \omega =\epsilon ={\sqrt {|\det g|}}\,\mathrm {d} x^{1}\wedge \cdots \wedge \mathrm {d} x^{n}} ここで det g はその座標系における計量テンソル g の行列式である。 n-次元ユークリッド空間に埋め込まれた二次元曲面を考えることで、体積要素の単純な例を考察することができる。この場合、体積要素は面積要素と呼ばれることもある。部分集合 U ⊂ R2 と写像 φ : U → R n {\displaystyle \varphi :U\to \mathbf {R} ^{n}} を考えることにより、Rn に埋め込まれた曲面を定義する。二次元では体積は面積であり、体積要素は曲面の任意の部分の面積を決定する方法を与える。したがって、体積要素は次の形式をとる。 f ( u 1 , u 2 ) d u 1 d u 2 {\displaystyle f(u_{1},u_{2})\,\mathrm {d} u_{1}\,\mathrm {d} u_{2}} これにより曲面上の集合 B の面積を積分を用いて以下のように計算できる。 Area ⁡ ( B ) = ∫ B f ( u 1 , u 2 ) d u 1 d u 2 {\displaystyle \operatorname {Area} (B)=\int _{B}f(u_{1},u_{2})\,\mathrm {d} u_{1}\,\mathrm {d} u_{2}} ここで、通常の意味での面積を与えるような体積要素を決定したい。写像のヤコブ行列は以下のように書ける。 λ i j = ∂ φ i ∂ u j {\displaystyle \lambda _{ij}={\frac {\partial \varphi _{i}}{\partial u_{j}}}} ここで添字 i は 1 から n を、j は 1 から 2 を走る。n-次元空間のユークリッド計量から、集合 U 上の計量 g = λTλ を計算でき、その行列要素は以下のように与えられる。 g i j = ∑ k = 1 n λ k i λ k j = ∑ k = 1 n ∂ φ k ∂ u i ∂ φ k ∂ u j {\displaystyle g_{ij}=\sum _{k=1}^{n}\lambda _{ki}\lambda _{kj}=\sum _{k=1}^{n}{\frac {\partial \varphi _{k}}{\partial u_{i}}}{\frac {\partial \varphi _{k}}{\partial u_{j}}}} この計量の行列式は次のようになる。 det g = | ∂ φ ∂ u 1 ∧ ∂ φ ∂ u 2 | 2 = det ( λ T λ ) {\displaystyle \det g=\left|{\frac {\partial \varphi }{\partial u_{1}}}\wedge {\frac {\partial \varphi }{\partial u_{2}}}\right|^{2}=\det(\lambda ^{\mathrm {T} }\lambda )} 正則曲面においては、この行列式はいたるところ非零、すなわちヤコブ行列のランクはいたるところで 2 である。 ここで、U 上の座標を微分同相写像 により変換し、座標 (u1, u2) が (v1, v2) に移るものとする。すなわち、 (u1, u2) = f(v1, v2) となる。この変換のヤコブ行列は次のように与えられる。 F i j = ∂ f i ∂ v j {\displaystyle F_{ij}={\frac {\partial f_{i}}{\partial v_{j}}}} この新しい座標では、 ∂ φ i ∂ v j = ∑ k = 1 2 ∂ φ i ∂ u k ∂ f k ∂ v j {\displaystyle {\frac {\partial \varphi _{i}}{\partial v_{j}}}=\sum _{k=1}^{2}{\frac {\partial \varphi _{i}}{\partial u_{k}}}{\frac {\partial f_{k}}{\partial v_{j}}}} となり、計量は以下のように変換される。 g ~ = F T g F {\displaystyle {\tilde {g}}=F^{\mathrm {T} }gF} ここで、 ~g は v 座標系に引き戻した計量である。この行列式は次のようになる。 det g ~ = det g ( det F ) 2 {\displaystyle \det {\tilde {g}}=\det g(\det F)^{2}} 以上から、体積要素が向きを保存する座標変換の下に不変であることがわかる。 二次元においては、体積は面積である。部分集合 B ⊂ U の面積は次のように積分で得られる。 Area ( B ) = ∬ B det g d u 1 d u 2 = ∬ B det g | det F | d v 1 d v 2 = ∬ B det g ~ d v 1 d v 2 {\displaystyle {\begin{aligned}{\mbox{Area}}(B)&=\iint _{B}{\sqrt {\det g}}\;\mathrm {d} u_{1}\;\mathrm {d} u_{2}\\&=\iint _{B}{\sqrt {\det g}}\;|\det F|\;\mathrm {d} v_{1}\;\mathrm {d} v_{2}\\&=\iint _{B}{\sqrt {\det {\tilde {g}}}}\;\mathrm {d} v_{1}\;\mathrm {d} v_{2}\end{aligned}}} したがって、どちらの座標系においても体積要素は同一の形式を持ち、すなわち座標変換で不変に保たれる。 以上の議論に二次元特有のものは一切ないので、任意の次元について成り立つ。 例として、半径 r の中心を原点とする R3 上の球を考える。この球は球面座標系を用いて次の写像で表現される。 ϕ ( u 1 , u 2 ) = ( r cos ⁡ u 1 sin ⁡ u 2 , r sin ⁡ u 1 sin ⁡ u 2 , r cos ⁡ u 2 ) {\displaystyle \phi (u_{1},u_{2})=(r\cos u_{1}\sin u_{2},r\sin u_{1}\sin u_{2},r\cos u_{2})} したがって、 g = ( r 2 sin 2 ⁡ u 2 0 0 r 2 ) {\displaystyle g={\begin{pmatrix}r^{2}\sin ^{2}u_{2}&0\\0&r^{2}\end{pmatrix}}} となり、面積要素は次のように得られる。 ω = det g d u 1 d u 2 = r 2 sin ⁡ u 2 d u 1 d u 2 {\displaystyle \omega ={\sqrt {\det g}}\;\mathrm {d} u_{1}\mathrm {d} u_{2}=r^{2}\sin u_{2}\,\mathrm {d} u_{1}\mathrm {d} u_{2}} 円柱座標変換#積分 面積分 体積分 ^ Carroll, Sean. Spacetime and Geometry. Addison Wesley, 2004, p. 90 Besse, Arthur L. (1987), Einstein manifolds, Ergebnisse der Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)], vol. 10, Berlin, New York: Springer-Verlag, pp. xii+510, ISBN 978-3-540-15279-8