Laplace作用素(Laplace operator)とは、$\mathbb{R}^n$ の開集合上の $C^2$ 級関数 $u$ に 2 階偏導関数の和 $\Delta u=\sum_j\partial^2u/\partial x_j^2$ を対応させる線形微分作用素であり、$\Delta u=0$ の解が調和関数である。平行移動と直交変換と交換し、2 階偏導関数の定数係数の 1 次結合で表される作用素のうち回転と交換するものはその定数倍に限る。$n\ge2$ なら $\Delta u(a)$ は小さな球面の上の平均と $u(a)$ の差の $2n/r^2$ 倍の極限に等しく、台がコンパクトな関数 $u$ について $\int(-\Delta u)\overline u\,dx=\int|\nabla u|^2dx\ge0$ が成り立つ。
Laplace 作用素(Laplace operator、Laplacian、ラプラシアン)とは、$\mathbb{R}^n$ の開集合上の 2 回連続微分可能な関数 $u$ に、2 階偏導関数の和
$$
\Delta u=\frac{\partial^2u}{\partial x_1^2}+\frac{\partial^2u}{\partial x_2^2}+\dots+\frac{\partial^2u}{\partial x_n^2}
$$
を対応させる 2 階の線形微分作用素である。勾配の発散 $\Delta u=\operatorname{div}(\nabla u)$ とも書ける。$\Delta u=0$ を満たす関数が 調和関数 であり、Laplace 作用素は Laplace 方程式・Poisson 方程式・熱方程式・波動方程式のいずれにも現れる。
Laplace 作用素の特徴は 2 つの見方にまとめられる。1 つは対称性で、平行移動と回転(直交変換)と交換し、2 階偏導関数の定数係数の 1 次結合で表される作用素のうち回転と交換するものは、Laplace 作用素の定数倍に限る。もう 1 つは平均との関係で、$\Delta u(a)$ は、$a$ のまわりの小さな球面の上の $u$ の平均が $u(a)$ からどれだけずれるかを測る。この記事では $\mathbb{R}^n$ の開集合の上の Laplace 作用素を扱い、多様体の微分形式の上の Laplace–de Rham作用素 や、グラフの Laplacian行列 との関係は最後の節で触れる。
以下、$\Omega\subset\mathbb{R}^n$($n\ge1$)を空でない開集合とし、関数は実数値または複素数値とする。$C^k(\Omega)$ は $\Omega$ 上の $k$ 回連続微分可能な関数の全体、$\nabla u=(\partial_1u,\dots,\partial_nu)$ は勾配、$\partial_j=\partial/\partial x_j$ である。
$u\in C^2(\Omega)$ に対して
$$
\Delta u:=\sum_{j=1}^n\partial_j^2u
$$
とおき、写像 $\Delta\colon C^2(\Omega)\to C(\Omega)$ を Laplace 作用素 という。$\nabla^2$ と書くこともある。
$\Delta$ は線形であり、ベクトル場 $w=(w_1,\dots,w_n)$ の発散 $\operatorname{div}w=\sum_j\partial_jw_j$ を使えば $\Delta u=\operatorname{div}(\nabla u)$ である。微分の順序を交換できるので、$u\in C^3(\Omega)$ なら $\partial_k\Delta u=\Delta\partial_ku$ である(ABR20 p. 2)。定義と名前は ABR20 p. 1 による。
符号の取り方には 2 つの流儀がある。解析学では上の $\Delta$ を使い、$-\Delta$ が「正の作用素」になる(prop-lap-positive)。微分幾何学では $-\sum_j\partial_j^2$ を Laplace 作用素と呼ぶことが多く、Laplace–de Rham作用素 は関数に対してこの符号のものになる。
3 は、平面波 $e^{i\xi\cdot x}$ が $\Delta$ の固有関数で固有値が $-|\xi|^2$ であることを示している。実部をとると、たとえば正方形 $(0,\pi)^2$ の上で $\sin(kx_1)\sin(lx_2)$($k,l$ は正の整数)は境界で $0$ になり、$\Delta$ をかけると $-(k^2+l^2)$ 倍になる。有界な領域の上で境界で $0$ になる固有関数の全体は Laplacianの固有関数 で扱う。
$u,v\in C^2(\Omega)$ について
$$
\Delta(uv)=u\,\Delta v+2\,\nabla u\cdot\nabla v+v\,\Delta u
$$
が成り立つ。ここで $\nabla u\cdot\nabla v=\sum_j\partial_ju\,\partial_jv$(複素共役はとらない)である。また $u$ が実数値で $\phi\in C^2(\mathbb{R})$ なら、$\Delta(\phi\circ u)=\phi''(u)\,|\nabla u|^2+\phi'(u)\,\Delta u$ である。
積の微分法を 2 回使うと $\partial_j^2(uv)=u\,\partial_j^2v+2\,\partial_ju\,\partial_jv+v\,\partial_j^2u$ であり、$j$ について足せば第 1 式である。合成関数の微分法から $\partial_j(\phi\circ u)=\phi'(u)\,\partial_ju$、$\partial_j^2(\phi\circ u)=\phi''(u)(\partial_ju)^2+\phi'(u)\,\partial_j^2u$ であり、$j$ について足せば第 2 式である。$\square$
第 1 式は ABR20 1.19(p. 13)にある。たとえば $u=v=x_1$ なら $\Delta(x_1^2)=0+2+0=2$ である。第 2 式で $\phi(s)=s^2$ とすると、実数値の調和関数 $u$ について $\Delta(u^2)=2|\nabla u|^2\ge0$ となる。
$u\in C^2(\Omega)$ とする。
1 と 2 は合成関数の微分法から $\partial_j^2u_y(x)=(\partial_j^2u)(x-y)$、$\partial_j^2u^{(r)}(x)=r^2(\partial_j^2u)(rx)$ となることによる。
3.$T=(t_{pj})$ とすると $(Tx)_p=\sum_jt_{pj}x_j$ なので、合成関数の微分法から
$$
\partial_j(u\circ T)=\sum_pt_{pj}\,(\partial_pu)\circ T,\qquad\partial_j^2(u\circ T)=\sum_{p,q}t_{pj}t_{qj}\,(\partial_q\partial_pu)\circ T
$$
である。$j$ について足すと、$\sum_jt_{pj}t_{qj}=(TT^{\mathsf T})_{pq}=\delta_{pq}$($T$ は直交行列)から $\Delta(u\circ T)=\sum_p(\partial_p^2u)\circ T=(\Delta u)\circ T$ である。$\square$
これは ABR20 の Invariance Properties の節(pp. 2–4)にある。3 は直交行列であることを本質的に使う:たとえば $T(x_1,x_2)=(2x_1,x_2)$ と調和関数 $u=x_1^2-x_2^2$ では $u\circ T=4x_1^2-x_2^2$、$\Delta(u\circ T)=6\ne0$ である。
逆に、回転との可換性は Laplace 作用素を特徴づける。
実数を成分とする対称行列 $A=(a_{jk})$ について、定数係数の 2 階の微分作用素を $L_Au:=\sum_{j,k}a_{jk}\,\partial_j\partial_ku$($u\in C^2(\mathbb{R}^n)$)とおく。すべての $u\in C^2(\mathbb{R}^n)$ とすべての直交行列 $T$ について $L_A(u\circ T)=(L_Au)\circ T$ が成り立つなら、ある実数 $c$ について $A=cI$、すなわち $L_A=c\,\Delta$ である。
段 1.thm-lap-invariance の 3 の証明と同じ計算で、$\partial_j\partial_k(u\circ T)=\sum_{p,q}t_{pj}t_{qk}(\partial_p\partial_qu)\circ T$ なので
$$
L_A(u\circ T)=\sum_{p,q}\Bigl(\sum_{j,k}t_{pj}a_{jk}t_{qk}\Bigr)(\partial_p\partial_qu)\circ T=\bigl(L_{TAT^{\mathsf T}}u\bigr)\circ T
$$
である($TAT^{\mathsf T}$ も対称行列)。
段 2.対称行列 $C=(c_{pq})$ は作用素 $L_C$ で決まる:$u=\frac12x_p^2$ なら $L_Cu=c_{pp}$、$p\ne q$ で $u=x_px_q$ なら $L_Cu=c_{pq}+c_{qp}=2c_{pq}$ である。仮定と段 1 から $L_{TAT^{\mathsf T}}u=L_Au$ がすべての $u$ で成り立つ($T$ は全単射)ので、$TAT^{\mathsf T}=A$、右から $T$ を掛けて $TA=AT$ である。
段 3.$e_j$ を標準基底とし、$e_j$ に垂直な超平面に関する鏡映 $R_j:=I-2e_je_j^{\mathsf T}$(直交行列)をとる。$R_jA=AR_j$ から $e_je_j^{\mathsf T}A=Ae_je_j^{\mathsf T}$ であり、両辺を $e_j$ に掛けると $Ae_j=a_{jj}e_j$ を得る。よって $A$ は対角行列である。次に $i\ne j$ について、$e_i$ と $e_j$ を入れ替える置換行列 $T$(直交行列)をとると、$TAT^{\mathsf T}=A$ の $(i,i)$ 成分を比べて $a_{jj}=a_{ii}$ を得る。したがって $A=cI$ である。$\square$
逆に $L=c\,\Delta$ が回転と交換することは thm-lap-invariance の 3 である。係数が定数であるという仮定は外せない:$|x|^2\Delta$ は $|Tx|=|x|$ からすべての回転と交換するが、$\Delta$ の定数倍ではない。
回転で不変な関数、つまり $|x|$ だけで決まる関数(動径関数)には、Laplace 作用素は 1 変数の微分作用素として働く。
$0\le a< b\le\infty$、$\phi\in C^2((a,b))$ とし、$A:=\{x\in\mathbb{R}^n\mid a< |x|< b\}$ 上で $u(x):=\phi(|x|)$ とおく。このとき $r=|x|$ として
$$
\Delta u(x)=\phi''(r)+\frac{n-1}r\,\phi'(r)=\frac1{r^{n-1}}\bigl(r^{n-1}\phi'(r)\bigr)'
$$
である。
$x\ne0$ で $\partial_j|x|=x_j/r$ なので $\partial_ju=\phi'(r)\,x_j/r$ であり、もう一度微分すると
$$
\partial_j^2u=\phi''(r)\frac{x_j^2}{r^2}+\phi'(r)\Bigl(\frac1r-\frac{x_j^2}{r^3}\Bigr)
$$
である。$\sum_jx_j^2=r^2$ を使って $j$ について足すと $\phi''(r)+\phi'(r)(n-1)/r$ になる。最後の等号は $\bigl(r^{n-1}\phi'\bigr)'=r^{n-1}\phi''+(n-1)r^{n-2}\phi'$ による。$\square$
$n\ge2$ とし、$A$ を prop-lap-radial の円環領域とする。$A$ 上の動径関数 $u(x)=\phi(|x|)$ が調和であるための必要十分条件は、定数 $c_1,c_2$ があって
$$
u(x)=\begin{cases}c_1+c_2|x|^{2-n}&(n\ge3)\\ c_1+c_2\log|x|&(n=2)\end{cases}
$$
となることである。とくに球 $B(0,b)$ の上の動径の調和関数は定数に限る。
$u$ が調和なら、$u$ は $C^2$ 級なので $\phi(r)=u(re_1)$ も $C^2$ 級であり、prop-lap-radial から $(r^{n-1}\phi')'=0$、すなわち $r^{n-1}\phi'(r)=C$(定数)である。$\phi'(r)=Cr^{1-n}$ を積分すると、$n\ge3$ なら $\phi=c_1+\frac{C}{2-n}r^{2-n}$、$n=2$ なら $\phi=c_1+C\log r$ である。逆にこの形の関数は、同じ式を逆にたどれば調和である。球 $B(0,b)$ の上の動径の調和関数は、$0< |x|< b$ で上の形をとり、原点で連続なので、$|x|\to0$ で発散する項の係数 $c_2$ は $0$ である。$\square$
$n=1$ では調和関数は 1 次関数で、動径関数 $\phi(|x|)$ が $0$ の両側で調和なら、それぞれの側で 1 次関数になる。$n\ge2$ の $|x|^{2-n}$ と $\log|x|$ を定数倍したものは Laplace 作用素の基本解と呼ばれる(Green関数)。
$n=2$ で $U(r,\theta):=u(r\cos\theta,r\sin\theta)$($r>0$)とおくと
$$
\Delta u=U_{rr}+\frac1rU_r+\frac1{r^2}U_{\theta\theta}
$$
である。たとえば $U=r^m\cos m\theta$ なら $m(m-1)r^{m-2}+mr^{m-2}-m^2r^{m-2}=0$ なので調和であり、これは $\operatorname{Re}(x_1+ix_2)^m$ である。
$c=\cos\theta$、$s=\sin\theta$ と書く。合成関数の微分法から $U_r=cu_{x_1}+su_{x_2}$、$U_{rr}=c^2u_{x_1x_1}+2csu_{x_1x_2}+s^2u_{x_2x_2}$、$U_\theta=r(-su_{x_1}+cu_{x_2})$、$U_{\theta\theta}=r^2(s^2u_{x_1x_1}-2csu_{x_1x_2}+c^2u_{x_2x_2})-r(cu_{x_1}+su_{x_2})$ である。$U_{rr}+U_{\theta\theta}/r^2$ は $u_{x_1x_1}+u_{x_2x_2}-U_r/r$ に等しいので、$U_r/r$ を足すと $\Delta u$ になる。
一般の $n$ では、$x=r\vartheta$($r>0$、$\vartheta\in S$)と書くと $\Delta u=\partial_r^2u+\frac{n-1}r\partial_ru+\frac1{r^2}\Delta_Su$ となる。$\Delta_S$ は単位球面の上の Laplace 作用素(Laplace–Beltrami 作用素)である(DV25 Theorem 2.5.2、p. 322。この記事では証明しない)。$\Delta_S$ の固有関数が 球面調和関数 である。
$n\ge2$ とし、$S$ を単位球面、$\sigma$ を $S$ 上の正規化された表面積測度($\sigma(S)=1$、直交変換で不変)とする(調和関数 の記事と同じ記号)。
$u\in C^2(\Omega)$、$a\in\Omega$ とすると
$$
\Delta u(a)=\lim_{r\to0+}\frac{2n}{r^2}\Bigl(\int_Su(a+r\zeta)\,d\sigma(\zeta)-u(a)\Bigr)
$$
が成り立つ。
球面上の積分.$\int_S\zeta_j\,d\sigma=0$、$\int_S\zeta_j\zeta_k\,d\sigma=\delta_{jk}/n$ を示す。第 $j$ 座標の符号を変える鏡映は直交変換で $\sigma$ を保ち、$\zeta_j$ と($k\ne j$ について)$\zeta_j\zeta_k$ の符号を変えるので、これらの積分は $0$ である。2 つの座標を入れ替える直交変換で $\int_S\zeta_j^2\,d\sigma$ はすべての $j$ で等しく、和は $\int_S|\zeta|^2\,d\sigma=1$ なので、どれも $1/n$ である。
Taylor の定理.$\overline B(a,r_0)\subset\Omega$ とする。$|h|\le r_0$ について、Taylorの定理 から
$$
u(a+h)=u(a)+\nabla u(a)\cdot h+\frac12\sum_{j,k}\partial_j\partial_ku(a)\,h_jh_k+R(h),\qquad|R(h)|\le\frac n2\,\omega(|h|)\,|h|^2
$$
が成り立つ。ここで $\omega(\rho):=\max_{j,k}\sup_{|y-a|\le\rho}|\partial_j\partial_ku(y)-\partial_j\partial_ku(a)|$ であり、2 階偏導関数の連続性から $\rho\to0$ で $\omega(\rho)\to0$ である(剰余を積分の形 $R(h)=\int_0^1(1-t)\sum_{j,k}\bigl(\partial_j\partial_ku(a+th)-\partial_j\partial_ku(a)\bigr)h_jh_k\,dt$ で書き($u$ が複素数値でもよい)、2 階偏導関数の差を $\omega$ で押さえ、$\sum_{j,k}|h_jh_k|\le n|h|^2$ を使う)。
$h=r\zeta$ を代入して $\sigma$ で積分すると、1 次の項は消え、2 次の項は $\frac{r^2}2\sum_j\partial_j^2u(a)\cdot\frac1n=\frac{r^2}{2n}\Delta u(a)$ になるので
$$
\Bigl|\int_Su(a+r\zeta)\,d\sigma(\zeta)-u(a)-\frac{r^2}{2n}\Delta u(a)\Bigr|\le\frac n2\,\omega(r)\,r^2
$$
である。両辺に $2n/r^2$ を掛けて $r\to0$ とすればよい。$\square$
$\Delta u(a)>0$ なら、小さい球面の上の $u$ の平均は $u(a)$ より大きい。$\Delta u=0$ なら平均と中心の値が(極限の意味で)つり合い、実際には調和関数の球面平均は中心の値にちょうど等しい(調和関数 の記事の定理「平均値の性質」)。この見方から、$\Delta u\ge0$ となる関数は劣平均値の性質をもつ 劣調和関数 になる。
$u\in C^2(\Omega)$、$a\in\Omega$、$h>0$ について
$$
\Delta_hu(a):=\sum_{j=1}^n\frac{u(a+he_j)+u(a-he_j)-2u(a)}{h^2}
$$
とおくと、$h\to0$ で $\Delta_hu(a)\to\Delta u(a)$ である。実部と虚部に分けて考えればよいので、$u$ は実数値としてよい。各 $j$ について、1 変数関数 $g(t):=u(a+te_j)$ に Taylor の定理を使うと $g(h)+g(-h)-2g(0)=h^2g''(\tau)$ となる $|\tau|\le h$ があり($g(\pm h)$ の 2 次の剰余 $\frac{h^2}2g''(\tau_\pm)$ の和に中間値の定理を使う)、$g''(\tau)=\partial_j^2u(a+\tau e_j)\to\partial_j^2u(a)$ だからである。
平面($n=2$)の $\Delta_h$ は、格子点 $a$ の値と上下左右の 4 点の値の差の和を $h^2$ で割ったもので、5 点差分と呼ばれる。格子 $h\mathbb{Z}^2$ をグラフとみると、$h^2\Delta_hu(a)=\sum_{b\sim a}\bigl(u(b)-u(a)\bigr)$($b$ は $a$ の隣の点)であり、これはグラフの Laplacian行列 $L=D-A$ を $u$ に掛けたものの符号を変えたものである。
$u,v\in C^2(\Omega)$ とし、$v$ の台は $\Omega$ のコンパクト集合とする。このとき
$$
\int_\Omega(\Delta u)\,v\,dx=-\int_\Omega\nabla u\cdot\nabla v\,dx=\int_\Omega u\,(\Delta v)\,dx
$$
が成り立つ。とくに $u\in C^2(\Omega)$ の台が $\Omega$ のコンパクト集合なら
$$
\int_\Omega(-\Delta u)\,\overline u\,dx=\int_\Omega|\nabla u|^2\,dx\ge0
$$
であり、等号は $u=0$ のときに限る。
$w:=v\,\partial_ju$ を $\Omega$ の外で $0$ と延ばすと、$\mathbb{R}^n$ 上の台がコンパクトな $C^1$ 級関数である。Fubiniの定理 で $x_j$ について先に積分すると、微積分学の基本定理 から $\int\partial_jw\,dx=0$ である。$\partial_jw=\partial_jv\,\partial_ju+v\,\partial_j^2u$ を $j$ について足せば第 1 の等号を得る。$u$ と $v$ の役割を入れ替えた同じ計算(今度は $u\,\partial_jv$ の台がコンパクト)で $\int u\,\Delta v=-\int\nabla u\cdot\nabla v$ となり、第 2 の等号を得る。後半は第 1 の等号で $v=\overline u$ とおけばよい($\nabla u\cdot\nabla\overline u=|\nabla u|^2$)。等号が成り立てば $\nabla u\equiv0$ で、$u$ を $0$ で延ばした台がコンパクトな関数は定数、すなわち $0$ である。$\square$
この性質から、境界で $0$ になる関数に $-\Delta$ を作用させた固有値問題 $-\Delta u=\lambda u$ の固有値は正になる(Laplacianの固有関数)。Fourier変換 の言葉では、$\partial_j$ が $i\xi_j$ を掛けることに対応するので、$-\Delta$ は $|\xi|^2$ を掛ける操作に対応する(ex-lap-basic の 3 の平面波がその固有関数である)。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| $T$ が直交行列(thm-lap-invariance の 3) | $T(x_1,x_2)=(2x_1,x_2)$、$u=x_1^2-x_2^2$ | $\Delta(u\circ T)=(\Delta u)\circ T$ |
| 係数が定数(thm-lap-rotation) | $L=\lvert x\rvert^2\Delta$ | 回転と交換する 2 階の作用素は $\Delta$ の定数倍 |
| $v$ の台がコンパクト(prop-lap-positive) | $\Omega=B(0,1)$、$u=v=\lvert x\rvert^2$ | $\int(-\Delta u)\,\overline u\,dx\ge0$ |
| $u$ が動径関数(cor-lap-radial-harmonic) | $u=x_1$($n\ge2$) | 調和関数は $c_1+c_2\lvert x\rvert^{2-n}$ または $c_1+c_2\log\lvert x\rvert$ の形 |
1 行目:$T$ は直交行列でない($\lvert Te_1\rvert=2$)。$u$ は調和だが $u\circ T=4x_1^2-x_2^2$ の Laplace 作用素は $8-2=6\ne0=(\Delta u)\circ T$ である。2 行目:$\lvert Tx\rvert=\lvert x\rvert$ から $\lvert x\rvert^2(\Delta(u\circ T))(x)=\lvert Tx\rvert^2(\Delta u)(Tx)$ で、回転と交換するが、係数行列 $\lvert x\rvert^2I$ は定数でない。3 行目:$u=\lvert x\rvert^2$ は $\mathbb{R}^n$ 上で $C^2$ 級だが、$B(0,1)$ の中のコンパクト集合に台が入らない。$-\Delta u=-2n$ なので $\int_{B(0,1)}(-2n)\lvert x\rvert^2\,dx< 0$ である。境界の項 $\int_{\partial B}u\,\partial_\nu u\,dS$ が残ることが原因である。4 行目:$x_1$ は調和だが、$x_1$ は $\lvert x\rvert$ だけでは決まらず、どの $c_1,c_2$ についても上の形にならない($x_1$ と $-x_1$ で値が違う)。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する