グリーン関数

グリーン関数（グリーンかんすう、英: Green's function）とは、微分方程式や偏微分方程式の解法の一つであるグリーン関数法に現れる関数である。グリーン関数法は、英国の数学者ジョージ・グリーンによって考案された。 物理学、数学、工学各分野において非常に重要な関数であり、広い用途で使用される。物理学におけるグリーン関数はプロパゲーター（伝播関数）とも呼ばれる。 J. A. Green により導入された組合せ論的関数のことをグリーン関数と呼ぶこともある。これはグリーン多項式とも呼ばれる。有限シュバレー群（オリジナルは有限体上の一般線型群）の既約表現を記述する数学的対象である。 下の偏微分方程式の（初期値）境界値問題を例に考える。 L y ( x ) = − f ( x ) x ∈ Ω y ( x ) = y 1 ¯ x ∈ Γ 1 ∂ y ( x ) ∂ n = y 2 ¯ x ∈ Γ 2 {\displaystyle {\begin{aligned}\mathrm {L} y({\boldsymbol {x}})&=-f({\boldsymbol {x}})&{\boldsymbol {x}}&\in \Omega \\y({\boldsymbol {x}})&={\bar {y_{1}}}&{\boldsymbol {x}}&\in \Gamma _{1}\\{\frac {\partial y({\boldsymbol {x}})}{\partial n}}&={\bar {y_{2}}}&{\boldsymbol {x}}&\in \Gamma _{2}\end{aligned}}} ここで、 L {\displaystyle \mathrm {L} } は微分作用素、 Ω {\displaystyle \Omega } は領域であり、領域の境界 Γ {\displaystyle \Gamma } は、 y 1 ¯ {\displaystyle {\bar {y_{1}}}} が規定されている境界 Γ 1 {\displaystyle \Gamma _{1}} と、 ∂ y / ∂ n {\displaystyle \partial y/\partial n} が規定されている境界 Γ 2 {\displaystyle \Gamma _{2}} からなり、 Γ 1 ∪ Γ 2 = Γ {\displaystyle \Gamma _{1}\cup \Gamma _{2}=\Gamma } 、 Γ 1 ∩ Γ 2 = ϕ {\displaystyle \Gamma _{1}\cap \Gamma _{2}=\phi } であるものとする。また、 n {\displaystyle n} は境界での外向き法線方向を示す。 上記の問題に対するグリーン関数 G ( x , x ′ ) {\displaystyle G({\boldsymbol {x}},{\boldsymbol {x'}})} とは次の条件を満たす関数のことである。 L G ( x , x ′ ) = − δ ( x − x ′ ) x ∈ Ω G ( x , x ′ ) = 0 x ∈ Γ 1 ∂ G ( x , x ′ ) ∂ n = 0 x ∈ Γ 2 {\displaystyle {\begin{aligned}\mathrm {L} G({\boldsymbol {x}},{\boldsymbol {x}}')&=-\delta ({\boldsymbol {x}}-{\boldsymbol {x}}')&{\boldsymbol {x}}&\in \Omega \\G({\boldsymbol {x}},{\boldsymbol {x}}')&=0&{\boldsymbol {x}}&\in \Gamma _{1}\\{\frac {\partial G({\boldsymbol {x}},{\boldsymbol {x}}')}{\partial n}}&=0&{\boldsymbol {x}}&\in \Gamma _{2}\end{aligned}}} ここに、x′ はソース点の位置を表す。 無限領域におけるグリーン関数を基本解という。 境界が単純（無限領域、半無限領域、無限平板領域など）でない場合にはグリーン関数を解析的に求めるのは大変困難である。 グリーン関数はもともと微分方程式の境界値問題に現れる関数である。物理学においても微分方程式を解くためにグリーン関数を用いることも多いが、量子物理学ではこれを拡張して使っている。つまり物理学においてグリーン関数は2通りの意味で扱われている。 境界値問題における微分方程式の主要解を意味し、与えられた全ての境界条件・初期条件を満足する。物理学では、微分方程式を直接解く代わりに、まず単純な点源問題の解であるグリーン関数を求めた後、重ね合わせの原理によって微分方程式の解をグリーン関数を用いて表す。 ある物理系を構成する個々の状態間の相関関数を与える関数として使われ、位置や時間などで指定されたある状態から他の状態への伝達(伝播)の特性を表す。詳細はプロパゲーターやグリーン関数 (多体理論)を参照。 電磁気学におけるポアソン方程式 Δ φ ( r ) = − ρ ( r ) {\displaystyle \Delta \varphi ({\boldsymbol {r}})=-\rho ({\boldsymbol {r}})} の解 φ ( r ) {\displaystyle \varphi ({\boldsymbol {r}})} を求めたい。この方程式の解として積分方程式 φ ( r ) = ∫ G ( r , r ′ ) ρ ( r ′ ) d r ′ {\displaystyle \varphi ({\boldsymbol {r}})=\int G({\boldsymbol {r}},{\boldsymbol {r'}})\rho ({\boldsymbol {r'}})d{\boldsymbol {r'}}} を仮定し、ポアソン方程式に代入するとグリーン関数 G ( r , r ′ ) {\displaystyle G({\boldsymbol {r}},{\boldsymbol {r'}})} の満たすべき式が得られる。 Δ G ( r , r ′ ) = − δ ( r − r ′ ) {\displaystyle \Delta G({\boldsymbol {r}},{\boldsymbol {r'}})=-\delta ({\boldsymbol {r}}-{\boldsymbol {r'}})} これを解くために両辺をフーリエ変換すると、 G ( r , r ′ ) {\displaystyle G({\boldsymbol {r}},{\boldsymbol {r'}})} のフーリエ変換 g ( k , r ′ ) = e − i k ⋅ r ′ k 2 {\displaystyle g({\boldsymbol {k}},{\boldsymbol {r'}})={\frac {e^{-i{\boldsymbol {k\cdot r'}}}}{{\boldsymbol {k}}^{2}}}} が得られる 。 これを逆フーリエ変換するとグリーン関数 G ( r , r ′ ) {\displaystyle G({\boldsymbol {r}},{\boldsymbol {r'}})} が求まる。 G ( r , r ′ ) = 1 4 π | r − r ′ | {\displaystyle G({\boldsymbol {r}},{\boldsymbol {r'}})={\frac {1}{4\pi |{\boldsymbol {r}}-{\boldsymbol {r'}}|}}} よってポアソン方程式の解は次のように求まる。 φ ( r ) = ∫ ρ ( r ′ ) 4 π | r − r ′ | d r ′ {\displaystyle \varphi ({\boldsymbol {r}})=\int {\frac {\rho ({\boldsymbol {r'}})}{4\pi |{\boldsymbol {r}}-{\boldsymbol {r'}}|}}d{\boldsymbol {r'}}} 以上のことから、位置 r {\displaystyle {\boldsymbol {r}}} の点電荷が別の位置 r ′ {\displaystyle {\boldsymbol {r'}}} に作る静電ポテンシャルを表したものがグリーン関数であり、これを重ね合わせたものが電荷分布 ρ ( r ) {\displaystyle \rho ({\boldsymbol {r}})} の作る静電ポテンシャル φ ( r ) {\displaystyle \varphi ({\boldsymbol {r}})} であることがわかる。 微分演算子を線型演算子 L と見て、微分方程式 Lφ = −ρ を解きたいとき、一種の逆演算子 L−1 を求めることができれば、φ = −L−1ρ というように微分方程式を解くことができる。これは線型代数における連立方程式において、係数行列の逆行列を求めることができれば連立方程式を解くことができることと対応している。このような L−1 をグリーン演算子 (Green's operator, Green operator) という。グリーン演算子を行列表示したときの行列要素をグリーン関数という。 このようにグリーン関数を抽象的な演算子と考えて取り扱うことには次のような利点がある。 微分演算子や積分演算子だけでなく、第二量子化のような抽象的な演算子を用いた理論に対してもそのまま用いられる。それは定常状態のシュレーディンガー方程式においてハミルトニアンを第二量子化における演算子で書かれていると考えるだけである。 複雑な関係式を簡潔に見通しよく書ける場合があり、一般的な性質の議論を見通し良く行える 例えば、方程式 ( − H ^ 0 + E ) | ϕ ⟩ = 0 {\displaystyle (-{\hat {H}}_{0}+E)|\phi \rangle =0} のグリーン演算子 ˆG0 が満たすべき方程式は ( − H ^ 0 + E ) G ^ 0 = − 1 {\displaystyle (-{\hat {H}}_{0}+E){\hat {G}}^{0}=-1} である。これを形式的に解くと G ^ 0 = − 1 E − H ^ 0 {\displaystyle {\hat {G}}^{0}=-{\frac {1}{E-{\hat {H}}_{0}}}} である。このグリーン演算子を具体的に計算するには ˆH0 の固有ベクトルを用いて G ^ 0 = − ∑ n 1 E − H ^ 0 | ϕ n ⟩ ⟨ ϕ n | = − ∑ n 1 E − E 0 | ϕ n ⟩ ⟨ ϕ n | {\displaystyle {\hat {G}}^{0}=-\sum _{n}{\frac {1}{E-{\hat {H}}_{0}}}|\phi _{n}\rangle \langle \phi _{n}|=-\sum _{n}{\frac {1}{E-E_{0}}}|\phi _{n}\rangle \langle \phi _{n}|} のように展開する。ただし E = E0 となるときは発散してしまう。それを避けるため分母をわずかに虚数軸方向にずらすことでこの問題は解消される G ^ ± 0 = ∑ n 1 E 0 − E ± i ϵ | ϕ n ⟩ ⟨ ϕ n | {\displaystyle {\hat {G}}_{\pm }^{0}=\sum _{n}{\frac {1}{E_{0}-E\pm i\epsilon }}|\phi _{n}\rangle \langle \phi _{n}|} たとえば ˆH0 = −Δ の場合のグリーン演算子の行列要素は、固有値を E = k2 として、次のように書ける。 G ± 0 ( r , r ′ ) = ⟨ r | G ^ ± 0 | r ′ ⟩ = 1 4 π | r − r ′ | e ± i k | r − r ′ | {\displaystyle G_{\pm }^{0}({\boldsymbol {r}},{\boldsymbol {r}}')=\langle {\boldsymbol {r}}|{\hat {G}}_{\pm }^{0}|{\boldsymbol {r}}'\rangle ={\frac {1}{4\pi |{\boldsymbol {r}}-{\boldsymbol {r}}'|}}e^{\pm ik|{\boldsymbol {r}}-{\boldsymbol {r}}'|}} ここで G + 0 {\displaystyle G_{+}^{0}} は外向き、 G − 0 {\displaystyle G_{-}^{0}} は内向きの球面波で、波が r' から r へ伝播する様子を示すものである。 散乱理論の形式論では、グリーン関数が用いられる。その基本方程式にもグリーン関数が含まれ、リップマン‐シュウィンガー方程式と呼ばれる。 | ψ ( ± ) ⟩ = | ϕ ⟩ + G 0 ^ ± V ^ | ψ ( ± ) ⟩ {\displaystyle |\psi ^{(\pm )}\rangle =|\phi \rangle +{\hat {G_{0}}}^{\pm }{\hat {V}}|\psi ^{(\pm )}\rangle \,} シュレーディンガー方程式を厳密に解く事は一般的に非常に困難な場合が多いが、近似に解く手法の一つとして摂動論がある。以下では摂動論におけるグリーン関数の形式理論について解説する。 系のハミルトニアン ˆH が無摂動項 ˆH0 と摂動項 ˆV の和で与えられた (ˆH = ˆH0 + ˆV) とする。無摂動ハミルトニアン ˆH0 に対して固有値方程式 H ^ 0 ϕ i ( 0 ) = E i ( 0 ) ϕ i ( 0 ) {\displaystyle {\hat {H}}_{0}\phi _{i}^{(0)}=E_{i}^{(0)}\phi _{i}^{(0)}} が成り立つ（例：ハートリー-フォック近似など）。 ω − ˆH0 を微分作用素として考えると非摂動グリーン関数 G (0)(ω)は 以下のように定義される（ここでデルタ関数 δ(x − x') は形式的に 1 とした）。 ( ω − H ^ 0 ) G ( 0 ) ( ω ) = − 1 {\displaystyle (\omega -{\hat {H}}_{0})G^{(0)}(\omega )=-1} 次に摂動ハミルトニアン ˆV で展開すると、 G ( ω ) = G ( 0 ) ( ω ) + G ( 0 ) ( ω ) V ^ G ( 0 ) ( ω ) + G ( 0 ) ( ω ) V ^ G ( 0 ) ( ω ) V ^ G ( 0 ) ( ω ) + ⋯ = G ( 0 ) ( ω ) + G ( 0 ) ( ω ) V ^ G ( ω ) {\displaystyle {\begin{aligned}G(\omega )&=G^{(0)}(\omega )+G^{(0)}(\omega ){\hat {V}}G^{(0)}(\omega )+G^{(0)}(\omega ){\hat {V}}G^{(0)}(\omega ){\hat {V}}G^{(0)}(\omega )+\dotsb \\&=G^{(0)}(\omega )+G^{(0)}(\omega ){\hat {V}}G(\omega )\end{aligned}}} この式の両辺に ω − ˆH0 を作用させ変形すると、摂動グリーン関数は次の関係を満たしていることがわかる。 ( ω − H ^ 0 − V ^ ) G ( ω ) = − 1 {\displaystyle (\omega -{\hat {H}}_{0}-{\hat {V}})G(\omega )=-1} また、この摂動グリーン関数が満たす関係式は ( H ^ 0 + V ^ ) ψ i = E i ψ i {\displaystyle ({\hat {H}}_{0}+{\hat {V}})\psi _{i}=E_{i}\psi _{i}} に対応している。 場の量子論や物性論においては、シュレーディンガー方程式に対するグリーン関数ではなくて、むしろ場の演算子に対する方程式に関連したものをグリーン関数と名付けて有効に用いている。それらの方程式は相互作用がない場合は、例えばスカラー場に対してクライン-ゴルドン方程式となるように、既に知られた方程式と同形のものになり、グリーン関数としても同じものとなる。しかし相互作用がある場合は方程式が非線形となり、摂動論的な扱いを除いて、古典的なグリーン関数の理論との対応を失う。 ^ Δ G ( r , r ′ ) = 1 ( 2 π ) 3 ∫ Δ g ( k , r ′ ) e i k ⋅ r d k = − 1 ( 2 π ) 3 ∫ g ( k , r ′ ) k 2 e i k ⋅ r d k {\displaystyle \Delta G({\boldsymbol {r}},{\boldsymbol {r'}})={\frac {1}{(2\pi )^{3}}}\int \Delta g({\boldsymbol {k,r'}})e^{i{\boldsymbol {k}}\cdot {\boldsymbol {r}}}d{\boldsymbol {k}}=-{\frac {1}{(2\pi )^{3}}}\int g({\boldsymbol {k,r'}})k^{2}e^{i{\boldsymbol {k}}\cdot {\boldsymbol {r}}}d{\boldsymbol {k}}} − δ ( r − r ′ ) = − 1 ( 2 π ) 3 ∫ e i k ⋅ ( r − r ′ ) d k {\displaystyle -\delta ({\boldsymbol {r}}-{\boldsymbol {r'}})=-{\frac {1}{(2\pi )^{3}}}\int e^{i{\boldsymbol {k}}\cdot ({\boldsymbol {r-r'}})}d{\boldsymbol {k}}} 両辺を比較すると、 g ( k , r ′ ) = e − i k ⋅ r ′ k 2 {\displaystyle g({\boldsymbol {k,r'}})={\frac {e^{-i{\boldsymbol {k}}\cdot {\boldsymbol {r'}}}}{{\boldsymbol {k}}^{2}}}} ^ 『物理学辞典』 培風館、1984年 ^ * 小泉義晴『微分方程式と量子統計力学のグリーン関数＜講義・演習＞』東海大学出版会、2010年。ISBN 978-4-486-01887-2。  ^ 今村勤『物理とグリーン関数』岩波書店、1976年。  ヘルムホルツ方程式 ラプラス方程式 境界要素法 インパルス応答 寺沢寛一:「自然科学者のための数学概論 増訂版」、岩波書店、ISBN 978-4000054805（1983年5月18日）# 第14章「境界値問題」。 小泉義晴：「現代工学のための量子物理学とグリーン関数 : 講義・演習ノート」，現代工学社， ISBN 978-4874721308（1987年1月25日）。 篠崎寿夫：「現代工学のための偏微分方程式とグリーン関数」，現代工学社，ISBN 978-4874721254 (1987年2月1日）。 篠崎寿夫：「現代工学のための常微分方程式とグリーン関数」，現代工学社，ISBN 978-4874721247 (1996年8月1日)。 今村勤: 「物理とグリーン関数」，岩波書店，ISBN 978-4000077194 (2016年2月18日)。 V.D. Seremet: ”Handbook of Green's Functions and Matrices”, WIT Press, ISBN 978-1-85312-933-9 (2002). 小泉義晴：「微分方程式と量子統計力学のグリーン関数―講義・演習」，東海大学出版会，ISBN 978-4486018872 (2010年11月)。 村上雅人:「なるほどグリーン関数」，海鳴社，ISBN 978-4875253549（2021年5月20日）。 亀高惟倫, 永井敦, 山岸弘幸:「グリーン関数」， 裳華房，ISBN 978-4785315979 (2022年11月25日)。