グリーン関数を重ね合わせて、任意のソースに従う微分方程式を解く方法を示すアニメーション。
点源に従う微分方程式の解がわかっていて、微分演算子が線形である場合、それらを重ね合わせて一般源の解を構築することができます

数学においてグリーン関数(またはグリーン関数[1])は、指定された初期条件または境界条件を持つ領域上で定義された 不同次線形微分演算子インパルス応答です。

これは、が線型微分演算子である場合、

  • グリーン関数は方程式 の解でありディラックのデルタ関数です
  • 初期値問題の解は畳み込み( )である

重ね合わせの原理により線形常微分方程式(ODE)が与えられた場合まず各sについてを解きソースがデルタ関数の和であるため、 Lの線形性により、解もグリーン関数の和であることがわかります

グリーン関数は、1820年代にこの概念を初めて提唱したイギリスの数学者 ジョージ・グリーンにちなんで名付けられました。現代の線形偏微分方程式の研究では、グリーン関数は主に基本解の観点から研究されています

多体理論においては、この用語は物理学、特に量子場の理論空気力学航空音響学電気力学、地震学統計場の理論においても用いられ、数学的な定義に当てはまらないものも含め、様々な種類の相関関数を指す。量子場の理論では、グリーン関数伝播関数の役割を果たす

定義と用途

ユークリッド空間の部分集合上の関数に作用する線型微分作用素L = L ( x )の点sにおけるグリーン関数G ( x , s )は、

ここでδはディラックのデルタ関数である。グリーン関数のこの性質は、次のような微分方程式を解くのに利用できる。

Lが自明でない場合、グリーン関数は一意ではありません。しかし実際には、対称性境界条件、および/またはその他の外部から課せられた基準の組み合わせによって、一意のグリーン関数が得られます。グリーン関数は、満たされる境界条件の種類に応じて、グリーン関数番号によって分類できます。グリーン関数は必ずしも実変数の関数ではありませんが、一般的には超関数の意味で理解されています

グリーン関数は波動方程式拡散方程式を解く際にも有用なツールである。[2] [3]量子力学において、ハミルトニアンのグリーン関数は状態密度の概念と重要な関連を持つ重要な概念である

物理学で用いられるグリーン関数は、通常、逆の符号で定義されます。つまり、 この定義は、ディラックのデルタ関数の 均一性 により、グリーン関数の特性を大きく変えるものではありません。

演算子が並進不変である場合、つまり、 がxに関して定数係数を持つ場合、グリーン関数は畳み込みカーネルと見なすことができます。つまり、 この場合、グリーン関数は線形時不変システム理論のインパルス応答 と同じです。

モチベーション

大まかに言えば、そのような関数G が演算子Lに対して見つかる場合グリーン関数の式 1にf ( s )を掛け、 sについて積分すると、次式が得られます。 演算子は線形であり、変数xにのみ作用し(積分変数sには作用しない)、演算子を積分の外側に取り出すと、次式が得られます。 これは、

方程式の解である

したがって、式1のグリーン関数と式2の右辺のソース項の知識があれば、関数u ( x )を得ることができます。このプロセスは、演算子Lの線形性に依存しています。

言い換えれば、式 2の解u ( x )は、式 3で与えられた積分によって決定できますf ( x )は既知ですが、 Gが既知でなければこの積分実行できません。ここでの問題は、式 1を満たすグリーン関数G を見つけることです。このため、グリーン関数は演算子Lの基本解と呼ばれることもあります。

すべての演算子がグリーン関数を持つわけではありません。グリーン関数はL右逆関数と考えることもできます。特定の演算子のグリーン関数を求めることの難しさに加え、式3の積分を評価するのは非常に難しい場合があります。しかし、この方法は理論的に正確な結果を与えます。

これは、ディラックのデルタ関数基底(fを 上に射影し、各射影上の解を重ね合わせたもの)に従ったfの展開として考えることができます。このような積分方程式はフレドホルム積分方程式として知られており、その研究がフレドホルム理論を構成します。

