ラプラス作用素

数学における**ラプラス作用素**（ラプラスさようそ、）あるいは**ラプラシアン**（）は、ユークリッド空間上の函数の勾配の発散として与えられる微分作用素である。記号では , , あるいは  で表されるのが普通である。函数  の点  におけるラプラシアン  は（次元に依存する定数の違いを除いて）点  を中心とする球面を半径が増大するように動かすときの  から得られる平均値になっている。直交座標系においては、ラプラシアンは各独立変数に関する函数の二階（非混合）偏導函数の和として与えられ、またほかに円筒座標系や球座標系などの座標系においても有用な表示を持つ。

ラプラス作用素の名称は、天体力学の研究に同作用素を最初に用いたフランス人数学者のピエール＝シモン・ド・ラプラス (1749–1827) に因んでいる。同作用素は与えられた重力ポテンシャルに適用すると質量密度の定数倍を与える。現在ではラプラス方程式と呼ばれる方程式  の解は調和函数と呼ばれ、自由空間において可能な重力場を表現するものである。

微分方程式においてラプラス作用素は電気ポテンシャル、重力ポテンシャル、熱や流体の拡散方程式、波の伝搬、量子力学といった、多くの物理現象を記述するのに現れる。ラプラシアンは、函数の勾配フローの流束密度を表す。

## 定義
ラプラス作用素は 次元ユークリッド空間上の函数  の勾配  の発散  として定義される二階の微分作用素である。つまり、 が二回微分可能実数値函数ならば  のラプラシアンは

:\Delta f \equiv \nabla^2 f := \nabla \cdot \nabla f 

で定義される。ただし、あとの記法は形式的に と書いたものである。あるいは同じことだが、 のラプラシアンは直交座標系  における**非混合**二階偏導函数の全てにわたる和

:\Delta f = \sum_{i=1}^n \frac {\partial^2 f}{\partial x^2_i}

としても書ける。二階の微分作用素として、ラプラス作用素は 級函数を  級の函数へ写す ()。つまり、式 1 (あるいは同値な 2) は作用素  を定める。あるいはより一般に任意の開集合  に対して作用素  を定める。

## 数学的特徴づけ
ラプラス作用素は、合同変換と可換である。すなわち、任意の級関数と任意の合同変換に対し、
: \Delta(\varphi(T(x)))=T(\Delta(\varphi(x)))
が成立する。

しかもラプラス作用素は、上記の性質を満たす非自明な微分演算子で最も簡単なものとして特徴づけることができる。これを説明する為、記号を導入する。を実数の集合とし、個の実数からなる組の集合をとする。と個の非負整数の組に対し、
:\frac{\partial}{\partial x^{\alpha}} := \frac{\partial^n}{\partial {x_1}^{\alpha_1} \dotsm \partial {x_n}^{\alpha_n}},
: |\alpha| := \alpha_1 + \dotsb + \alpha_n
と表記する。微分演算子
:D:=\sum_{\alpha\colon |\alpha| \le k} a_{\alpha} \frac{\partial}{\partial x^{\alpha}}
が任意の級関数と向きを保つ任意の合同変換に対し、
: D(\varphi(T(x)))=T(D(\varphi(x)))
が成立していたとする。このとき、実数係数の1変数多項式 p(X) = \sum_m u_m X^m が存在し、
: D = p(\Delta) = \sum_m u_m \Delta^m
が成立する。

よってラプラス作用素は、合同変換に対して不変な微分演算子の中で、自明なもの（＝恒等的にを対応させる微分演算子）を除けば最も簡単なものである。

## 動機付け
### 拡散
拡散の物理理論において、ラプラス作用素は（ラプラス方程式を通じて）平衡の数学的記述に自然に現れる。具体的に、 が化学濃度のような適当な量の平衡密度であるとき、 の滑らかな境界を持つ領域  を通る流束が、 に流入も漏出も無いとすれば、 であるから
: \int_{\partial V} \nabla u \cdot \boldsymbol{n}\, dS = 0
と書ける。ただし、 は領域  の境界に対して外側を向く単位法ベクトルである。発散定理により
:\int_V \operatorname{div} \nabla u\, dV = \int_{\partial V} \nabla u \cdot \boldsymbol{n}\, dS = 0
は領域  が滑らかな境界を持つ限りにおいて成り立つから、これにより
: \operatorname{div} \nabla u = \Delta u = 0
が導かれる。方程式の左辺はラプラス作用素である。ラプラス作用素それ自身は拡散方程式によって記述されるような、化学濃度の流入や漏出を表す点を含む非平衡拡散に対する物理的解釈を持つ。

