§E15.30調和関数

最終更新

球の Poisson 表示は、内部の値を境界値の正の重み付き平均として表す。この表示によって境界値から解を構成することができるが、解の内部での変動や特異点の有無にも制約が生じる。Laplace 方程式を満たす調和関数は楕円型方程式の解の最も基本的なクラスの一つであり、例えばu(x)=x12−x22u(x)=x_1^2-x_2^2は各点で周囲の球面上の平均値に一致する。境界値と内部の値を結ぶ平均の構造は、非負の解の値を比較し、孤立点の近くでの増大から特異点の除去可能性を判定するための基礎となる。本記事では、調和関数の平均値の性質と、その性質に基づく評価や特異点の扱いについて解説する。

1 平均値による特徴づけ

定義 1.1 (調和関数).n≥1n\geq1、Ω⊆Rn\Omega\subseteq\R^nを開集合とする。関数u∈C2(Ω;C)u\in C^2(\Omega;\C)が

Δu=∑j=1n∂xj2u=0\Delta u=\sum_{j=1}^n\partial_{x_j}^2u=0

を満たすとき、uuをΩ\Omega上の 調和関数 (harmonic function) という。実数値関数についても同じ定義を用いる。

定理 1.2 (平均値の性質).n≥2n\geq2、Ω⊆Rn\Omega\subseteq\R^nを開集合、u∈C(Ω;C)u\in C(\Omega;\C)とする。σn−1=∣Sn−1∣\sigma_{n-1}=|S^{n-1}|とおく。このとき、次の条件は同値である。

  1. uuは調和関数である。
  2. Br(x)‾⊆Ω\overline{B_r(x)}\subseteq\Omegaを満たす任意のx∈Ωx\in\Omegaとr>0r>0に対して、 u(x)=1σn−1rn−1∫∂Br(x)u(y) dσ(y)u(x)=\frac{1}{\sigma_{n-1}r^{n-1}}\int_{\partial B_r(x)}u(y)\,d\sigma(y) が成り立つ。
  3. Br(x)‾⊆Ω\overline{B_r(x)}\subseteq\Omegaを満たす任意のx∈Ωx\in\Omegaとr>0r>0に対して、 u(x)=1∣Br∣∫Br(x)u(y) dyu(x)=\frac{1}{|B_r|}\int_{B_r(x)}u(y)\,dy が成り立つ。

いずれかの条件を満たすuuはC∞(Ω)C^\infty(\Omega)に属する。

証明.(1)⇒\Rightarrow(2)を示す。Br(x)‾⊆Ω\overline{B_r(x)}\subseteq\Omegaとする。§E15.29 定理 5.1を平行移動して適用し、境界値u∣∂Br(x)u|_{\partial B_r(x)}の Poisson 積分h∈C∞(Br(x))∩C(Br(x)‾)h\in C^\infty(B_r(x))\cap C(\overline{B_r(x)})をとる。u−hu-hの実部と虚部はC2(Br(x))∩C(Br(x)‾)C^2(B_r(x))\cap C(\overline{B_r(x)})に属し、球内で調和的で、球面上で零である。−Δ-\Deltaに対する§E15.18 系 1.2をそれぞれの関数とその負に適用すると、実部と虚部はともに零となる。従ってu=hu=hであり、uuはBr(x)B_r(x)内部で滑らかである。中心における核は1/(σn−1rn−1)1/(\sigma_{n-1}r^{n-1})であるから、球面平均の等式を得る。

(2)⇒\Rightarrow(3)を示す。極座標積分により

∫Br(x)u(y) dy=∫0r∫∂Bt(x)u(y) dσ(y) dt=σn−1rnnu(x)=∣Br∣u(x)\int_{B_r(x)}u(y)\,dy =\int_0^r\int_{\partial B_t(x)}u(y)\,d\sigma(y)\,dt =\frac{\sigma_{n-1}r^n}{n}u(x)=|B_r|u(x)

が成り立つ。

(3)⇒\Rightarrow(2)を示す。球平均の等式をrrで微分すると、極座標積分とuuの連続性により

∫∂Br(x)u(y) dσ(y)=ddr(∣Br∣u(x))=σn−1rn−1u(x)\int_{\partial B_r(x)}u(y)\,d\sigma(y) =\frac{d}{dr}\bigl(|B_r|u(x)\bigr) =\sigma_{n-1}r^{n-1}u(x)

