Greenova funkce

Greenova  funkce je funkce používaná k řešení lineárních nehomogenních diferenciálních rovnic s okrajovými podmínkami (nehomogenní okrajová úloha ). Pojmenován po anglickém matematikovi George Greenovi , který tuto teorii poprvé rozvinul ve 30. letech 19. století.

Greenovy funkce jsou užitečné v elektrostatice - pro řešení Poissonovy rovnice ; v teorii kondenzované hmoty  umožňují řešit rovnici difúze (a s ní splývající rovnici tepla); v kvantové mechanice  je Greenova funkce Hamiltoniána jednou z klíčových funkcí a souvisí s hustotou stavů. Greenovy funkce používané v těchto oblastech jsou velmi podobné, protože difúzní rovnice a Schrödingerova rovnice jsou si v určitém smyslu podobné. Všechny oblasti matematické a teoretické fyziky , kde jsou Greenovy funkce extrémně užitečné, je možná obtížné dokonce vyjmenovat. Pomáhají nacházet stacionární i nestacionární řešení, a to i za různých okrajových podmínek.

V částicové fyzice a statistické fyzice se Greenovy funkce používají jako propagátory ve Feynmanových diagramech (a výraz „Greenova funkce“ se často obecně aplikuje na funkci korelace v kvantové teorii pole ). Greenova funkce je široce používána v aplikacích teorie rozptylu ve fyzice pevných látek ( rentgenová difrakce , výpočty elektronových spekter kovových materiálů).

Definice a použití

Greenova funkce lineárního diferenciálního operátora působícího na zobecněné funkce na podmnožině euklidovského prostoru v bodě je libovolné řešení rovnice

,

kde  je funkce Dirac delta . Tato vlastnost Greenovy funkce může být použita k řešení diferenciální rovnice tvaru

,

Greenova funkce je inverzní operátor k , takže je často symbolicky označována jako .

Pokud je jádro operátoru netriviální, pak Greenova funkce není jedinečná. V praxi však použití principu symetrie, okrajových podmínek nebo jiných doplňkových podmínek umožňuje určit konkrétní Greenovu funkci. Obecně řečeno, Greenova funkce není obyčejná, ale zobecněná funkce , to znamená, že může vypadnout z třídy obyčejných funkcí, například mít vlastnosti tvaru delta funkce nebo jejích derivátů.

Greenova funkce je také užitečným nástrojem pro řešení vlnové rovnice, difúzní rovnice a kvantově mechanických rovnic, kde hraje zásadní roli Greenova funkce Hamiltonova operátoru a souvisí s hustotou stavů . Ve fyzice je Greenova funkce obvykle definována s opačným znaménkem:

,

který výrazně nemění jeho vlastnosti.

Pokud je operátor translačně invariantní , to znamená, že má konstantní koeficienty vzhledem k , pak lze Greenovu funkci zvolit jako konvoluční operátor .

.

V tomto případě se shoduje s impulsní přechodovou funkcí z teorie lineárních stacionárních systémů .

Poznámka

Někdy, když nehomogenní rovnice obsahuje na pravé straně konstantní koeficient, to znamená, že má tvar

.

V tomto případě se řešení původní nehomogenní rovnice s libovolnou funkcí na pravé straně zapíše jako

.
  1. Je zřejmé, že rozdíl v definici funkce Greena popsané v této části od definice uvedené v článku výše se netýká podstaty věci, ale pouze preferované formy zápisu

Greenova funkce Sturm-Liouvilleova operátoru (jednorozměrný případ)

Prohlášení o problému

Nechť je operátor Sturm  - Liouville , lineární diferenciální operátor tvaru:

,

a nechť  je operátor okrajové podmínky:

Greenova věta

Dovolit být  spojitá funkce na intervalu . Předpokládejme také, že úkol

je regulární, to znamená, že existuje pouze triviální řešení homogenního problému.

Pak je tu unikátní řešení vyhovující systému

,

což je dáno výrazem

,

kde  je funkce Greena, která splňuje následující požadavky (jsou to také vlastnosti funkce Greena):

  1. kontinuální v a .
  2. Pro , .
  3. Pro , .
  4. Derivační skok: .
  5. Symetrické: .