非同次境界値問題を解くためのグリーン関数

数学におけるグリーン関数の主な用途は、非同次境界値問題を解くことです。現代の理論物理学では、グリーン関数はファインマン図における伝播関数としても用いられることが多く、グリーン関数という用語は、あらゆる相関関数を指す場合にも用いられます

フレームワーク

シュトゥルム・リウヴィル作用素、すなわち形式が成り立つ線型微分作用素と し、をベクトル値境界条件作用素 とする。

を における連続関数とするさらに、問題が「正則」であるとする 。すなわち、すべてのxに対して の唯一の解はであるとする[a]

定理

を満たす 解は 1 つだけ存在し 、それは によって与えられます 。 ここで、は次の条件を満たすグリーン関数です。

  1. はおよびにおいて連続です
  2. のために  
  3. のために  
  4. 派生的な「ジャンプ」: . 
  5. 対称性: . 

グリーン関数の発達と発達遅延

グリーン関数は必ずしも一意ではない。なぜなら、同次方程式の任意の解をあるグリーン関数に加えると、別のグリーン関数が得られるからである。したがって、同次方程式が非自明な解を持つ場合、複数のグリーン関数が存在する。ある種の境界値問題初期値問題では、 に対してのみ零でないグリーン関数を求める必要がある。この場合、その解は遅延グリーン関数と呼ばれることがある。[4]同様に、 に対してのみ零でないグリーン関数は、高度グリーン関数と呼ばれる。[5]このような場合、2つのグリーン関数の任意の線形結合も有効なグリーン関数となる。高度グリーン関数と遅延グリーン関数はどちらも片側グリーン関数と呼ばれ、定義域内のすべての に対して零でないグリーン関数は両側グリーン関数と呼ばれる。[6]

高度な、あるいは遅延したという用語は、変数 x が時間に対応する場合に特に有用である。このような場合、遅延グリーン関数を用いて得られる解は過去の発生源のみに依存し因果的であるのに対し、高度なグリーン関数を用いて得られる解は未来の発生源のみに依存し非因果的である。これらの問題では、因果的解が物理的に重要な場合が多い。しかし、高度なグリーン関数は、境界データから発生源を求める逆問題の解を求める際に有用である。高度なグリーン関数と遅延グリーン関数は、特に非同次電磁波方程式の解の解析においてよく用いられる[7]

グリーン関数の発見

ユニット

グリーン関数の形は一義的に決まるわけではありませんが、次元解析を行ってグリーン関数の単位を求めることは、他の方法で求めたグリーン関数の妥当性確認において重要な意味を持ちます。定義式を簡単に調べると、 の単位は の単位だけでなく、位置ベクトルと が要素となっている空間の数と単位にも依存することがわかります。このことから、次の関係式が導かれます。 ここでは「 の物理的単位[さらに説明が必要]と定義され、 は空間(または時空)の体積要素です

たとえば、および時間が唯一の変数である場合、次のようになります。 ダランベール演算子および空間が 3 次元である 場合、次のようになります。

固有値展開