を得る。

(2)⇒\Rightarrow(1)を示す。非負の放射対称関数η∈Cc∞(B1)\eta\in C_c^\infty(B_1)を∫η=1\int\eta=1となるようにとり、ηε(y)=ε−nη(y/ε)\eta_\varepsilon(y)=\varepsilon^{-n}\eta(y/\varepsilon)とおく。B2ε(a)‾⊆Ω\overline{B_{2\varepsilon}(a)}\subseteq\Omegaとする。各x∈Bε(a)x\in B_\varepsilon(a)に対して、球面平均の等式を半径について積分すると

∫ηε(y)u(x−y) dy=u(x)∫ηε(y) dy=u(x)\int\eta_\varepsilon(y)u(x-y)\,dy=u(x)\int\eta_\varepsilon(y)\,dy=u(x)

を得る。B2ε(a)B_{2\varepsilon}(a)の近傍で11となるχ∈Cc∞(Ω)\chi\in C_c^\infty(\Omega)をとり、χu\chi uを零で全空間へ延長する。Bε(a)B_\varepsilon(a)上でu=ηε∗(χu)u=\eta_\varepsilon*(\chi u)であるから、§E15.2 命題 4.4によりuuは滑らかである。

§E4.5 系 3.1をuuの実部と虚部へ適用する。球面上での座標の奇対称性と回転対称性により、

1σn−1∫Sn−1ωi dσ=0,1σn−1∫Sn−1ωiωj dσ=δijn\frac{1}{\sigma_{n-1}}\int_{S^{n-1}}\omega_i\,d\sigma=0, \qquad \frac{1}{\sigma_{n-1}}\int_{S^{n-1}}\omega_i\omega_j\,d\sigma=\frac{\delta_{ij}}{n}

である。従って、球面平均と Taylor 展開から

0=1σn−1∫Sn−1(u(x+rω)−u(x)) dσ(ω)=r22nΔu(x)+o(r2)0=\frac{1}{\sigma_{n-1}}\int_{S^{n-1}}\bigl(u(x+r\omega)-u(x)\bigr)\,d\sigma(\omega) =\frac{r^2}{2n}\Delta u(x)+o(r^2)

を得る。r2r^2で割ってr↓0r\downarrow0とすればΔu(x)=0\Delta u(x)=0である。▨

系 1.3.n≥2n\geq2、Ω⊆Rn\Omega\subseteq\R^nを開集合とし、T∈D′(Ω)T\in\mathcal D'(\Omega)がΔT=0\Delta T=0を満たすとする。このとき、TTは一意の滑らかな調和関数uuによって与えられ、uuは定理 1.2 (2)と定理 1.2 (3)の等式を満たす。

証明.§E15.21 系 2.3により、TTは滑らかな調和関数uuによって与えられる。定理 1.2をuuへ適用する。二つの連続な代表はほとんど至る所で一致し、連続性により至る所で一致する。▨

例 1.4.n≥2n\geq2とし、u(x)=x12−x22u(x)=x_1^2-x_2^2、v(x)=∣x∣2v(x)=|x|^2とおく。このとき、Δu=2−2=0\Delta u=2-2=0である一方、Δv=2n\Delta v=2nである。球Br(0)B_r(0)上では各座標の積分が零であり、回転対称性と極座標積分から

1∣Br∣∫Br(0)yj2 dy=1n∣Br∣∫Br(0)∣y∣2 dy=r2n+2\frac{1}{|B_r|}\int_{B_r(0)}y_j^2\,dy =\frac{1}{n|B_r|}\int_{B_r(0)}|y|^2\,dy =\frac{r^2}{n+2}

となる。従って任意のx∈Rnx\in\R^nに対して

1∣Br∣∫Br(x)u(y) dy=u(x),1∣Br∣∫Br(x)v(y) dy=v(x)+nr2n+2\frac{1}{|B_r|}\int_{B_r(x)}u(y)\,dy=u(x), \qquad \frac{1}{|B_r|}\int_{B_r(x)}v(y)\,dy=v(x)+\frac{nr^2}{n+2}

である。

命題 1.5.n≥1n\geq1とm≥0m\geq0を整数、r>0r>0を実数とし、h∈L2(Br(x0);C)h\in L^2(B_r(x_0);\C)を調和関数とする。このとき、n,mn,mだけに依存する定数Cn,mC_{n,m}が存在し、

