ラプラス作用素

数学におけるラプラス作用素（ラプラスさようそ、英: Laplace operator）あるいはラプラシアン（英: Laplacian）は、ユークリッド空間上の函数の勾配の発散として与えられる微分作用素である。記号では ∇·∇, ∇2, あるいは ∆ で表されるのが普通である。函数 f の点 p におけるラプラシアン ∆f(p) は（次元に依存する定数の違いを除いて）点 p を中心とする球面を半径が増大するように動かすときの f(p) から得られる平均値になっている。直交座標系においては、ラプラシアンは各独立変数に関する函数の二階（非混合）偏導函数の和として与えられ、またほかに円筒座標系や球座標系などの座標系においても有用な表示を持つ。 ラプラス作用素の名称は、天体力学の研究に同作用素を最初に用いたフランス人数学者のピエール＝シモン・ド・ラプラス (1749–1827) に因んでいる。同作用素は与えられた重力ポテンシャルに適用すると質量密度の定数倍を与える。現在ではラプラス方程式と呼ばれる方程式 ∆f = 0 の解は調和函数と呼ばれ、自由空間において可能な重力場を表現するものである。 微分方程式においてラプラス作用素は電気ポテンシャル、重力ポテンシャル、熱や流体の拡散方程式、波の伝搬、量子力学といった、多くの物理現象を記述するのに現れる。ラプラシアンは、函数の勾配フローの流束密度を表す。 ラプラス作用素はn 次元ユークリッド空間上の函数 f の勾配 ∇f の発散 ∇· として定義される二階の微分作用素である。つまり、f が二回微分可能実数値函数ならば f のラプラシアンは Δ f ≡ ∇ 2 f := ∇ ⋅ ∇ f {\displaystyle \Delta f\equiv \nabla ^{2}f:=\nabla \cdot \nabla f} (1) で定義される。ただし、あとの記法は形式的に ∇ = ( ∂ / ∂ x 1 , … , ∂ / ∂ x n ) {\displaystyle \nabla =(\partial /\partial x_{1},\dotsc ,\partial /\partial x_{n})} と書いたものである。あるいは同じことだが、f のラプラシアンは直交座標系 xi における非混合二階偏導函数の全てにわたる和 Δ f = ∑ i = 1 n ∂ 2 f ∂ x i 2 {\displaystyle \Delta f=\sum _{i=1}^{n}{\frac {\partial ^{2}f}{\partial x_{i}^{2}}}} (2) としても書ける。二階の微分作用素として、ラプラス作用素はCk 級函数を Ck − 2 級の函数へ写す (k ≥ 2)。つまり、式 1 (あるいは同値な 2) は作用素 ∆: Ck(Rn) → Ck − 2(Rn) を定める。あるいはより一般に任意の開集合 Ω に対して作用素 ∆: Ck(Ω) → Ck − 2(Ω) を定める。 ラプラス作用素は、合同変換と可換である。すなわち、任意のC∞級関数φ : Rn → Rと任意の合同変換Tに対し、 Δ ( φ ( T ( x ) ) ) = T ( Δ ( φ ( x ) ) ) {\displaystyle \Delta (\varphi (T(x)))=T(\Delta (\varphi (x)))} が成立する。 しかもラプラス作用素は、上記の性質を満たす非自明な微分演算子で最も簡単なものとして特徴づけることができる。これを説明する為、記号を導入する。Rを実数の集合とし、n個の実数からなる組の集合をRnとする。x=(x1,…,xn)∈Rnとn個の非負整数の組α=(α1,…,αn)に対し、 ∂ ∂ x α := ∂ n ∂ x 1 α 1 ⋯ ∂ x n α n , {\displaystyle {\partial \over \partial x^{\alpha }}:={\partial ^{n} \over \partial x_{1}{}^{\alpha _{1}}\cdots \partial x_{n}{}^{\alpha _{n}}},} | α | := α 1 + ⋯ + α n {\displaystyle |\alpha |:=\alpha _{1}+\cdots +\alpha _{n}} と表記する。微分演算子 D := ∑ α : | α | ≤ k a α ∂ ∂ x α {\displaystyle D:=\sum _{\alpha :|\alpha |\leq k}a_{\alpha }{\partial \over \partial x^{\alpha }}} が任意のC∞級関数φ : Rn → Rと向きを保つ任意の合同変換Tに対し、 D ( φ ( T ( x ) ) ) = T ( D ( φ ( x ) ) ) {\displaystyle D(\varphi (T(x)))=T(D(\varphi (x)))} が成立していたとする。このとき、実数係数の1変数多項式p(X)=Σm umXmが存在し、 D = p ( Δ ) = ∑ m u m Δ m {\displaystyle D=p(\Delta )=\sum _{m}u_{m}\Delta ^{m}} が成立する。 よってラプラス作用素は、合同変換に対して不変な微分演算子の中で、自明なもの（＝恒等的に0を対応させる微分演算子）を除けば最も簡単なものである。 拡散の物理理論において、ラプラス作用素は（ラプラス方程式を通じて）平衡の数学的記述に自然に現れる。具体的に、u が化学濃度のような適当な量の平衡密度であるとき、u の滑らかな境界を持つ領域 V を通る流束が、V に流入も漏出も無いとすれば、0 であるから ∫ ∂ V ∇ u ⋅ n d S = 0 {\displaystyle \int _{\partial V}\nabla u\cdot {\boldsymbol {n}}\,dS=0} と書ける。ただし、 n {\textstyle {\boldsymbol {n}}} は領域 V の境界に対して外側を向く単位法ベクトルである。発散定理により ∫ V div ⁡ ∇ u d V = ∫ ∂ V ∇ u ⋅ n d S = 0 {\displaystyle \int _{V}\operatorname {div} \nabla u\,dV=\int _{\partial V}\nabla u\cdot {\boldsymbol {n}}\,dS=0} は領域 V が滑らかな境界を持つ限りにおいて成り立つから、これにより div ⁡ ∇ u = Δ u = 0 {\displaystyle \operatorname {div} \nabla u=\Delta u=0} が導かれる。方程式の左辺はラプラス作用素である。ラプラス作用素それ自身は拡散方程式によって記述されるような、化学濃度の流入や漏出を表す点を含む非平衡拡散に対する物理的解釈を持つ。 φ が電荷分布 q に付随した電位を記述するものとすると、電荷分布自身は φ のラプラシアンとして q = Δ φ {\displaystyle q=\Delta \varphi } (1) で与えられる。これはガウスの法則の帰結である。実際、V が任意の滑らかな領域ならば、電場 E {\textstyle {\boldsymbol {E}}} の電束に関するガウスの法則により、（単位当たりの）電荷は ∫ ∂ V E ⋅ n d S = ∫ ∂ V ∇ φ ⋅ n d S = ∫ V q d V {\displaystyle \int _{\partial V}{\boldsymbol {E}}\cdot {\boldsymbol {n}}\,dS=\int _{\partial V}\nabla \varphi \cdot {\boldsymbol {n}}\,dS=\int _{V}q\,dV} になる。ただし、最初の等号は静電場は静電位の勾配に等しいという事実を用いた。発散定理により、 ∫ V Δ φ d V = ∫ V q d V {\displaystyle \int _{V}\Delta \varphi \,dV=\int _{V}q\,dV} が成り立ち、これは任意の領域 V に対して成り立つことから (1) を得る。 同じ説明によって、重力ポテンシャルのラプラシアンが質量分布となることが導かれる。電荷や質量の分布が与えられていてそれらに付随するポテンシャルは未知ということはよくあることである。適当な境界条件の下でポテンシャル函数を求めるということは、ポワソン方程式を解くことに同じである。 物理学においてラプラス作用素が現れる別な理由は、領域 U における方程式 ∆f = 0 の解はディリクレエネルギー汎函数を停留させる函数 E ( f ) := 1 2 ∫ U ‖ ∇ f ‖ 2 d x {\displaystyle E(f):={\frac {1}{2}}\int _{U}\Vert \nabla f\Vert ^{2}\,dx} となることである。これを見るために f: U → R は函数で、函数 u: U → R は U の境界上で消えていると仮定する。このとき d d ε | ε = 0 E ( f + ε u ) = ∫ U ∇ f ⋅ ∇ u d x = − ∫ U u Δ f d x {\displaystyle {\frac {d}{d\varepsilon }}{\Bigg |}_{\varepsilon =0}E(f+\varepsilon u)=\int _{U}\nabla f\cdot \nabla u\,dx=-\int _{U}u\Delta f\,dx} が成り立つ（ただし、最後の等号はグリーンの第一恒等式（英語版）を用いた）。この計算により、∆f = 0 ならば E は f の周りで停留する。逆に E が f の周りで停留するならば変分法の基本補題（英語版） により ∆f = 0 である。 二次元のラプラス作用素は x, y を xy-平面上の標準直交座標として Δ f := ∂ 2 f ∂ x 2 + ∂ 2 f ∂ y 2 {\displaystyle \Delta f:={\frac {\partial ^{2}f}{\partial x^{2}}}+{\frac {\partial ^{2}f}{\partial y^{2}}}} で与えられる。 極座標 Δ f = 1 r ∂ ∂ r ( r ∂ f ∂ r ) + 1 r 2 ∂ 2 f ∂ θ 2 = 1 r ∂ f ∂ r + ∂ 2 f ∂ r 2 + 1 r 2 ∂ 2 f ∂ θ 2 . {\displaystyle {\begin{aligned}\Delta f&={1 \over r}{\partial \over \partial r}\left(r{\partial f \over \partial r}\right)+{1 \over r^{2}}{\partial ^{2}f \over \partial \theta ^{2}}\\&={1 \over r}{\partial f \over \partial r}+{\partial ^{2}f \over \partial r^{2}}+{1 \over r^{2}}{\partial ^{2}f \over \partial \theta ^{2}}.\end{aligned}}} 三次元では様々な座標系がラプラシアンを記述するために広く用いられる。 直交座標系 Δ f = ∂ 2 f ∂ x 2 + ∂ 2 f ∂ y 2 + ∂ 2 f ∂ z 2 . {\displaystyle \Delta f={\frac {\partial ^{2}f}{\partial x^{2}}}+{\frac {\partial ^{2}f}{\partial y^{2}}}+{\frac {\partial ^{2}f}{\partial z^{2}}}.} 円筒座標系 Δ f = 1 ρ ∂ ∂ ρ ( ρ ∂ f ∂ ρ ) + 1 ρ 2 ∂ 2 f ∂ φ 2 + ∂ 2 f ∂ z 2 . {\displaystyle \Delta f={1 \over \rho }{\partial \over \partial \rho }\left(\rho {\partial f \over \partial \rho }\right)+{1 \over \rho ^{2}}{\partial ^{2}f \over \partial \varphi ^{2}}+{\partial ^{2}f \over \partial z^{2}}.} 球面座標系 Δ f = 1 r 2 ∂ ∂ r ( r 2 ∂ f ∂ r ) + 1 r 2 sin ⁡ θ ∂ ∂ θ ( sin ⁡ θ ∂ f ∂ θ ) + 1 r 2 sin 2 ⁡ θ ∂ 2 f ∂ φ 2 . {\displaystyle \Delta f={1 \over r^{2}}{\partial \over \partial r}\left(r^{2}{\partial f \over \partial r}\right)+{1 \over r^{2}\sin \theta }{\partial \over \partial \theta }\left(\sin \theta {\partial f \over \partial \theta }\right)+{1 \over r^{2}\sin ^{2}\theta }{\partial ^{2}f \over \partial \varphi ^{2}}.} 一般の曲線座標系（英語版） ( ξ 1 , ξ 2 , ξ 3 ) {\displaystyle (\xi ^{1},\xi ^{2},\xi ^{3})} ∇ 2 = ∇ ξ m ⋅ ∇ ξ n ∂ 2 ∂ ξ m ∂ ξ n + ∇ 2 ξ m ∂ ∂ ξ m , {\displaystyle \nabla ^{2}=\nabla \xi ^{m}\cdot \nabla \xi ^{n}{\partial ^{2} \over \partial \xi ^{m}\partial \xi ^{n}}+\nabla ^{2}\xi ^{m}{\partial \over \partial \xi ^{m}},} ここでアインシュタインの和の規約を用いた。 N 次元球座標系において、r を正の実数をとる半径、θ は単位球面 SN−1 の元として、パラメータ表示 x = rθ ∈ RN をすれば Δ f = ∂ 2 f ∂ r 2 + N − 1 r ∂ f ∂ r + 1 r 2 Δ S N − 1 f {\displaystyle \Delta f={\frac {\partial ^{2}f}{\partial r^{2}}}+{\frac {N-1}{r}}{\frac {\partial f}{\partial r}}+{\frac {1}{r^{2}}}\Delta _{S^{N-1}}f} と書ける。ただし、∆SN−1 は球ラプラシアンとも呼ばれる (N−1)-次元球面上のラプラス＝ベルトラミ作用素である。二つの球対称微分項は 1 r N − 1 ∂ ∂ r ( r N − 1 ∂ f ∂ r ) {\displaystyle {\frac {1}{r^{N-1}}}{\frac {\partial }{\partial r}}{\Bigl (}r^{N-1}{\frac {\partial f}{\partial r}}{\Bigr )}} と書いても同じことである。一つの帰結として、SN−1 ⊂ RN 上で定義される函数の球ラプラシアンは R N ∖ { 0 } {\textstyle \mathbf {R} ^{N}\backslash \{0\}} へ延長した函数の通常のラプラシアンとして計算することができて、それは半直線に沿って定数（つまり、斉零次の斉次函数）になる。 ラプラス作用素のスペクトルは、対応する固有函数 f が − Δ f = λ f {\displaystyle -\Delta f=\lambda f} を満たすようにできる固有値 −λ の全てからなる[要検証 – ノート]。上の式はヘルムホルツ方程式と呼ばれるものである。 Ω を Rn の有界領域とすれば、ラプラス作用素の固有函数全体はヒルベルト空間 L2(Ω) の正規直交基底を成す。この結果は本質的にはコンパクト自己随伴作用素に関するスペクトル定理をラプラス作用素の逆作用素（これはポワンカレ不等式およびコンドラコフ埋蔵定理（英語版）によってコンパクト）に適用することにより従う。固有函数が無限回微分可能函数であることも示せる。この結果はより一般に、任意の境界付きコンパクトリーマン多様体上のラプラス＝ベルトラム作用素について成り立ち、また実際に有界領域上滑らかな係数を持つ任意の楕円型作用素に対するディリクレ固有値問題についても正しい。Ω が超球面であるときの、ラプラス作用素の固有函数は球面調和函数と呼ばれる。 ラプラス作用素の概念は、リーマン多様体上で定義されたラプラス＝ベルトラミ作用素（英語版）と呼ばれる楕円型作用素に一般化することができる。同様にダランベール作用素は擬リーマン多様体上の双曲型作用素に一般化される。ラプラス＝ベルトラミ作用素を函数に適用すれば、その函数のヘッセ行列のトレース Δ f = tr ⁡ ( H ( f ) ) {\displaystyle \Delta f=\operatorname {tr} (H(f))} が得られる。ただし、トレースは計量テンソルの逆に関して取るものとする。ラプラス＝ベルトラミ作用素を同様の式でテンソル場に作用する作用素（これもまたラプラス＝ベルトラミ作用素と呼ばれる）に一般化することができる。 ラプラス作用素の別な一般化として、擬リーマン多様体上で定義される外微分を用いた「幾何学者のラプラシアン」と呼ばれる Δ f = d ∗ d f {\displaystyle \Delta f=d^{*}df} を考えることもできる。ここで d∗は余微分で、ホッジ双対を使って書くこともできる。これが上に述べた「解析学者のラプラシアン」とは異なるものであることには注意すべきである。そのことは大域解析学の論文を読むときには常に気を付けねばならない。 より一般に、微分形式に対して定義される「ホッジ」ラプラシアン α は Δ α = d ∗ d α + d d ∗ α {\displaystyle \Delta \alpha =d^{*}d\alpha +dd^{*}\alpha } と書ける。これはまたラプラス＝ドラーム作用素（英語版）とも呼ばれ、ヴァイツェンベック不等式（英語版）によってラプラス＝ベルトラミ作用素と関係する。 ラプラシアンを適当な仕方によって非ユークリッド空間に一般化することができて、それには楕円型、双曲型、超双曲型（英語版）などが可能である。 ミンコフスキー空間におけるラプラス＝ベルトラミ作用素はダランベール作用素 ◻ = 1 c 2 ∂ 2 ∂ t 2 − ∂ 2 ∂ x 2 − ∂ 2 ∂ y 2 − ∂ 2 ∂ z 2 {\displaystyle \square ={\frac {1}{c^{2}}}{\partial ^{2} \over \partial t^{2}}-{\partial ^{2} \over \partial x^{2}}-{\partial ^{2} \over \partial y^{2}}-{\partial ^{2} \over \partial z^{2}}} となる。これは考える空間上の等長変換群のもとで不変な微分作用素であるという意味においてラプラス作用素の一般化となるものであり、時間不変函数へ制限する限りにおいてはラプラス作用素に帰着される。ここでは計量の符号を作用素の空間成分に関して負符号を許すようにしてあることに注意（高エネルギー粒子物理学ではこう仮定するのが普通）。ダランベール作用素は波動方程式に現れる微分作用素であるという理由で波動作用素と呼ばれることもある。これはまたクライン＝ゴルドン方程式（質量の無い場合には波動方程式に帰着される）の成分でもある。 計量における余分な因子 c は、物理学において空間と時間を異なる単位で測っている場合に必要となるものである（例えば同様のことは x-方向をメートルで y-方向をセンチメートルで測ったりするような場合にも出てくる）。実際、理論物理学では方程式を簡単にする目的で、自然単位系などの単位系のもと c = 1 として扱うのがふつうである。 ベクトルラプラス作用素（英語版）: ベクトル場に対する一般化 微分幾何におけるラプラス作用素（英語版） 離散ラプラス作用素（英語版）: 本項の連続的なラプラシアンに対して、有限差分化した類似対応物。 ラプラシアン行列