微分演算子 L完全な固有ベクトル の集合Ψn ( x ) (つまり、 LΨn = λnΨnなる関数Ψnとスカラーλn集合を許容する場合、これらの固有ベクトル固有値からグリーン関数 構築することができます

「完全」とは、関数の集合{Ψn }次の完全性関係を満たすことを意味する。

すると次のことが成り立つ。

ここで、 は複素共役を表します。

この方程式の各辺に演算子Lを適用すると、想定されていた完全性関係が得られます。

上記の形式で書かれたグリーン関数と、その固有ベクトルによって形成される関数空間との関係についての一般的な研究は、フレドホルム理論として知られています。

グリーン関数を求める方法は他にもいくつかあり、像法変数分離法、ラプラス変換法などがあります。[8]

ヴロンスキアンによる表現

を で定義される一般の線型2階微分作用素とする

とが一緒になって同次問題に対する線型独立な解の基底を形成すると仮定する。グリーン関数 の同次境界条件が与えられれば、を要求することで構築することができる。これらの条件を満たすグリーン関数は、 の連続性およびその導関数「ジャンプ」とともに、次のように書くことができる。

ここで、 はおよび のWronskian行列式として知られています。これはやや限定的なケースですが、Wronskian は、境界導関数に関する条件 (ノイマン条件) や、単一境界上の関数とその正規導関数に関する一対の条件 (コーシー条件) を含む、片側 (前進/後進) グリーン関数を必要とする他の境界値問題にも頻繁に現れます

グリーン関数の組み合わせ

微分演算子を として因数分解できる場合、 のグリーン関数はおよびのグリーン関数から構築できます 上記の恒等式はを の右逆演算子の表現とみなすことで直ちに得られます。これはによって定義される可逆線形演算子 がその行列要素 によって表現される方法に似ています

導関数 のスカラー多項式である微分演算子については、さらに恒等式が成り立ちます代数 の基本定理は、が自身と可換であるという事実と組み合わさって、多項式が因数分解されて、次の形式になることを保証します。 ここで、 の零点ですをとの両方についてフーリエ変換すると次のようになります。 次に、部分分数分解 を使用して分数を和に分割してから、と空間に戻すフーリエ変換を行うことができます。 この処理により、グリーン関数の積分とその和を関連付ける恒等式が得られます。 たとえば、 の場合、のグリーン関数の 1 つの形式は次のようになります。 提示された例は解析的に扱いやすいものですが、積分が自明でない場合 (たとえば、 が多項式の演算子である場合) に機能する処理を示しています。

グリーン関数の表

次の表は、頻繁に現れる微分演算子のグリーン関数の概要を示しています。ここで、、ヘビサイドステップ関数ベッセル関数第1種修正ベッセル関数、は第2種修正ベッセル関数です[9]最初の列に時間( t )が現れている場合は、遅延(因果)グリーン関数がリストされています。

ラプラシアンのグリーン関数

ラプラシアンを含む線形微分演算子のグリーン関数は、グリーン恒等式の 2 番目を使用して簡単に使用できます

グリーンの定理を導くには、発散定理ガウスの定理とも呼ばれる)から始める。

をガウスの法則に代入します

∇演算子の積の法則 を計算して適用する。

これを発散定理に当てはめるとグリーン定理が得られる。

線型微分作用素Lラプラシアン, ∇ 2であり、そのラプラシアンに対するグリーン関数Gが存在すると仮定する。グリーン関数の定義的性質は依然として成り立つ。

グリーンの2番目の恒等式を導入し、グリーンの恒等式を見てみましょうそして、

この式を用いることで、ノイマン境界条件またはディリクレ境界条件のいずれかの下で、ラプラス方程式 2 φ ( x ) = 0またはポアソン方程式 2 φ ( x ) = − ρ ( x )を解くことができます。言い換えれば、 (1) φ ( x )の値が体積の境界面上で指定されている場合(ディリクレ境界条件)、または(2) φ ( x )の法線微分が境界面上で指定されている場合(ノイマン境界条件)、体積内のあらゆる場所でφ ( x )を解くことができます。

問題が領域内の φ ( x )を解くことであると仮定する。すると、ディラックのデルタ関数の定義性により、積分 は単純にφ ( x )となり、以下の式が得られる。

この形式は調和関数のよく知られた特性すなわち境界面上の値または正規導関数が既知であれば、体積内の関数の値はどこでも既知であるという特性を表現します。

静電気学ではφ ( x )は電位ρ ( x )電荷 密度、法線微分は電場の法線成分として 解釈されます。