∥Dmh∥L∞(Br/2(x0))≤Cn,mr−m−n/2∥h∥L2(Br(x0))\|D^m h\|_{L^\infty(B_{r/2}(x_0))} \leq C_{n,m}r^{-m-n/2}\|h\|_{L^2(B_r(x_0))}

が成り立つ。ここでDmhD^m hは∣α∣=m|\alpha|=mの偏導関数の族であり、そのノルムは各成分のノルムの最大値とする。m=0m=0ではD0h=hD^0h=hとする。

証明. 非負の放射対称関数η∈Cc∞(B1/2)\eta\in C_c^\infty(B_{1/2})を∫η=1\int\eta=1となるように固定し、ηr(y)=r−nη(y/r)\eta_r(y)=r^{-n}\eta(y/r)とおく。n≥2n\geq2では定理 1.2 (2)と極座標積分により、各x∈Br/2(x0)x\in B_{r/2}(x_0)に対して

h(x)=∫Br(x0)ηr(x−y)h(y) dyh(x)=\int_{B_r(x_0)}\eta_r(x-y)h(y)\,dy

が成り立つ。n=1n=1ではh′′=0h''=0よりhhはアフィン関数であり、ηr\eta_rの偶対称性と積分の正規化から同じ等式を得る。

hhの球外への零延長とηr\eta_rの畳み込みに§E15.2 命題 4.4を適用し、∣α∣=m|\alpha|=mとすると、Cauchy–Schwarz の不等式により

∣∂αh(x)∣≤∥∂αηr∥L2(Rn)∥h∥L2(Br(x0))=r−m−n/2∥∂αη∥L2(Rn)∥h∥L2(Br(x0))|\partial^\alpha h(x)| \leq\|\partial^\alpha\eta_r\|_{L^2(\R^n)}\|h\|_{L^2(B_r(x_0))} =r^{-m-n/2}\|\partial^\alpha\eta\|_{L^2(\R^n)}\|h\|_{L^2(B_r(x_0))}

となる。有限個の多重指数について最大値をとれば結論を得る。▨

2 Harnack の不等式と Liouville の定理

定理 2.1 (Harnack の不等式).n≥2n\geq2とする。このとき、uuがBR(0)⊆RnB_R(0)\subseteq\R^n上の非負の実数値調和関数ならば、任意の0<r<R0<r<Rに対して

sup⁡Br(0)u≤(R+rR−r)ninf⁡Br(0)u\sup_{B_r(0)}u\leq\left(\frac{R+r}{R-r}\right)^n\inf_{B_r(0)}u

が成り立つ。また、Ω⊆Rn\Omega\subseteq\R^nが連結な開集合、K⊆ΩK\subseteq\Omegaが空でないコンパクト集合ならば、K,ΩK,\Omegaだけに依存する定数CK,Ω≥1C_{K,\Omega}\geq1が存在し、Ω\Omega上の任意の非負の実数値調和関数uuに対して

sup⁡Ku≤CK,Ωinf⁡Ku\sup_Ku\leq C_{K,\Omega}\inf_Ku

が成り立つ。

証明.r<ρ<Rr<\rho<Rとする。§E15.29 定理 5.1により、uuはBρ(0)B_\rho(0)上でその境界値の Poisson 積分に一致する。∣x∣=s≤r|x|=s\leq r、∣ζ∣=ρ|\zeta|=\rhoに対して、ρ−s≤∣x−ζ∣≤ρ+s\rho-s\leq|x-\zeta|\leq\rho+sより

ρn−2(ρ−r)(ρ+r)n−1≤Pρ(x,ζ)Pρ(0,ζ)≤ρn−2(ρ+r)(ρ−r)n−1\frac{\rho^{n-2}(\rho-r)}{(\rho+r)^{n-1}} \leq\frac{P_\rho(x,\zeta)}{P_\rho(0,\zeta)} \leq\frac{\rho^{n-2}(\rho+r)}{(\rho-r)^{n-1}}

となる。境界値u(ζ)u(\zeta)は非負であるから、積分して

