Mathematische Architektur

Die Schubart Master-Formel ist ein Doppelintegral über Energie und Winkel, dessen Integrand ein Produkt aus fünf Faktoren ist. Dieses Kapitel analysiert die formale Struktur dieser Beziehung: die Integrationsdomänen, die Dimensionen der beteiligten Größen und die Art und Weise, wie die einzelnen Faktoren zusammenwirken. Ziel ist eine vollständige mathematische Beschreibung des Rahmenwerks, unabhängig von der physikalischen Interpretation der einzelnen Terme.

Die erweiterte physikalische Darstellung

Die kanonische Schubart Master-Formel ist die kompakte Volumenform $P(t) = \eta \cdot \int_V \Phi_{\text{eff}}(\mathbf{r},t) \cdot \sigma_{\text{eff}}(E)\,dV$. Um ihre mathematische Architektur vollständig explizit zu machen, arbeitet dieses Kapitel mit ihrer erweiterten physikalischen Darstellung:

$$P(t) = \iint \Phi(E,\theta,t)\,\sigma(E)\,\Delta p\,v_{\text{ph}}\,\eta_{\text{conv}}\;dE\,d\theta$$

Der Integrand ist ein Produkt, in dem jeder Faktor eine wohldefinierte Rolle spielt. Die Flussdichte $\Phi(E,\theta,t)$ trägt die gesamte Abhängigkeit von den Integrationsvariablen $E$ und $\theta$ sowie von der Zeit. Der Streuquerschnitt $\sigma(E)$ hängt allein von der Energie ab. Die drei verbleibenden Faktoren — Impulsübertrag $\Delta p$, Phononengeschwindigkeit $v_{\text{ph}}$ und Wandlungswirkungsgrad $\eta_{\text{conv}}$ — werden im Grundmodell als über die Integrationsdomäne konstant behandelt, können aber in erweiterten Formulierungen selbst von $E$ und $\theta$ abhängen, wie in Zukünftige Erweiterungen diskutiert.

Es ist aufschlussreich, den Integranden explizit als differenzielle Leistung zu lesen. Schreibt man

$$\frac{\partial^2 P}{\partial E\,\partial \theta} = \Phi(E,\theta,t)\,\sigma(E)\,\Delta p\,v_{\text{ph}}\,\eta_{\text{conv}},$$

so wird deutlich, dass die Master-Formel die zweifach differenzielle Leistungsdichte im Energie-Winkel-Raum über die zulässige Domäne aufsummiert. Diese Lesart trennt sauber, was von den Integrationsvariablen abhängt (der linke Teil $\Phi\,\sigma$) und was im Grundmodell als konstanter Vorfaktor herausgezogen werden kann (das Produkt $\Delta p\,v_{\text{ph}}\,\eta_{\text{conv}}$). Sie ist zugleich der Ausgangspunkt jeder Diskretisierung: Ein numerisches Schema ersetzt die kontinuierliche Domäne durch ein Gitter und die Integration durch eine gewichtete Summe über die Gitterzellen.

Die Integrationsdomänen

Das Integral erstreckt sich über zwei Variablen. Die Energieintegration läuft über den relevanten Energiebereich des einfallenden Feldes, formal von einer unteren Schwelle bis zu einer oberen Grenze. Außerhalb dieses Bereichs verschwindet entweder der Fluss $\Phi$ oder der Streuquerschnitt $\sigma$, sodass keine Beiträge entstehen. Die Winkelintegration läuft über den relevanten Raumwinkelbereich, aus dem Fluss einfallen kann.

Die Wahl der Integrationsgrenzen ist Teil der Modellspezifikation und hängt vom betrachteten Anwendungsfall ab. Das Rahmenwerk selbst schreibt keine festen Grenzen vor; es verlangt lediglich, dass die Domänen so gewählt werden, dass alle physikalisch relevanten Beiträge erfasst und keine unphysikalischen Bereiche einbezogen werden. Die formale Behandlung der Grenzen erfolgt in Systemgrenzen.

Dimensionsanalyse

Ein zentraler Konsistenztest jeder physikalischen Gleichung ist die Dimensionsanalyse. Das Endergebnis $P(t)$ muss die Einheit Watt, also $\text{kg}\,\text{m}^2\,\text{s}^{-3}$, tragen. Man prüft dies, indem man die Dimensionen der Faktoren zusammensetzt.

Die Flussdichte $\Phi(E,\theta,t)$ ist pro Fläche, Zeit, Raumwinkel und Energie definiert. Der Streuquerschnitt $\sigma(E)$ trägt die Dimension einer Fläche, also $\text{m}^2$. Der Impulsübertrag $\Delta p$ hat die Dimension $\text{kg}\,\text{m}\,\text{s}^{-1}$, die Phononengeschwindigkeit $v_{\text{ph}}$ die Dimension $\text{m}\,\text{s}^{-1}$, und der Wandlungswirkungsgrad $\eta_{\text{conv}}$ ist dimensionslos. Die Integration über $dE$ und $d\theta$ trägt die Dimensionen Energie und Raumwinkel bei. In der Zusammensetzung heben sich Fläche, Energie und Raumwinkel gegen die entsprechenden inversen Dimensionen der Flussdichte weg, sodass am Ende die Leistungsdimension verbleibt. Diese Konsistenz ist eine notwendige Bedingung dafür, dass das Rahmenwerk überhaupt eine wohldefinierte Leistung liefert; die vollständige Buchführung findet sich in der Notationsreferenz.

Die multiplikative Struktur