### ポテンシャルに付随する密度
 が電荷分布  に付随した電位を記述するものとすると、電荷分布自身は  のラプラシアンとして

で与えられる。これはガウスの法則の帰結である。実際、 が任意の滑らかな領域ならば、電場  の電束に関するガウスの法則により、（単位当たりの）電荷は
: \int_{\partial V} \boldsymbol{E}\cdot \boldsymbol{n}\, dS = \int_{\partial V} \nabla \varphi \cdot \boldsymbol{n}\, dS = \int_V q\,dV
になる。ただし、最初の等号は静電場は静電位の勾配に等しいという事実を用いた。発散定理により、
: \int_V \Delta\varphi\,dV = \int_V q\, dV
が成り立ち、これは任意の領域  に対して成り立つことから () を得る。

同じ説明によって、重力ポテンシャルのラプラシアンがとなることが導かれる。電荷や質量の分布が与えられていてそれらに付随するポテンシャルは未知ということはよくあることである。適当な境界条件の下でポテンシャル函数を求めるということは、ポワソン方程式を解くことに同じである。

### エネルギー最小化
物理学においてラプラス作用素が現れる別な理由は、領域  における方程式  の解はディリクレエネルギー・汎函数を停留させる函数
: E(f) := \frac{1}{2} \int_U \lVert \nabla f \rVert^2 \,dx
となることである。これを見るために  は函数で、函数  は  の境界上で消えていると仮定する。このとき
: 
 \frac{d}{d\varepsilon} E(f+\varepsilon u) \Bigg|_{\varepsilon = 0}
= \int_U \nabla f \cdot \nabla u \, dx
= -\int_U u \Delta f\, dx 

が成り立つ（ただし、最後の等号はグリーンの恒等式を用いた）。この計算により、 ならば  は  の周りで停留する。逆に  が  の周りで停留するならば により  である。

## 各種座標表示
### 二次元
二次元のラプラス作用素は  を -平面上の標準直交座標として
:\Delta f := \frac{\partial^2f}{\partial x^2} + \frac{\partial^2 f}{\partial y^2}
で与えられる。

; 極座標
: \begin{align}
 \Delta f &= \frac{1}{r} \frac{\partial}{\partial r}
  \left( r \frac{\partial f}{\partial r} \right) 
+ \frac{1}{r^2} \frac{\partial^2 f}{\partial \theta^2}\\
&= \frac{1}{r} \frac{\partial f}{\partial r} 
+ \frac{\partial^2 f}{\partial r^2}
+ \frac{1}{r^2} \frac{\partial^2 f}{\partial \theta^2}.
\end{align} 


### 三次元

三次元では様々な座標系がラプラシアンを記述するために広く用いられる。

; 直交座標系
: 
\Delta f = \frac{\partial^2 f}{\partial x^2} + \frac{\partial^2 f}{\partial y^2} + \frac{\partial^2 f}{\partial z^2}.

; 円筒座標系
:  \Delta f 
= \frac{1}{\rho} \frac{\partial}{\partial \rho}
  \left(\rho \frac{\partial f}{\partial \rho} \right) 
+ \frac{1}{\rho^2} \frac{\partial^2 f}{\partial \varphi^2}
+ \frac{\partial^2 f}{\partial z^2 }. 

; 球面座標系
:  \Delta f 
= \frac{1}{r^2} \frac{\partial}{\partial r}
  \left(r^2 \frac{\partial f}{\partial r} \right) 
+ \frac{1}{r^2 \sin \theta} \frac{\partial}{\partial \theta}
  \left(\sin \theta \frac{\partial f}{\partial \theta} \right) 
+ \frac{1}{r^2 \sin^2 \theta} \frac{\partial^2 f}{\partial \varphi^2}.