ρn−2(ρ−r)(ρ+r)n−1u(0)≤u(x)≤ρn−2(ρ+r)(ρ−r)n−1u(0)\frac{\rho^{n-2}(\rho-r)}{(\rho+r)^{n-1}}u(0) \leq u(x)\leq \frac{\rho^{n-2}(\rho+r)}{(\rho-r)^{n-1}}u(0)

を得る。任意のx,y∈Br(0)x,y\in B_r(0)に対して両側の評価を組み合わせると

u(x)≤(ρ+rρ−r)nu(y)u(x)\leq\left(\frac{\rho+r}{\rho-r}\right)^n u(y)

である。ρ↑R\rho\uparrow Rとしてから上限と下限をとる。

各a∈Ωa\in\Omegaに対し、B2sa(a)‾⊆Ω\overline{B_{2s_a}(a)}\subseteq\Omegaとなるsa>0s_a>0を一つ固定する。球Bsa(a)B_{s_a}(a)内の任意の二点の値は、既に示した不等式により係数3n3^nで比較される。基点a0∈Ωa_0\in\Omegaを固定し、球族{Bsa(a):a∈Ω}\{B_{s_a}(a):a\in\Omega\}の有限個の球で、隣り合う球が交わる鎖によってa0a_0と結ばれる点の集合をEEとする。EEは空でなく開集合である。x∉Ex\notin Eならば、Bsx(x)B_{s_x}(x)はEEと交わらない。交われば鎖へBsx(x)B_{s_x}(x)を付け加えてx∈Ex\in Eとなるためである。従ってΩ∖E\Omega\setminus Eも開集合であり、連結性からE=ΩE=\Omegaである。

KKを有限個の球Bsaj(aj)B_{s_{a_j}}(a_j)で覆い、各球をa0a_0と結ぶ有限の鎖を選ぶ。交点を介して球内の比較を反復すると、有限な定数A≥1A\geq1が存在して、任意のx∈Kx\in Kに対し

A−1u(a0)≤u(x)≤Au(a0)A^{-1}u(a_0)\leq u(x)\leq Au(a_0)

となる。従ってu(x)≤A2u(y)u(x)\leq A^2u(y)が任意のx,y∈Kx,y\in Kに対して成り立ち、CK,Ω=A2C_{K,\Omega}=A^2とすればよい。▨

系 2.2 (Liouville の定理).n≥1n\geq1とする。このとき、Rn\R^n上の実数値調和関数uuが上に有界または下に有界ならば、uuは定数である。Rn\R^n上の複素数値調和関数vvがsup⁡Rn∣v∣<∞\sup_{\R^n}|v|<\inftyを満たす場合も、vvは定数である。

証明.n≥2n\geq2とし、u≥Mu\geq Mとする。w=u−Mw=u-Mは非負の調和関数である。x,y∈Rnx,y\in\R^nを固定し、x,y∈Br(0)x,y\in B_r(0)となるr>0r>0をとる。任意のR>rR>rに対して定理 2.1より

w(x)≤(R+rR−r)nw(y)w(x)\leq\left(\frac{R+r}{R-r}\right)^nw(y)

である。R→∞R\to\inftyとすればw(x)≤w(y)w(x)\leq w(y)となる。x,yx,yを交換すれば等号を得る。上に有界な場合は−u-uに同じ議論を適用する。n=1n=1ではu′′=0u''=0よりuuはアフィン関数であり、一側有界性からその傾きは零である。複素数値の場合はvvの実部と虚部がともに有界な実数値調和関数なので、それぞれ定数である。▨

例 2.3.n≥2n\geq2とし、u(x)=x1u(x)=x_1とおく。このとき、Δu=0\Delta u=0であるが、任意のr>0r>0に対しsup⁡Br(0)u=r\sup_{B_r(0)}u=r、inf⁡Br(0)u=−r\inf_{B_r(0)}u=-rである。従って、非負性を外した Harnack 不等式は成り立たない。またuuはRn\R^n全体で調和的で非定数であるから、一側有界性を外した Liouville の結論も成り立たない。

3 孤立特異点の除去

定理 3.1 (除去可能特異点定理).n≥2n\geq2、Ω⊆Rn\Omega\subseteq\R^nを開集合、a∈Ωa\in\Omegaとする。このとき、u ⁣:Ω∖{a}→Cu\colon\Omega\setminus\{a\}\to\Cが調和関数であり、x→ax\to aのとき

