Lippmann- Schwinger- Gleichung

Quantenmechanikvorlesung von Prof. Dr. E. Schöll, PhD

65px|Kein GFDL

Der Artikel Lippmann- Schwinger- Gleichung basiert auf der Vorlesungsmitschrift von Franz- Josef Schmitt des 6.Kapitels (Abschnitt 1) der Quantenmechanikvorlesung von Prof. Dr. E. Schöll, PhD.

   |}}

{{#set:Urheber=Prof. Dr. E. Schöll, PhD|Inhaltstyp=Script|Kapitel=6|Abschnitt=1}} Kategorie:Quantenmechanik __SHOWFACTBOX__

Man betrachte Teilchen, die in Wechselwirkung stehen, jedoch keine gebundenen Zustände miteinander einnehmen:

Der Hamiltonoperator{{#set:Fachbegriff=Hamiltonoperator|Index=Hamiltonoperator}} kann geschrieben werden als:

$\hat{H}={\hat{H}}^{(0)}+{\hat{H}}^{(1)}$

Dabei bezeichne ${\hat{H}}^{(0)}$ die kinetische Energie{{#set:Fachbegriff=kinetische Energie|Index=kinetische Energie}} und ${\hat{H}}^{(1)}$ die Wechselwirkungsenergie{{#set:Fachbegriff=Wechselwirkungsenergie|Index=Wechselwirkungsenergie}}.

stationäre Streuung

Im Falle stationärer Streuung{{#set:Fachbegriff=stationärer Streuung|Index=stationärer Streuung}} erhalten wir:

$\hat{H}|\mathrm{\Psi }⟩=E|\mathrm{\Psi }⟩.$

$|\mathrm{\Psi }⟩$ beschreibt ein am Anfang einlaufendes Teilchen (ohne Wechselwirkung), die anschließende Streuung und schließlich wieder auseinanderlaufende Teilchen:

BILD WW: Streuung

Stationär bedeutet hier: Das Gleichgewicht hat sich bereits eingeregelt. Der Prozess ist zeitlich stationär, weil jede Veränderung an einem Teilchenzustand 1 durch ein nachrückendes Teilchen, dessen Zustand sich in den des ersten (Zustand 1) begibt, aufgefüllt wird.

Die Schrödingergleichung{{#set:Fachbegriff=Schrödingergleichung|Index=Schrödingergleichung}} lautet:

$(E-{\hat{H}}_0)|\mathrm{\Psi }⟩={\hat{H}}^{(1)}|\mathrm{\Psi }⟩$

Erster Schritt bei derartigen Problemen: Isolation der Störung!

Die formale Lösung kann angegeben werden mittels:

$\begin{array}{rl}& |\mathrm{\Psi }⟩=|\mathrm{\Phi }⟩+\frac{1}{(E-{\hat{H}}_0)}{\hat{H}}^{(1)}|\mathrm{\Psi }⟩\\ & \frac{1}{(E-{\hat{H}}_0)}:={(E-{\hat{H}}_0)}^{-1}\\ \end{array}$

Die Division zwischen 1 und dem Operator der linken Seite ist dabei als Ausführung der inversen Operation{{#set:Fachbegriff=inversen Operation|Index=inversen Operation}} zu verstehen!

$|\mathrm{\Phi }⟩$

ist eine beliebige Lösung der wechselwirkungsfreien Gleichung

$({\hat{H}}_0-E)|\mathrm{\Phi }⟩=0$

Beweis

$\begin{array}{rl}& (E-{\hat{H}}_0)|\mathrm{\Psi }⟩=(E-{\hat{H}}_0)|\mathrm{\Phi }⟩+(E-{\hat{H}}_0)\frac{1}{(E-{\hat{H}}_0)}{\hat{H}}^{(1)}|\mathrm{\Psi }⟩\\ & (E-{\hat{H}}_0)\frac{1}{(E-{\hat{H}}_0)}:=1\\ & \Rightarrow (E-{\hat{H}}_0)|\mathrm{\Psi }⟩={\hat{H}}^{(1)}|\mathrm{\Psi }⟩\iff (E-{\hat{H}}_0)|\mathrm{\Phi }⟩=0\\ \end{array}$

Bemerkung

Die Gleichung $|\mathrm{\Psi }⟩=|\mathrm{\Phi }⟩+\frac{1}{(E-{\hat{H}}_0)}{\hat{H}}^{(1)}|\mathrm{\Psi }⟩$

ist dabei eine Integralgleichung, beispielsweise in der Ortsdarstellung:

$⟨\overline{r}|\mathrm{\Psi }⟩=⟨\overline{r}|\mathrm{\Phi }⟩+\underset{}{\overset{}{\int }}\underset{}{\overset{}{\int }}⟨\overline{r}|\frac{1}{(E-{\hat{H}}_0)}|\overline{r}\stackrel{´}{}⟩⟨\overline{r}\stackrel{´}{}|{\hat{H}}^{(1)}|\overline{r}\stackrel{´}{}\stackrel{´}{}⟩⟨\overline{r}\stackrel{´}{}\stackrel{´}{}|\mathrm{\Psi }⟩d^{3}r\stackrel{´}{}\stackrel{´}{}{d}^{3}r\stackrel{´}{}$

Berechnung des inversen Operators

$\frac{1}{(E-{\hat{H}}_0)}$

Hier: Greenscher Operator{{#set:Fachbegriff=Greenscher Operator|Index=Greenscher Operator}}, sogenannte Resolvente{{#set:Fachbegriff=Resolvente|Index=Resolvente}} (auch: Residuum!) der Schrödingergleichung.

Methode: Transformation auf Impulsdarstellung (Fourier- Transformation) und komplexe Integration.

Aber: die Lösung ist nicht eindeutig, je nach Wahl des Integrationsweges ergeben sich unterschiedliche Egebnisse (Integrationsweg in der komplexen Ebene). Dementsprechend ergeben sich verschiedene Randbedingungen

Die Festlegung erfolgt durch Addition eines kleinen komplexen Terms $i\epsilon$. Am Schluss kann man dann $\epsilon \to 0$ gehen lassen.

Damit ergibt sich als Lippmann- Schwinger- Gleichung{{#set:Fachbegriff=Lippmann- Schwinger- Gleichung|Index=Lippmann- Schwinger- Gleichung}}

$|{\mathrm{\Psi }}^{(+)}⟩=|\mathrm{\Phi }⟩+\frac{1}{(E-{\hat{H}}_0+i\epsilon )}{\hat{H}}^{(1)}|{\mathrm{\Psi }}^{(+)}⟩$

Wesentlicher Vorteil zur Schrödingergleichung: Die Lippmann- Schwinger- Gleichung ist die im Vergleich zur Schrödingergleichung komplexe Erweiterung mit reeller Polstellenfreiheit!

Mit

auslaufender Welle

$|{\mathrm{\Psi }}^{(+)}⟩$

Streuwelle

$\frac{1}{(E-{\hat{H}}_0+i\epsilon )}{\hat{H}}^{(1)}|{\mathrm{\Psi }}^{(+)}⟩$ und

einlaufender Welle

$|\mathrm{\Phi }⟩$ (Lösung des ungestörten Problems)

Die auslaufende Welle $|{\mathrm{\Psi }}^{(+)}⟩$ ist die Summe aus einlaufender Welle und Streuwelle!

Greensche Funktion des freien Teilchens

(= Ortsdartellung des Greenschen Operators)

$G_{+}(\overline{r},\overline{r}\stackrel{´}{}):=\frac{\hslash }{2m}⟨\overline{r}|\frac{1}{(E-{\hat{H}}_0+i\epsilon )}|\overline{r}\stackrel{´}{}⟩$

Dabei werden zwei "Einsen" eingeschoben und wir gewinnen:

$G_{+}(\overline{r},\overline{r}\stackrel{´}{})=\frac{\hslash }{2m}\underset{}{\overset{}{\int }}{d}^{3}q\underset{}{\overset{}{\int }}{d}^{3}q\stackrel{´}{}⟨\overline{r}|\overline{q}⟩⟨\overline{q}|\frac{1}{(E-{\hat{H}}_0+i\epsilon )}|\overline{q}\stackrel{´}{}⟩⟨\overline{q}\stackrel{´}{}|\overline{r}\stackrel{´}{}⟩$

Also: Das Problem wird nach der auslaufenden Welle aufgelöst. Zur Polstellenfreiheit erweitert man komplex.

Dann isoliert man den Greenschen Operator und führt mit diesem eine Fouriertransormation{{#set:Fachbegriff=Fouriertransormation|Index=Fouriertransormation}} durch!

Der obige Einschub einer Basis ist noch keine Fouriertransformation. Wir befinden uns dann immer noch im Ortstraum!

Dabei bezeichnen $\overline{q},\overline{q}\stackrel{´}{}$ die Wellenvektoren mit der entsprechenden Impulsdarstellung $\hslash \overline{q}$.

Für ein freies Teilchen, für das der Hamiltonian direkt angegeben werden kann ${\hat{H}}_0=\frac{{\hat{\overline{p}}}^2}{2m}$ gilt:

$⟨\overline{q}|{\hat{H}}_0|\overline{q}\stackrel{´}{}⟩=\frac{{\hslash }^2{\overline{q}}^2}{2m}\delta (\overline{q}-\overline{q}\stackrel{´}{})$

Somit also

$⟨\overline{q}|\frac{1}{E-{\hat{H}}_0+i\epsilon }|\overline{q}\stackrel{´}{}⟩=\frac{\delta (\overline{q}-\overline{q}\stackrel{´}{})}{E-\frac{{\hslash }^2{\overline{q}}^2}{2m}+i\epsilon }$

Asymptotisch gelte für das einlaufende Teilchen als Anfangsbedingung sozusagen

$\overline{p}=\hslash \overline{k}\Rightarrow E=\frac{{\hslash }^2{\overline{k}}^2}{2m}$

$\begin{array}{rl}& \Rightarrow ⟨\overline{q}|\frac{1}{E-{\hat{H}}_0+i\epsilon }|\overline{q}\stackrel{´}{}⟩=\frac{2m}{{\hslash }^2}\frac{\delta (\overline{q}-\overline{q}\stackrel{´}{})}{{\overline{k}}^2-{\overline{q}}^2+i\eta }=:\frac{2m}{{\hslash }^2}{\tilde{G}}_{+}(\overline{q})\delta (\overline{q}-\overline{q}\stackrel{´}{})\\ & \eta =\frac{2m}{{\hslash }^2}\epsilon \\ & ⟨\overline{r}|\overline{q}⟩=\frac{1}{{\left(2\pi \right)}^{\frac{3}{2}}}{e}^{i\overline{q}\overline{r}}\\ \end{array}$

Damit folgt dann als angekündigte Fouriertrafo:

$\begin{array}{rl}& G_{+}(\overline{r},\overline{r}\stackrel{´}{})=\frac{1}{{\left(2\pi \right)}^3}\underset{}{\overset{}{\int }}{d}^{3}q{\tilde{G}}_{+}(\overline{q})e^{i\overline{q}(\overline{r}-\overline{r}\stackrel{´}{})}\\ & {\tilde{G}}_{+}(\overline{q})=\frac{1}{{\overline{k}}^2-{\overline{q}}^2+i\eta }\\ \end{array}$

()

Also: Wir führen die Fouriertransformation durch und gewinnen als Fouriertransformierte ${\tilde{G}}_{+}(\overline{q})=\frac{1}{{\overline{k}}^2-{\overline{q}}^2+i\eta }$.

Die Rücktransformation liefert die gesuchte Greensfunktion $G_{+}(\overline{r},\overline{r}\stackrel{´}{})$, die mittels Residuensatz aus der bekannten ${\tilde{G}}_{+}(\overline{q})=\frac{1}{{\overline{k}}^2-{\overline{q}}^2+i\eta }$ durch Fouriertrafo gewonnen werden kann!

$G_{+}(\overline{r},\overline{r}\stackrel{´}{})$ hängt also nur von $(\overline{r}-\overline{r}\stackrel{´}{})$ ab!

Berechnung von $G_{+}(\overline{r}-\overline{r}\stackrel{´}{}):=G_{+}(\overline{R})$ in Polarkoordinaten $\overline{q}$ erfolgt mittels Residuensatz{{#set:Fachbegriff=Residuensatz|Index=Residuensatz}}

$G_{+}(\overline{R})=\frac{1}{{\left(2\pi \right)}^3}\underset{}{\overset{}{\int }}{d}^{3}q\frac{1}{{\overline{k}}^2-{\overline{q}}^2+i\eta }{e}^{i\overline{q}\overline{R}}$

Dabei lege man $\overline{R}=\overline{r}-\overline{r}\stackrel{´}{}$ entlang der z- Achse, so dass zwischen $\overline{R}$ und $\overline{q}$ gerade der Winkel $\vartheta$ liegt:

$\begin{array}{rl}& G_{+}(\overline{R})=\frac{1}{{\left(2\pi \right)}^3}\underset{0}{\overset{\mathrm{\infty }}{\int }}dq\underset{-1}{\overset{1}{\int }}d\cos \vartheta \underset{0}{\overset{2\pi }{\int }}d\varphi \frac{q^2}{{\overline{k}}^2-{\overline{q}}^2+i\eta }{e}^{iqR\cos \vartheta }\\ & G_{+}(\overline{R})=\frac{1}{4{\pi }^{2}iqR}\underset{0}{\overset{\mathrm{\infty }}{\int }}dq{q}^2\frac{e^{iqR}-e^{-iqR}}{q({\overline{k}}^2-{\overline{q}}^2+i\eta )}\\ \end{array}$

Dies kann man leicht weiter zusammenfassen, indem im zweiten Term einfach q durch -q ersetzt wird:

$G_{+}(\overline{R})=\frac{1}{4{\pi }^{2}iR}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dqq\frac{e^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }$

Die Integration erfolgt mittels Residuensatz in der komplexen q- Ebene:

Dazu ist es nötig, die komplexe Zahl q in Polarkoordinaten umzuschreiben:

$\begin{array}{rl}& q=\rho \cdot e^{i\mathrm{\Phi }}\\ & 0\le \mathrm{\Phi }\le \pi \\ & dq=\rho \cdot e^{i\mathrm{\Phi }}id\mathrm{\Phi }\\ \end{array}$

Skizzenhaft:

Da die Integration im Unendlichen (Halbkreisbogen) verschwindet kann man das Integral von Minus bis Plus Unendlich auch gleich als Ringintegral schreiben. Wesentlich ist dann: dass es nur Beiträge aus den Polstellen der Funktion gibt. demnach müssen diese gesucht werden:

Die Pole des Integranden:

$\begin{array}{rl}& \frac{1}{{\overline{k}}^2-{\overline{q}}^2+i\eta }\\ & q_{1/2}=\pm \sqrt{{\overline{k}}^2+i\eta }\approx (k+\frac{i\eta }{2k})\\ \end{array}$

Genau genommen muss noch gezeigt werden, dass das Integral über den Kreisbogen für Radius gegen Unendlich verschwindet:

$\begin{array}{c}\lim \\ \rho \to \mathrm{\infty }\\ \end{array}\underset{}{\oint }dqq\frac{e^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dqq\frac{e^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }+\begin{array}{c}\lim \\ \rho \to \mathrm{\infty }\\ \end{array}\underset{0}{\overset{\pi }{\int }}d\mathrm{\Phi }i{e}^{2i\mathrm{\Phi }}\frac{{\rho }^2{e}^{i\rho R\cos \mathrm{\Phi }}{e^{-}}^{\rho R\sin \mathrm{\Phi }}}{{\overline{k}}^2-{\rho }^2{e}^{2i\mathrm{\Phi }}+i\eta }$

Aber:

$\begin{array}{rl}& \begin{array}{c}\lim \\ \rho \to \mathrm{\infty }\\ \end{array}\underset{0}{\overset{\pi }{\int }}d\mathrm{\Phi }i{e}^{2i\mathrm{\Phi }}\frac{{\rho }^2{e}^{i\rho R\cos \mathrm{\Phi }}{e^{-}}^{\rho R\sin \mathrm{\Phi }}}{{\overline{k}}^2-{\rho }^2{e}^{2i\mathrm{\Phi }}+i\eta }=0\\ & da\\ & \begin{array}{c}\lim \\ \rho \to \mathrm{\infty }\\ \end{array}{e^{-}}^{\rho R\sin \mathrm{\Phi }}=0\\ & \Rightarrow \begin{array}{c}\lim \\ \rho \to \mathrm{\infty }\\ \end{array}\underset{}{\oint }dqq\frac{e^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dqq\frac{e^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }\\ \end{array}$

Mittels Residuensatz ergibt sich dann

$\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dqq\frac{e^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }=2\pi i{\left(RES\left(\frac{q{e}^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }\right)\right)|}_{q=q_1}$

Dies ist der Beitrag vom oberen Integrationsweg, weshalb das Residuum an q=q1 ausgewertet werden muss. Ebenso hätte man den unteren Integrationsweg nehmen können. Dann wäre das Residuum eben an q=q2 auszuwerten gewesen. Da wir die Polstellen des Arguments des Residuums gefunden haben können wir umschreiben:

$\begin{array}{rl}& RES{\frac{q{e}^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }|}_{q=q_1}=RES{\frac{q{e}^{iqR}}{(\sqrt{{\overline{k}}^2+i\eta }-q)(\sqrt{{\overline{k}}^2+i\eta }+q)}|}_{q_1\equiv \sqrt{{\overline{k}}^2+i\eta }}\\ & \sqrt{{\overline{k}}^2+i\eta }=q_1=-q_2\\ & RES{\frac{q{e}^{iqR}}{(\sqrt{{\overline{k}}^2+i\eta }-q)(\sqrt{{\overline{k}}^2+i\eta }+q)}|}_{q_1\equiv \sqrt{{\overline{k}}^2+i\eta }}=\begin{array}{c}\lim \\ q->q1\\ \end{array}\frac{(q-q_1)q{e}^{iqR}}{(q_1-q)(q-q_2)}=\frac{q_1{e}^{i{q}_{1}R}}{(q_1-q_2)}\\ & =-\frac{e^{i\sqrt{{\overline{k}}^2+i\eta }R}}{2}\\ & \begin{array}{c}\lim \\ \eta \to 0\\ \end{array}-\frac{e^{i\sqrt{{\overline{k}}^2+i\eta }R}}{2}=-\frac{e^{ikR}}{2}\\ \end{array}$

Also hat man ein Ergebnis für $G_{+}(\overline{R})=\frac{1}{4{\pi }^{2}iR}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dqq\frac{e^{iqR}}{{\overline{k}}^2-{\overline{q}}^2+i\eta }$, man erhält

$G_{+}(\overline{R})=\frac{1}{4{\pi }^{2}iR}2\pi iRES{|}_{q_1}=\frac{-e^{ikR}}{4\pi R}$

()

Wesentlich: $G_{+}(\overline{R})=G_{+}(\overline{r}-\overline{r}\stackrel{´}{})$ erfüllt die Differenzialgleichung für die Greensche Funktion:

$(\mathrm{\Delta }+k^2)G_{+}(\overline{r}-\overline{r}\stackrel{´}{})=\delta (\overline{r}-\overline{r}\stackrel{´}{})$

Denn:

$\begin{array}{rl}& \delta (\overline{r}-\overline{r}\stackrel{´}{})=⟨\overline{r}|\overline{r}\stackrel{´}{}⟩=G_{+}(\overline{r},\overline{r}\stackrel{´}{})=⟨\overline{r}|(E-{\hat{H}}_0+i\epsilon )\frac{1}{(E-{\hat{H}}_0+i\epsilon )}|\overline{r}\stackrel{´}{}⟩\\ & \cong ⟨\overline{r}|(\frac{{\hslash }^2{k}^2}{2m}-\frac{{\hat{p}}^2}{2m})\frac{1}{(E-{\hat{H}}_0+i\epsilon )}|\overline{r}\stackrel{´}{}⟩=\frac{{\hslash }^2}{2m}(k^2+\mathrm{\Delta })⟨\overline{r}|\frac{1}{(E-{\hat{H}}_0+i\epsilon )}|\overline{r}\stackrel{´}{}⟩\\ \end{array}$

Ortsdarstellung der Lippmann- Schwinger- Gleichung

$\begin{array}{rl}⟨\overline{r}|{\mathrm{\Psi }}^{(+)}⟩& =⟨\overline{r}|\mathrm{\Phi }⟩+\underset{}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}⟨\overline{r}|\frac{1}{(E-{\hat{H}}_0+i\epsilon )}|\overline{r}\stackrel{´}{}⟩⟨\overline{r}\stackrel{´}{}|{\hat{H}}^1|{\mathrm{\Psi }}^{(+)}⟩\\ & =⟨\overline{r}|\mathrm{\Phi }⟩+\frac{2m}{{\hslash }^2}\underset{}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}{G}_{+}(\overline{r}-\overline{r}\stackrel{´}{})⟨\overline{r}\stackrel{´}{}|{\hat{H}}^1|{\mathrm{\Psi }}^{(+)}⟩=e^{i\overline{k}\overline{r}}+\frac{2m}{{\hslash }^2}\underset{}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\frac{-e^{ik\left|(\overline{r}-\overline{r}\stackrel{´}{})\right|}}{4\pi \left|(\overline{r}-\overline{r}\stackrel{´}{})\right|}⟨\overline{r}\stackrel{´}{}|{\hat{H}}^1|{\mathrm{\Psi }}^{(+)}⟩\\ \end{array}$

Mit der

durchlaufenden freien Welle

$⟨\overline{r}|\mathrm{\Phi }⟩=e^{i\overline{k}\overline{r}}$ und der

Streuwelle

$\frac{2m}{{\hslash }^2}\underset{}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}{G}_{+}(\overline{r}-\overline{r}\stackrel{´}{})⟨\overline{r}\stackrel{´}{}|{\hat{H}}^1|{\mathrm{\Psi }}^{(+)}⟩=+\frac{2m}{{\hslash }^2}\underset{}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\frac{-e^{ik\left|(\overline{r}-\overline{r}\stackrel{´}{})\right|}}{4\pi \left|(\overline{r}-\overline{r}\stackrel{´}{})\right|}⟨\overline{r}\stackrel{´}{}|{\hat{H}}^1|{\mathrm{\Psi }}^{(+)}⟩$

Zusammenfassung

Aus der Schrödingergleichung

Die Schrödingergleichung lautet:

$(E-{\hat{H}}_0)|\mathrm{\Psi }⟩={\hat{H}}^{(1)}|\mathrm{\Psi }⟩$

Mit dem linearen Differentialoperator $(E-{\hat{H}}_0)|\mathrm{\Psi }⟩$ und der Inhomogenität ${\hat{H}}^{(1)}|\mathrm{\Psi }⟩$

kann man formal lösen:

$\begin{array}{rl}& |{\mathrm{\Psi }}^{(+)}⟩=|\mathrm{\Phi }⟩+\frac{1}{(E-{\hat{H}}_0+i\epsilon )}{\hat{H}}^{(1)}|{\mathrm{\Psi }}^{(+)}⟩\\ & \frac{1}{(E-{\hat{H}}_0)}:={(E-{\hat{H}}_0)}^{-1}\\ \end{array}$

eine Form der Lippmann- Schwinger- Gleichung mit auslaufender Welle $|{\mathrm{\Psi }}^{(+)}⟩$ Greenschen Operator (auch sogenannte RESOLVENTE)$\frac{1}{(E-{\hat{H}}_0+i\epsilon )}$ und durchlaufender Welle (freie einlaufende Lösung) $|\mathrm{\Phi }⟩$

Die Berechnung der Greenschen Funktion des freien Teilchens:

Als Operator: ${\hat{G}}_{+}:=\frac{1}{(E-{\hat{H}}_0+i\epsilon )}$ erfüllt $(E-{\hat{H}}_0){\hat{G}}_{+}=1$

Übergang in die Impulsdarstellung:

$⟨\overline{q}|{\hat{G}}_{+}|\overline{q}\stackrel{´}{}⟩=\frac{2m}{\hslash }{\hat{G}}_{+}(\overline{q})\delta (\overline{q}-\overline{q}\stackrel{´}{})$

Mit

${\hat{G}}_{+}(\overline{q}):=\frac{1}{k^2-q^2+i\eta }$

Mittels Fouriertrafo erfolgt der Übergang in die Ortsdarstellung:

${\hat{G}}_{+}(\overline{r}-\overline{r}\stackrel{´}{}):=\frac{{\hslash }^2}{2m}⟨\overline{r}|{\hat{G}}_{+}(\overline{q})|\overline{r}\stackrel{´}{}⟩=-\frac{e^{ik|\overline{r}-\overline{r}\stackrel{´}{}|}}{4\pi |\overline{r}-\overline{r}\stackrel{´}{}|}$

Dieser erfüllt dann eine Relation des Impulsoperators in Ortsdarstellung (orts- Differeniationsrelation):

$(\mathrm{\Delta }+k^2){\hat{G}}_{+}(\overline{r}-\overline{r}\stackrel{´}{})=\delta (\overline{r}-\overline{r}\stackrel{´}{})$

(dies ist die skalare Helmholtzgleichung!)

Potenzialstreuungen

${\hat{H}}^{(1)}$ sei ein Potenzial, das die Wechselwirkung mit einem schweren Teilchen als Streuzentrum{{#set:Fachbegriff=Streuzentrum|Index=Streuzentrum}}(Target) beschreibt. Allgemein: Beschreibung im Schwerpunktsystem

Hier kann man als Schwerpunktsystem näherungsweise den Schwerpunkt des schweren Teilchen annehmen

In Ortsdarstellung schreiben wir:

$\begin{array}{rl}& ⟨\overline{r}\stackrel{´}{}|{\hat{H}}^{(1)}|{\mathrm{\Psi }}^{(+)}⟩=V(\overline{r}\stackrel{´}{}){\mathrm{\Psi }}^{(+)}(\overline{r}\stackrel{´}{})\\ & \Rightarrow {\mathrm{\Psi }}^{(+)}(\overline{r})=e^{i\overline{k}\overline{r}}-\frac{2m}{{\hslash }^2}\underset{}{\overset{}{\int }}{d}^{3}r\stackrel{´}{}\frac{e^{i\overline{k}|\overline{r}-\overline{r}\stackrel{´}{}|}}{4\pi |\overline{r}-\overline{r}\stackrel{´}{}|}V(\overline{r}\stackrel{´}{}){\mathrm{\Psi }}^{(+)}(\overline{r}\stackrel{´}{})\\ \end{array}$

Dies ist die Lippmann Schwingergleichung für eine Potenzialstreuung.

Diese Gleichung ist völlig äquivalent zur Schrödingergeleichung mit Randbedingungen.

Als Randbedingungen sind in der Streutheorie prinzipiell die asymptotischen Wege für r gegen Plus oder Minus UNENDLICH zu verstehen.