ディリクレ境界値問題を解く問題であれば、グリーン関数はxまたはx′が境界面上にあるときにG ( x , x ′)がゼロになるように選択する必要があります。こうすることで、面積分の2つの項のうち1つだけが残ります。ノイマン境界値問題を解く問題であれば、境界面上で法線微分がゼロになるようにグリーン関数を選択するのが論理的に思えるかもしれません。しかし、グリーン関数を定義する微分方程式にガウスの定理を適用すると、 G ( x , x ′) の法線微分は境界面上でゼロにはならないことがわかります。なぜなら、法線微分は境界面上で1に積分されなければならないからです。[10]

正規微分が最も単純な形をとるのは定数、すなわち1/ Sです。ここでSは表面の表面積です。解の表面積項は、 表面上の電位の平均値 です。この値は一般には知られていませんが、多くの場合、目標は電位そのものではなく、電位の勾配によって与えられる電場を求めることであるため、重要ではありません。

境界条件がない場合、ラプラス方程式のグリーン関数(3変数ラプラス方程式のグリーン関数)は

境界面が無限大であると仮定し、この式をグリーン関数に代入すると、最終的に電荷密度に関する電位の標準的な式が得られる。

グリーン関数番号が X11 である 次の問題のグリーン関数を求めます。

第一段階:手元の線形演算子のグリーン関数は、次の解として定義される。

ならばデルタ関数はゼロとなり、一般解は

の場合、 における境界条件は、 かつ であること意味します

の場合、 における境界条件は次を意味する 。

同様の理由から、 の式も省略されます。

これまでの結果をまとめると次のようになります。

ステップ 2:次のタスクは、 とを決定することです

グリーン関数の連続性を確保することは

定義微分方程式(すなわち式 *)を から まで積分し、 がゼロに近づくときの極限をとることで、一階微分における適切な不連続性を確保できます。残りの項は構成上連続となるため、二階微分のみを積分することに注意してください。

2つの(不)連続方程式を解く

したがって、この問題のグリーン関数は次のようになります。