u(x)={o(∣x−a∣2−n),n≥3,o(log⁡(1/∣x−a∣)),n=2u(x)= \begin{cases} o(|x-a|^{2-n}),&n\geq3,\\ o(\log(1/|x-a|)),&n=2 \end{cases}

を満たすならば、uuはΩ\Omega上の調和関数に一意に延長される。特に、aaのある穿孔近傍で有界な調和関数の特異点は除去可能である。

証明. 実数値の場合を考え、BR(a)‾⊆Ω\overline{B_R(a)}\subseteq\OmegaとなるR>0R>0をとる。§E15.29 定理 5.1により、境界値u∣∂BR(a)u|_{\partial B_R(a)}の調和拡張h∈C∞(BR(a))∩C(BR(a)‾)h\in C^\infty(B_R(a))\cap C(\overline{B_R(a)})が存在する。w=u−hw=u-hとおくと、wwは穿孔球上で調和的であり、外側の球面上で零である。hhは有界であるから、wwもuuと同じ little-o 条件を満たす。

§E15.6 定理 2.2の Newton 核を用い、σn−1=∣Sn−1∣\sigma_{n-1}=|S^{n-1}|、0<s≤R0<s\leq Rに対して

q(s)=Φn(se1)−Φn(Re1)={s2−n−R2−n(n−2)σn−1,n≥3,12πlog⁡Rs,n=2q(s)=\Phi_n(se_1)-\Phi_n(Re_1) =\begin{cases} \displaystyle\frac{s^{2-n}-R^{2-n}}{(n-2)\sigma_{n-1}},&n\geq3,\\[4pt] \displaystyle\frac{1}{2\pi}\log\frac{R}{s},&n=2 \end{cases}

とおく。q(∣x−a∣)q(|x-a|)は穿孔球上で調和的で、0<∣x−a∣<R0<|x-a|<Rでは正である。任意のε>0\varepsilon>0に対し、little-o 条件から、十分小さい任意のδ>0\delta>0について

∣w(x)∣≤εq(δ)(∣x−a∣=δ)|w(x)|\leq\varepsilon q(\delta)\qquad(|x-a|=\delta)

が成り立つ。滑らかな有界領域BR(a)∖Bδ(a)‾B_R(a)\setminus\overline{B_\delta(a)}上で−Δ-\Deltaに§E15.18 系 1.2を適用すると、内外両方の境界で成立する±w≤εq(∣x−a∣)\pm w\leq\varepsilon q(|x-a|)が領域内部でも成立する。

0<∣x−a∣<R0<|x-a|<Rを固定してδ↓0\delta\downarrow0とすれば、∣w(x)∣≤εq(∣x−a∣)|w(x)|\leq\varepsilon q(|x-a|)を得る。次にε↓0\varepsilon\downarrow0とするとw(x)=0w(x)=0である。従ってu=hu=hが穿孔球上で成立し、u(a)=h(a)u(a)=h(a)と定めれば調和拡張を得る。複素数値の場合は実部と虚部に同じ議論を適用する。二つの拡張はaa以外で一致するため、連続性からaaでも一致する。n≥3n\geq3では∣x−a∣2−n→∞|x-a|^{2-n}\to\infty、n=2n=2ではlog⁡(1/∣x−a∣)→∞\log(1/|x-a|)\to\inftyであるから、有界な関数は各次元の little-o 条件を満たす。▨

例 3.2.n≥2n\geq2とする。このとき、§E15.6 定理 2.2の基本解u(x)=Φn(x−a)u(x)=\Phi_n(x-a)はRn∖{a}\R^n\setminus\{a\}上で調和的であり、x→ax\to aのとき、n≥3n\geq3ではO(∣x−a∣2−n)O(|x-a|^{2-n})、n=2n=2ではO(log⁡(1/∣x−a∣))O(\log(1/|x-a|))である。しかしu(x)→+∞u(x)\to+\inftyであるから、連続に延長することはできない。従って定理 3.1の little-o 条件を big-O 条件へ弱めることはできない。

Φn\Phi_nは局所可積分であるが、§E15.6 定理 2.2を平行移動すると、超関数の意味で

Δu=−δa\Delta u=-\delta_a

となる。穿孔領域で調和的かつ特異点の近傍で局所可積分であることだけから、全領域で超関数的に調和的であるとは結論することができない。

前提記事