Nalezení Greenovy funkce

Jako řada přes vlastní funkce operátora

Pokud je množina vlastních vektorů ( vlastních funkcí ) diferenciálního operátoru

(to je množina takových funkcí , že pro každou existuje číslo , které )

je kompletní, pak lze zkonstruovat Greenovu funkci pomocí vlastních vektorů a vlastních hodnot .

Úplnost systému funkcí znamená naplnění vztahu

.

Dá se to ukázat

.

Ve skutečnosti, jednáme -li s tímto součtem jako operátor, získáme delta funkci (kvůli vztahu úplnosti).

(Nadčárka, , označuje komplexní konjugaci ; jsou-li  reálné funkce , lze ji vynechat).

Pro parabolické rovnice

Tepelná rovnice , Schrödingerova rovnice a difúzní rovnice mohou být reprezentovány jako parciální diferenciální rovnice :

 

 

 

 

(jeden)

kde  je hermitovský operátor , jsou prostorové souřadnice

 - teplota ,.

 je vlnová funkce , .

 je koncentrace látky, .

Vlastní funkce operátoru tvoří úplný ortonormální systém a splňují rovnici

.

Předpokládejme, že řešení rovnice (1) lze znázornit jako:

 

 

 

 

(2)

Dosazením navrženého tvaru řešení do rovnice (1) získáme:

.

Takto:

.

Tato rovnice musí platit pro všechna m. Dostaneme rovnici:

,

kde

.

Proto lze řešení původní rovnice (1) reprezentovat jako:

.

Uvážíme-li řady (2) rovnoměrně konvergentní, můžeme zjistit, že:

,

kde  je prvek hlasitosti.

Z tohoto vzorce vyplývá:

Pokud je tedy dán počáteční stav, pak

Tato rovnice může být zapsána v pohodlnější formě:

,

kde:

.

Tento výraz se nazývá Greenova funkce pro rovnici (1).

Greenova funkce pro Laplaciana

Greenovu funkci pro Laplaciana lze odvodit z Greenovy věty .

Abychom získali Greenovu větu, začněme Gaussovým zákonem :

.

Přijímáme a dosazujeme Gaussův zákon. Pojďme vypočítat a použít řetězové pravidlo pro operátora :

.

Dosazením výsledku do Gaussovy věty dostaneme Greenovu větu:

.

Za předpokladu, že náš lineární diferenciální operátor je Laplaciánský , , a že pro něj máme Greenovu funkci . Definici Greenovy funkce lze v tomto případě zapsat takto:

.

Vložili jsme Greenovu větu. Pak dostaneme:

.

Pomocí výrazu můžeme řešit Laplaceovu rovnici ( ) a Poissonovu rovnici ( ) s Neumannovými nebo Dirichletovými okrajovými podmínkami. Jinými slovy, můžeme najít řešení všude uvnitř dané oblasti, pokud (1) je dána hodnota na hranici této oblasti ( Dirichletovy okrajové podmínky ), nebo (2) je dána normální derivace na hranici této oblasti ( Neumannovy okrajové podmínky).

Zajímejme se o řešení uvnitř domény. V tomto případě se integrál zjednoduší na kvůli hlavní vlastnosti funkce delta a máme:

.

Tento vzorec vyjadřuje známou vlastnost harmonických funkcí , která spočívá v tom, že pokud je známa hodnota normální derivace na hranici oblasti, pak jsou všechny hodnoty funkce v jakémkoli vnitřním bodě této oblasti také známý.

V elektrostatice se jím rozumí elektrostatický potenciál , hustota elektrického náboje a normální derivace jako normální složka elektrického pole.

Při řešení Dirichletovy okrajové úlohy se volí Greenova funkce ve tvaru . Tato funkce zmizí, když nebo je na rozhraní; a naopak, při řešení Neumannovy okrajové úlohy bychom měli zvolit Greenovu funkci tak, aby její normální derivace zmizela na povrchu. V integrálu tak po povrchu zůstává pouze jeden ze dvou členů.

Při absenci okrajových podmínek má Greenova funkce pro Laplaciána tvar:

.