Die Faktoren treten multiplikativ auf, was bedeutet, dass das Verschwinden eines einzigen Faktors das gesamte Produkt zum Verschwinden bringt. Diese Struktur kodiert die Bedingungskette aus Nichtgleichgewichts-Energiewandlung: Ohne Fluss, ohne Kopplung, ohne gerichteten Impuls, ohne Transport oder ohne Gleichrichtung gibt es keine Leistung. Jeder Faktor ist notwendig, keiner allein hinreichend.

Diese multiplikative Verkettung lässt sich als Faktorisierung des Integranden auffassen, in der jeder Faktor eine eigene physikalische Frage beantwortet: $\Phi$ fragt „wie viel fällt ein?“, $\sigma(E)$ „wie stark koppelt es?“, $\Delta p$ „wie gerichtet ist der Übertrag?“, $v_{\text{ph}}$ „wie schnell wird er abtransportiert?“ und $\eta_{\text{conv}}$ „welcher Anteil bleibt nutzbar?“. Weil diese Fragen konzeptionell unabhängig sind, lässt sich jede für sich modellieren und messen, ohne die übrigen festzulegen. Genau diese Trennbarkeit macht die Faktorisierung zum organisierenden Prinzip des gesamten Rahmenwerks — sie erlaubt es, konkurrierende physikalische Modelle allein über die funktionale Form einzelner Faktoren zu vergleichen, während die Integralstruktur unverändert bleibt.

Zugleich erlaubt die multiplikative Form eine einfache Sensitivitätsanalyse. Da die Leistung in jedem Faktor linear ist, wenn die übrigen konstant gehalten werden, überträgt sich eine relative Änderung eines Faktors unmittelbar in eine gleich große relative Änderung der Leistung. Diese Eigenschaft vereinfacht die numerische Untersuchung erheblich.

Sensitivität und logarithmische Ableitung

Die multiplikative Faktorisierung erlaubt eine besonders einfache Sensitivitätsanalyse, die für die numerische und experimentelle Praxis von großem Nutzen ist. Weil die Leistung ein Produkt der Faktoren ist, ist ihre relative Änderung gleich der Summe der relativen Änderungen der Faktoren. Formal ausgedrückt gilt für die logarithmische Ableitung

$$\frac{d\ln P}{P} = \frac{d\Phi}{\Phi} + \frac{d\sigma}{\sigma} + \frac{d(\Delta p)}{\Delta p} + \frac{dv_{\text{ph}}}{v_{\text{ph}}} + \frac{d\eta_{\text{conv}}}{\eta_{\text{conv}}},$$

sofern die Faktoren als konstant über die Domäne behandelt werden. Jeder Faktor trägt also mit einem Sensitivitätskoeffizienten von exakt eins zur relativen Änderung der Leistung bei. Eine Verdopplung eines beliebigen Faktors verdoppelt die Leistung; eine Halbierung halbiert sie. Diese Eigenschaft macht die multiplikative Struktur nicht nur konzeptionell klar, sondern auch numerisch bequem: Unsicherheiten in den Eingangsgrößen propagieren additiv in der logarithmischen Darstellung, und die dominierende Unsicherheitsquelle lässt sich unmittelbar identifizieren. Für die praktische Fehleranalyse bedeutet dies, dass die Gesamtunsicherheit von $P(t)$ von demjenigen Faktor bestimmt wird, dessen relative Unsicherheit am größten ist — ein Ergebnis, das ohne die Faktorisierung nur mit erheblichem Aufwand zu gewinnen wäre.

Linearität und Trennbarkeit

Die Gleichung ist linear im Fluss $\Phi$, was die in Mehrkanaliger Energiefluss beschriebene Superposition ermöglicht. Sie ist außerdem in der Energievariable separierbar, sofern $\Phi$ selbst faktorisiert. Nimmt man etwa an, dass sich die Flussdichte als Produkt einer Energie- und einer Winkelverteilung schreiben lässt, $\Phi(E,\theta,t) = f(E)\,g(\theta)\,h(t)$, so zerfällt das Doppelintegral in ein Produkt eindimensionaler Integrale:

$$P(t) = \Delta p\,v_{\text{ph}}\,\eta_{\text{conv}}\,h(t)\left(\int f(E)\,\sigma(E)\,dE\right)\left(\int g(\theta)\,d\theta\right).$$

Diese Separation reduziert den Rechenaufwand erheblich und macht die Zeitabhängigkeit von $P(t)$ im faktorisierten Fall vollständig durch $h(t)$ bestimmt. Sie gilt allerdings nur, wenn die Verteilung tatsächlich faktorisiert; korreliert die Winkelverteilung mit der Energie, bleibt das volle Doppelintegral erforderlich. Diese strukturellen Eigenschaften sind rein mathematischer Natur und gelten unabhängig davon, welcher physikalische Mechanismus den einzelnen Faktoren zugrunde liegt. Sie machen die Master-Formel zu einem Rahmen, der sich sowohl analytisch als auch numerisch handhaben lässt, und bilden die Grundlage für die in Numerische Simulation beschriebene Diskretisierung. Die vollständige Herleitung der Integralform findet sich in vollständige Herleitung, die eigenständige Dimensionskontrolle in Dimensionsanalyse.

Weiterführende Literatur

• G. B. Arfken, H. J. Weber, F. E. Harris: Mathematical Methods for Physicists. 7. Auflage, Academic Press, 2013. Mehrdimensionale Integration, Separierbarkeit und Dimensionsbetrachtungen.
• J. R. Taylor: Scattering Theory: The Quantum Theory of Nonrelativistic Collisions. Wiley, New York 1972 (Nachdruck Dover 2006). Formale Rolle des Streuquerschnitts im Integranden.
• N. W. Ashcroft, N. D. Mermin: Solid State Physics. Holt, Rinehart and Winston, New York 1976. Transportintegrale über Energie und Winkel im Festkörperkontext.