楕円型方程式の弱解は、方程式を部分積分した形でしか満たさない「ゆるい」解です。たとえば板の温度分布を決める方程式で、ゆるい意味でしか方程式を満たさない解でも、内部では滑らかに変わると言えるのでしょうか。この記事は、その問いに答える道具として、差分商による L 2 L^2 L 2 の評価と、積分で Hölder 連続性をとらえる Campanato の方法を順に見ます。読み終えると、弱解が 2 回微分できて方程式が各点で成り立つ条件と、その証明の仕組みが分かります。
係数が C 1 C^1 C 1 級なら、L u = f ∈ L 2 Lu=f\in L^2 Lu = f ∈ L 2 の弱解 u ∈ H 1 u\in H^1 u ∈ H 1 は内部で H 2 H^2 H 2 に入ります(内部 H 2 H^2 H 2 正則性 )。証明は、方程式に差分商を施して得られるエネルギー評価と、「差分商が一様に有界なら弱微分が存在する」という補題です。
凸領域では、u = 0 u=0 u = 0 を境界値とする関数について ∥ D 2 u ∥ L 2 ≤ ∥ Δ u ∥ L 2 \|D^2u\|_{L^2}\le\|\Delta u\|_{L^2} ∥ D 2 u ∥ L 2 ≤ ∥Δ u ∥ L 2 が成り立ちます。部分積分の境界項が境界の曲がり方(凸性)で符号を持つからです(Grisvard の評価 の核心)。
関数が Hölder 連続(2 点の値の差が距離の α \alpha α 乗で抑えられること)であることは、球の上での平均からのずれ ∫ B ρ ∣ g − g ρ ∣ 2 \int_{B_\rho}|g-g_\rho|^2 ∫ B ρ ∣ g − g ρ ∣ 2 が ρ n + 2 α \rho^{n+2\alpha} ρ n + 2 α で減衰することと同値です(Campanato の定理 )。調和関数 との比較と反復補題でこの減衰を示すと、Δ u = div F \Delta u=\operatorname{div}F Δ u = div F 、F ∈ C 0 , α F\in C^{0,\alpha} F ∈ C 0 , α なら ∇ u ∈ C 0 , α \nabla u\in C^{0,\alpha} ∇ u ∈ C 0 , α 、したがって Δ u = f ∈ C 0 , α \Delta u=f\in C^{0,\alpha} Δ u = f ∈ C 0 , α なら D 2 u ∈ C 0 , α D^2u\in C^{0,\alpha} D 2 u ∈ C 0 , α です(Schauder 評価 )。
c ≥ 0 c\ge0 c ≥ 0 の発散形作用素では、境界で u ≤ 0 u\le0 u ≤ 0 (u + ∈ H 0 1 u^+\in H^1_0 u + ∈ H 0 1 )かつ L u ≤ 0 Lu\le0 Lu ≤ 0 の弱解は u ≤ 0 u\le0 u ≤ 0 です(弱最大値原理 )。
− Δ + V -\Delta+V − Δ + V が正の(優)解 ψ \psi ψ をもてば、二次形式は E E E 以上です(Allegretto–Piepenbrink の定理 )。基底状態変換 φ = ψ g \varphi=\psi g φ = ψ g (正の解 ψ \psi ψ で割って、二次形式を非負の項と残りの項に書き直す変換)による恒等式から出ます。Hardy の不等式はその一例です。
この記事で使う事前知識
弱解 : 方程式を試験関数との積分に直した形で満たす解のこと。→ 弱解
Sobolev 空間と弱微分 : 微分を積分の部分積分で定義し直し、2 乗可積分な弱微分をもつ関数の空間のこと。→ Sobolev 空間
調和関数 : Δ u = 0 \Delta u=0 Δ u = 0 を満たす関数のこと。球の上の平均が中心の値に等しくなる。→ 調和関数
弱収束 : すべての元との内積が収束する、という弱い意味の収束。→ 弱収束
凸関数 : グラフの弦が常にグラフの上側にある関数。2 階微分の行列は半正定値になる。→ 凸関数
一様楕円型 : 係数行列の固有値が定数 θ > 0 \theta\gt0 θ > 0 以上で下から抑えられる、という方程式の条件(第 3 節で使う)。
前章では、Hilbert 空間の方法で弱解 u ∈ H 0 1 u\in H^1_0 u ∈ H 0 1 を作りました。しかし弱解は、方程式 L u = f Lu=f Lu = f を部分積分した形でしか満たしていません。2 階導関数が存在して方程式がほとんど至るところ成り立つか、さらに古典解になるかは、別に示す必要があります。これが正則性 の問題です。
正則性は応用の上でも欠かせません。Dirichlet ラプラシアンの固有関数が滑らかであること、自己共役作用素の定義域が H 2 H^2 H 2 と一致すること(第 6 章や Schrödinger 作用素の理論で使う)、非線形方程式で弱解から古典解に戻ることは、すべてこの章の評価に依存します。
正則性の理論には 2 つの系統があります。
L 2 L^2 L 2 理論 (第 3・4 節): エネルギー評価に差分商を組み合わせ、f ∈ H k f\in H^k f ∈ H k なら u ∈ H k + 2 u\in H^{k+2} u ∈ H k + 2 を示します。Hilbert 空間の方法と相性がよい。
Hölder 理論 (第 5 節): f ∈ C k , α f\in C^{k,\alpha} f ∈ C k , α なら u ∈ C k + 2 , α u\in C^{k+2,\alpha} u ∈ C k + 2 , α を示します(Schauder 評価)。f ∈ C 0 f\in C^0 f ∈ C 0 では u ∈ C 2 u\in C^2 u ∈ C 2 とは限らないので、Hölder 空間が「ちょうど 2 階分」の得をする正しい尺度です。Campanato の方法は、Hölder 連続性を積分の言葉で言い換えることで、L 2 L^2 L 2 のエネルギー評価だけから Hölder 評価を引き出します。
以下 Ω ⊂ R n \Omega\subset\mathbb{R}^n Ω ⊂ R n は開集合、関数は実数値です。V ⋐ W V\Subset W V ⋐ W は V ˉ \bar V V ˉ がコンパクトで W W W に含まれることを表します。e k e_k e k を第 k k k 単位ベクトル、h ≠ 0 h\ne0 h = 0 について差分商 を
D k h u ( x ) = u ( x + h e k ) − u ( x ) h D_k^hu(x)=\frac{u(x+he_k)-u(x)}h D k h u ( x ) = h u ( x + h e k ) − u ( x )
と書きます(dist ( x , ∂ Ω ) > ∣ h ∣ \operatorname{dist}(x,\partial\Omega)\gt|h| dist ( x , ∂ Ω ) > ∣ h ∣ で定義される)。v h ( x ) = v ( x + h e k ) v^h(x)=v(x+he_k) v h ( x ) = v ( x + h e k ) とおくと、積の差分商は D k h ( v w ) = v h D k h w + ( D k h v ) w D_k^h(vw)=v^hD_k^hw+(D_k^hv)\,w D k h ( v w ) = v h D k h w + ( D k h v ) w で、平行移動の不変性から、片方の台がコンパクトなら差分の部分積分 ∫ v D k h w d x = − ∫ ( D k − h v ) w d x \int v\,D_k^hw\,dx=-\int(D_k^{-h}v)\,w\,dx ∫ v D k h w d x = − ∫ ( D k − h v ) w d x が成り立ちます。
補題 2.1 (局所的な軟化 ) ✓
u ∈ H 1 ( Ω ) u\in H^1(\Omega) u ∈ H 1 ( Ω ) 、W ⋐ Ω W\Subset\Omega W ⋐ Ω とする。ε < dist ( W , ∂ Ω ) \varepsilon\lt\operatorname{dist}(W,\partial\Omega) ε < dist ( W , ∂ Ω ) なら W W W 上で ∂ k ( ρ ε ∗ u ) = ρ ε ∗ ∂ k u \partial_k(\rho_\varepsilon*u)=\rho_\varepsilon*\partial_ku ∂ k ( ρ ε ∗ u ) = ρ ε ∗ ∂ k u で、ε → 0 \varepsilon\to0 ε → 0 のとき ρ ε ∗ u → u \rho_\varepsilon*u\to u ρ ε ∗ u → u 、∇ ( ρ ε ∗ u ) → ∇ u \nabla(\rho_\varepsilon*u)\to\nabla u ∇ ( ρ ε ∗ u ) → ∇ u が L 2 ( W ) L^2(W) L 2 ( W ) で成り立つ。特に、台がコンパクトな u ∈ H 1 ( Ω ) u\in H^1(\Omega) u ∈ H 1 ( Ω ) は H 0 1 ( Ω ) H^1_0(\Omega) H 0 1 ( Ω ) に属する。
証明(補題 2.1) x ∈ W x\in W x ∈ W なら y ↦ ρ ε ( x − y ) y\mapsto\rho_\varepsilon(x-y) y ↦ ρ ε ( x − y ) は C c ∞ ( Ω ) C_c^\infty(\Omega) C c ∞ ( Ω ) に属するので、弱微分の定義(定義 2.1[Sobolev 空間] )から
∂ k ( ρ ε ∗ u ) ( x ) = ∫ ∂ x k ρ ε ( x − y ) u ( y ) d y = − ∫ ∂ y k [ ρ ε ( x − y ) ] u ( y ) d y = ∫ ρ ε ( x − y ) ∂ k u ( y ) d y \partial_k(\rho_\varepsilon*u)(x)=\int\partial_{x_k}\rho_\varepsilon(x-y)\,u(y)\,dy=-\int\partial_{y_k}\bigl[\rho_\varepsilon(x-y)\bigr]u(y)\,dy=\int\rho_\varepsilon(x-y)\,\partial_ku(y)\,dy ∂ k ( ρ ε ∗ u ) ( x ) = ∫ ∂ x k ρ ε ( x − y ) u ( y ) d y = − ∫ ∂ y k [ ρ ε ( x − y ) ] u ( y ) d y = ∫ ρ ε ( x − y ) ∂ k u ( y ) d y です。u u u 、∂ k u \partial_ku ∂ k u を Ω \Omega Ω の外で 0 0 0 とおいた L 2 ( R n ) L^2(\mathbb{R}^n) L 2 ( R n ) の関数に 補題 2.2[Schwartz 空間と緩増加超関数] の 1 を使えば収束が出ます。台が K ⋐ Ω K\Subset\Omega K ⋐ Ω なら、ε \varepsilon ε が小さいとき ρ ε ∗ u ∈ C c ∞ ( Ω ) \rho_\varepsilon*u\in C_c^\infty(\Omega) ρ ε ∗ u ∈ C c ∞ ( Ω ) で、W W W を K K K の近傍にとれば H 1 ( Ω ) H^1(\Omega) H 1 ( Ω ) で u u u に収束します。
∎
補題 2.2 (差分商と弱微分 ) ✓
u ∈ H 1 ( Ω ) u\in H^1(\Omega) u ∈ H 1 ( Ω ) 、V ⋐ Ω V\Subset\Omega V ⋐ Ω 、0 < ∣ h ∣ < dist ( V , ∂ Ω ) 0\lt|h|\lt\operatorname{dist}(V,\partial\Omega) 0 < ∣ h ∣ < dist ( V , ∂ Ω ) なら ∥ D k h u ∥ L 2 ( V ) ≤ ∥ ∂ k u ∥ L 2 ( Ω ) \|D_k^hu\|_{L^2(V)}\le\|\partial_ku\|_{L^2(\Omega)} ∥ D k h u ∥ L 2 ( V ) ≤ ∥ ∂ k u ∥ L 2 ( Ω ) 。
u ∈ L 2 ( Ω ) u\in L^2(\Omega) u ∈ L 2 ( Ω ) 、V ⋐ Ω V\Subset\Omega V ⋐ Ω とし、ある C C C と δ > 0 \delta\gt0 δ > 0 について、0 < ∣ h ∣ < δ 0\lt|h|\lt\delta 0 < ∣ h ∣ < δ で ∥ D k h u ∥ L 2 ( V ) ≤ C \|D_k^hu\|_{L^2(V)}\le C ∥ D k h u ∥ L 2 ( V ) ≤ C とする。このとき V V V 上で弱微分 ∂ k u ∈ L 2 ( V ) \partial_ku\in L^2(V) ∂ k u ∈ L 2 ( V ) が存在し、∥ ∂ k u ∥ L 2 ( V ) ≤ C \|\partial_ku\|_{L^2(V)}\le C ∥ ∂ k u ∥ L 2 ( V ) ≤ C 。
証明(補題 2.2)
まず u u u が V V V の ∣ h ∣ |h| ∣ h ∣ 近傍で C 1 C^1 C 1 なら、D k h u ( x ) = ∫ 0 1 ∂ k u ( x + t h e k ) d t D_k^hu(x)=\int_0^1\partial_ku(x+the_k)\,dt D k h u ( x ) = ∫ 0 1 ∂ k u ( x + t h e k ) d t なので、Cauchy–Schwarz の不等式と Fubini の定理から ∫ V ∣ D k h u ∣ 2 ≤ ∫ 0 1 ∫ V ∣ ∂ k u ( x + t h e k ) ∣ 2 d x d t ≤ ∥ ∂ k u ∥ L 2 ( Ω ) 2 \int_V|D_k^hu|^2\le\int_0^1\int_V|\partial_ku(x+the_k)|^2dx\,dt\le\|\partial_ku\|_{L^2(\Omega)}^2 ∫ V ∣ D k h u ∣ 2 ≤ ∫ 0 1 ∫ V ∣ ∂ k u ( x + t h e k ) ∣ 2 d x d t ≤ ∥ ∂ k u ∥ L 2 ( Ω ) 2 です。一般の u u u は、補題 2.1 で ρ ε ∗ u \rho_\varepsilon*u ρ ε ∗ u に置き換えて ε → 0 \varepsilon\to0 ε → 0 とします。
h j → 0 h_j\to0 h j → 0 を ∥ D k h j u ∥ L 2 ( V ) ≤ C \|D_k^{h_j}u\|_{L^2(V)}\le C ∥ D k h j u ∥ L 2 ( V ) ≤ C となるようにとると、定理 5.5[Hilbert 空間の基礎] により部分列で D k h j u ⇀ w D_k^{h_j}u\rightharpoonup w D k h j u ⇀ w (L 2 ( V ) L^2(V) L 2 ( V ) で弱収束、∥ w ∥ ≤ C \|w\|\le C ∥ w ∥ ≤ C )です。φ ∈ C c ∞ ( V ) \varphi\in C_c^\infty(V) φ ∈ C c ∞ ( V ) について、差分の部分積分で ∫ u D k − h j φ = − ∫ ( D k h j u ) φ \int u\,D_k^{-h_j}\varphi=-\int(D_k^{h_j}u)\varphi ∫ u D k − h j φ = − ∫ ( D k h j u ) φ で、D k − h j φ → ∂ k φ D_k^{-h_j}\varphi\to\partial_k\varphi D k − h j φ → ∂ k φ は一様収束(台は一定のコンパクト集合)なので、極限で ∫ u ∂ k φ = − ∫ w φ \int u\,\partial_k\varphi=-\int w\varphi ∫ u ∂ k φ = − ∫ w φ です。よって w = ∂ k u w=\partial_ku w = ∂ k u です。
∎
前章と同じく L u = − ∑ ∂ j ( a i j ∂ i u ) + ∑ b i ∂ i u + c u Lu=-\sum\partial_j(a^{ij}\partial_iu)+\sum b^i\partial_iu+cu Lu = − ∑ ∂ j ( a ij ∂ i u ) + ∑ b i ∂ i u + c u 、B B B をその双線形形式とし、L L L は一様楕円型(定数 θ \theta θ )とします。
定理 3.1 (内部 H² 正則性 ) ✓
a i j ∈ C 1 ( Ω ˉ ) a^{ij}\in C^1(\bar\Omega) a ij ∈ C 1 ( Ω ˉ ) 、b i , c ∈ L ∞ ( Ω ) b^i,c\in L^\infty(\Omega) b i , c ∈ L ∞ ( Ω ) 、f ∈ L 2 ( Ω ) f\in L^2(\Omega) f ∈ L 2 ( Ω ) とし、u ∈ H 1 ( Ω ) u\in H^1(\Omega) u ∈ H 1 ( Ω ) がすべての v ∈ H 0 1 ( Ω ) v\in H^1_0(\Omega) v ∈ H 0 1 ( Ω ) で B [ u , v ] = ⟨ f , v ⟩ B[u,v]=\langle f,v\rangle B [ u , v ] = ⟨ f , v ⟩ を満たすとする(u u u の境界値は問わない)。このとき V ⋐ Ω V\Subset\Omega V ⋐ Ω について u ∈ H 2 ( V ) u\in H^2(V) u ∈ H 2 ( V ) で、V , Ω V,\Omega V , Ω と係数だけで決まる C C C について
∥ u ∥ H 2 ( V ) ≤ C ( ∥ f ∥ L 2 ( Ω ) + ∥ u ∥ H 1 ( Ω ) ) . \|u\|_{H^2(V)}\le C\bigl(\|f\|_{L^2(\Omega)}+\|u\|_{H^1(\Omega)}\bigr). ∥ u ∥ H 2 ( V ) ≤ C ( ∥ f ∥ L 2 ( Ω ) + ∥ u ∥ H 1 ( Ω ) ) . さらに L u = f Lu=f Lu = f が V V V 上ほとんど至るところ成り立つ。
証明の方針は次のとおりです。弱解はまだ 2 階微分をもつか分からないので、微分の代わりに差分商を方程式に施し、差分商が h h h によらず一様に有界であることを示します。すると前節の補題の 2 により弱微分が存在します。試験関数に切り落とし関数 ζ 2 \zeta^2 ζ 2 を掛けるのは、台を W W W の内側に収めて、部分積分で領域の境界の項が出ないようにするためです。
証明(定理 3.1) V ⋐ W ⋐ Ω V\Subset W\Subset\Omega V ⋐ W ⋐ Ω をとり、ζ ∈ C c ∞ ( W ) \zeta\in C_c^\infty(W) ζ ∈ C c ∞ ( W ) 、0 ≤ ζ ≤ 1 0\le\zeta\le1 0 ≤ ζ ≤ 1 、V V V 上 ζ = 1 \zeta=1 ζ = 1 とします(補題 2.3[Schwartz 空間と緩増加超関数] )。1 階と 0 階の項を右辺に移し、f ~ = f − ∑ b i ∂ i u − c u ∈ L 2 \tilde f=f-\sum b^i\partial_iu-cu\in L^2 f ~ = f − ∑ b i ∂ i u − c u ∈ L 2 とおくと、∑ i , j ∫ a i j ∂ i u ∂ j v = ∫ f ~ v \sum_{i,j}\int a^{ij}\partial_iu\,\partial_jv=\int\tilde fv ∑ i , j ∫ a ij ∂ i u ∂ j v = ∫ f ~ v (v ∈ H 0 1 v\in H^1_0 v ∈ H 0 1 )です。
∣ h ∣ |h| ∣ h ∣ を小さくとり、v = − D k − h ( ζ 2 D k h u ) v=-D_k^{-h}(\zeta^2D_k^hu) v = − D k − h ( ζ 2 D k h u ) とおきます。ζ 2 D k h u \zeta^2D_k^hu ζ 2 D k h u は台がコンパクトな H 1 H^1 H 1 関数なので、v ∈ H 0 1 ( Ω ) v\in H^1_0(\Omega) v ∈ H 0 1 ( Ω ) です(補題 2.1 )。
左辺。 差分の部分積分と積の差分商の公式から
∑ ∫ a i j ∂ i u ∂ j v = ∑ ∫ D k h ( a i j ∂ i u ) ∂ j ( ζ 2 D k h u ) = A 1 + A 2 , \sum\int a^{ij}\partial_iu\,\partial_jv=\sum\int D_k^h(a^{ij}\partial_iu)\,\partial_j(\zeta^2D_k^hu)=A_1+A_2, ∑ ∫ a ij ∂ i u ∂ j v = ∑ ∫ D k h ( a ij ∂ i u ) ∂ j ( ζ 2 D k h u ) = A 1 + A 2 , A 1 = ∑ ∫ ζ 2 a i j , h D k h ∂ i u D k h ∂ j u , A 2 = ∑ ∫ [ ( D k h a i j ) ∂ i u ζ 2 D k h ∂ j u + D k h ( a i j ∂ i u ) 2 ζ ∂ j ζ D k h u ] A_1=\sum\int\zeta^2a^{ij,h}\,D_k^h\partial_iu\,D_k^h\partial_ju,\qquad A_2=\sum\int\Bigl[(D_k^ha^{ij})\partial_iu\,\zeta^2D_k^h\partial_ju+D_k^h(a^{ij}\partial_iu)\,2\zeta\partial_j\zeta\,D_k^hu\Bigr] A 1 = ∑ ∫ ζ 2 a ij , h D k h ∂ i u D k h ∂ j u , A 2 = ∑ ∫ [ ( D k h a ij ) ∂ i u ζ 2 D k h ∂ j u + D k h ( a ij ∂ i u ) 2 ζ ∂ j ζ D k h u ] です(a i j , h = a i j ( ⋅ + h e k ) a^{ij,h}=a^{ij}(\cdot+he_k) a ij , h = a ij ( ⋅ + h e k ) )。楕円性から A 1 ≥ θ ∫ ζ 2 ∣ D k h ∇ u ∣ 2 A_1\ge\theta\int\zeta^2|D_k^h\nabla u|^2 A 1 ≥ θ ∫ ζ 2 ∣ D k h ∇ u ∣ 2 です。∣ D k h a i j ∣ ≤ ∥ ∇ a i j ∥ ∞ |D_k^ha^{ij}|\le\|\nabla a^{ij}\|_\infty ∣ D k h a ij ∣ ≤ ∥∇ a ij ∥ ∞ と 補題 2.2 の 1(∫ W ∣ D k h u ∣ 2 ≤ ∥ ∇ u ∥ 2 \int_W|D_k^hu|^2\le\|\nabla u\|^2 ∫ W ∣ D k h u ∣ 2 ≤ ∥∇ u ∥ 2 )を使い、a b ≤ ε a 2 + 1 4 ε b 2 ab\le\varepsilon a^2+\frac1{4\varepsilon}b^2 ab ≤ ε a 2 + 4 ε 1 b 2 を a = ζ ∣ D k h ∇ u ∣ a=\zeta|D_k^h\nabla u| a = ζ ∣ D k h ∇ u ∣ に対して使うと、∣ A 2 ∣ ≤ θ 2 ∫ ζ 2 ∣ D k h ∇ u ∣ 2 + C ∥ u ∥ H 1 ( Ω ) 2 |A_2|\le\frac\theta2\int\zeta^2|D_k^h\nabla u|^2+C\|u\|_{H^1(\Omega)}^2 ∣ A 2 ∣ ≤ 2 θ ∫ ζ 2 ∣ D k h ∇ u ∣ 2 + C ∥ u ∥ H 1 ( Ω ) 2 です。
右辺。 補題 2.2 の 1 を ζ 2 D k h u \zeta^2D_k^hu ζ 2 D k h u に使うと ∥ v ∥ ≤ ∥ ∂ k ( ζ 2 D k h u ) ∥ ≤ C ( ∥ ζ D k h ∇ u ∥ + ∥ u ∥ H 1 ) \|v\|\le\|\partial_k(\zeta^2D_k^hu)\|\le C\bigl(\|\zeta D_k^h\nabla u\|+\|u\|_{H^1}\bigr) ∥ v ∥ ≤ ∥ ∂ k ( ζ 2 D k h u ) ∥ ≤ C ( ∥ ζ D k h ∇ u ∥ + ∥ u ∥ H 1 ) なので、
∣ ∫ f ~ v ∣ ≤ θ 4 ∫ ζ 2 ∣ D k h ∇ u ∣ 2 + C ( ∥ f ∥ 2 + ∥ u ∥ H 1 ( Ω ) 2 ) \Bigl|\int\tilde fv\Bigr|\le\frac\theta4\int\zeta^2|D_k^h\nabla u|^2+C\bigl(\|f\|^2+\|u\|_{H^1(\Omega)}^2\bigr) ∫ f ~ v ≤ 4 θ ∫ ζ 2 ∣ D k h ∇ u ∣ 2 + C ( ∥ f ∥ 2 + ∥ u ∥ H 1 ( Ω ) 2 ) です(∥ f ~ ∥ ≤ ∥ f ∥ + C ∥ u ∥ H 1 \|\tilde f\|\le\|f\|+C\|u\|_{H^1} ∥ f ~ ∥ ≤ ∥ f ∥ + C ∥ u ∥ H 1 )。
結論。 合わせて ∫ V ∣ D k h ∇ u ∣ 2 ≤ ∫ ζ 2 ∣ D k h ∇ u ∣ 2 ≤ C ( ∥ f ∥ 2 + ∥ u ∥ H 1 ( Ω ) 2 ) \int_V|D_k^h\nabla u|^2\le\int\zeta^2|D_k^h\nabla u|^2\le C(\|f\|^2+\|u\|_{H^1(\Omega)}^2) ∫ V ∣ D k h ∇ u ∣ 2 ≤ ∫ ζ 2 ∣ D k h ∇ u ∣ 2 ≤ C ( ∥ f ∥ 2 + ∥ u ∥ H 1 ( Ω ) 2 ) が、小さなすべての h h h と k k k で成り立ちます。補題 2.2 の 2 を ∂ i u \partial_iu ∂ i u に使えば ∂ k ∂ i u ∈ L 2 ( V ) \partial_k\partial_iu\in L^2(V) ∂ k ∂ i u ∈ L 2 ( V ) と評価が出ます。最後に、v ∈ C c ∞ ( V ) v\in C_c^\infty(V) v ∈ C c ∞ ( V ) について弱形式を部分積分し直す(a i j ∂ i u ∈ H 1 ( V ) a^{ij}\partial_iu\in H^1(V) a ij ∂ i u ∈ H 1 ( V ) なので積の微分が使える)と ∫ ( L u − f ) v = 0 \int(Lu-f)v=0 ∫ ( Lu − f ) v = 0 で、補題 2.2[Schwartz 空間と緩増加超関数] の 3 から L u = f Lu=f Lu = f a.e. です。
∎
同じ議論を、方程式を x k x_k x k で微分したもの(∂ k u \partial_ku ∂ k u の満たす方程式)に繰り返すと、a i j , b i , c ∈ C m + 1 a^{ij},b^i,c\in C^{m+1} a ij , b i , c ∈ C m + 1 、f ∈ H m f\in H^m f ∈ H m なら u ∈ H l o c m + 2 u\in H^{m+2}_{\mathrm{loc}} u ∈ H loc m + 2 が帰納的に出ます。係数と f f f が C ∞ C^\infty C ∞ なら、Sobolev 埋め込み(定理 4.1[Sobolev 空間] )と合わせて u ∈ C ∞ ( Ω ) u\in C^\infty(\Omega) u ∈ C ∞ ( Ω ) です。とくに Dirichlet ラプラシアンの固有関数は Ω \Omega Ω の内部で C ∞ C^\infty C ∞ です。
境界まで含めた H 2 H^2 H 2 正則性には境界の滑らかさが要ります。境界を局所的に平らにして接線方向にだけ差分商を使い、法線方向の 2 階導関数は方程式から求めるのが標準的な方法です(Evans 6.3.2 節)。この記事では証明しません。次節では、凸領域に特有の別の機構を見ます。
x′ p = 0 ν = −e_n Ω x_n = ψ(x′)(凸) 凸領域の境界の点 p での座標。境界は凸関数 x_n = ψ(x′) のグラフで、領域はその上側にある。D²ψ(0) のトレースが境界項の符号を決める
定理 4.1 (凸領域の H² 評価(Grisvard) ) ✓
Ω \Omega Ω を滑らかな境界をもつ有界な凸領域とする。u ∈ C 3 ( Ω ˉ ) u\in C^3(\bar\Omega) u ∈ C 3 ( Ω ˉ ) が ∂ Ω \partial\Omega ∂ Ω 上で u = 0 u=0 u = 0 を満たすなら
∫ Ω ∣ D 2 u ∣ 2 d x ≤ ∫ Ω ( Δ u ) 2 d x , ∣ D 2 u ∣ 2 = ∑ i , j ( ∂ i ∂ j u ) 2 . \int_\Omega|D^2u|^2dx\le\int_\Omega(\Delta u)^2dx,\qquad|D^2u|^2=\sum_{i,j}(\partial_i\partial_ju)^2 . ∫ Ω ∣ D 2 u ∣ 2 d x ≤ ∫ Ω ( Δ u ) 2 d x , ∣ D 2 u ∣ 2 = i , j ∑ ( ∂ i ∂ j u ) 2 .
証明の方針は、( Δ u ) 2 − ∣ D 2 u ∣ 2 (\Delta u)^2-|D^2u|^2 ( Δ u ) 2 − ∣ D 2 u ∣ 2 をあるベクトル場の発散として書き、発散定理で境界の積分に直して、その符号を調べることです。符号を決めるのが領域の凸性で、境界の曲がり方が Δ u \Delta u Δ u と D 2 u D^2u D 2 u の差になって現れます。
証明(定理 4.1) 恒等式。 u i = ∂ i u u_i=\partial_iu u i = ∂ i u と書き、ベクトル場 Y i = u i Δ u − ∑ j u j ∂ j u i Y_i=u_i\,\Delta u-\sum_ju_j\,\partial_ju_i Y i = u i Δ u − ∑ j u j ∂ j u i を考えます。u ∈ C 3 u\in C^3 u ∈ C 3 なので Y ∈ C 1 ( Ω ˉ ) Y\in C^1(\bar\Omega) Y ∈ C 1 ( Ω ˉ ) で、3 階導関数の項が打ち消し合って
div Y = ∑ i , j ( ∂ i u i ∂ j u j − ∂ i u j ∂ j u i ) = ( Δ u ) 2 − ∣ D 2 u ∣ 2 \operatorname{div}Y=\sum_{i,j}\bigl(\partial_iu_i\,\partial_ju_j-\partial_iu_j\,\partial_ju_i\bigr)=(\Delta u)^2-|D^2u|^2 div Y = i , j ∑ ( ∂ i u i ∂ j u j − ∂ i u j ∂ j u i ) = ( Δ u ) 2 − ∣ D 2 u ∣ 2 です(∂ i u j = ∂ j u i \partial_iu_j=\partial_ju_i ∂ i u j = ∂ j u i )。発散定理(補題 2.1[ラプラス方程式と調和関数] )から
∫ Ω ( Δ u ) 2 − ∫ Ω ∣ D 2 u ∣ 2 = ∫ ∂ Ω ( ∂ ν u Δ u − ν ⋅ D 2 u ∇ u ) d S . \int_\Omega(\Delta u)^2-\int_\Omega|D^2u|^2=\int_{\partial\Omega}\Bigl(\partial_\nu u\,\Delta u-\nu\cdot D^2u\,\nabla u\Bigr)dS . ∫ Ω ( Δ u ) 2 − ∫ Ω ∣ D 2 u ∣ 2 = ∫ ∂ Ω ( ∂ ν u Δ u − ν ⋅ D 2 u ∇ u ) d S . 境界項の計算。 ∂ Ω \partial\Omega ∂ Ω 上で u = 0 u=0 u = 0 なので ∇ u \nabla u ∇ u は法線方向を向き、∇ u = ( ∂ ν u ) ν \nabla u=(\partial_\nu u)\nu ∇ u = ( ∂ ν u ) ν です。よって被積分関数は ∂ ν u ( Δ u − ∂ ν 2 u ) \partial_\nu u\,\bigl(\Delta u-\partial_\nu^2u\bigr) ∂ ν u ( Δ u − ∂ ν 2 u ) (∂ ν 2 u = ν ⋅ D 2 u ν \partial_\nu^2u=\nu\cdot D^2u\,\nu ∂ ν 2 u = ν ⋅ D 2 u ν )です。境界の点 p p p を固定し、回転と平行移動で p = 0 p=0 p = 0 、ν ( p ) = − e n \nu(p)=-e_n ν ( p ) = − e n とし、p p p の近くで Ω = { x n > ψ ( x ′ ) } \Omega=\{x_n\gt\psi(x')\} Ω = { x n > ψ ( x ′ )} (ψ ( 0 ) = 0 \psi(0)=0 ψ ( 0 ) = 0 、∇ ψ ( 0 ) = 0 \nabla\psi(0)=0 ∇ ψ ( 0 ) = 0 )と表します。Ω \Omega Ω が凸なので ψ \psi ψ は凸関数で、D 2 ψ ( 0 ) D^2\psi(0) D 2 ψ ( 0 ) は半正定値です。u ( x ′ , ψ ( x ′ ) ) = 0 u(x',\psi(x'))=0 u ( x ′ , ψ ( x ′ )) = 0 を x α x_\alpha x α (α < n \alpha\lt n α < n )で 2 回微分して x ′ = 0 x'=0 x ′ = 0 を代入すると
∂ α 2 u ( 0 ) + ∂ n u ( 0 ) ∂ α 2 ψ ( 0 ) = 0 \partial_\alpha^2u(0)+\partial_nu(0)\,\partial_\alpha^2\psi(0)=0 ∂ α 2 u ( 0 ) + ∂ n u ( 0 ) ∂ α 2 ψ ( 0 ) = 0 です(∇ ψ ( 0 ) = 0 \nabla\psi(0)=0 ∇ ψ ( 0 ) = 0 により、1 階微分を含む項は消える)。α \alpha α について足すと Δ u − ∂ n 2 u = − ∂ n u tr D 2 ψ ( 0 ) \Delta u-\partial_n^2u=-\partial_nu\,\operatorname{tr}D^2\psi(0) Δ u − ∂ n 2 u = − ∂ n u tr D 2 ψ ( 0 ) で、∂ ν u = − ∂ n u \partial_\nu u=-\partial_nu ∂ ν u = − ∂ n u 、∂ ν 2 u = ∂ n 2 u \partial_\nu^2u=\partial_n^2u ∂ ν 2 u = ∂ n 2 u なので、被積分関数は p p p で
∂ ν u ( Δ u − ∂ ν 2 u ) = ( ∂ ν u ) 2 tr D 2 ψ ( 0 ) ≥ 0 \partial_\nu u\bigl(\Delta u-\partial_\nu^2u\bigr)=(\partial_\nu u)^2\operatorname{tr}D^2\psi(0)\ge0 ∂ ν u ( Δ u − ∂ ν 2 u ) = ( ∂ ν u ) 2 tr D 2 ψ ( 0 ) ≥ 0 です。Δ u \Delta u Δ u 、∂ ν u \partial_\nu u ∂ ν u 、ν ⋅ D 2 u ∇ u \nu\cdot D^2u\,\nabla u ν ⋅ D 2 u ∇ u は回転で不変なので、この計算は元の座標でも正しく、境界項は 0 0 0 以上です。
∎
tr D 2 ψ ( 0 ) \operatorname{tr}D^2\psi(0) tr D 2 ψ ( 0 ) は p p p での平均曲率の n − 1 n-1 n − 1 倍です。境界項は「境界が内側に曲がっているほど D 2 u D^2u D 2 u を Δ u \Delta u Δ u より小さくする」ことを言っています。凸でない領域(たとえば内角が π \pi π を超える角をもつ多角形)では、− Δ u = f ∈ L 2 -\Delta u=f\in L^2 − Δ u = f ∈ L 2 の解が H 2 H^2 H 2 に入らないことがあります。
定義 5.1 (Hölder 空間 ) ↗
U ⊂ R n U\subset\mathbb{R}^n U ⊂ R n を開集合、0 < α ≤ 1 0\lt\alpha\le1 0 < α ≤ 1 、k ≥ 0 k\ge0 k ≥ 0 を整数とする。U U U 上の関数 g g g の Hölder 半ノルム を
[ g ] α ; U = sup x ≠ y ∈ U ∣ g ( x ) − g ( y ) ∣ ∣ x − y ∣ α [g]_{\alpha;U}=\sup_{x\ne y\in U}\frac{|g(x)-g(y)|}{|x-y|^\alpha} [ g ] α ; U = x = y ∈ U sup ∣ x − y ∣ α ∣ g ( x ) − g ( y ) ∣ と置く。g ∈ C k ( U ) g\in C^k(U) g ∈ C k ( U ) で、多重指数 ∣ β ∣ ≤ k |\beta|\le k ∣ β ∣ ≤ k のすべての偏導関数 ∂ β g \partial^\beta g ∂ β g が U U U 上有界かつ一様連続(したがって U ˉ \bar U U ˉ へ連続に延びる)であり、∣ β ∣ = k |\beta|=k ∣ β ∣ = k について [ ∂ β g ] α ; U < ∞ [\partial^\beta g]_{\alpha;U}\lt\infty [ ∂ β g ] α ; U < ∞ となるもの全体を C k , α ( U ˉ ) C^{k,\alpha}(\bar U) C k , α ( U ˉ ) と書き、Hölder 空間 という。ノルムは
∥ g ∥ C k , α ( U ˉ ) = ∑ ∣ β ∣ ≤ k sup U ∣ ∂ β g ∣ + ∑ ∣ β ∣ = k [ ∂ β g ] α ; U \|g\|_{C^{k,\alpha}(\bar U)}=\sum_{|\beta|\le k}\sup_U|\partial^\beta g|+\sum_{|\beta|=k}[\partial^\beta g]_{\alpha;U} ∥ g ∥ C k , α ( U ˉ ) = ∣ β ∣ ≤ k ∑ U sup ∣ ∂ β g ∣ + ∣ β ∣ = k ∑ [ ∂ β g ] α ; U とする。α = 1 \alpha=1 α = 1 の C 0 , 1 ( U ˉ ) C^{0,1}(\bar U) C 0 , 1 ( U ˉ ) は有界な Lipschitz 関数の空間である。U U U の各点がある開球 B B B (B ˉ ⊂ U \bar B\subset U B ˉ ⊂ U )を近傍にもち、g ∣ B ∈ C k , α ( B ˉ ) g|_B\in C^{k,\alpha}(\bar B) g ∣ B ∈ C k , α ( B ˉ ) となるとき、g ∈ C k , α ( U ) g\in C^{k,\alpha}(U) g ∈ C k , α ( U ) (局所 Hölder 空間)と書く。
C k , α C^{k,\alpha} C k , α は C k C^k C k と C k + 1 C^{k+1} C k + 1 の間を指数 α \alpha α で細かく刻む尺度です。f ∈ C k , α f\in C^{k,\alpha} f ∈ C k , α から u ∈ C k + 2 , α u\in C^{k+2,\alpha} u ∈ C k + 2 , α を導くのが、この節の Schauder 評価の考え方です。
Hölder 半ノルム [ g ] α ; U [g]_{\alpha;U} [ g ] α ; U は 定義 5.1 のとおりです。g x , ρ = 1 ω n ρ n ∫ B ( x , ρ ) g g_{x,\rho}=\frac1{\omega_n\rho^n}\int_{B(x,\rho)}g g x , ρ = ω n ρ n 1 ∫ B ( x , ρ ) g を球平均とします。値がベクトルの場合は成分ごとに考えます。次の事実を何度も使います: 定数 c c c のうち ∫ B ∣ g − c ∣ 2 \int_{B}|g-c|^2 ∫ B ∣ g − c ∣ 2 を最小にするのは c = g B c=g_B c = g B (球 B B B 上の平均)です(c c c についての 2 次式の平方完成)。特に ∫ B ( x , ρ ) ∣ g − g x , ρ ∣ 2 \int_{B(x,\rho)}|g-g_{x,\rho}|^2 ∫ B ( x , ρ ) ∣ g − g x , ρ ∣ 2 は ρ \rho ρ について単調増加です。
定理 5.2 (Campanato の定理 ) ✓
0 < α ≤ 1 0\lt\alpha\le1 0 < α ≤ 1 、R > 0 R\gt0 R > 0 、g ∈ L 2 ( B ( x 0 , 2 R ) ) g\in L^2(B(x_0,2R)) g ∈ L 2 ( B ( x 0 , 2 R )) とし、ある M M M について
∫ B ( x , ρ ) ∣ g − g x , ρ ∣ 2 d y ≤ M 2 ρ n + 2 α ( x ∈ B ( x 0 , R ) , 0 < ρ ≤ R ) \int_{B(x,\rho)}|g-g_{x,\rho}|^2dy\le M^2\rho^{n+2\alpha}\qquad(x\in B(x_0,R),\ 0\lt\rho\le R) ∫ B ( x , ρ ) ∣ g − g x , ρ ∣ 2 d y ≤ M 2 ρ n + 2 α ( x ∈ B ( x 0 , R ) , 0 < ρ ≤ R ) が成り立つとする。このとき g g g は B ( x 0 , R ) B(x_0,R) B ( x 0 , R ) 上で連続な代表元をもち、[ g ] α ; B ( x 0 , R ) ≤ C ( n , α ) M [g]_{\alpha;B(x_0,R)}\le C(n,\alpha)\,M [ g ] α ; B ( x 0 , R ) ≤ C ( n , α ) M である。
この定理のねらいは、積分の評価から点ごとの評価を取り戻すことです。Schauder 評価で実際に手に入るのは球の上の積分の評価ですが、示したいのは Hölder 連続性(各点の値の差の評価)です。証明は、半径を半分ずつ縮めた球の平均の差を足し合わせて、平均が各点で収束することを示し、その極限の差を評価する流れです。
証明(定理 5.2) 平均の極限。 x ∈ B ( x 0 , R ) x\in B(x_0,R) x ∈ B ( x 0 , R ) 、ρ / 2 ≤ r ≤ ρ ≤ R \rho/2\le r\le\rho\le R ρ /2 ≤ r ≤ ρ ≤ R なら、Cauchy–Schwarz の不等式で
∣ g x , r − g x , ρ ∣ 2 ≤ 1 ω n r n ∫ B ( x , r ) ∣ g − g x , ρ ∣ 2 ≤ 2 n ω n ρ n M 2 ρ n + 2 α |g_{x,r}-g_{x,\rho}|^2\le\frac1{\omega_nr^n}\int_{B(x,r)}|g-g_{x,\rho}|^2\le\frac{2^n}{\omega_n\rho^n}M^2\rho^{n+2\alpha} ∣ g x , r − g x , ρ ∣ 2 ≤ ω n r n 1 ∫ B ( x , r ) ∣ g − g x , ρ ∣ 2 ≤ ω n ρ n 2 n M 2 ρ n + 2 α なので ∣ g x , r − g x , ρ ∣ ≤ c 1 M ρ α |g_{x,r}-g_{x,\rho}|\le c_1M\rho^\alpha ∣ g x , r − g x , ρ ∣ ≤ c 1 M ρ α (c 1 = ( 2 n / ω n ) 1 / 2 c_1=(2^n/\omega_n)^{1/2} c 1 = ( 2 n / ω n ) 1/2 )です。ρ j = 2 − j ρ \rho_j=2^{-j}\rho ρ j = 2 − j ρ について足し合わせると、r ≤ ρ r\le\rho r ≤ ρ のとき ∣ g x , r − g x , ρ ∣ ≤ c 1 M ρ α ∑ j ≥ 0 2 − j α = c 2 M ρ α |g_{x,r}-g_{x,\rho}|\le c_1M\rho^\alpha\sum_{j\ge0}2^{-j\alpha}=c_2M\rho^\alpha ∣ g x , r − g x , ρ ∣ ≤ c 1 M ρ α ∑ j ≥ 0 2 − j α = c 2 M ρ α です。(ρ j + 1 < r ≤ ρ j \rho_{j+1}\lt r\le\rho_j ρ j + 1 < r ≤ ρ j となる j j j をとり、g x , r g_{x,r} g x , r と g x , ρ j g_{x,\rho_j} g x , ρ j 、g x , ρ i g_{x,\rho_i} g x , ρ i と g x , ρ i + 1 g_{x,\rho_{i+1}} g x , ρ i + 1 の差をそれぞれ評価する。)これを ρ = max ( r , r ′ ) \rho=\max(r,r') ρ = max ( r , r ′ ) に使えば ∣ g x , r − g x , r ′ ∣ ≤ 2 c 2 M max ( r , r ′ ) α |g_{x,r}-g_{x,r'}|\le2c_2M\max(r,r')^\alpha ∣ g x , r − g x , r ′ ∣ ≤ 2 c 2 M max ( r , r ′ ) α なので、極限 g ~ ( x ) = lim r → 0 g x , r \tilde g(x)=\lim_{r\to0}g_{x,r} g ~ ( x ) = lim r → 0 g x , r が存在し、∣ g ~ ( x ) − g x , ρ ∣ ≤ c 2 M ρ α |\tilde g(x)-g_{x,\rho}|\le c_2M\rho^\alpha ∣ g ~ ( x ) − g x , ρ ∣ ≤ c 2 M ρ α です。
代表元であること。 上の評価から g ⋅ , r → g ~ g_{\cdot,r}\to\tilde g g ⋅ , r → g ~ は B ( x 0 , R ) B(x_0,R) B ( x 0 , R ) 上一様です。一方 g ⋅ , r = χ r ∗ g g_{\cdot,r}=\chi_r*g g ⋅ , r = χ r ∗ g (χ r = 1 B ( 0 , r ) / ( ω n r n ) \chi_r=\mathbf 1_{B(0,r)}/(\omega_nr^n) χ r = 1 B ( 0 , r ) / ( ω n r n ) )で、補題 2.2[Schwartz 空間と緩増加超関数] の 1 の証明は ρ ε \rho_\varepsilon ρ ε を非負で積分 1 1 1 、台が B ˉ ( 0 , r ) \bar B(0,r) B ˉ ( 0 , r ) の χ r \chi_r χ r に替えてもそのまま通るので、χ r ∗ g → g \chi_r*g\to g χ r ∗ g → g は L 2 ( B ( x 0 , R ) ) L^2(B(x_0,R)) L 2 ( B ( x 0 , R )) で成り立ちます。よって g ~ = g \tilde g=g g ~ = g a.e. です。
Hölder 評価。 x , y ∈ B ( x 0 , R ) x,y\in B(x_0,R) x , y ∈ B ( x 0 , R ) 、ρ = ∣ x − y ∣ ≤ R / 2 \rho=|x-y|\le R/2 ρ = ∣ x − y ∣ ≤ R /2 とします。B ( x , ρ ) ⊂ B ( y , 2 ρ ) B(x,\rho)\subset B(y,2\rho) B ( x , ρ ) ⊂ B ( y , 2 ρ ) なので
∣ g x , ρ − g y , 2 ρ ∣ 2 ≤ 1 ω n ρ n ∫ B ( y , 2 ρ ) ∣ g − g y , 2 ρ ∣ 2 ≤ 2 n + 2 α ω n M 2 ρ 2 α |g_{x,\rho}-g_{y,2\rho}|^2\le\frac1{\omega_n\rho^n}\int_{B(y,2\rho)}|g-g_{y,2\rho}|^2\le\frac{2^{n+2\alpha}}{\omega_n}M^2\rho^{2\alpha} ∣ g x , ρ − g y , 2 ρ ∣ 2 ≤ ω n ρ n 1 ∫ B ( y , 2 ρ ) ∣ g − g y , 2 ρ ∣ 2 ≤ ω n 2 n + 2 α M 2 ρ 2 α で、∣ g ~ ( x ) − g ~ ( y ) ∣ ≤ ∣ g ~ ( x ) − g x , ρ ∣ + ∣ g x , ρ − g y , 2 ρ ∣ + ∣ g y , 2 ρ − g ~ ( y ) ∣ ≤ c 3 M ρ α |\tilde g(x)-\tilde g(y)|\le|\tilde g(x)-g_{x,\rho}|+|g_{x,\rho}-g_{y,2\rho}|+|g_{y,2\rho}-\tilde g(y)|\le c_3M\rho^\alpha ∣ g ~ ( x ) − g ~ ( y ) ∣ ≤ ∣ g ~ ( x ) − g x , ρ ∣ + ∣ g x , ρ − g y , 2 ρ ∣ + ∣ g y , 2 ρ − g ~ ( y ) ∣ ≤ c 3 M ρ α です。∣ x − y ∣ > R / 2 |x-y|\gt R/2 ∣ x − y ∣ > R /2 なら、線分 [ x , y ] ⊂ B ( x 0 , R ) [x,y]\subset B(x_0,R) [ x , y ] ⊂ B ( x 0 , R ) を長さ R / 2 R/2 R /2 以下の 4 つ以下の区間に分けて足せば、4 c 3 M ∣ x − y ∣ α 4c_3M|x-y|^\alpha 4 c 3 M ∣ x − y ∣ α 以下です。
∎
逆に [ g ] α ≤ M [g]_{\alpha}\le M [ g ] α ≤ M なら ∫ B ( x , ρ ) ∣ g − g x , ρ ∣ 2 ≤ ∫ B ( x , ρ ) ∣ g − g ( x ) ∣ 2 ≤ ω n M 2 ρ n + 2 α \int_{B(x,\rho)}|g-g_{x,\rho}|^2\le\int_{B(x,\rho)}|g-g(x)|^2\le\omega_nM^2\rho^{n+2\alpha} ∫ B ( x , ρ ) ∣ g − g x , ρ ∣ 2 ≤ ∫ B ( x , ρ ) ∣ g − g ( x ) ∣ 2 ≤ ω n M 2 ρ n + 2 α なので、この条件は Hölder 連続性と同値です。積分の条件は、エネルギー評価と直接結びつくのが利点です。
補題 5.3 (Weyl の補題 ) ✓
v ∈ H 1 ( B ) v\in H^1(B) v ∈ H 1 ( B ) (B B B は球)がすべての φ ∈ C c ∞ ( B ) \varphi\in C_c^\infty(B) φ ∈ C c ∞ ( B ) で ∫ B ∇ v ⋅ ∇ φ = 0 \int_B\nabla v\cdot\nabla\varphi=0 ∫ B ∇ v ⋅ ∇ φ = 0 を満たすなら、v v v は B B B 上の調和関数(に a.e. で等しい)である。
証明(補題 5.3) v ε = ρ ε ∗ v v_\varepsilon=\rho_\varepsilon*v v ε = ρ ε ∗ v は B ε = { x ∈ B ∣ dist ( x , ∂ B ) > ε } B_\varepsilon=\{x\in B\mid\operatorname{dist}(x,\partial B)\gt\varepsilon\} B ε = { x ∈ B ∣ dist ( x , ∂ B ) > ε } で C ∞ C^\infty C ∞ で、補題 2.1 の計算と同様に Δ v ε ( x ) = − ∫ ∇ y [ ρ ε ( x − y ) ] ⋅ ∇ v ( y ) d y = 0 \Delta v_\varepsilon(x)=-\int\nabla_y\bigl[\rho_\varepsilon(x-y)\bigr]\cdot\nabla v(y)\,dy=0 Δ v ε ( x ) = − ∫ ∇ y [ ρ ε ( x − y ) ] ⋅ ∇ v ( y ) d y = 0 です(ρ ε ( x − ⋅ ) ∈ C c ∞ ( B ) \rho_\varepsilon(x-\cdot)\in C_c^\infty(B) ρ ε ( x − ⋅ ) ∈ C c ∞ ( B ) )。K ⋐ B K\Subset B K ⋐ B と r < dist ( K , ∂ B ) / 2 r\lt\operatorname{dist}(K,\partial B)/2 r < dist ( K , ∂ B ) /2 をとると、ε , δ < r \varepsilon,\delta\lt r ε , δ < r で v ε − v δ v_\varepsilon-v_\delta v ε − v δ は K K K の r r r 近傍で調和なので、命題 3.4[ラプラス方程式と調和関数] の第 1 式と Cauchy–Schwarz の不等式から sup K ∣ v ε − v δ ∣ ≤ C ∥ v ε − v δ ∥ L 2 ( B ) → 0 \sup_K|v_\varepsilon-v_\delta|\le C\|v_\varepsilon-v_\delta\|_{L^2(B)}\to0 sup K ∣ v ε − v δ ∣ ≤ C ∥ v ε − v δ ∥ L 2 ( B ) → 0 です。よって v ε v_\varepsilon v ε は K K K 上一様に連続関数 v ˉ \bar v v ˉ に収束し、L 2 L^2 L 2 の極限 v v v と a.e. で一致します。v ˉ \bar v v ˉ は平均値の性質をもつ関数の一様極限なので平均値の性質をもち、定理 3.3[ラプラス方程式と調和関数] により調和です。
∎
補題 5.4 (調和関数の勾配の減衰 ) ✓
v v v が B ( x 0 , R ) B(x_0,R) B ( x 0 , R ) で調和で ∇ v ∈ L 2 ( B ( x 0 , R ) ) \nabla v\in L^2(B(x_0,R)) ∇ v ∈ L 2 ( B ( x 0 , R )) なら、0 < ρ ≤ R 0\lt\rho\le R 0 < ρ ≤ R について
∫ B ( x 0 , ρ ) ∣ ∇ v − ( ∇ v ) x 0 , ρ ∣ 2 ≤ C ( n ) ( ρ R ) n + 2 ∫ B ( x 0 , R ) ∣ ∇ v − ( ∇ v ) x 0 , R ∣ 2 . \int_{B(x_0,\rho)}|\nabla v-(\nabla v)_{x_0,\rho}|^2\le C(n)\Bigl(\frac\rho R\Bigr)^{n+2}\int_{B(x_0,R)}|\nabla v-(\nabla v)_{x_0,R}|^2 . ∫ B ( x 0 , ρ ) ∣∇ v − ( ∇ v ) x 0 , ρ ∣ 2 ≤ C ( n ) ( R ρ ) n + 2 ∫ B ( x 0 , R ) ∣∇ v − ( ∇ v ) x 0 , R ∣ 2 .
証明(補題 5.4) ρ ≥ R / 4 \rho\ge R/4 ρ ≥ R /4 なら、球平均の最小性と積分範囲の包含から C = 4 n + 2 C=4^{n+2} C = 4 n + 2 で成り立ちます。ρ < R / 4 \rho\lt R/4 ρ < R /4 とし、h = ∂ k v − c h=\partial_kv-c h = ∂ k v − c (c = ( ∂ k v ) x 0 , R c=(\partial_kv)_{x_0,R} c = ( ∂ k v ) x 0 , R )とおくと h h h は調和です。z ∈ B ( x 0 , R / 2 ) z\in B(x_0,R/2) z ∈ B ( x 0 , R /2 ) なら 命題 3.4[ラプラス方程式と調和関数] と Cauchy–Schwarz の不等式で ∣ h ( z ) ∣ ≤ ( ω n ( R / 2 ) n ) − 1 / 2 ∥ h ∥ L 2 ( B ( x 0 , R ) ) |h(z)|\le(\omega_n(R/2)^n)^{-1/2}\|h\|_{L^2(B(x_0,R))} ∣ h ( z ) ∣ ≤ ( ω n ( R /2 ) n ) − 1/2 ∥ h ∥ L 2 ( B ( x 0 , R )) 、x ∈ B ( x 0 , R / 4 ) x\in B(x_0,R/4) x ∈ B ( x 0 , R /4 ) なら同じ命題の第 2 式(半径 R / 4 R/4 R /4 )で ∣ ∇ h ( x ) ∣ ≤ 4 n 3 / 2 R sup B ( x 0 , R / 2 ) ∣ h ∣ |\nabla h(x)|\le\frac{4n^{3/2}}R\sup_{B(x_0,R/2)}|h| ∣∇ h ( x ) ∣ ≤ R 4 n 3/2 sup B ( x 0 , R /2 ) ∣ h ∣ です。合わせて sup B ( x 0 , R / 4 ) ∣ ∇ h ∣ 2 ≤ C R − n − 2 ∥ h ∥ L 2 ( B ( x 0 , R ) ) 2 \sup_{B(x_0,R/4)}|\nabla h|^2\le CR^{-n-2}\|h\|_{L^2(B(x_0,R))}^2 sup B ( x 0 , R /4 ) ∣∇ h ∣ 2 ≤ C R − n − 2 ∥ h ∥ L 2 ( B ( x 0 , R )) 2 です。x ∈ B ( x 0 , ρ ) x\in B(x_0,\rho) x ∈ B ( x 0 , ρ ) では ∣ h ( x ) − h ( x 0 ) ∣ ≤ ρ sup ∣ ∇ h ∣ |h(x)-h(x_0)|\le\rho\sup|\nabla h| ∣ h ( x ) − h ( x 0 ) ∣ ≤ ρ sup ∣∇ h ∣ なので
∫ B ( x 0 , ρ ) ∣ ∂ k v − ( ∂ k v ) x 0 , ρ ∣ 2 ≤ ∫ B ( x 0 , ρ ) ∣ h − h ( x 0 ) ∣ 2 ≤ ω n ρ n + 2 sup B ( x 0 , R / 4 ) ∣ ∇ h ∣ 2 ≤ C ( ρ R ) n + 2 ∫ B ( x 0 , R ) ∣ h ∣ 2 \int_{B(x_0,\rho)}|\partial_kv-(\partial_kv)_{x_0,\rho}|^2\le\int_{B(x_0,\rho)}|h-h(x_0)|^2\le\omega_n\rho^{n+2}\sup_{B(x_0,R/4)}|\nabla h|^2\le C\Bigl(\frac\rho R\Bigr)^{n+2}\int_{B(x_0,R)}|h|^2 ∫ B ( x 0 , ρ ) ∣ ∂ k v − ( ∂ k v ) x 0 , ρ ∣ 2 ≤ ∫ B ( x 0 , ρ ) ∣ h − h ( x 0 ) ∣ 2 ≤ ω n ρ n + 2 B ( x 0 , R /4 ) sup ∣∇ h ∣ 2 ≤ C ( R ρ ) n + 2 ∫ B ( x 0 , R ) ∣ h ∣ 2 です。k k k について足せば主張です。
∎
補題 5.5 (反復補題 ) ✓
ϕ \phi ϕ を ( 0 , R 0 ] (0,R_0] ( 0 , R 0 ] 上の非負の単調増加関数とし、定数 A ≥ 1 A\ge1 A ≥ 1 、B ≥ 0 B\ge0 B ≥ 0 、0 < b < a 0\lt b\lt a 0 < b < a について
ϕ ( ρ ) ≤ A ( ρ R ) a ϕ ( R ) + B R b ( 0 < ρ ≤ R ≤ R 0 ) \phi(\rho)\le A\Bigl(\frac\rho R\Bigr)^a\phi(R)+BR^b\qquad(0\lt\rho\le R\le R_0) ϕ ( ρ ) ≤ A ( R ρ ) a ϕ ( R ) + B R b ( 0 < ρ ≤ R ≤ R 0 ) が成り立つとする。このとき C = C ( A , a , b ) C=C(A,a,b) C = C ( A , a , b ) があって、ϕ ( ρ ) ≤ C [ ( ρ / R ) b ϕ ( R ) + B ρ b ] \phi(\rho)\le C\bigl[(\rho/R)^b\phi(R)+B\rho^b\bigr] ϕ ( ρ ) ≤ C [ ( ρ / R ) b ϕ ( R ) + B ρ b ] (0 < ρ ≤ R ≤ R 0 0\lt\rho\le R\le R_0 0 < ρ ≤ R ≤ R 0 )。
この補題は、Schauder 評価の証明の最後の段で使います。そこで得られる不等式は、半径 R R R から ρ \rho ρ へ縮めるたびに減衰が ( ρ / R ) a (\rho/R)^a ( ρ / R ) a だけ効き、誤差が R b R^b R b だけ加わる形です。補題はこれを繰り返して、小さな半径でも ρ b \rho^b ρ b の減衰(Campanato の条件)にまとめ直します。証明は、一定の比 τ \tau τ ずつ半径を縮めて不等式を使い続けるものです。
証明(補題 5.5) c = ( a + b ) / 2 c=(a+b)/2 c = ( a + b ) /2 とし、τ ∈ ( 0 , 1 ) \tau\in(0,1) τ ∈ ( 0 , 1 ) を A τ a ≤ τ c A\tau^a\le\tau^c A τ a ≤ τ c となるようにとります(τ a − c ≤ 1 / A \tau^{a-c}\le1/A τ a − c ≤ 1/ A )。仮定を ρ = τ R \rho=\tau R ρ = τ R に使うと ϕ ( τ R ) ≤ τ c ϕ ( R ) + B R b \phi(\tau R)\le\tau^c\phi(R)+BR^b ϕ ( τ R ) ≤ τ c ϕ ( R ) + B R b で、これを R , τ R , τ 2 R , … R,\tau R,\tau^2R,\ldots R , τ R , τ 2 R , … に繰り返すと
ϕ ( τ k R ) ≤ τ k c ϕ ( R ) + B R b ∑ j = 0 k − 1 τ c ( k − 1 − j ) τ j b ≤ τ k b ϕ ( R ) + B R b τ ( k − 1 ) b 1 − τ c − b \phi(\tau^kR)\le\tau^{kc}\phi(R)+BR^b\sum_{j=0}^{k-1}\tau^{c(k-1-j)}\tau^{jb}\le\tau^{kb}\phi(R)+\frac{BR^b\tau^{(k-1)b}}{1-\tau^{c-b}} ϕ ( τ k R ) ≤ τ k c ϕ ( R ) + B R b j = 0 ∑ k − 1 τ c ( k − 1 − j ) τ j b ≤ τ k b ϕ ( R ) + 1 − τ c − b B R b τ ( k − 1 ) b です。τ k + 1 R < ρ ≤ τ k R \tau^{k+1}R\lt\rho\le\tau^kR τ k + 1 R < ρ ≤ τ k R なら、単調性と τ k b ≤ τ − b ( ρ / R ) b \tau^{kb}\le\tau^{-b}(\rho/R)^b τ k b ≤ τ − b ( ρ / R ) b から ϕ ( ρ ) ≤ ϕ ( τ k R ) ≤ τ − b ( ρ / R ) b ϕ ( R ) + τ − 2 b 1 − τ c − b B ρ b \phi(\rho)\le\phi(\tau^kR)\le\tau^{-b}(\rho/R)^b\phi(R)+\frac{\tau^{-2b}}{1-\tau^{c-b}}B\rho^b ϕ ( ρ ) ≤ ϕ ( τ k R ) ≤ τ − b ( ρ / R ) b ϕ ( R ) + 1 − τ c − b τ − 2 b B ρ b です。
∎
定理 5.6 (Schauder 評価(Campanato の方法) ) ✓
0 < α < 1 0\lt\alpha\lt1 0 < α < 1 、B 1 = B ( 0 , 1 ) B_1=B(0,1) B 1 = B ( 0 , 1 ) とし、F ∈ C 0 , α ( B ˉ 1 ; R n ) F\in C^{0,\alpha}(\bar B_1;\mathbb{R}^n) F ∈ C 0 , α ( B ˉ 1 ; R n ) 、u ∈ H 1 ( B 1 ) u\in H^1(B_1) u ∈ H 1 ( B 1 ) がすべての φ ∈ C c ∞ ( B 1 ) \varphi\in C_c^\infty(B_1) φ ∈ C c ∞ ( B 1 ) で
∫ B 1 ∇ u ⋅ ∇ φ d x = ∫ B 1 F ⋅ ∇ φ d x \int_{B_1}\nabla u\cdot\nabla\varphi\,dx=\int_{B_1}F\cdot\nabla\varphi\,dx ∫ B 1 ∇ u ⋅ ∇ φ d x = ∫ B 1 F ⋅ ∇ φ d x を満たす(Δ u = div F \Delta u=\operatorname{div}F Δ u = div F の弱解)とする。このとき ∇ u \nabla u ∇ u は B 1 / 2 B_{1/2} B 1/2 上で連続な代表元をもち、
[ ∇ u ] α ; B 1 / 2 ≤ C ( n , α ) ( ∥ ∇ u ∥ L 2 ( B 1 ) + [ F ] α ; B 1 ) . [\nabla u]_{\alpha;B_{1/2}}\le C(n,\alpha)\bigl(\|\nabla u\|_{L^2(B_1)}+[F]_{\alpha;B_1}\bigr). [ ∇ u ] α ; B 1/2 ≤ C ( n , α ) ( ∥∇ u ∥ L 2 ( B 1 ) + [ F ] α ; B 1 ) .
x₀ B(x₀, ρ) B(x₀, R) u = v + w v : 調和(F を F(x₀) に凍結) → 勾配の振れが (ρ/R)^(n+2) で減衰 w : 誤差(F の Hölder 連続性が効く) → 大きさは R^(n+2α) 程度
凍結の考え方。小さな球 B(x₀, ρ) を含む球 B(x₀, R) の中で、u を調和関数 v と誤差 w の和に分ける。w は右辺 F − F(x₀) だけを受け持つ
証明の方針は、各点 x 0 x_0 x 0 の近くで右辺 F F F を定数に置き換えた方程式(解は調和関数)と比べることです。u u u の勾配の「平均からのずれ」が、小さい球ほどどう減るかを、調和関数の減衰と誤差の評価に分けて調べます。
証明(定理 5.6) x 0 ∈ B 1 / 2 x_0\in B_{1/2} x 0 ∈ B 1/2 、0 < R ≤ 1 / 2 0\lt R\le1/2 0 < R ≤ 1/2 を固定し、B R = B ( x 0 , R ) B_R=B(x_0,R) B R = B ( x 0 , R ) と書きます。
分解。 H 0 1 ( B R ) H^1_0(B_R) H 0 1 ( B R ) 上の汎関数 φ ↦ ∫ B R ( F − F ( x 0 ) ) ⋅ ∇ φ \varphi\mapsto\int_{B_R}(F-F(x_0))\cdot\nabla\varphi φ ↦ ∫ B R ( F − F ( x 0 )) ⋅ ∇ φ は有界なので、Lax–Milgram の定理(定理 3.2[楕円型方程式の弱解と変分法] 、補題 2.1[楕円型方程式の弱解と変分法] により ⟨ ∇ ⋅ , ∇ ⋅ ⟩ \langle\nabla\cdot,\nabla\cdot\rangle ⟨ ∇ ⋅ , ∇ ⋅ ⟩ は H 0 1 H^1_0 H 0 1 の同値な内積)で w ∈ H 0 1 ( B R ) w\in H^1_0(B_R) w ∈ H 0 1 ( B R ) を
∫ B R ∇ w ⋅ ∇ φ = ∫ B R ( F − F ( x 0 ) ) ⋅ ∇ φ ( φ ∈ H 0 1 ( B R ) ) \int_{B_R}\nabla w\cdot\nabla\varphi=\int_{B_R}(F-F(x_0))\cdot\nabla\varphi\qquad(\varphi\in H^1_0(B_R)) ∫ B R ∇ w ⋅ ∇ φ = ∫ B R ( F − F ( x 0 )) ⋅ ∇ φ ( φ ∈ H 0 1 ( B R )) となるようにとれます。φ = w \varphi=w φ = w として Cauchy–Schwarz の不等式を使うと
∫ B R ∣ ∇ w ∣ 2 ≤ ∫ B R ∣ F − F ( x 0 ) ∣ 2 ≤ ω n [ F ] α 2 R n + 2 α \int_{B_R}|\nabla w|^2\le\int_{B_R}|F-F(x_0)|^2\le\omega_n[F]_\alpha^2R^{n+2\alpha} ∫ B R ∣∇ w ∣ 2 ≤ ∫ B R ∣ F − F ( x 0 ) ∣ 2 ≤ ω n [ F ] α 2 R n + 2 α です。v = u − w v=u-w v = u − w は、φ ∈ C c ∞ ( B R ) \varphi\in C_c^\infty(B_R) φ ∈ C c ∞ ( B R ) について ∫ ∇ v ⋅ ∇ φ = ∫ F ⋅ ∇ φ − ∫ ( F − F ( x 0 ) ) ⋅ ∇ φ = F ( x 0 ) ⋅ ∫ ∇ φ = 0 \int\nabla v\cdot\nabla\varphi=\int F\cdot\nabla\varphi-\int(F-F(x_0))\cdot\nabla\varphi=F(x_0)\cdot\int\nabla\varphi=0 ∫ ∇ v ⋅ ∇ φ = ∫ F ⋅ ∇ φ − ∫ ( F − F ( x 0 )) ⋅ ∇ φ = F ( x 0 ) ⋅ ∫ ∇ φ = 0 を満たすので、補題 5.3 により B R B_R B R で調和です。
減衰の不等式。 ϕ ( ρ ) = ∫ B ( x 0 , ρ ) ∣ ∇ u − ( ∇ u ) x 0 , ρ ∣ 2 \phi(\rho)=\int_{B(x_0,\rho)}|\nabla u-(\nabla u)_{x_0,\rho}|^2 ϕ ( ρ ) = ∫ B ( x 0 , ρ ) ∣∇ u − ( ∇ u ) x 0 , ρ ∣ 2 とおきます。∇ u = ∇ v + ∇ w \nabla u=\nabla v+\nabla w ∇ u = ∇ v + ∇ w と球平均の最小性、補題 5.4 から、ρ ≤ R \rho\le R ρ ≤ R で
ϕ ( ρ ) ≤ 2 ∫ B ρ ∣ ∇ v − ( ∇ v ) ρ ∣ 2 + 2 ∫ B ρ ∣ ∇ w ∣ 2 ≤ 2 C ( ρ R ) n + 2 ∫ B R ∣ ∇ v − ( ∇ v ) R ∣ 2 + 2 ∫ B R ∣ ∇ w ∣ 2 \phi(\rho)\le2\int_{B_\rho}|\nabla v-(\nabla v)_\rho|^2+2\int_{B_\rho}|\nabla w|^2\le2C\Bigl(\frac\rho R\Bigr)^{n+2}\int_{B_R}|\nabla v-(\nabla v)_R|^2+2\int_{B_R}|\nabla w|^2 ϕ ( ρ ) ≤ 2 ∫ B ρ ∣∇ v − ( ∇ v ) ρ ∣ 2 + 2 ∫ B ρ ∣∇ w ∣ 2 ≤ 2 C ( R ρ ) n + 2 ∫ B R ∣∇ v − ( ∇ v ) R ∣ 2 + 2 ∫ B R ∣∇ w ∣ 2 で、∫ B R ∣ ∇ v − ( ∇ v ) R ∣ 2 ≤ ∫ B R ∣ ∇ v − ( ∇ u ) R ∣ 2 ≤ 2 ϕ ( R ) + 2 ∫ B R ∣ ∇ w ∣ 2 \int_{B_R}|\nabla v-(\nabla v)_R|^2\le\int_{B_R}|\nabla v-(\nabla u)_R|^2\le2\phi(R)+2\int_{B_R}|\nabla w|^2 ∫ B R ∣∇ v − ( ∇ v ) R ∣ 2 ≤ ∫ B R ∣∇ v − ( ∇ u ) R ∣ 2 ≤ 2 ϕ ( R ) + 2 ∫ B R ∣∇ w ∣ 2 を代入すると
ϕ ( ρ ) ≤ 4 C ( ρ R ) n + 2 ϕ ( R ) + C ′ [ F ] α 2 R n + 2 α \phi(\rho)\le4C\Bigl(\frac\rho R\Bigr)^{n+2}\phi(R)+C'[F]_\alpha^2R^{n+2\alpha} ϕ ( ρ ) ≤ 4 C ( R ρ ) n + 2 ϕ ( R ) + C ′ [ F ] α 2 R n + 2 α です。
反復と結論。 補題 5.5 (a = n + 2 a=n+2 a = n + 2 、b = n + 2 α b=n+2\alpha b = n + 2 α 、R 0 = 1 / 2 R_0=1/2 R 0 = 1/2 )から ϕ ( ρ ) ≤ C [ ( ρ / R ) n + 2 α ϕ ( R ) + [ F ] α 2 ρ n + 2 α ] \phi(\rho)\le C\bigl[(\rho/R)^{n+2\alpha}\phi(R)+[F]_\alpha^2\rho^{n+2\alpha}\bigr] ϕ ( ρ ) ≤ C [ ( ρ / R ) n + 2 α ϕ ( R ) + [ F ] α 2 ρ n + 2 α ] で、R = 1 / 2 R=1/2 R = 1/2 とおき ϕ ( 1 / 2 ) ≤ ∫ B 1 ∣ ∇ u ∣ 2 \phi(1/2)\le\int_{B_1}|\nabla u|^2 ϕ ( 1/2 ) ≤ ∫ B 1 ∣∇ u ∣ 2 を使うと、全ての x 0 ∈ B 1 / 2 x_0\in B_{1/2} x 0 ∈ B 1/2 、ρ ≤ 1 / 2 \rho\le1/2 ρ ≤ 1/2 で
∫ B ( x 0 , ρ ) ∣ ∇ u − ( ∇ u ) x 0 , ρ ∣ 2 ≤ C ( ∥ ∇ u ∥ L 2 ( B 1 ) 2 + [ F ] α 2 ) ρ n + 2 α \int_{B(x_0,\rho)}|\nabla u-(\nabla u)_{x_0,\rho}|^2\le C\bigl(\|\nabla u\|_{L^2(B_1)}^2+[F]_\alpha^2\bigr)\rho^{n+2\alpha} ∫ B ( x 0 , ρ ) ∣∇ u − ( ∇ u ) x 0 , ρ ∣ 2 ≤ C ( ∥∇ u ∥ L 2 ( B 1 ) 2 + [ F ] α 2 ) ρ n + 2 α です。定理 5.2 (中心 0 0 0 、R = 1 / 2 R=1/2 R = 1/2 、B ( 0 , 1 ) ⊂ B 1 B(0,1)\subset B_1 B ( 0 , 1 ) ⊂ B 1 )を ∇ u \nabla u ∇ u の各成分に使えば主張です。
∎
証明の構造は「凍結」です。各点 x 0 x_0 x 0 の近くで、右辺を定数 F ( x 0 ) F(x_0) F ( x 0 ) に凍結した方程式の解(調和関数 v v v )は ρ n + 2 \rho^{n+2} ρ n + 2 という最良の減衰をもち、凍結の誤差 w w w は F F F の Hölder 連続性のおかげで R n + 2 α R^{n+2\alpha} R n + 2 α しか寄与しません。α < 1 \alpha\lt1 α < 1 が必要なのは反復補題の b < a b\lt a b < a のためで、α = 1 \alpha=1 α = 1 (Lipschitz)では結論が成り立たないことが知られています。
系 5.7 (Poisson 方程式の Schauder 評価 ) ✓
0 < α < 1 0\lt\alpha\lt1 0 < α < 1 、f ∈ C 0 , α ( B ˉ 1 ) f\in C^{0,\alpha}(\bar B_1) f ∈ C 0 , α ( B ˉ 1 ) とし、u ∈ H 1 ( B 1 ) u\in H^1(B_1) u ∈ H 1 ( B 1 ) が − Δ u = f -\Delta u=f − Δ u = f の弱解(∫ ∇ u ⋅ ∇ φ = ∫ f φ \int\nabla u\cdot\nabla\varphi=\int f\varphi ∫ ∇ u ⋅ ∇ φ = ∫ f φ 、φ ∈ C c ∞ ( B 1 ) \varphi\in C_c^\infty(B_1) φ ∈ C c ∞ ( B 1 ) )であるとする。このとき 2 階の弱微分 D 2 u D^2u D 2 u は B 1 / 4 B_{1/4} B 1/4 上で α \alpha α -Hölder 連続な代表元をもち、[ D 2 u ] α ; B 1 / 4 ≤ C ( n , α ) ( ∥ f ∥ C 0 , α ( B ˉ 1 ) + ∥ u ∥ H 1 ( B 1 ) ) [D^2u]_{\alpha;B_{1/4}}\le C(n,\alpha)\bigl(\|f\|_{C^{0,\alpha}(\bar B_1)}+\|u\|_{H^1(B_1)}\bigr) [ D 2 u ] α ; B 1/4 ≤ C ( n , α ) ( ∥ f ∥ C 0 , α ( B ˉ 1 ) + ∥ u ∥ H 1 ( B 1 ) ) 。ここで ∥ f ∥ C 0 , α = sup ∣ f ∣ + [ f ] α \|f\|_{C^{0,\alpha}}=\sup|f|+[f]_\alpha ∥ f ∥ C 0 , α = sup ∣ f ∣ + [ f ] α である。
証明(系 5.7) 定理 3.1 (a i j = δ i j a^{ij}=\delta_{ij} a ij = δ ij 、b = c = 0 b=c=0 b = c = 0 、Ω = B 1 \Omega=B_1 Ω = B 1 、V = B 3 / 4 V=B_{3/4} V = B 3/4 )から u ∈ H 2 ( B 3 / 4 ) u\in H^2(B_{3/4}) u ∈ H 2 ( B 3/4 ) 、∥ D 2 u ∥ L 2 ( B 3 / 4 ) ≤ C ( ∥ f ∥ L 2 + ∥ u ∥ H 1 ) \|D^2u\|_{L^2(B_{3/4})}\le C(\|f\|_{L^2}+\|u\|_{H^1}) ∥ D 2 u ∥ L 2 ( B 3/4 ) ≤ C ( ∥ f ∥ L 2 + ∥ u ∥ H 1 ) です。U k = ∂ k u ∈ H 1 ( B 3 / 4 ) U_k=\partial_ku\in H^1(B_{3/4}) U k = ∂ k u ∈ H 1 ( B 3/4 ) について、φ ∈ C c ∞ ( B 3 / 4 ) \varphi\in C_c^\infty(B_{3/4}) φ ∈ C c ∞ ( B 3/4 ) とすると、弱微分の定義と弱形式(試験関数 ∂ k φ \partial_k\varphi ∂ k φ )から
∫ ∇ U k ⋅ ∇ φ = − ∫ ∇ u ⋅ ∇ ∂ k φ = − ∫ f ∂ k φ = ∫ ( − f e k ) ⋅ ∇ φ \int\nabla U_k\cdot\nabla\varphi=-\int\nabla u\cdot\nabla\partial_k\varphi=-\int f\,\partial_k\varphi=\int(-f\,e_k)\cdot\nabla\varphi ∫ ∇ U k ⋅ ∇ φ = − ∫ ∇ u ⋅ ∇ ∂ k φ = − ∫ f ∂ k φ = ∫ ( − f e k ) ⋅ ∇ φ です。U k U_k U k は F = − f e k ∈ C 0 , α F=-fe_k\in C^{0,\alpha} F = − f e k ∈ C 0 , α について 定理 5.6 の方程式を B 3 / 4 B_{3/4} B 3/4 で満たすので、相似変換 x ↦ 3 4 x x\mapsto\frac34x x ↦ 4 3 x で B 1 B_1 B 1 に移して定理を使えば、∇ U k \nabla U_k ∇ U k は B 3 / 8 ⊃ B 1 / 4 B_{3/8}\supset B_{1/4} B 3/8 ⊃ B 1/4 上で Hölder 連続で、評価が出ます。
∎
2 階の弱微分が連続なら u u u は古典的な意味で C 2 C^2 C 2 で、− Δ u = f -\Delta u=f − Δ u = f が各点で成り立ちます。第 1 章で f ∈ C c 2 f\in C_c^2 f ∈ C c 2 と仮定して Φ ∗ f ∈ C 2 \Phi*f\in C^2 Φ ∗ f ∈ C 2 を示しましたが、この系により f f f が Hölder 連続なら十分です。変数係数 a i j ∈ C 0 , α a^{ij}\in C^{0,\alpha} a ij ∈ C 0 , α の場合も、係数を x 0 x_0 x 0 で凍結して定数係数の方程式(線形変換で Laplace 方程式に帰着)と比べる同じ方法で ∇ u ∈ C 0 , α \nabla u\in C^{0,\alpha} ∇ u ∈ C 0 , α が示されます。ただし凍結の誤差 ∫ B R ∣ ( a ( x ) − a ( x 0 ) ) ∇ u ∣ 2 ≤ C R 2 α ∫ B R ∣ ∇ u ∣ 2 \int_{B_R}|(a(x)-a(x_0))\nabla u|^2\le CR^{2\alpha}\int_{B_R}|\nabla u|^2 ∫ B R ∣ ( a ( x ) − a ( x 0 )) ∇ u ∣ 2 ≤ C R 2 α ∫ B R ∣∇ u ∣ 2 を扱うために、先に ∫ B ρ ∣ ∇ u ∣ 2 ≤ C ρ n − 2 + 2 β \int_{B_\rho}|\nabla u|^2\le C\rho^{n-2+2\beta} ∫ B ρ ∣∇ u ∣ 2 ≤ C ρ n − 2 + 2 β の形の評価を同じ反復で得る必要があります(Giaquinta–Martinazzi の第 5 章)。
補題 6.1 (正の部分の弱微分 ) ✓
u ∈ H 1 ( Ω ) u\in H^1(\Omega) u ∈ H 1 ( Ω ) なら u + = max ( u , 0 ) ∈ H 1 ( Ω ) u^+=\max(u,0)\in H^1(\Omega) u + = max ( u , 0 ) ∈ H 1 ( Ω ) で、∇ u + = 1 { u > 0 } ∇ u \nabla u^+=\mathbf 1_{\{u\gt0\}}\nabla u ∇ u + = 1 { u > 0 } ∇ u a.e. である。さらに u ∈ H 0 1 ( Ω ) u\in H^1_0(\Omega) u ∈ H 0 1 ( Ω ) なら u + ∈ H 0 1 ( Ω ) u^+\in H^1_0(\Omega) u + ∈ H 0 1 ( Ω ) である。
証明(補題 6.1) 連鎖律。 G ∈ C 1 ( R ) G\in C^1(\mathbb{R}) G ∈ C 1 ( R ) 、G ′ G' G ′ は有界、G ( 0 ) = 0 G(0)=0 G ( 0 ) = 0 とすると、G ( u ) ∈ H 1 ( Ω ) G(u)\in H^1(\Omega) G ( u ) ∈ H 1 ( Ω ) で ∇ G ( u ) = G ′ ( u ) ∇ u \nabla G(u)=G'(u)\nabla u ∇ G ( u ) = G ′ ( u ) ∇ u です。実際、W ⋐ Ω W\Subset\Omega W ⋐ Ω 上で 補題 2.1 の u k = ρ 1 / k ∗ u u_k=\rho_{1/k}*u u k = ρ 1/ k ∗ u をとり、L 2 ( W ) L^2(W) L 2 ( W ) で u k → u u_k\to u u k → u 、∇ u k → ∇ u \nabla u_k\to\nabla u ∇ u k → ∇ u 、さらに部分列で a.e. 収束とします(∥ u k j − u ∥ L 2 ( W ) < 2 − j \|u_{k_j}-u\|_{L^2(W)}\lt2^{-j} ∥ u k j − u ∥ L 2 ( W ) < 2 − j となる部分列を取ると、単調収束定理とミンコフスキーの不等式から ∥ ∑ j ∣ u k j − u ∣ ∥ L 2 ( W ) ≤ ∑ j 2 − j < ∞ \|\sum_j|u_{k_j}-u|\|_{L^2(W)}\le\sum_j2^{-j}\lt\infty ∥ ∑ j ∣ u k j − u ∣ ∥ L 2 ( W ) ≤ ∑ j 2 − j < ∞ なので ∑ j ∣ u k j − u ∣ \sum_j|u_{k_j}-u| ∑ j ∣ u k j − u ∣ はほとんど至るところ有限で、したがって u k j → u u_{k_j}\to u u k j → u a.e. です)。φ ∈ C c ∞ ( W ) \varphi\in C_c^\infty(W) φ ∈ C c ∞ ( W ) について古典的に ∫ G ( u k ) ∂ i φ = − ∫ G ′ ( u k ) ∂ i u k φ \int G(u_k)\partial_i\varphi=-\int G'(u_k)\partial_iu_k\,\varphi ∫ G ( u k ) ∂ i φ = − ∫ G ′ ( u k ) ∂ i u k φ で、∣ G ( u k ) − G ( u ) ∣ ≤ sup ∣ G ′ ∣ ∣ u k − u ∣ |G(u_k)-G(u)|\le\sup|G'|\,|u_k-u| ∣ G ( u k ) − G ( u ) ∣ ≤ sup ∣ G ′ ∣ ∣ u k − u ∣ 、∣ G ′ ( u k ) ∂ i u k − G ′ ( u ) ∂ i u ∣ ≤ sup ∣ G ′ ∣ ∣ ∂ i u k − ∂ i u ∣ + ∣ G ′ ( u k ) − G ′ ( u ) ∣ ∣ ∂ i u ∣ |G'(u_k)\partial_iu_k-G'(u)\partial_iu|\le\sup|G'|\,|\partial_iu_k-\partial_iu|+|G'(u_k)-G'(u)||\partial_iu| ∣ G ′ ( u k ) ∂ i u k − G ′ ( u ) ∂ i u ∣ ≤ sup ∣ G ′ ∣ ∣ ∂ i u k − ∂ i u ∣ + ∣ G ′ ( u k ) − G ′ ( u ) ∣∣ ∂ i u ∣ で、最後の項は優収束定理で L 1 ( W ) L^1(W) L 1 ( W ) で 0 0 0 に収束します。極限で ∫ G ( u ) ∂ i φ = − ∫ G ′ ( u ) ∂ i u φ \int G(u)\partial_i\varphi=-\int G'(u)\partial_iu\,\varphi ∫ G ( u ) ∂ i φ = − ∫ G ′ ( u ) ∂ i u φ です。
正の部分。 G ε ( t ) = t 2 + ε 2 − ε G_\varepsilon(t)=\sqrt{t^2+\varepsilon^2}-\varepsilon G ε ( t ) = t 2 + ε 2 − ε (t > 0 t\gt0 t > 0 )、G ε ( t ) = 0 G_\varepsilon(t)=0 G ε ( t ) = 0 (t ≤ 0 t\le0 t ≤ 0 )は連鎖律の条件を満たし、ε → 0 \varepsilon\to0 ε → 0 で G ε ( u ) → u + G_\varepsilon(u)\to u^+ G ε ( u ) → u + 、G ε ′ ( u ) ∇ u → 1 { u > 0 } ∇ u G_\varepsilon'(u)\nabla u\to\mathbf 1_{\{u\gt0\}}\nabla u G ε ′ ( u ) ∇ u → 1 { u > 0 } ∇ u が各点で、優関数 ∣ u ∣ |u| ∣ u ∣ 、∣ ∇ u ∣ |\nabla u| ∣∇ u ∣ つきで成り立ちます。弱微分の等式で極限をとれば前半です。
H 0 1 H^1_0 H 0 1 の場合。 φ k ∈ C c ∞ ( Ω ) \varphi_k\in C_c^\infty(\Omega) φ k ∈ C c ∞ ( Ω ) を H 1 H^1 H 1 で u u u に収束するようにとると、φ k + \varphi_k^+ φ k + は台がコンパクトな H 1 H^1 H 1 関数なので 補題 2.1 から H 0 1 H^1_0 H 0 1 に属し、前半から ∥ φ k + ∥ H 1 ≤ ∥ φ k ∥ H 1 \|\varphi_k^+\|_{H^1}\le\|\varphi_k\|_{H^1} ∥ φ k + ∥ H 1 ≤ ∥ φ k ∥ H 1 は有界です。H 0 1 H^1_0 H 0 1 で弱収束する部分列の極限を w w w とすると、∣ φ k + − u + ∣ ≤ ∣ φ k − u ∣ |\varphi_k^+-u^+|\le|\varphi_k-u| ∣ φ k + − u + ∣ ≤ ∣ φ k − u ∣ から φ k + → u + \varphi_k^+\to u^+ φ k + → u + (L 2 L^2 L 2 )なので、w = u + w=u^+ w = u + です。よって u + ∈ H 0 1 ( Ω ) u^+\in H^1_0(\Omega) u + ∈ H 0 1 ( Ω ) です。
∎
定理 6.2 (弱最大値原理 ) ✓
Ω \Omega Ω を有界な開集合、L = − ∑ ∂ j ( a i j ∂ i ⋅ ) + c L=-\sum\partial_j(a^{ij}\partial_i\,\cdot)+c L = − ∑ ∂ j ( a ij ∂ i ⋅ ) + c を一様楕円型(a i j , c ∈ L ∞ a^{ij},c\in L^\infty a ij , c ∈ L ∞ )で c ≥ 0 c\ge0 c ≥ 0 とする。u ∈ H 1 ( Ω ) u\in H^1(\Omega) u ∈ H 1 ( Ω ) が、v ≥ 0 v\ge0 v ≥ 0 を満たすすべての v ∈ H 0 1 ( Ω ) v\in H^1_0(\Omega) v ∈ H 0 1 ( Ω ) について B [ u , v ] ≤ 0 B[u,v]\le0 B [ u , v ] ≤ 0 (弱い意味で L u ≤ 0 Lu\le0 Lu ≤ 0 )を満たし、u + ∈ H 0 1 ( Ω ) u^+\in H^1_0(\Omega) u + ∈ H 0 1 ( Ω ) (弱い意味で境界上 u ≤ 0 u\le0 u ≤ 0 )であるとする。このとき Ω \Omega Ω 上 a.e. で u ≤ 0 u\le0 u ≤ 0 である。
証明(定理 6.2) v = u + ≥ 0 v=u^+\ge0 v = u + ≥ 0 を代入します。補題 6.1 から ∂ i u ∂ j u + = ∂ i u + ∂ j u + \partial_iu\,\partial_ju^+=\partial_iu^+\,\partial_ju^+ ∂ i u ∂ j u + = ∂ i u + ∂ j u + 、u u + = ( u + ) 2 uu^+=(u^+)^2 u u + = ( u + ) 2 なので
0 ≥ B [ u , u + ] = ∫ ( ∑ a i j ∂ i u + ∂ j u + + c ( u + ) 2 ) ≥ θ ∥ ∇ u + ∥ 2 0\ge B[u,u^+]=\int\Bigl(\sum a^{ij}\partial_iu^+\partial_ju^++c(u^+)^2\Bigr)\ge\theta\|\nabla u^+\|^2 0 ≥ B [ u , u + ] = ∫ ( ∑ a ij ∂ i u + ∂ j u + + c ( u + ) 2 ) ≥ θ ∥∇ u + ∥ 2 です。よって ∇ u + = 0 \nabla u^+=0 ∇ u + = 0 で、u + ∈ H 0 1 u^+\in H^1_0 u + ∈ H 0 1 なので Poincaré の不等式(補題 2.1[楕円型方程式の弱解と変分法] )から u + = 0 u^+=0 u + = 0 です。
∎
M ≥ 0 M\ge0 M ≥ 0 について ( u − M ) + ∈ H 0 1 (u-M)^+\in H^1_0 ( u − M ) + ∈ H 0 1 なら、c ≥ 0 c\ge0 c ≥ 0 から L ( u − M ) = L u − c M ≤ 0 L(u-M)=Lu-cM\le0 L ( u − M ) = Lu − c M ≤ 0 なので u ≤ M u\le M u ≤ M です。これが古典的な sup Ω u ≤ sup ∂ Ω u + \sup_\Omega u\le\sup_{\partial\Omega}u^+ sup Ω u ≤ sup ∂ Ω u + の弱い版です。− Δ u = f ≥ 0 -\Delta u=f\ge0 − Δ u = f ≥ 0 の弱解 u ∈ H 0 1 u\in H^1_0 u ∈ H 0 1 は、− u -u − u に定理を使って u ≥ 0 u\ge0 u ≥ 0 です(正値性の保存)。1 階の項 b i ∂ i u b^i\partial_iu b i ∂ i u がある場合も同じ結論が成り立ちますが、試験関数を工夫する必要があります(Gilbarg–Trudinger, Theorem 8.1)。
Schrödinger 作用素 − Δ + V -\Delta+V − Δ + V が下に有界かどうかを、正の解の存在で判定するのがこの節の定理です。Ω ⊂ R n \Omega\subset\mathbb{R}^n Ω ⊂ R n を開集合、V V V を Ω \Omega Ω 上の実数値連続関数とし、φ ∈ C c ∞ ( Ω ) \varphi\in C_c^\infty(\Omega) φ ∈ C c ∞ ( Ω ) について h [ φ ] = ∫ Ω ( ∣ ∇ φ ∣ 2 + V φ 2 ) d x h[\varphi]=\int_\Omega\bigl(|\nabla\varphi|^2+V\varphi^2\bigr)dx h [ φ ] = ∫ Ω ( ∣∇ φ ∣ 2 + V φ 2 ) d x とおきます。
定理 7.1 (Allegretto–Piepenbrink の定理 ) ✓
E ∈ R E\in\mathbb{R} E ∈ R とし、ψ ∈ C 2 ( Ω ) \psi\in C^2(\Omega) ψ ∈ C 2 ( Ω ) が Ω \Omega Ω 上で ψ > 0 \psi\gt0 ψ > 0 、− Δ ψ + V ψ ≥ E ψ -\Delta\psi+V\psi\ge E\psi − Δ ψ + V ψ ≥ E ψ を満たすとする。このとき、すべての φ ∈ C c ∞ ( Ω ) \varphi\in C_c^\infty(\Omega) φ ∈ C c ∞ ( Ω ) で
h [ φ ] = ∫ Ω ψ 2 ∣ ∇ φ ψ ∣ 2 d x + ∫ Ω φ 2 ψ ( − Δ ψ + V ψ ) d x ≥ E ∫ Ω φ 2 d x h[\varphi]=\int_\Omega\psi^2\Bigl|\nabla\frac\varphi\psi\Bigr|^2dx+\int_\Omega\frac{\varphi^2}\psi\bigl(-\Delta\psi+V\psi\bigr)dx\ge E\int_\Omega\varphi^2dx h [ φ ] = ∫ Ω ψ 2 ∇ ψ φ 2 d x + ∫ Ω ψ φ 2 ( − Δ ψ + V ψ ) d x ≥ E ∫ Ω φ 2 d x が成り立つ。
証明(定理 7.1) g = φ / ψ ∈ C c 2 ( Ω ) g=\varphi/\psi\in C_c^2(\Omega) g = φ / ψ ∈ C c 2 ( Ω ) とおくと ∇ φ = ψ ∇ g + g ∇ ψ \nabla\varphi=\psi\nabla g+g\nabla\psi ∇ φ = ψ ∇ g + g ∇ ψ なので
∣ ∇ φ ∣ 2 = ψ 2 ∣ ∇ g ∣ 2 + ψ ∇ ψ ⋅ ∇ ( g 2 ) + g 2 ∣ ∇ ψ ∣ 2 |\nabla\varphi|^2=\psi^2|\nabla g|^2+\psi\nabla\psi\cdot\nabla(g^2)+g^2|\nabla\psi|^2 ∣∇ φ ∣ 2 = ψ 2 ∣∇ g ∣ 2 + ψ ∇ ψ ⋅ ∇ ( g 2 ) + g 2 ∣∇ ψ ∣ 2 です。g 2 g^2 g 2 の台はコンパクトなので、中央の項を部分積分すると ∫ ψ ∇ ψ ⋅ ∇ ( g 2 ) = − ∫ g 2 ( ∣ ∇ ψ ∣ 2 + ψ Δ ψ ) \int\psi\nabla\psi\cdot\nabla(g^2)=-\int g^2\bigl(|\nabla\psi|^2+\psi\Delta\psi\bigr) ∫ ψ ∇ ψ ⋅ ∇ ( g 2 ) = − ∫ g 2 ( ∣∇ ψ ∣ 2 + ψ Δ ψ ) で、∣ ∇ ψ ∣ 2 |\nabla\psi|^2 ∣∇ ψ ∣ 2 の項が打ち消し合って ∫ ∣ ∇ φ ∣ 2 = ∫ ψ 2 ∣ ∇ g ∣ 2 − ∫ g 2 ψ Δ ψ \int|\nabla\varphi|^2=\int\psi^2|\nabla g|^2-\int g^2\psi\Delta\psi ∫ ∣∇ φ ∣ 2 = ∫ ψ 2 ∣∇ g ∣ 2 − ∫ g 2 ψ Δ ψ です。∫ V φ 2 = ∫ g 2 ψ ⋅ V ψ \int V\varphi^2=\int g^2\psi\cdot V\psi ∫ V φ 2 = ∫ g 2 ψ ⋅ V ψ を足せば等式で、第 1 項は非負、第 2 項は仮定から E ∫ g 2 ψ 2 = E ∫ φ 2 E\int g^2\psi^2=E\int\varphi^2 E ∫ g 2 ψ 2 = E ∫ φ 2 以上です。
∎
等式の部分は基底状態変換 と呼ばれ、ψ \psi ψ が固有関数(− Δ ψ + V ψ = E ψ -\Delta\psi+V\psi=E\psi − Δ ψ + V ψ = E ψ )なら h [ φ ] − E ∥ φ ∥ 2 = ∫ ψ 2 ∣ ∇ ( φ / ψ ) ∣ 2 h[\varphi]-E\|\varphi\|^2=\int\psi^2|\nabla(\varphi/\psi)|^2 h [ φ ] − E ∥ φ ∥ 2 = ∫ ψ 2 ∣∇ ( φ / ψ ) ∣ 2 です。したがって、ψ \psi ψ 自身が二次形式の(閉包の)定義域に入る正の固有関数なら、その固有値 E E E は h [ ψ ] / ∥ ψ ∥ 2 = E h[\psi]/\|\psi\|^2=E h [ ψ ] /∥ ψ ∥ 2 = E から二次形式 h [ φ ] / ∥ φ ∥ 2 h[\varphi]/\|\varphi\|^2 h [ φ ] /∥ φ ∥ 2 の下限以下であり、定理と合わせて下限そのものです。下に有界な二次形式からは Friedrichs 拡大(定理 6.1[自己共役性の判定] 。下に有界な二次形式を、下限を保ったまま自己共役作用素 に結びつける標準的な方法)で − Δ + V ≥ E -\Delta+V\ge E − Δ + V ≥ E の自己共役作用素が定まります。
例 7.2 (Hardy の不等式 )
n ≥ 3 n\ge3 n ≥ 3 、Ω = R n ∖ { 0 } \Omega=\mathbb{R}^n\setminus\{0\} Ω = R n ∖ { 0 } 、ψ ( x ) = ∣ x ∣ − ( n − 2 ) / 2 \psi(x)=|x|^{-(n-2)/2} ψ ( x ) = ∣ x ∣ − ( n − 2 ) /2 とします。Δ ∣ x ∣ a = a ( a + n − 2 ) ∣ x ∣ a − 2 \Delta|x|^a=a(a+n-2)|x|^{a-2} Δ∣ x ∣ a = a ( a + n − 2 ) ∣ x ∣ a − 2 に a = − n − 2 2 a=-\frac{n-2}2 a = − 2 n − 2 を入れると − Δ ψ = ( n − 2 ) 2 4 ∣ x ∣ − 2 ψ -\Delta\psi=\frac{(n-2)^2}4|x|^{-2}\psi − Δ ψ = 4 ( n − 2 ) 2 ∣ x ∣ − 2 ψ です。V = − ( n − 2 ) 2 4 ∣ x ∣ 2 V=-\frac{(n-2)^2}{4|x|^2} V = − 4∣ x ∣ 2 ( n − 2 ) 2 、E = 0 E=0 E = 0 として 定理 7.1 を使うと
∫ R n ∣ ∇ φ ∣ 2 d x ≥ ( n − 2 ) 2 4 ∫ R n φ 2 ∣ x ∣ 2 d x ( φ ∈ C c ∞ ( R n ∖ { 0 } ) ) \int_{\mathbb{R}^n}|\nabla\varphi|^2dx\ge\frac{(n-2)^2}4\int_{\mathbb{R}^n}\frac{\varphi^2}{|x|^2}dx\qquad(\varphi\in C_c^\infty(\mathbb{R}^n\setminus\{0\})) ∫ R n ∣∇ φ ∣ 2 d x ≥ 4 ( n − 2 ) 2 ∫ R n ∣ x ∣ 2 φ 2 d x ( φ ∈ C c ∞ ( R n ∖ { 0 })) です。定数 ( n − 2 ) 2 4 \frac{(n-2)^2}4 4 ( n − 2 ) 2 は最良であることが知られています(ψ \psi ψ を切り落として近似すると比が定数に近づく)。ψ \psi ψ は原点の近くでちょうど ∇ ψ ∉ L 2 \nabla\psi\notin L^2 ∇ ψ ∈ / L 2 となる境目の関数です。
演習 8.1 易
φ ∈ C c ∞ ( ( 0 , π ) ) \varphi\in C_c^\infty((0,\pi)) φ ∈ C c ∞ (( 0 , π )) なら ∫ 0 π φ ′ 2 d x ≥ ∫ 0 π φ 2 d x \int_0^\pi\varphi'^2dx\ge\int_0^\pi\varphi^2dx ∫ 0 π φ ′2 d x ≥ ∫ 0 π φ 2 d x であることを、定理 7.1 を使って示せ。
解答 n = 1 n=1 n = 1 、Ω = ( 0 , π ) \Omega=(0,\pi) Ω = ( 0 , π ) 、V = 0 V=0 V = 0 、ψ = sin x > 0 \psi=\sin x\gt0 ψ = sin x > 0 とすると − ψ ′ ′ = ψ -\psi''=\psi − ψ ′′ = ψ なので、E = 1 E=1 E = 1 として定理から ∫ φ ′ 2 ≥ ∫ φ 2 \int\varphi'^2\ge\int\varphi^2 ∫ φ ′2 ≥ ∫ φ 2 です。基底状態変換の等式から、等号に近づくのは φ / sin \varphi/\sin φ / sin がほぼ定数のときで、定数 1 1 1 は最良です(− Δ D -\Delta_D − Δ D の最低固有値が 1 1 1 )。
演習 8.2 標準
n = 2 n=2 n = 2 、Ω \Omega Ω を単位円板、u = 1 − ∣ x ∣ 2 u=1-|x|^2 u = 1 − ∣ x ∣ 2 とする。∫ Ω ( Δ u ) 2 − ∫ Ω ∣ D 2 u ∣ 2 \int_\Omega(\Delta u)^2-\int_\Omega|D^2u|^2 ∫ Ω ( Δ u ) 2 − ∫ Ω ∣ D 2 u ∣ 2 と、定理 4.1 の証明の境界項 ∫ ∂ Ω ( ∂ ν u ) 2 tr D 2 ψ d S \int_{\partial\Omega}(\partial_\nu u)^2\operatorname{tr}D^2\psi\,dS ∫ ∂ Ω ( ∂ ν u ) 2 tr D 2 ψ d S をそれぞれ計算し、一致することを確かめよ。
解答 D 2 u = − 2 I D^2u=-2I D 2 u = − 2 I なので ∣ D 2 u ∣ 2 = 8 |D^2u|^2=8 ∣ D 2 u ∣ 2 = 8 、( Δ u ) 2 = 16 (\Delta u)^2=16 ( Δ u ) 2 = 16 で、左辺は 8 π 8\pi 8 π です。境界では ∂ ν u = − 2 \partial_\nu u=-2 ∂ ν u = − 2 、円の曲率は 1 1 1 なので tr D 2 ψ = 1 \operatorname{tr}D^2\psi=1 tr D 2 ψ = 1 (n − 1 = 1 n-1=1 n − 1 = 1 個の主曲率の和)で、右辺は 4 ⋅ 2 π = 8 π 4\cdot2\pi=8\pi 4 ⋅ 2 π = 8 π です。
演習 8.3 標準
定理 5.6 で F F F が定数なら、u u u は B 1 B_1 B 1 で調和(に a.e. で等しい)であることを示し、評価が [ ∇ u ] α ; B 1 / 2 ≤ C ∥ ∇ u ∥ L 2 ( B 1 ) [\nabla u]_{\alpha;B_{1/2}}\le C\|\nabla u\|_{L^2(B_1)} [ ∇ u ] α ; B 1/2 ≤ C ∥∇ u ∥ L 2 ( B 1 ) になることを、命題 3.4[ラプラス方程式と調和関数] から直接確かめよ。
解答 F F F が定数なら ∫ F ⋅ ∇ φ = F ⋅ ∫ ∇ φ = 0 \int F\cdot\nabla\varphi=F\cdot\int\nabla\varphi=0 ∫ F ⋅ ∇ φ = F ⋅ ∫ ∇ φ = 0 なので、補題 5.3 から u u u は調和です。h = ∂ k u − ( ∂ k u ) B 1 h=\partial_ku-(\partial_ku)_{B_1} h = ∂ k u − ( ∂ k u ) B 1 は調和で、補題 5.4 の証明と同じ評価で sup B 1 / 2 ∣ ∇ h ∣ ≤ C ∥ h ∥ L 2 ( B 1 ) ≤ C ∥ ∇ u ∥ L 2 ( B 1 ) \sup_{B_{1/2}}|\nabla h|\le C\|h\|_{L^2(B_1)}\le C\|\nabla u\|_{L^2(B_1)} sup B 1/2 ∣∇ h ∣ ≤ C ∥ h ∥ L 2 ( B 1 ) ≤ C ∥∇ u ∥ L 2 ( B 1 ) です。B 1 / 2 B_{1/2} B 1/2 は凸なので平均値の不等式から ∣ h ( x ) − h ( y ) ∣ ≤ sup ∣ ∇ h ∣ ∣ x − y ∣ ≤ sup ∣ ∇ h ∣ ∣ x − y ∣ α |h(x)-h(y)|\le\sup|\nabla h|\,|x-y|\le\sup|\nabla h|\,|x-y|^\alpha ∣ h ( x ) − h ( y ) ∣ ≤ sup ∣∇ h ∣ ∣ x − y ∣ ≤ sup ∣∇ h ∣ ∣ x − y ∣ α (∣ x − y ∣ ≤ 1 |x-y|\le1 ∣ x − y ∣ ≤ 1 )で、主張の評価が出ます。調和関数では Hölder どころか Lipschitz(実際には C ∞ C^\infty C ∞ )評価が成り立ちます。
L. C. Evans, Partial Differential Equations , 2nd ed., American Mathematical Society, 2010 — Chapter 5.8.2(差分商)、Chapter 6.3(正則性)、Chapter 6.4(最大値原理)。
D. Gilbarg, N. S. Trudinger, Elliptic Partial Differential Equations of Second Order , Springer, 2001(Classics in Mathematics 版)— Chapters 7–8(弱微分の連鎖律・弱最大値原理)。
M. Giaquinta, L. Martinazzi, An Introduction to the Regularity Theory for Elliptic Systems, Harmonic Maps and Minimal Graphs , 2nd ed., Edizioni della Normale, 2012 — Chapters 4–5(Campanato 空間と Schauder 評価)。
P. Grisvard, Elliptic Problems in Nonsmooth Domains , Pitman, 1985 — Chapter 3(凸領域の H 2 H^2 H 2 正則性)。
H. L. Cycon, R. G. Froese, W. Kirsch, B. Simon, Schrödinger Operators with Application to Quantum Mechanics and Global Geometry , Springer, 1987 — Section 2.4(Allegretto–Piepenbrink の定理)。
執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 4 回(Claude Code)・最終 2026-10-12