; 一般の 
: \nabla^2 = \nabla \xi^m \cdot \nabla \xi^n \frac{\partial^2}{\partial \xi^m \partial \xi^n} + \nabla^2 \xi^m \frac{\partial}{\partial \xi^m }, 
: ここでアインシュタインの和の規約を用いた。

### 一般次元
 次元球座標系において、 を正の実数をとる半径、 は単位球面 の元として、パラメータ表示  をすれば

:  \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


と書ける。ただし、 は球ラプラシアンとも呼ばれる (*N*&minus;1)-次元球面上のラプラス＝ベルトラミ作用素である。二つの球対称微分項は
:\frac{1}{r^{N-1}} \frac{\partial}{\partial r} \biggl(r^{N-1} \frac{\partial f}{\partial r} \biggr)
と書いても同じことである。一つの帰結として、 上で定義される函数の球ラプラシアンは  へ延長した函数の通常のラプラシアンとして計算することができて、それは半直線に沿って定数（つまり、斉零次の斉次函数）になる。

## スペクトル論

ラプラス作用素のスペクトルは、対応する固有函数  が
:-\Delta f = \lambda f
を満たすようにできる固有値  の全てからなる。上の式はヘルムホルツ方程式と呼ばれるものである。  を  の有界領域とすれば、ラプラス作用素の固有函数全体はヒルベルト空間 ''L''<sup>2</sup>(&Omega;) の正規直交基底を成す。この結果は本質的にはコンパクト自己随伴作用素に関するスペクトル定理をラプラス作用素の逆作用素（これはポワンカレ不等式およびレリッヒ＝コンドラショフの定理によってコンパクト）に適用することにより従う。固有函数が無限回微分可能函数であることも示せる。この結果はより一般に、任意の境界付きコンパクトリーマン多様体上のラプラス＝ベルトラム作用素について成り立ち、また実際に有界領域上滑らかな係数を持つ任意の楕円型作用素に対するディリクレ固有値問題についても正しい。 が超球面であるときの、ラプラス作用素の固有函数は球面調和函数と呼ばれる。

## 一般化
### ダランベール作用素
ラプラシアンを適当な仕方によって非ユークリッド空間に一般化することができて、それには楕円型、双曲型、超双曲型などが可能である。

ミンコフスキー空間におけるラプラス作用素はダランベール作用素（ダランベルシアン）
:
となる。これは考える空間上の等長写像群 (:en:Isometry group) のもとで不変な微分作用素であるという意味においてラプラス作用素の一般化となるものであり、時間不変函数へ制限する限りにおいてはラプラス作用素に帰着される。ここでは計量の符号を作用素の空間成分に関して負符号を許すようにしてあることに注意（高エネルギー素粒子物理学ではこう仮定するのが普通）。ダランベール作用素は波動方程式に現れる微分作用素であるという理由で波動作用素と呼ばれることもある。これはまたクライン＝ゴルドン方程式（質量の無い場合には波動方程式に帰着される）の成分でもある。

計量における余分な因子  は、物理学において空間と時間を異なる単位で測っている場合に必要となるものである（例えば同様のことは -方向をメートルで -方向をセンチメートルで測ったりするような場合にも出てくる）。実際、理論物理学では方程式を簡単にする目的で、自然単位系などの単位系のもと  として扱うのがふつうである。

### リーマン多様体上のラプラス作用素

ラプラス作用素の概念は、勾配  と発散  をリーマン多様体  上に一般化する事で定義でき、 上定義された関数  に対して
: \Delta f := \operatorname{div}\operatorname{grad}f 
 = - \frac{1}{\sqrt{|\operatorname{det}g|\,}}
     \frac{\partial}{\partial x^i} 
     \left( \sqrt}
ここでは曲率テンソルであり、(\iota_{e_j}R(e_i,e_j) \lrcorner\alpha)(X_1,\ldots,X_{n-1}) =\alpha(R(e_i,e_j)e_j,X_1,\ldots,X_{n-1})である。

上記の公式を**ヴァイツェンベック・ボホナーの公式**（）あるいは**ヴァイツェンベックの公式**（）という。

## 関連項目
* : ベクトル場に対する一般化
* : 本項の連続的なラプラシアンに対して、有限差分化した類似対応物。
* ラプラシアン行列

## 脚注

## 参考文献
*.
* .
* .
* .
* 

## 外部リンク
*
*