Uvažujeme-li hraniční plochu jako nekonečně velkou a dosadíme-li do tohoto výrazu Greenovu funkci, dojdeme k podobnému vyjádření pro elektrický potenciál ve smyslu hustoty elektrického náboje .

.

Příklad

(Tento příklad slouží jako ilustrace k odstavci Greenova funkce operátoru Sturm-Liouville (jednorozměrný případ) a zde popsané úvahy ilustrují body věty z odpovídajícího odstavce, odkazy na body, které jsou uvedeny v text níže).

Zadán úkol

; .

Najděte Greenovu funkci.

První krok: Greenova funkce v tomto případě musí být podle definice řešením rovnice

 

 

 

 

(3)

kde dva tahy označují druhou derivaci s ohledem na .

Pro , kde je -funkce rovna nule, je tato rovnice redukována na homogenní (bod 2 zmíněné věty):

,

to znamená, že pro všechny body kromě , Greenova funkce bude řešením takové homogenní rovnice.

Obecné řešení takové rovnice

,

kde a  jsou konstanty (nezávisí na ).

Všude tedy musí mít přesně tento tvar, kromě bodu , navíc nalevo a napravo od něj mohou mít (a budou) koeficienty a různé hodnoty.

Na Greenovu funkci klademe okrajové podmínky, které se shodují s okrajovými podmínkami původního problému (bod 3 věty zmíněné v úvodní poznámce). Greenova funkce s takto uloženými okrajovými podmínkami je výhodná, protože řešení konstruovaná součtem nebo integrací takových Greenových funkcí tyto okrajové podmínky automaticky splní.

Z levé okrajové podmínky:  - kladené na Greenovu funkci vidíme, že pro obecné řešení musí být koeficient nulový, tj.

.

Stejně tak z pravé okrajové podmínky:  - získáme koeficient rovný nule , tedy pro

.

V důsledku toho, vezmeme-li v úvahu, že koeficienty a obecně řečeno mohou záviset na , můžeme napsat:

Druhý krok:

Musíme definovat a .

Integrací dvojnásobku levé a pravé strany rovnice (3) s delta funkcí na pravé straně vidíme, že Greenova funkce musí být spojitá (bod 1 zmíněné věty), a tedy podmínka pro shodu řešení a :

.

Integrací levé a pravé části stejné rovnice od do získáme podmínku pro skok první derivace (bod 4 věty) a pomocí ní získáme:

.

Pomocí Cramerova pravidla nebo jednoduše uhodnutím řešení soustavy těchto dvou rovnic to dostaneme

.

Tyto výrazy splňují podmínku bodu 5 věty.

Pak Greenova funkce problému:

,

které lze napsat jako

Tabulka s Greenovými funkcemi

Tato tabulka uvádí Greenovy funkce pro běžně se vyskytující diferenciální operátory, kde , ,  je Heavisideova funkce ,  je Besselova funkce , je modifikovaná Besselova funkce prvního druhu a  je modifikovaná Besselova funkce druhého druhu . [2] Kde se v prvním sloupci objevuje čas ( t ) a jsou zobrazeny kauzální Greenovy funkce .

Diferenciální operátor L Greenova funkce G Příklad aplikace
, Harmonický oscilátor
, Poissonova rovnice
, Poissonova rovnice
stacionární 3D Schrödingerova rovnice pro volnou částici
v prostoru s rozměry Potenciální Yukawa , Propagátor
1D vlnová rovnice
2D vlnová rovnice
3D vlnová rovnice
1D difúzní rovnice
2D difúzní rovnice
3D difúzní rovnice
1D Klein-Gordonova rovnice
2D Klein-Gordonova rovnice
3D Klein-Gordonova rovnice
telegrafní rovnice
2D relativistická rovnice tepla
3D relativistická rovnice tepla

Další příklady

Viz také

Poznámky

  1. Li Tsung-dao Matematické metody ve fyzice. - M.: Mir, 1965. - str. 200
  2. Některé příklady jsou převzaty z Schulz, Hermann: Physik mit Bleistift. Frankfurt nad Mohanem: Deutsch, 2001. ISBN 3-8171-1661-6 (německy)

Literatura