その他の例

  • n = 1とし、 R全体を部分集合とする。L を とするこのとき、ヘヴィサイドステップ関数Θ( xx 0 )はLのx 0におけるグリーン関数となる
  • n = 2し、部分集合を1/4平面{( x , y ) : x , y ≥ 0}とし、Lをラプラシアンとする。また、 x = 0ディリクレ境界条件y = 0でノイマン境界条件が課されると仮定する。このとき、X10Y20グリーン関数は
  • 、これら3つすべてが実数の元であるとします。すると、区間 で積分可能な - 階導関数を持つ任意の関数に対して、 が成り立ちます上式のグリーン関数 は一意ではありません。を に加えると、がすべての に対してを満たす場合(例えば の場合、この式はどのように変形されますかまた、上式をを中心とするテイラー級数の形式と比較してください

参照

脚注

  1. ^ 専門用語では、「正則」とは、同次問題( )に対して自明な( )のみが存在することを意味します

参考文献

  1. ^ Wright, MCM (2006-10-01). 「グリーン関数かグリーンの関数か?」Nature Physics . 2 (10): 646. Bibcode :2006NatPh...2..646W. doi :10.1038/nphys411. ISSN  1745-2473.
  2. ^ Khorrami, Mohammad (2021). 「異なる次元における波動方程式のグリーン関数」.ミシュコルツ数学ノート. 22 (2): 721. doi : 10.18514/mmn.2021.2922 . ISSN  1787-2405.
  3. ^ Yu, Xiao; Lan, Kunquan; Wu, Jianhong (2021). 「グリーン関数、線形2階微分方程式、および1次元拡散移流モデル」 .応用数学研究. 147 (1): 319– 362. doi :10.1111/sapm.12384. ISSN  1467-9590.
  4. ^ Chu, Yi-Zen; Starkman, Glenn D. (2011-11-10) 「宇宙論と重力物理学のための摂動時空における遅延グリーン関数」、Physical Review D84 (12) 124020、arXiv : 1108.1825Bibcode :2011PhRvD..84l4020C、doi :10.1103/PhysRevD.84.124020
  5. ^ Soares de Castro, Antonio (2024-09-01). 「一次微分方程式のための遅延グリーン関数と高度グリーン関数」 . Physica Scripta . 99 (9): 096101. Bibcode :2024PhyS...99i6101S. doi :10.1088/1402-4896/ad6bc7. ISSN  0031-8949.
  6. ^ カマルゴ、ルーベンス・デ・フィゲイレド;キアッキオ、アリー・オロジンボ。デ・オリベイラ、エドムンド・カペラス (2013-03-06)。 「グリーンの片面関数と両面関数」。境界値の問題2013 (1): 45.土井: 10.1186/1687-2770-2013-45hdl : 11449/113526ISSN  1687-2770。
  7. ^ ジャクソン、ジョン・デヴィッド (2009).古典電気力学(第 3 版、[Nachdr.] 版)。ニューヨーク州ホーボーケン:ワイリー。ISBN 978-0-471-30932-1
  8. ^ Cole, KD; Beck, JV; Haji-Sheikh, A.; Litkouhi, B. (2011). 「グリーン関数を求める方法」.グリーン関数を用いた熱伝導. Taylor and Francis. pp.  101– 148. ISBN 978-1-4398-1354-6
  9. ^いくつかの例は Schulz、Hermann (2001)から引用。 Physik mit Bleistift: das Analytische Handwerkszeug des Naturwissenschaftlers (4. Aufl ed.)。フランクフルト・アム・マイン: ドイツ語。ISBN 978-3-8171-1661-4
  10. ^ ジャクソン、ジョン・デイビッド (1998年8月14日).古典電気力学. ジョン・ワイリー・アンド・サンズ. p. 39.

引用文献

  • Bayin, SS (2006).科学と工学における数学的手法. Wiley. 第18章と第19章.
  • アイゲス、レナード(1972年)『古典電磁場』ニューヨーク:ドーバー出版、ISBN 0-486-63947-9
    第 5 章では、グリーン関数を使用して静電気の境界値問題を解決する方法について非常に分かりやすく説明しています。
  • ポリアニン, AD; ザイツェフ, VF (2003). 『常微分方程式の厳密解ハンドブック(第2版)』 ボカラトン, フロリダ州: チャップマン・アンド・ホール/CRCプレス. ISBN 1-58488-297-2
  • バートン、ガブリエル (1989). 『グリーン関数と伝播の要素:ポテンシャル、拡散、波動』 . オックスフォード科学出版. オックスフォード:ニューヨーク:クラレンドン・プレス;オックスフォード大学出版局. ISBN 978-0-19-851988-1
    手順を説明したグリーン関数の教科書。
  • ポリアニン, AD (2002).エンジニアと科学者のための線形偏微分方程式ハンドブック. フロリダ州ボカラトン: チャップマン・アンド・ホール/CRCプレス. ISBN 1-58488-299-9
  • マシューズ、ジョン; ウォーカー、ロバート・L. (1970). 『物理学の数学的方法』(第2版). ニューヨーク: WAベンジャミン. ISBN 0-8053-7002-1
  • Folland, GB フーリエ解析とその応用数学シリーズ Wadsworth and Brooks/Cole.
  • グリーン、G (1828). 『電気と磁気の理論への数学的解析の応用に関する試論』ノッティンガム、イギリス: T. ホイールハウス、10-12ページ。
  • Faryad and, M.; Lakhtakia, A. (2018). 電磁気学における無限空間二項グリーン関数. ロンドン, イギリス / サンラファエル, カリフォルニア: IoP Science (イギリス) / Morgan and Claypool (アメリカ).書誌コード:2018idgf.book.....F. ISBN 978-1-68174-557-2
  • Şeremet, VD (2003).グリーン関数と行列ハンドブック. サウサンプトン: WIT Press. ISBN 978-1-85312-933-9