Samstag, 1. August 2026

Multiphysik, Materialdesign und Abschluss – vom atomistischen Mechanismus zum technischen Bauteil




Titelbild: Eigene Illustration der Skalenbrücke vom Elektron zum Bauteil sowie zentraler Kontinuumsgleichungen für Stofftransport, Wärme, Mechanik und Elektrostatik.

Materialien funktionieren nie durch einen einzelnen Mechanismus. Die elektrische Leitfähigkeit hängt von Elektronenstruktur und Defekten ab, die Diffusion von atomaren Sprungbarrieren, die Lebensdauer von Rissen und Grenzflächen, die Bauteiltemperatur von Wärmeentstehung und Kühlung. In einem technischen System wirken diese Prozesse gleichzeitig und beeinflussen sich gegenseitig.

Der Abschluss von Physikalische Chemie IV verfolgt daher eine übergeordnete Frage: Wie werden atomistische Mechanismen in verlässliche Vorhersagen für ein reales Bauteil übersetzt? Dazu müssen Modelle über viele Größen- und Zeitskalen gekoppelt, experimentell geprüft und mit Unsicherheiten versehen werden.

Einordnung: Die Teile 1 bis 19 entwickelten Kristallstruktur, Phononen, Defekte, Elektronen, Halbleiter, Magnetismus, Nanomaterialien, Oberflächen, Katalyse, Korrosion, Photovoltaik und Batteriematerialien. Teil 20 verbindet diese Themen zu einer gemeinsamen Methodik aus Multiskalenmodellierung, Multiphysik, Materialdesign und Lebenszyklusdenken.

Lernziele

Nach diesem Beitrag kannst du:
  • räumliche und zeitliche Skalen eines Materialproblems strukturieren,
  • Elektronenstruktur-, atomistische, mesoskopische und Kontinuumsmodelle einordnen,
  • Homogenisierung und repräsentative Volumenelemente erklären,
  • Stoff-, Ladungs-, Wärme- und Impulsbilanzen formulieren,
  • chemische, thermische und mechanische Rückkopplungen erkennen,
  • dimensionslose Kennzahlen für Regimeabschätzungen verwenden,
  • Phasenfeld- und Schädigungsmodelle einordnen,
  • Verifikation, Validierung und Kalibrierung unterscheiden,
  • Sensitivitäts- und Unsicherheitsanalysen durchführen,
  • Surrogatmodelle und digitale Zwillinge kritisch bewerten,
  • mehrzielige Materialoptimierung und Pareto-Fronten erklären,
  • Lebensdauer, Herstellbarkeit und Nachhaltigkeit in das Design integrieren.

1. Warum Skalenkopplung notwendig ist

Ein einzelnes Materialproblem kann Größenordnungen von Ångström bis Meter und von Femtosekunden bis Jahrzehnte umfassen. Elektronen relaxieren extrem schnell, Diffusion in Festkörpern kann dagegen Stunden oder Jahre benötigen. Ein direktes Modell aller Atome eines technischen Bauteils ist weder rechnerisch noch konzeptionell sinnvoll.

Skalenkopplung bedeutet deshalb, relevante Information gezielt von einer Beschreibungsebene zur nächsten zu übertragen.

2. Elektronenstruktur

Elektronenstrukturmethoden liefern unter anderem:

  • Bandstrukturen und Zustandsdichten,
  • Defektbildungsenergien,
  • Adsorptionsenergien,
  • Reaktionspfade und Aktivierungsbarrieren,
  • magnetische und dielektrische Eigenschaften.

Die typische räumliche Skala liegt bei wenigen Ångström bis Nanometern. Die Ergebnisse hängen von Austausch-Korrelations-Näherung, Zellgröße, Randbedingungen und Temperaturbehandlung ab.

3. Molekulardynamik

Molekulardynamik integriert klassische oder quantenmechanisch erzeugte Kräfte über die Zeit. Sie beschreibt:

  • Schwingungen und lokale Struktur,
  • Diffusion und Transport,
  • mechanische Antwort,
  • Grenzflächenbewegung,
  • thermische Fluktuationen.

Die zugänglichen Zeiten sind begrenzt. Seltene Ereignisse mit hohen Barrieren werden in einer kurzen Simulation möglicherweise nie beobachtet.

4. Kinetische Monte-Carlo-Methoden

Wenn die relevanten Ereignisse bekannt sind, können Sprungraten verwendet werden:

kjjexp(−Ea,j/kBT)

Kinetisches Monte Carlo überspringt die schnellen Schwingungen und simuliert direkt die Abfolge seltener Zustandswechsel. Die Methode ist effizient, solange der Ereigniskatalog vollständig ist.

5. Mesoskopische Modelle

Auf der Mesoskala werden einzelne Atome nicht mehr aufgelöst. Stattdessen verwendet man Felder oder statistische Objekte:

  • Phasenanteile,
  • Kornorientierungen,
  • Porenstrukturen,
  • Versetzungsdichten,
  • lokale Zusammensetzungen,
  • Schädigungsvariablen.

Diese Ebene verbindet atomistische Materialparameter mit beobachtbaren Mikrostrukturen.

6. Kontinuumsmechanik und Transportmodelle

Auf Bauteilskala werden Größen als kontinuierliche Felder behandelt. Dazu gehören Temperatur T, Konzentration c, Potential φ, Verschiebung u und Spannung σ.

Die zugrunde liegenden Bilanzgleichungen gelten nur, wenn das betrachtete Volumen groß genug gegenüber der Mikrostruktur ist.

7. Repräsentatives Volumenelement

Ein repräsentatives Volumenelement, kurz RVE, muss groß genug sein, um die statistische Mikrostruktur abzubilden, aber klein genug, um als lokaler Materialpunkt im Bauteilmodell zu dienen.

Seine Größe hängt von der Zielgröße ab. Für Wärmeleitung kann ein kleineres Volumen repräsentativ sein als für Rissausbreitung.

8. Homogenisierung

Homogenisierung ersetzt eine heterogene Mikrostruktur durch effektive Eigenschaften:

J̄=−Keff∇X̄

Die effektive Leitfähigkeit Keff ist im Allgemeinen ein Tensor und hängt von Volumenanteilen, Form, Orientierung, Konnektivität und Grenzflächenwiderständen ab.

9. Voigt- und Reuss-Grenzen

Für zwei Phasen mit Volumenanteil φ:

Kparallel=φK₁+(1−φ)K₂
1/Kseriell=φ/K₁+(1−φ)/K₂

Diese Grenzfälle zeigen, wie stark Mikrostruktur und Flussrichtung die effektive Eigenschaft verändern können.

10. Skalenübergang nach oben

Von kleiner zu großer Skala werden beispielsweise übertragen:

  • freie Energien,
  • Diffusionskoeffizienten,
  • Reaktionskonstanten,
  • Elastizitätstensoren,
  • Grenzflächenenergien,
  • effektive Leitfähigkeiten.

Jeder Parameter benötigt eine klare Definition und einen Gültigkeitsbereich.

11. Rückkopplung nach unten

Bauteilmodelle liefern lokale Zustände wie Temperatur, Druck, Stromdichte oder Spannung. Diese Größen können atomistische Barrieren und Defektpopulationen verändern.

Skalenkopplung ist daher häufig bidirektional.

12. Hierarchische und simultane Kopplung

  • Hierarchisch: kleinere Skala liefert Parameter, danach wird das größere Modell unabhängig gerechnet.
  • Simultan: mehrere Skalen werden während der Simulation gegenseitig aktualisiert.

Simultane Kopplung ist leistungsfähig, aber numerisch aufwendig und schwerer zu validieren.

13. Modellreduktion

Ein gutes Modell enthält nicht möglichst viel Physik, sondern die für die Zielgröße notwendige Physik. Modellreduktion kann erfolgen durch:

  • Zeitskalentrennung,
  • Quasigleichgewicht,
  • stationäre Näherung,
  • Homogenisierung,
  • reduzierte Basisverfahren,
  • datenbasierte Surrogate.

14. Gültigkeitsbereich

Jede Näherung besitzt Bedingungen. Eine konstante Diffusivität kann bei kleiner Konzentrationsänderung sinnvoll sein, aber in einem Zweiphasengebiet scheitern. Ein lineares Elastizitätsmodell gilt nicht bei plastischer Verformung oder großen Rissen.



Abbildung 1. Jede Skala besitzt eigene Zustandsgrößen und numerische Methoden. Verlässliche Skalenkopplung überträgt nur physikalisch identifizierbare und experimentell prüfbare Information.

15. Allgemeine Bilanzgleichung

Viele Kontinuumsmodelle lassen sich als Bilanz formulieren:

Akkumulation=−Divergenz des Flusses+Quelle−Senke

Für eine extensive Größe mit Dichte y und Fluss J:

∂y/∂t+∇·J=s

Diese Struktur gilt für Masse, Ladung, Energie und Impuls. Die physikalische Modellierung steckt in der Flussbeziehung und im Quellterm.

16. Stofftransport

Für eine neutrale Spezies mit Fick-Diffusion und Reaktion:

∂c/∂t=∇·(D∇c)+R(c,T,φ)

D kann von Konzentration, Temperatur, Spannung und Mikrostruktur abhängen. Bei anisotropen Kristallen ist D ein Tensor.

17. Transport geladener Spezies

Für ein Ion kommen Diffusion und Migration zusammen:

Ji=−Di∇ci−ziuiFci∇φ

In konzentrierten Elektrolyten müssen Aktivitäten, gekoppelte Flüsse und lokale Elektroneutralität berücksichtigt werden.

18. Poisson-Gleichung

∇·(ε∇φ)=−ρe

Sie koppelt Ladungsverteilung und elektrisches Potential. In makroskopisch elektroneutralen Bereichen wird häufig eine vereinfachte Stromerhaltung verwendet:

∇·i=0

19. Wärmeleitung

ρcp∂T/∂t=∇·(k∇T)+Q̇

Q̇ kann ohmische Wärme, Reaktionsenthalpie, Phasenumwandlungswärme, plastische Dissipation oder Strahlungsabsorption enthalten.

20. Mechanisches Gleichgewicht

∇·σ+b=ρ∂²u/∂t²

Bei langsamen Prozessen wird die Trägheit häufig vernachlässigt:

∇·σ+b=0

Die Materialgleichung verbindet Spannung und Dehnung.

21. Chemische Dehnung

Einlagerung, Defektbildung oder Phasenumwandlung verändern das Gittervolumen:

ε=εelchemthpl
εchem=β(c−c₀)

Räumliche Konzentrationsgradienten erzeugen dadurch Spannungsgradienten und Rissrisiken.

22. Thermische Dehnung

εthT(T−T₀)

Unterschiedliche Ausdehnungskoeffizienten in Verbundmaterialien erzeugen Grenzflächenspannungen. Temperaturzyklen können auch ohne chemische Reaktion zu Delamination führen.

23. Reaktions-Wärme-Kopplung

Eine Reaktion kann Wärme erzeugen, während ihre Geschwindigkeit nach Arrhenius mit der Temperatur wächst:

k(T)=k₀exp(−Ea/RT)

Damit entsteht eine positive Rückkopplung. Thermisches Durchgehen, Hotspots in Katalysatoren und lokale Batterieüberhitzung sind Beispiele.

24. Transport-Struktur-Kopplung

Poren, Risse und Phasengrenzen verändern Leitfähigkeit und Diffusionswege. Umgekehrt kann Transport neue Phasen, Poren oder Risse erzeugen. Eine konstante effektive Eigenschaft ist dann unzureichend.

25. Phasenfeldmethode

Ein Ordnungsparameter η beschreibt kontinuierlich verschiedene Phasen oder Kornorientierungen. Eine freie Energie lautet schematisch:

F=∫[f(c,η,T)+κ/2|∇η|²]dV

Der Gradiententerm erzeugt eine endliche Grenzflächenbreite und Grenzflächenenergie.

26. Allen–Cahn- und Cahn–Hilliard-Dynamik

Für einen nicht erhaltenen Ordnungsparameter:

∂η/∂t=−LδF/δη

Für eine erhaltene Zusammensetzung:

∂c/∂t=∇·[M∇(δF/δc)]

Diese Gleichungen beschreiben Kornwachstum, Entmischung, Ausscheidung und bewegte Phasengrenzen.

27. Bruch- und Schädigungsmodelle

Risse können explizit durch diskrete Geometrie oder diffus durch eine Schädigungsvariable d beschrieben werden. Phasenfeldbruch ersetzt die scharfe Rissfläche durch eine endliche Übergangszone.

Die Kopplung an chemische Dehnung und Transport ermöglicht Modelle für Elektrodenrisse, Korrosionsrisse und thermische Delamination.

28. Randbedingungen

Ein Gleichungssystem ist ohne Rand- und Anfangsbedingungen unvollständig. Typische Bedingungen sind:

  • vorgegebene Konzentration oder Fluss,
  • vorgegebenes Potential oder Strom,
  • Temperatur oder Wärmeübergang,
  • Verschiebung oder Oberflächenlast,
  • periodische Bedingungen eines RVE.

Falsche Randbedingungen können stärker wirken als kleine Parameterfehler.

29. Dimensionsanalyse

Durch Skalierung werden Variablen dimensionslos:

x*=x/L, t*=t/t₀, c*=c/c₀

Dadurch erscheinen Verhältnisse charakteristischer Zeiten und Kräfte. Diese Kennzahlen zeigen dominante Regime und erleichtern den Vergleich verschiedener Bauteilgrößen.

30. Fourier-Zahl

Fo=αt/L²

Fo misst, wie weit thermische Diffusion während der Zeit t fortgeschritten ist. Bei Fo≫1 hat sich Temperatur über die Länge L weitgehend ausgeglichen.

31. Biot-Zahl

Bi=hL/k

Bi vergleicht äußeren Wärmeübergang mit innerer Wärmeleitung. Bei Bi≪0,1 kann ein Körper häufig als räumlich isotherm behandelt werden.

32. Péclet-Zahl

Pe=vL/D

Pe vergleicht konvektiven mit diffusivem Stofftransport. Große Werte erzeugen dünne Grenzschichten und starke Richtungsabhängigkeit.

33. Damköhler-Zahl

Da=kL²/D

Da vergleicht Reaktion und Diffusion. Kleine Werte bedeuten schnelle Durchmischung gegenüber Reaktion, große Werte starke Konzentrationsgradienten.

34. Thiele-Modul

φ=L√(k/D)

Der Thiele-Modul ist eng mit Da verwandt. Er wird besonders für Reaktion in porösen Körnern und Elektroden verwendet.

35. Mechanische Kennzahlen

Für gekoppelte chemisch-mechanische Probleme sind Verhältnisse nützlich, etwa:

Πσ=EβΔc/σkrit

Πσ vergleicht eine chemisch erzeugte Spannungsskala mit einer kritischen Festigkeit. Werte über eins signalisieren mögliches Versagen.

36. Regimekarten

Werden mehrere Kennzahlen gegeneinander aufgetragen, entstehen Regimekarten. Sie zeigen beispielsweise Bereiche:

  • kinetischer Kontrolle,
  • Diffusionskontrolle,
  • thermischer Stabilität,
  • mechanischen Versagens,
  • gekoppelter Instabilität.

Regimekarten sind oft informativer als eine einzelne Simulation.



Abbildung 2. Multiphysik entsteht durch gegenseitige Quellterme und zustandsabhängige Materialparameter. Dimensionslose Kennzahlen zeigen, welche Prozesse in einem gegebenen Größen-, Zeit- und Belastungsbereich dominieren.

37. Numerische Diskretisierung

Kontinuumsgleichungen werden räumlich und zeitlich diskretisiert. Wichtige Verfahren sind:

  • Finite-Differenzen-Methode,
  • Finite-Volumen-Methode,
  • Finite-Elemente-Methode,
  • Spektralmethoden,
  • Partikelmethoden.

Die Wahl hängt von Geometrie, Erhaltungseigenschaften, Nichtlinearität und gewünschten lokalen Größen ab.

38. Netzunabhängigkeit

Eine numerische Lösung muss gegen Netzverfeinerung geprüft werden. Wird ein Riss, eine Grenzschicht oder eine Phasengrenze nicht räumlich aufgelöst, kann das Ergebnis scheinbar konvergieren und dennoch physikalisch falsch sein.

39. Zeitintegration

Explizite Verfahren sind einfach, aber oft durch Stabilitätsbedingungen begrenzt. Implizite Verfahren erlauben größere Zeitschritte, benötigen jedoch nichtlineare Gleichungslöser.

Bei stark getrennten Zeitskalen entstehen steife Systeme.

40. Verifikation

Verifikation beantwortet:

Lösen wir die gewählten Gleichungen numerisch richtig?

Werkzeuge sind analytische Speziallösungen, Herstellungsfälle, Konvergenzraten, Bilanzfehler und Vergleich unabhängiger Codes.

41. Validierung

Validierung beantwortet:

Beschreiben die Gleichungen die reale Zielgröße ausreichend?

Dafür werden unabhängige Experimente mit kontrollierten Randbedingungen und Messunsicherheiten benötigt.

42. Kalibrierung

Kalibrierung passt unbekannte Parameter an Daten an. Sie ist keine Validierung. Ein Modell kann Daten mit falscher Physik fitten, wenn genügend flexible Parameter vorhanden sind.

43. Parameteridentifizierbarkeit

Ein Parameter ist nur identifizierbar, wenn die Messdaten empfindlich und spezifisch auf ihn reagieren. Stark korrelierte Parameter können viele gleich gute Fits erzeugen.

Mehr Messgrößen und gezielt variierte Bedingungen verbessern die Identifizierbarkeit.

44. Lokale Sensitivität

Si=∂y/∂pi

Eine dimensionslose Sensitivität lautet:

Si*=(pi/y)∂y/∂pi

Sie zeigt die lokale Reaktion der Zielgröße y auf einen Parameter pi.

45. Globale Sensitivität

Globale Verfahren variieren alle Parameter über ihre zulässigen Bereiche. Sie erfassen Nichtlinearität, Wechselwirkungen und Regimewechsel.

Varianzbasierte Indizes, Screening-Verfahren und Stichprobenmethoden sind typische Ansätze.

46. Unsicherheitsquellen

UnsicherheitsartBeispielBehandlung
MessunsicherheitTemperatur, Strom, SchichtdickeKalibrierung und Fehlerfortpflanzung
ParameterstreuungKorngröße, DefektdichteWahrscheinlichkeitsverteilung
Modellformfehlende Reaktion oder falsche RandbedingungModellvergleich und Diskrepanzterm
Numerischer FehlerNetz und ZeitschrittKonvergenzanalyse
Extrapolationneue Temperatur oder ZusammensetzungGültigkeitsprüfung und konservative Intervalle

47. Lineare Unsicherheitsfortpflanzung

Für unabhängige Größen pi:

σy²≈Σi(∂y/∂pi)²σp_i²

Die Näherung gilt bei kleinen Unsicherheiten und annähernd linearem Verhalten.

48. Monte-Carlo-Unsicherheitsanalyse

Parameter werden aus ihren Verteilungen gezogen, das Modell wiederholt gerechnet und die Zielgröße statistisch ausgewertet. Das Verfahren ist allgemein, aber bei teuren Multiphysikmodellen rechenintensiv.

49. Bayes’sche Aktualisierung

Vorwissen und Messdaten werden kombiniert:

Posterior ∝ Likelihood × Prior

Das Ergebnis ist keine einzelne optimale Zahl, sondern eine Verteilung plausibler Parameter.

50. Surrogatmodelle

Ein Surrogat approximiert ein teures Modell durch eine schnell auswertbare Funktion. Beispiele sind Polynome, Gauß-Prozesse, neuronale Netze und reduzierte Ordnungsmodelle.

Ein Surrogat ist nur innerhalb des abgedeckten Parameterraums vertrauenswürdig. Extrapolation muss erkannt und gekennzeichnet werden.

51. Physikinformierte Modelle

Datenbasierte Modelle können Bilanzgleichungen, Symmetrien, Erhaltungssätze und Dimensionsanalyse als Struktur verwenden. Dadurch sinkt der Datenbedarf und die Vorhersage wird physikalisch konsistenter.

Physikinformiert bedeutet jedoch nicht automatisch validiert. Falsche Randbedingungen oder unvollständige Physik können weiterhin systematische Fehler erzeugen.

52. Digitaler Zwilling

Ein digitaler Zwilling verbindet ein laufendes physikalisches System mit einem aktualisierbaren Modell. Sensordaten korrigieren Zustände und Parameter, während das Modell nicht direkt messbare Größen und zukünftige Entwicklung schätzt.

Ein digitaler Zwilling benötigt:

  • ein beobachtbares Zustandsmodell,
  • zuverlässige Sensoren,
  • Unsicherheitsbewertung,
  • regelmäßige Zustandsaktualisierung,
  • klare Entscheidungskriterien.

53. Inverses Materialdesign

Vorwärtsmodelle berechnen Eigenschaften aus Struktur und Zusammensetzung. Inverses Design startet mit Zielgrößen und sucht passende Materialien oder Mikrostrukturen:

Struktur → Eigenschaft → Funktion
Ziel-Funktion → zulässige Strukturen

Die inverse Abbildung ist meist nicht eindeutig. Viele Strukturen können ähnliche Zielwerte besitzen.

54. Designraum

Ein Materialdesign kann Variablen enthalten wie:

  • Elementzusammensetzung,
  • Dotierung und Defektkonzentration,
  • Partikelgröße und -form,
  • Porosität und Tortuosität,
  • Schichtdicke,
  • Grenzflächenchemie,
  • Prozesstemperatur und Zeit.

Der Designraum wächst kombinatorisch. Physikalische Regeln und Herstellungsgrenzen reduzieren ihn.

55. Mehrzieloptimierung

Materialien werden selten für nur eine Größe optimiert. Eine Batterieelektrode soll gleichzeitig hohe Energie, Leistung, Lebensdauer und Sicherheit besitzen. Ein Katalysator soll aktiv, selektiv, stabil und ressourceneffizient sein.

Ein gewichtetes Ziel kann lauten:

J=w₁f₁+w₂f₂+…+wnfn

Gewichte drücken Entscheidungen aus und sind keine Naturkonstanten.

56. Pareto-Front

Ein Entwurf ist Pareto-optimal, wenn kein Ziel verbessert werden kann, ohne mindestens ein anderes zu verschlechtern. Die Pareto-Front macht Zielkonflikte sichtbar, bevor eine konkrete Gewichtung gewählt wird.

57. Nebenbedingungen

Optimierung muss harte Grenzen berücksichtigen:

  • thermodynamische Stabilität,
  • mechanische Festigkeit,
  • toxikologische und rechtliche Vorgaben,
  • Rohstoffverfügbarkeit,
  • Prozessfenster,
  • Kosten und Skalierbarkeit.

58. Aktives Lernen

Aktives Lernen wählt gezielt das nächste Experiment oder die nächste Simulation aus. Bevorzugt werden Punkte mit hohem erwarteten Informationsgewinn oder großer möglicher Verbesserung.

So wird nicht der gesamte Designraum gleichmäßig abgesucht.

59. Design of Experiments

Statistische Versuchsplanung variiert mehrere Faktoren systematisch. Sie trennt Haupteffekte und Wechselwirkungen und vermeidet die ineffiziente Änderung nur eines Parameters pro Versuch.

60. Lebensdauerprognose

Lebensdauer ist häufig das Ergebnis gekoppelter, zeitabhängiger Schädigung. Ein allgemeines Zustandsmodell kann lauten:

dX/dt=f(X,T,Last,Umgebung)

X kann Risslänge, Grenzschichtdicke, Defektpopulation oder aktive Fläche darstellen.

61. Arrhenius-Beschleunigung

Für einen einzelnen thermisch aktivierten Mechanismus:

AF=exp[Ea/R(1/Tuse−1/Ttest)]

Der Beschleunigungsfaktor ist nur gültig, wenn bei Test- und Einsatztemperatur derselbe Mechanismus dominiert.

62. Mechanismenwechsel

Zu aggressive Prüfbedingungen können neue Reaktionen, Phasen oder Transportregime erzeugen. Dann ist die beschleunigte Alterung nicht mehr repräsentativ.

Eine lineare Extrapolation über einen Mechanismenwechsel ist unzulässig.

63. Zustands- und Gesundheitsgrößen

Für technische Systeme werden interne Größen oft zu Zustandsindikatoren verdichtet:

  • State of Charge,
  • State of Health,
  • Restlebensdauer,
  • Schädigungsindex,
  • Sicherheitsmarge.

Diese Größen sind modellabhängig und müssen mit beobachtbaren Messsignalen verbunden sein.

64. Nachhaltiges Materialdesign

Eine hohe Funktionsleistung kann durch kritische Rohstoffe, energieintensive Herstellung oder schwieriges Recycling erkauft werden. Nachhaltiges Design bewertet daher den gesamten Lebensweg.

65. Lebenszyklusanalyse

Eine Lebenszyklusanalyse umfasst typischerweise:

  1. Rohstoffgewinnung,
  2. Material- und Bauteilherstellung,
  3. Nutzungsphase,
  4. Wartung und Ersatz,
  5. Recycling oder Entsorgung.

Systemgrenze und funktionelle Einheit bestimmen das Ergebnis wesentlich.

66. Funktionelle Einheit

Materialien müssen auf dieselbe Funktion bezogen werden. Für einen Speicher kann die geeignete Einheit eine über die Lebensdauer abgegebene Kilowattstunde sein, nicht nur ein Kilogramm Material.

67. Zirkularität und Design for Recycling

Recyclingfähigkeit beginnt bei der Material- und Bauteilarchitektur. Trennbare Schichten, bekannte Additive, reversible Verbindungen und geringe Stoffvielfalt erleichtern Rückgewinnung.

Ein extrem komplexer Hochleistungsverbund kann am Lebensende ungünstig sein.

68. Kritische Rohstoffe

Rohstoffkritikalität umfasst nicht nur geologische Häufigkeit, sondern auch:

  • geografische Konzentration,
  • Koproduktion mit anderen Metallen,
  • politische und soziale Risiken,
  • Substituierbarkeit,
  • Recyclingquote,
  • Nachfragewachstum.

69. Herstellbarkeit

Ein Material kann im Labor hervorragend funktionieren und dennoch industriell scheitern. Wichtige Fragen sind:

  • Sind Reinheit und Toleranzen realistisch?
  • Ist der Prozess schnell und skalierbar?
  • Entstehen toxische oder flüchtige Stoffe?
  • Wie empfindlich ist die Leistung gegenüber Prozessschwankungen?
  • Lässt sich Qualität zerstörungsfrei prüfen?

70. Robust Design

Ein robuster Entwurf behält seine Funktion trotz unvermeidbarer Streuung. Optimiert wird nicht nur der Mittelwert, sondern auch die Empfindlichkeit:

Ziel=hohe Leistung−λ·Streuung

Ein geringfügig weniger leistungsfähiges, aber tolerantes Material kann technisch überlegen sein.

71. Fallstudie Batterie

Eine Batterie koppelt Elektrochemie, Festkörperdiffusion, Elektrolyttransport, Wärme und Mechanik. Hohe Ladeleistung erhöht Konzentrationsgradienten, Polarisation und Wärme. Diese steigern mechanische Spannungen und Nebenreaktionen, wodurch Widerstand und Wärme weiter zunehmen.

Das optimale Partikeldesign ist daher kein reines Diffusionsproblem.

72. Fallstudie Photovoltaik

Bei einer Solarzelle koppeln Lichtausbreitung, Ladungstransport, Rekombination und Wärme. Dünnere Schichten verringern Transportwege, können aber Absorption reduzieren. Stärkere Texturierung erhöht den optischen Weg, kann jedoch Oberflächenrekombination oder mechanische Rauheit fördern.

73. Fallstudie Katalysator

Ein poröser Katalysator koppelt Adsorption, Reaktion, Stofftransport und Wärme. Sehr aktive Zentren können nutzlos sein, wenn Reaktanden die Poren nicht erreichen oder ein Hotspot die Selektivität verändert.

74. Fallstudie Korrosion

Korrosion verbindet Potentialfelder, Ionentransport, Filmwachstum, lokale Chemie und Mechanik. Ein Passivfilm kann thermodynamisch stabil sein und dennoch durch Spannung oder Chloridanreicherung lokal versagen.

75. Gesamtübersicht Physikalische Chemie IV

BlockZentrale FragenTypische Modelle
Kristallstruktur und DefekteWie ist Materie geordnet und wie weicht sie von Idealität ab?Gitter, reziproker Raum, Defektgleichgewichte
Elektronische EigenschaftenWie entstehen Leitfähigkeit, Halbleitung, Magnetismus und Optik?Bänder, Zustandsdichte, Trägerstatistik
Oberflächen und NanomaterialienWie verändern Größe und Grenzflächen die freie Energie?Kapillarität, Adsorption, Quantenbegrenzung
Reaktionen und DegradationWie laufen Oberflächenreaktionen und Materialabbau ab?Mikrokinetik, Passivierung, Transport
Energieumwandlung und SpeicherungWie werden Licht und chemische Potentiale technisch genutzt?Solarzellen- und Batteriemodelle
Multiphysik und DesignWie werden Mechanismen zu robusten Bauteilen?Skalenkopplung, Optimierung, Unsicherheit


Abbildung 3. Materialentwicklung ist ein iterativer Regelkreis. Modelle wählen Experimente, Messungen aktualisieren Modelle, und technische Entscheidungen berücksichtigen Leistung, Streuung, Lebensdauer, Fertigung und Umweltwirkung.

76. Übungsaufgaben

Aufgabe 1: Diffusionszeit

Eine Spezies diffundiert über L=50,0 µm mit D=2,00×10−14 m² s−1. Berechne die charakteristische Zeit t=L²/D.

Aufgabe 2: Fourier-Zahl

Für ein Bauteil mit L=5,00 mm, thermischer Diffusivität α=1,20×10−6 m² s−1 und t=600 s: Berechne Fo=αt/L² und bewerte den Temperaturausgleich.

Aufgabe 3: Thermische Spannung

Ein vollständig behinderter Werkstoff besitzt E=150 GPa, αT=12,0×10−6 K−1, ν=0,300 und ΔT=80,0 K. Verwende:

σth=EαTΔT/(1−ν)

Aufgabe 4: Damköhler-Zahl

Für k=0,800 s−1, L=2,00 mm und D=5,00×10−6 m² s−1: Berechne Da=kL²/D.

Aufgabe 5: Effektive Leitfähigkeit

Ein zweiphasiger Verbund enthält φ=0,350 einer Phase mit K₁=100 S m−1 und 0,650 einer Phase mit K₂=1,00 S m−1. Berechne Parallel- und Seriengrenze.

Aufgabe 6: Unsicherheitsfortpflanzung

Für y=ab/c gelten a=2,00±0,10, b=5,00±0,20 und c=4,00±0,08. Die Größen seien unabhängig. Berechne y und σy.

Aufgabe 7: Arrhenius-Beschleunigung

Ein Alterungsmechanismus besitzt Ea=0,750 eV. Berechne den Beschleunigungsfaktor zwischen 25 °C Einsatz- und 80 °C Prüftemperatur.

Aufgabe 8: Gewichtete Mehrzielbewertung

Die Gewichte für Leistung, Lebensdauer und Nachhaltigkeit lauten 0,40, 0,35 und 0,25. Entwurf A besitzt normierte Werte 0,85, 0,65 und 0,90; Entwurf B 0,75, 0,88 und 0,70. Berechne beide gewichteten Scores.

Aufgabe 9: Adiabatische Temperaturerhöhung

Ein Bauteil setzt kurzfristig 12,0 kJ Wärme frei. Masse m=0,800 kg, Wärmekapazität cp=900 J kg−1 K−1. Berechne ΔT=Q/(mcp).

Aufgabe 10: Lebenszyklusbezogene Emission

Die Herstellung eines Speichers verursacht 75,0 kg CO₂-Äquivalente pro kWh Nennkapazität. Er erreicht 2500 Zyklen bei 80 Prozent Entladetiefe und 92 Prozent Energie-Wirkungsgrad. Berechne die herstellungsbezogene Emission pro abgegebener kWh.

77. Vollständige Lösungen

Lösung 1
t=(50,0×10−6)²/(2,00×10−14)
t=1,25×105 s=34,72 h

Die lange Zeit zeigt, warum kleine Diffusionswege für hohe Leistung entscheidend sind.

Lösung 2
Fo=(1,20×10−6)(600)/(5,00×10−3
Fo=28,8

Fo ist deutlich größer als eins. Über die betrachtete Länge hat thermische Diffusion genügend Zeit für einen weitgehenden Temperaturausgleich.

Lösung 3
σth=(150×109)(12,0×10−6)(80,0)/(1−0,300)
σth=2,057×108 Pa=205,7 MPa

Eine reale Struktur kann durch plastische Verformung, Kriechen oder Rissbildung entspannen.

Lösung 4
Da=(0,800)(2,00×10−3)²/(5,00×10−6)
Da=0,640

Reaktion und Diffusion liegen in ähnlicher Größenordnung. Eine reine Grenzfallnäherung ist daher nicht sicher.

Lösung 5

Parallelgrenze:

Kparallel=(0,350)(100)+(0,650)(1,00)=35,65 S m−1

Seriengrenze:

Kseriell=1/[0,350/100+0,650/1,00]
Kseriell=1,530 S m−1

Die enorme Spanne zeigt die Bedeutung der Phasenkonnektivität.

Lösung 6
y=(2,00)(5,00)/4,00=2,50

Relative Unsicherheit:

σy/y=√[(0,10/2,00)²+(0,20/5,00)²+(0,08/4,00)²]
σy/y=0,06708
σy=0,1677

Damit lautet das Ergebnis näherungsweise y=2,50±0,17.

Lösung 7

0,750 eV pro Teilchen entspricht 72,36 kJ mol−1.

AF=exp[(72364/8,314)(1/298,15−1/353,15)]
AF=94,28

Die Prüfung läuft nach dem Einmechanismus-Modell rund 94-mal schneller. Das Ergebnis ist unbrauchbar, falls bei 80 °C ein anderer Mechanismus einsetzt.

Lösung 8
SA=0,40(0,85)+0,35(0,65)+0,25(0,90)=0,7925
SB=0,40(0,75)+0,35(0,88)+0,25(0,70)=0,7830

Mit diesen Gewichten liegt A knapp vorn. Andere Gewichte können die Entscheidung umkehren.

Lösung 9
ΔT=12000/[(0,800)(900)]
ΔT=16,67 K

Die Rechnung vernachlässigt Wärmeabfuhr und Temperaturabhängigkeit der Wärmekapazität.

Lösung 10

Über die Lebensdauer abgegebene Energie pro kWh Nennkapazität:

Elife=(2500)(0,80)(0,92)=1840 kWh

Herstellungsbezogene Emission:

75,0 kg/1840 kWh=0,04076 kg kWh−1
40,76 g CO₂e kWh−1

Recycling, Strommix, Wartung und Kapazitätsabnahme sind in dieser vereinfachten Rechnung nicht enthalten.

78. Zusammenfassung

  • Materialfunktion entsteht aus Mechanismen über viele Längen- und Zeitskalen.
  • Elektronenstruktur liefert Energien, Barrieren und elektronische Eigenschaften.
  • Atomistische Methoden beschreiben Dynamik, Diffusion und lokale Struktur.
  • Mesoskopische Modelle bilden Phasen, Körner, Poren und Grenzflächen ab.
  • Kontinuumsmodelle lösen Bilanzgleichungen auf Bauteilskala.
  • Homogenisierung ersetzt Mikrostruktur durch effektive Tensoren.
  • Ein RVE muss für die jeweilige Zielgröße repräsentativ sein.
  • Skalenkopplung kann hierarchisch oder simultan erfolgen.
  • Stoff-, Ladungs-, Wärme- und Impulsbilanzen besitzen eine gemeinsame Struktur.
  • Chemische und thermische Dehnung koppeln Transport und Mechanik.
  • Reaktionswärme kann positive thermische Rückkopplung erzeugen.
  • Phasenfeldmodelle beschreiben diffuse Grenzflächen und Morphologieentwicklung.
  • Dimensionslose Kennzahlen identifizieren dominante Regime.
  • Fo charakterisiert thermischen Ausgleich, Pe Konvektion, Da Reaktion gegen Diffusion.
  • Numerische Ergebnisse benötigen Netz- und Zeitschrittkonvergenz.
  • Verifikation prüft die Numerik, Validierung die Realitätsnähe.
  • Kalibrierung ersetzt keine unabhängige Validierung.
  • Sensitivität zeigt, welche Parameter die Zielgröße kontrollieren.
  • Unsicherheit umfasst Messung, Parameter, Modellform und Numerik.
  • Surrogatmodelle sind schnell, aber nur im abgedeckten Bereich vertrauenswürdig.
  • Digitale Zwillinge verbinden Sensoren, Zustandsmodelle und Prognose.
  • Inverses Design sucht Strukturen für vorgegebene Funktionen.
  • Mehrzieloptimierung macht Zielkonflikte sichtbar.
  • Pareto-Fronten trennen physikalische Möglichkeiten von subjektiver Gewichtung.
  • Aktives Lernen wählt informative Experimente.
  • Beschleunigte Alterung ist nur bei gleichem Mechanismus extrapolierbar.
  • Nachhaltigkeit benötigt eine definierte funktionelle Einheit.
  • Robustes Design minimiert Empfindlichkeit gegenüber Streuung.
  • Herstellbarkeit und Recycling müssen früh in den Designraum eingehen.
  • Physikalische Chemie verbindet atomistische Erklärung und technische Vorhersage.

79. Häufig gestellte Fragen

Warum kann nicht einfach das detaillierteste Modell verwendet werden?

Das detaillierteste Modell ist oft zu teuer, enthält unbekannte Parameter und liefert nicht automatisch bessere Bauteilvorhersagen. Entscheidend ist die zur Zielgröße passende Auflösung.

Was ist der Unterschied zwischen Multiskalen- und Multiphysikmodell?

Multiskala koppelt verschiedene Größen- oder Zeitebenen. Multiphysik koppelt verschiedene Felder wie Wärme, Transport und Mechanik. Ein Modell kann beides zugleich sein.

Wann ist Homogenisierung unzulässig?

Wenn die Mikrostruktur nicht klein gegenüber dem Bauteil ist, lokale Hotspots entscheidend sind oder kein repräsentatives Volumen existiert.

Ist ein Modell mit guter Datenanpassung automatisch richtig?

Nein. Flexible Modelle können falsche Mechanismen kompensieren. Unabhängige Validierung und physikalische Plausibilität bleiben notwendig.

Warum sind dimensionslose Zahlen hilfreich?

Sie vergleichen konkurrierende Prozesse unabhängig von Einheiten und zeigen Grenzfälle, Skalierungsgesetze und dominante Regime.

Kann maschinelles Lernen physikalische Modelle ersetzen?

Es kann sie beschleunigen oder ergänzen. Außerhalb der Datenbasis und bei neuen Mechanismen benötigt es jedoch physikalische Struktur und Unsicherheitskontrolle.

Was macht einen digitalen Zwilling aus?

Nicht nur eine Simulation, sondern die regelmäßige Verbindung mit Messdaten, Zustandsaktualisierung und konkreter Entscheidungsfunktion.

Warum ist eine Pareto-Front besser als ein einziger Score?

Sie zeigt zunächst die objektiven Zielkonflikte. Ein Score versteckt diese durch vorab gewählte Gewichte.

Warum kann ein beschleunigter Test irreführen?

Hohe Temperatur, Spannung oder Feuchte können einen anderen Schädigungsmechanismus aktivieren als im realen Einsatz.

Was ist die zentrale Botschaft des Kurses?

Makroskopische Materialeigenschaften entstehen aus Struktur, Energie, Transport und Nichtgleichgewicht. Verlässliches Design verlangt die kontrollierte Verbindung dieser Ebenen.

80. Abschluss von Physikalische Chemie IV

Mit Teil 20 ist der vierte Kursteil abgeschlossen. Der Weg führte von idealen Kristallgittern über Defekte, Elektronen, Phononen, Oberflächen und Nanostrukturen bis zu Katalysatoren, Korrosionssystemen, Solarzellen und Batterien.

Die verbindende Idee lautet: Ein Material ist kein statischer Stoffwertträger. Es ist ein dynamisches System aus Zuständen, Flüssen, Grenzflächen und Rückkopplungen. Seine technische Funktion entsteht erst unter realen Randbedingungen, und seine Lebensdauer ist Teil dieser Funktion.

Literatur und weiterführende Quellen

  1. J. Newman & K. E. Thomas-Alyea: Electrochemical Systems. Wiley.
  2. W. W. Mullins und moderne Fachliteratur zu Phasenfeld- und Mikrostrukturmodellen.
  3. T. J. R. Hughes: The Finite Element Method sowie weiterführende Literatur zur numerischen Kontinuumsmechanik.
  4. G. E. P. Box, W. G. Hunter & J. S. Hunter: Statistics for Experimenters. Wiley.
  5. Fachliteratur zu Multiskalenmodellierung, Unsicherheitsquantifizierung, Materialinformatik und Lebenszyklusanalyse.

Didaktischer Hinweis: Kontinuumsfelder, effektive Eigenschaften, dimensionslose Kennzahlen, lineare Unsicherheitsfortpflanzung, Arrhenius-Beschleunigung und gewichtete Zielgrößen sind kontrollierte Modelle. Reale Materialsysteme können nichtstationäre Mikrostrukturen, Mechanismenwechsel, korrelierte Parameter, diskrete Defekte, stochastische Risse und gesellschaftlich definierte Zielkonflikte besitzen.

Bild- und Nutzungsnachweis: Das Titelbild und die Abbildungen 1–3 wurden eigens für diesen Blogartikel mit selbst erstelltem Programmcode erzeugt. Es wurden keine Grafiken aus Lehrbüchern oder anderen Publikationen übernommen. Die Abbildungen sind für die Verwendung in diesem monetarisierbaren Blog vorgesehen.

Moderne Batterie- und Festelektrolytmaterialien – Interkalation, Phasengrenzen und Degradation




Titelbild: Eigene Illustration einer elektrochemischen Zelle sowie zentraler Beziehungen für Zellspannung, theoretische Kapazität, Ionentransport und Überspannungen.

Eine wiederaufladbare Batterie speichert elektrische Energie als chemische freie Energie. Beim Laden wird die Zellreaktion gegen ihre spontane Richtung betrieben; beim Entladen läuft sie kontrolliert zurück. Moderne Batteriematerialien müssen dabei weit mehr leisten als eine hohe theoretische Kapazität: Ionen und Elektronen müssen transportiert, Phasenänderungen reversibel durchlaufen, Grenzflächen stabilisiert und mechanische Spannungen über viele Zyklen aufgenommen werden.

Die Zellleistung entsteht aus einem gekoppelten System:

  • thermodynamische Stabilität und Zellspannung,
  • Einlagerung oder Umwandlung in den Elektroden,
  • Diffusion in Festkörpern und Poren,
  • Ladungsübertragung an Grenzflächen,
  • Ionenleitung im Elektrolyten,
  • elektronische Leitung im Elektrodennetzwerk,
  • Wärmeentwicklung und Wärmeabfuhr,
  • chemische, mechanische und elektrochemische Degradation.
Einordnung: Frühere Teile behandelten chemisches Potential, Phasengleichgewichte, Festkörperdefekte, Diffusion, Elektrodenkinetik und Materialtransport. Teil 19 verbindet diese Grundlagen zur Physikalischen Chemie moderner Batteriesysteme. Teil 20 schließt den Kurs mit Multiphysik, Materialdesign und der Kopplung atomistischer Mechanismen an technische Bauteile ab.

Lernziele

Nach diesem Beitrag kannst du:
  • Zellspannung aus Gibbs-Energie und chemischen Potentialen ableiten,
  • Interkalation, Legierung und Konversionsreaktion unterscheiden,
  • Spannungsprofile aus Festlösungs- und Zweiphasengleichgewichten erklären,
  • theoretische Kapazität und spezifische Energie berechnen,
  • C-Rate, Diffusionszeit und Leistungsbegrenzung einordnen,
  • Ionenleitfähigkeit und Überführungszahl interpretieren,
  • Flüssig-, Polymer-, Sulfid- und Oxidfestelektrolyte vergleichen,
  • Butler–Volmer-, ohmische und Konzentrationsüberspannung unterscheiden,
  • SEI und CEI als kinetisch stabilisierende Grenzphasen erklären,
  • Graphit, Silizium, Lithium-Metall und wichtige Kathodenklassen vergleichen,
  • Dendriten, Rissbildung und Kontaktverlust physikalisch einordnen,
  • Degradationsmodi, Sicherheit und diagnostische Methoden bewerten.

1. Elektrochemische Zellreaktion

Eine vollständige Zelle überträgt Ionen durch den Elektrolyten und Elektronen durch den äußeren Stromkreis. Für eine allgemeine Zellreaktion:

ΣiνiAi=0

ist die molare Gibbs-Energieänderung:

ΔrG=Σiνiμi

Die Zellspannung ist daher unmittelbar mit chemischen Potentialen verknüpft.

2. Reversible Zellspannung

Erev=−ΔrG/(nF)

n ist die Zahl übertragener Elektronen pro stöchiometrischem Reaktionsumsatz. Eine negative Gibbs-Energie der Entladung erzeugt eine positive Zellspannung.

3. Nernst-Gleichung der Gesamtzelle

E=E°−RT/(nF)lnQ

Q ist der Reaktionsquotient. In Batterien hängen Aktivitäten häufig nicht nur von gelösten Spezies, sondern von der Zusammensetzung fester Einlagerungsphasen ab.

4. Elektrochemisches Potential

Für eine geladene Spezies i:

μ̃ii+zi

Im Gleichgewicht ist das elektrochemische Potential der mobilen Spezies über jede leitfähige Phase konstant. Spannungsabfälle entstehen dort, wo elektrochemische Potentiale zwischen Elektroden und Kontakten unterschiedlich sind.

5. Lithium als Referenz

Viele Batteriereaktionen werden relativ zum Paar Li+/Li angegeben. Für eine Einlagerungselektrode:

Host+xLi++xe ⇌ LixHost

entspricht die Elektrodenspannung gegenüber Lithium-Metall dem Unterschied des Lithium-Chemischen Potentials im Host und im Metall.

6. Spannung und Zusammensetzung

Für eine Einlagerungsverbindung LixHost gilt näherungsweise:

E(x)=−[μLi,Host(x)−μLi,Metall]/F

Die Steigung eines Spannungsprofils spiegelt daher direkt die Zusammensetzungsabhängigkeit des chemischen Potentials wider.

7. Interkalation

Bei einer Interkalation werden Ionen in ein bestehendes Wirtsgitter eingelagert, ohne dass dessen Grundgerüst vollständig zerstört wird. Geeignete Wirte besitzen:

  • freie Zwischengitter- oder Schichtplätze,
  • niedrige Migrationsbarrieren,
  • redoxaktive Gitterbausteine,
  • begrenzte Volumenänderung,
  • ausreichende elektronische Leitfähigkeit oder Leitadditive.

8. Festlösung

In einer einphasigen Festlösung ändert sich x kontinuierlich. Ein ideales Mischungsmodell lautet:

G(x)=xG1+(1−x)G0+RT[xlnx+(1−x)ln(1−x)]

Reale Systeme enthalten zusätzlich Wechselwirkungs- und elastische Beiträge. Die Spannung verläuft bei einer Festlösung gewöhnlich geneigt.

9. Reguläre Lösung

Gmix=RT[xlnx+(1−x)ln(1−x)]+Ωx(1−x)

Ein positiver Wechselwirkungsparameter Ω kann eine Mischungslücke erzeugen. Unterhalb einer kritischen Temperatur trennt sich das Material in zwei Zusammensetzungen.

10. Zweiphasengebiet

Koexistieren zwei Phasen α und β, werden ihre Gleichgewichtszusammensetzungen durch eine gemeinsame Tangente an G(x) bestimmt. Im idealen Gleichgewicht bleibt das chemische Potential während der Umwandlung konstant.

Dadurch entsteht ein Spannungsplateau.

11. Phasengrenzen

Eine Zweiphasenreaktion benötigt Keimbildung und Bewegung von Phasengrenzen. Grenzflächenenergie, elastische Kohärenzspannung und Partikelgröße können das Gleichgewicht verschieben und Hysterese erzeugen.

12. Spannungs-Hysterese

Die Ladespannung liegt oft über der Entladespannung. Beiträge stammen aus:

  • Ladungsübertragungsüberspannung,
  • ohmischem Widerstand,
  • Diffusionsgradienten,
  • Keimbildungsbarrieren,
  • elastischer Energie und Reibung bewegter Phasengrenzen,
  • irreversiblen Nebenreaktionen.

13. Interkalation, Legierung und Konversion

ReaktionstypStrukturprinziptypische Stärketypische Herausforderung
InterkalationIoneneinbau in erhaltenes Wirtsskelettgute Reversibilitätbegrenzte Zahl verfügbarer Plätze
LegierungBildung Li-reicher Metall- oder Halbmetallphasensehr hohe Kapazitätgroße Volumenänderung
KonversionZerfall in neue Phasen, etwa Metall und Li₂Ohohe Elektronenzahlgroße Hysterese und langsame Rekonstruktion

14. Theoretische spezifische Kapazität

Für n übertragene Elektronen pro Formeleinheit mit molarer Masse M:

Qth=nF/(3,6M)

Ist M in g mol−1 eingesetzt, ergibt sich Qth in mAh g−1.

15. Praktische Kapazität

Die praktische Kapazität ist kleiner als die theoretische, weil:

  • nicht alle Plätze zugänglich sind,
  • Spannungsgrenzen Nebenreaktionen vermeiden müssen,
  • Diffusion und Leitfähigkeit die Ausnutzung begrenzen,
  • Leitadditive, Binder und Stromableiter keine aktive Kapazität liefern.

16. Spezifische Energie

Für eine reale Spannungs-Kapazitäts-Kurve:

w=∫V(Q)dQ/m

Bei angenähert konstanter mittlerer Spannung:

w≈QspezVmittel

Elektrodenwerte dürfen nicht direkt mit Zellwerten verwechselt werden, weil Elektrolyt, Separator, Gehäuse und inaktive Komponenten Masse hinzufügen.

17. Energiedichte und Leistungsdichte

Energiedichte beschreibt gespeicherte Energie pro Masse oder Volumen. Leistungsdichte beschreibt Energieabgabe pro Zeit. Hohe Energie verlangt oft dicke, hoch beladene Elektroden; hohe Leistung bevorzugt kurze Wege und große aktive Fläche.

18. C-Rate

Eine Rate von 1C entlädt die Nennkapazität theoretisch in einer Stunde. Für eine Kapazität CN in Ah:

I=CRateCN

Bei 0,5C beträgt die ideale Entladezeit zwei Stunden, bei 2C eine halbe Stunde.

19. Coulomb-Wirkungsgrad

ηC=QEntladung/QLadung

Ein geringfügig unter eins liegender Coulomb-Wirkungsgrad akkumuliert über viele Zyklen zu erheblichem Kapazitätsverlust.

20. Energie-Wirkungsgrad

ηE=EnergieEntladung/EnergieLadung

Spannungshysterese und Widerstände senken den Energie-Wirkungsgrad auch dann, wenn der Coulomb-Wirkungsgrad hoch ist.



Abbildung 1. Die Zellspannung folgt aus einer Differenz chemischer Potentiale. Eine kontinuierliche Festlösung liefert eine geneigte Spannungskurve; ein Zweiphasengebiet ein thermodynamisches Plateau.

21. Festkörperdiffusion

Die Einlagerung erfordert Ionentransport vom Partikelrand in das aktive Material. Für eine Konzentration c:

J=−D∇c
∂c/∂t=D∇²c

Die charakteristische Diffusionszeit skaliert näherungsweise wie:

tdiff≈L²/D

Eine Halbierung der Partikelgröße verkürzt diese Zeit um etwa den Faktor vier.

22. Chemische Diffusivität

In konzentrierten Festlösungen wird der Fluss durch den Gradienten des chemischen Potentials angetrieben. Eine chemische Diffusivität kann geschrieben werden als:

D̃=D*·∂ln a/∂ln c

Der thermodynamische Faktor kann D̃ stark von der Tracer- oder Sprungdiffusivität unterscheiden. Nahe Phasengrenzen kann die einfache Diffusionsbeschreibung versagen.

23. Anisotrope Diffusion

Schichtoxide besitzen oft schnelle zweidimensionale Diffusionswege, Olivinstrukturen eher eindimensionale Kanäle. Defekte, Antisite-Besetzung und Korngrenzen können Wege blockieren oder verbinden.

Partikelform und Kristallorientierung beeinflussen deshalb die Leistungsfähigkeit.

24. Galvanostatische intermittierende Titration

Bei GITT wird ein kleiner Strompuls angelegt und anschließend entspannt. Aus transienter Spannungsänderung, Gleichgewichtsschritt und Geometrie kann eine scheinbare chemische Diffusivität abgeschätzt werden.

Die Auswertung setzt unter anderem eindimensionale Diffusion, kleine Zusammensetzungsänderung und vernachlässigbare Nebenreaktionen voraus.

25. Potentiostatische intermittierende Titration

PITT verwendet kleine Potentialsprünge und misst den abklingenden Strom. Sie ist besonders nützlich zur Untersuchung von Phasenumwandlungen und Diffusionsregimen, benötigt aber ebenfalls kontrollierte Geometrie und Grenzflächenkinetik.

26. Poröse Elektroden

Technische Elektroden enthalten aktive Partikel, Leitadditiv, Binder und Elektrolytporen. Elektronen fließen durch die Feststoffmatrix, Ionen durch die Porenphase. Die lokale Reaktionsrate verteilt sich über die Elektrodendicke.

Eine hohe Materialkapazität ist wertlos, wenn das elektronische oder ionische Netzwerk nicht durchgängig ist.

27. Nernst–Planck-Transport

Für eine geladene Spezies i in verdünnter Form:

Ji=−Di∇ci−ziuiFci∇φ+civ

Die Terme beschreiben Diffusion, Migration im elektrischen Feld und Konvektion. In konzentrierten Batterieelektrolyten sind Aktivitäten, gekoppelte Flüsse und Stefan–Maxwell-Beschreibungen angemessener.

28. Ladungsübertragung

Die Reaktion an einer Elektrolyt-Elektroden-Grenze wird häufig durch Butler–Volmer beschrieben:

i=i₀{exp[αaFη/(RT)]−exp[−αcFη/(RT)]}

i₀ hängt von Konzentrationen, Bedeckung, Oberflächenstruktur und Temperatur ab.

29. Aktivierungsüberspannung

Für einen symmetrischen Ein-Elektronen-Prozess mit αac=0,5:

ηakt=2RT/F·asinh[i/(2i₀)]

Eine große Austauschstromdichte reduziert die notwendige Aktivierungsüberspannung.

30. Ohmische Überspannung

ηohm=iR

Beiträge stammen aus Elektrolyt, Stromableitern, Kontaktstellen, Leitadditivnetzwerk und Grenzschichten. Flächenspezifische Widerstände werden häufig in Ω cm² angegeben.

31. Konzentrationsüberspannung

Bei hoher Rate entstehen Konzentrationsgradienten in Elektrolyt und Festkörper. Das lokale Gleichgewichtspotential verschiebt sich, und aktive Bereiche können verarmen oder übersättigt werden.

Die gemessene Zellspannung kann daher deutlich von der Gleichgewichtsspannung abweichen.

32. Polarisation der Gesamtzelle

VZelle=Eeq−ηakt−ηohm−ηdiff

Beim Laden kehren sich die Vorzeichen entsprechend um. Hysterese enthält reversible und irreversible Beiträge.

33. Flüssige Elektrolyte

Konventionelle Lithium-Ionen-Zellen verwenden Lithiumsalze in organischen Carbonaten oder verwandten Lösungsmitteln. Vorteile sind:

  • hohe Raumtemperaturleitfähigkeit,
  • gute Benetzung poröser Elektroden,
  • geringer Kontaktwiderstand.

Nachteile sind Brennbarkeit, flüchtige Komponenten und begrenzte elektrochemische Stabilität.

34. Elektrochemisches Stabilitätsfenster

Das thermodynamische Fenster liegt zwischen Reduktions- und Oxidationsgrenze des Elektrolyten. Reale Zellen können außerhalb dieses Fensters arbeiten, wenn kinetisch passivierende Grenzschichten entstehen.

Eine hohe gemessene Zersetzungsspannung allein beweist kein breites thermodynamisches Fenster.

35. Ionenleitfähigkeit

σ=Σinizi²q²Di/(kBT)

Die Nernst–Einstein-Beziehung gilt ideal für unabhängige Ladungsträger. Korrelationen, Ionenpaare und kollektive Bewegung können die reale Leitfähigkeit verändern.

36. Arrhenius- und VTF-Verhalten

In kristallinen Festelektrolyten ist häufig:

σT=σ₀exp(−Ea/kBT)

Bei polymeren oder glasigen Elektrolyten koppelt die Ionenbewegung oft an Segment- oder Strukturrelaxationen. Dann wird häufig eine Vogel–Tammann–Fulcher-Form verwendet.

37. Überführungszahl

t+=iKation/igesamt

In Flüssigelektrolyten tragen Anionen häufig einen erheblichen Stromanteil. Ein kleines t+ fördert Konzentrationspolarisation. Ein idealer Lithium-Festelektrolyt nähert sich tLi+=1.

38. Polymer-Festelektrolyte

Polymere solvatisieren Lithiumionen über polare Gruppen. Vorteile sind Flexibilität und guter Grenzflächenkontakt. Herausforderungen sind:

  • begrenzte Raumtemperaturleitfähigkeit,
  • Kristallisation des Polymers,
  • kleine Kationenüberführungszahl bei mobilen Anionen,
  • mechanische und oxidative Stabilität.

39. Oxidische Festelektrolyte

Granate, NASICON-artige Phasen und andere Oxide können hohe chemische und mechanische Stabilität besitzen. Typische Schwierigkeiten sind:

  • spröde Keramik,
  • hohe Sintertemperatur,
  • großer Grenzflächenwiderstand,
  • Poren und schlechte Festkörperbenetzung.

40. Sulfidische Festelektrolyte

Sulfide besitzen weiche, polarisierbare Gitter und können sehr hohe Lithiumleitfähigkeiten erreichen. Sie lassen sich teilweise kalt verdichten und kontaktieren Elektroden gut.

Empfindlichkeit gegenüber Feuchte, mögliche H₂S-Bildung und begrenzte elektrochemische Stabilität erfordern kontrollierte Verarbeitung und Schutzschichten.

41. Halogenid-Festelektrolyte

Halogenidische Leiter können gute Hochvoltstabilität und günstige Kathodenkontakte bieten. Reduktive Stabilität gegenüber Lithium-Metall und Rohstoff- beziehungsweise Prozessfragen müssen materialspezifisch gelöst werden.

42. Raumladung an Festkörpergrenzen

Treffen zwei Festkörper mit unterschiedlichen Defektchemischen Potentialen aufeinander, verteilen sich mobile Ionen und elektronische Defekte neu. Eine Raumladungszone kann entstehen.

Sie kann die lokale Ionenleitfähigkeit erhöhen oder senken, ist aber nur ein möglicher Beitrag neben chemischer Reaktion, Poren und mechanischem Kontakt.

43. Flächenspezifischer Elektrolytwiderstand

RA=L/σ

Dünne Elektrolyte senken den Widerstand, erhöhen aber Anforderungen an Defektfreiheit, mechanische Stabilität und Kurzschlussfestigkeit.

44. Grenzflächenreaktionsschichten

Ein Festelektrolyt kann an der Elektrode thermodynamisch instabil sein. Die Zersetzung kann:

  • eine ionisch leitende, elektronisch isolierende Passivschicht bilden,
  • eine elektronisch leitende und fortschreitend wachsende Mischphase erzeugen,
  • mechanisch instabile oder poröse Produkte bilden.

Nicht jede Zersetzung ist daher gleich schädlich.

45. Solid Electrolyte Interphase

Die SEI entsteht meist durch Elektrolytreduktion an der negativen Elektrode. Eine gute SEI ist:

  • elektronisch isolierend,
  • lithiumionenleitend,
  • chemisch stabil,
  • mechanisch anpassungsfähig,
  • möglichst dünn und gleichmäßig.

Ihre Bildung verbraucht zyklierbares Lithium und erklärt einen Teil des irreversiblen Erstzyklusverlusts.

46. Cathode Electrolyte Interphase

Die CEI entsteht durch oxidative Reaktionen an der positiven Elektrode. Sie kann weitere Elektrolytoxidation hemmen, aber zugleich Widerstand, Gasentwicklung und Übergangsmetallauflösung beeinflussen.

47. Lithium-Metall-Anode

Lithium-Metall besitzt eine sehr hohe theoretische spezifische Kapazität und das niedrigste elektrochemische Potential üblicher Anoden. Herausforderungen sind:

  • ungleiche Abscheidung,
  • SEI-Bruch und Neubildung,
  • totes Lithium,
  • Volumenänderung der gesamten Anode,
  • interne Kurzschlüsse.

48. Dendriten und Filamente

Lokale Stromdichteerhöhungen, Konzentrationsverarmung, Defekte und mechanische Heterogenität können nadel- oder filamentartige Metallabscheidung fördern.

Ein hoher Elastizitätsmodul allein garantiert keine Unterdrückung. Grenzflächenkontakt, Elektronenleckage, Bruchmechanik und Defektpfade sind ebenso entscheidend.

49. Kritische Stromdichte

Die kritische Stromdichte ist die höchste Rate, bei der ein bestimmter Aufbau unter definierten Bedingungen stabil zyklieren kann. Sie hängt von:

  • Temperatur,
  • Druck,
  • Elektrolytdicke,
  • Kapazität pro Halbzyklus,
  • Oberflächenvorbehandlung,
  • vorhandenen Defekten

ab und ist daher keine universelle Materialkonstante.



Abbildung 2. Ionenleitung im Festelektrolyten wird durch Defekte, Korngrenzen und Raumladung beeinflusst. Die Zellpolarisation enthält zusätzlich Ladungsübertragung, ohmischen Widerstand, Konzentrationsgradienten und Grenzschichten.

50. Kathodenmaterialien

Eine positive Elektrode soll hohe Spannung, hohe reversible Kapazität, schnelle Diffusion, geringe Volumenänderung und thermische Stabilität verbinden. Diese Anforderungen stehen teilweise im Konflikt.

51. Schichtoxide

Schichtoxide wie LiMO₂ ermöglichen zweidimensionale Lithiumdiffusion. Nickel erhöht häufig die Kapazität, kann aber strukturelle und thermische Stabilität verschlechtern. Kobalt verbessert Leitfähigkeit und Struktur, ist jedoch teuer und ressourcenkritisch. Mangan kann Stabilität fördern, besitzt aber eigene Auflösungs- und Ordnungsprobleme.

52. Lithium-Eisenphosphat

LiFePO₄ besitzt eine robuste Olivinstruktur und ein ausgeprägtes Spannungsplateau. Vorteile sind thermische Stabilität und lange Lebensdauer. Nachteile sind niedrigere Spannung und geringe intrinsische elektronische Leitfähigkeit, die durch Kohlenstoffbeschichtung und kleine Partikel kompensiert wird.

53. Spinellmaterialien

Spinelle bieten dreidimensionale Diffusionswege. Hochvoltspinelle können hohe Leistung liefern, erfordern jedoch oxidationsstabile Elektrolyte und kontrollierte Grenzflächenchemie.

54. Schwefelkathoden

Schwefel besitzt eine hohe theoretische Kapazität. Die Reaktion verläuft über Polysulfide bis Li₂S. Lösliche Zwischenprodukte können zwischen den Elektroden pendeln, aktive Masse verlagern und die Coulomb-Effizienz senken.

55. Graphitanode

Graphit interkaliert Lithium schrittweise bis ungefähr LiC₆. Die geringe Spannung erhöht die Zellenergie, liegt aber nahe der Lithiumabscheidung. Bei niedriger Temperatur, hoher Rate oder hoher Ladezustandsgrenze kann Lithium-Plating auftreten.

56. Siliziumanode

Silizium besitzt eine sehr hohe theoretische Kapazität, erfährt aber enorme Volumenänderungen. Folgen sind:

  • Partikelrisse,
  • Verlust elektrischen Kontakts,
  • ständiger SEI-Neuaufbau,
  • Elektrodenschwellung.

Nanostrukturierung, elastische Binder und Verbundkonzepte vermindern, beseitigen diese Probleme aber nicht vollständig.

57. Natrium-Ionen- und andere Systeme

Natrium-Ionen-Batterien verwenden ähnliche Prinzipien, aber andere Wirtsstrukturen und Elektrolyte. Natrium ist größer und besitzt ein anderes Standardpotential. Weitere Systeme nutzen Magnesium, Zink, Aluminium oder redoxaktive organische Materialien.

Mehrwertige Ionen bieten hohe Ladung pro Ion, diffundieren aber häufig langsamer und wechselwirken stärker mit Wirtsgittern und Elektrolyten.

58. Mechanische Spannungen in Partikeln

Zusammensetzungsänderungen erzeugen chemische Dehnung. Bei räumlichen Konzentrationsgradienten entstehen Spannungen:

εchem=β(c−c₀)

Elastische Einschränkung kann Risse, plastische Verformung oder Phasengrenzenbewegung beeinflussen.

59. Partikelriss und Kontaktverlust

Risse vergrößern die Oberfläche, legen frisches Material frei und unterbrechen Leitpfade. Elektrolyt dringt in neue Flächen ein, SEI oder CEI wächst weiter, und aktive Fragmente können elektrisch isoliert werden.

60. Verlust zyklierbaren Lithiums

Loss of Lithium Inventory, kurz LLI, entsteht durch:

  • SEI- und CEI-Wachstum,
  • totes Lithium,
  • irreversible Nebenprodukte,
  • Lithiumbindung in elektrisch isolierten Bereichen.

Die Elektroden können strukturell intakt erscheinen, während der Zell-Ladezustandsbereich dennoch schrumpft.

61. Verlust aktiven Materials

Loss of Active Material, kurz LAM, umfasst:

  • Rissbildung und elektrische Isolation,
  • Auflösung und Migration von Übergangsmetallen,
  • irreversible Phasenänderungen,
  • Porenkollaps oder Delamination.

62. Impedanzanstieg

Wachsende Grenzschichten, schlechter Kontakt, Elektrolytverlust und Porenverstopfung erhöhen die Polarisation. Eine Zelle kann deshalb unter geringer Rate noch hohe Kapazität zeigen, bei hoher Rate aber stark einbrechen.

63. Kalender- und Zyklenalterung

  • Kalenderalterung: zeitabhängige Reaktionen bei Lagerung, beeinflusst von Temperatur und Ladezustand.
  • Zyklenalterung: zusätzliche Schäden durch Strom, Konzentrationsgradienten, Volumenänderung und Spannungswechsel.

Beide Mechanismen überlagern sich im Betrieb.

64. Wärmeentstehung

Eine vereinfachte Wärmebilanz enthält:

Q̇=I(V−Eeq)−IT∂Eeq/∂T+Q̇Nebenreaktionen

Der erste Term ist irreversible Polarisation, der zweite reversible entropische Wärme. Nebenreaktionen können bei Missbrauch dominieren.

65. Thermisches Durchgehen

Steigt die Wärmeproduktion schneller als die Wärmeabfuhr, kann positive Rückkopplung entstehen:

Temperatur↑ → Reaktionsrate↑ → Wärme↑ → Temperatur↑

Elektrolytzersetzung, SEI-Zerfall, Reaktion geladener Elektroden und Sauerstofffreisetzung können sich gegenseitig beschleunigen.

66. Sicherheitsdesign

Schutzmaßnahmen umfassen:

  • Separatorabschaltung,
  • Druckentlastung,
  • Stromunterbrechung,
  • Temperaturüberwachung,
  • balancierte Ladezustände,
  • flammenhemmende oder weniger flüchtige Elektrolyte,
  • thermisch stabile Elektrodenmaterialien.


Abbildung 3. Moderne Batterien kombinieren unterschiedliche Elektroden- und Elektrolytklassen. Alterung entsteht durch Lithiumverlust, aktives Material, Impedanz, Risse und Dendriten; Sicherheitsversagen koppelt Defekte, Nebenreaktionen, Gasbildung und Wärme.

67. Diagnostische Methoden

MethodePrimäre InformationWichtige Einschränkung
galvanostatisches ZyklierenKapazität, Effizienz, Rate und AlterungMechanismen sind überlagert
zyklische VoltammetrieRedoxlagen und PolarisationPoren- und Transporteffekte beeinflussen Peaks
GITT/PITTGleichgewichtsspannung und scheinbare DiffusionModellgeometrie und Relaxationszeit
ImpedanzspektroskopieBulk-, Grenzflächen- und DiffusionsbeiträgeErsatzschaltbilder sind nicht eindeutig
operando XRDPhasen und Gitterparameteramorphe und nanoskalige Phasen schwer erfassbar
Röntgenabsorptionsspektroskopielokale Struktur und Oxidationszustandräumliche Mittelung
Neutronenmethodenleichte Elemente und Lithiumverteilungaufwendige Probenumgebung
ElektronenmikroskopieRisse, Grenzflächen und NanostrukturPräparation und Strahlschäden
KalorimetrieWärmefreisetzung und SicherheitsreaktionenProbenmaßstab und Wärmeübertragung
Post-mortem-AnalytikProdukte und SchadensbilderEntnahme verändert reaktive Zustände

68. Typische Fehler

FehlerKorrektur
theoretische Elektrodenkapazität als Zellkapazität ausgebenGegenelektrode, inaktive Masse und Zellbalance berücksichtigen.
Spannungsplateau ausschließlich als langsame Kinetik interpretierenZweiphasengleichgewicht und gemeinsames chemisches Potential prüfen.
C-Rate ohne Bezugs-Nennkapazität vergleichenKapazitätsdefinition, Spannungslimits und Temperatur nennen.
hohe Leitfähigkeit mit vollständiger Elektrolytstabilität gleichsetzenelektronische Leitfähigkeit, Reaktionsprodukte und Grenzflächen prüfen.
elektrochemisches Fenster nur aus linearem Sweep ableitenStromnachweisgrenze, Elektrodenmaterial und Kinetik berücksichtigen.
GITT-Diffusionskoeffizient als universelle Materialkonstante ansehenZusammensetzung, Partikelgeometrie und Phasenregime dokumentieren.
SEI ausschließlich als schädliche Schicht behandelnSie verbraucht Lithium, ermöglicht aber oft erst kinetische Stabilität.
Festelektrolyt als automatisch dendritensicher ansehenDefekte, Elektronenleitung, Kontakt und Bruchmechanik untersuchen.
hohe Coulomb-Effizienz mit hohem Energie-Wirkungsgrad gleichsetzenSpannungshysterese und Polarisation einbeziehen.
Kapazitätsverlust nur LLI zuschreibenLAM, Impedanz und Verschiebung der Elektrodenfenster trennen.
beschleunigte Alterung ohne Mechanismenvergleich extrapolierenTemperatur-, Ladezustands- und Ratenabhängigkeit validieren.
eine sichere Einzelkomponente mit einer sicheren Zelle gleichsetzenGesamtarchitektur, Wärmewege und Fehlerketten prüfen.

69. Übungsaufgaben

Aufgabe 1: Zellspannung aus der Gibbs-Energie

Die molare Gibbs-Energie der Entladung beträgt ΔrG=−350 kJ mol−1. Pro Reaktionsumsatz wird ein Elektron übertragen.

Berechne:

E=−ΔrG/(nF)

Aufgabe 2: Theoretische Kapazität von LiFePO₄

LiFePO₄ überträgt ideal ein Elektron pro Formeleinheit. Die molare Masse beträgt 157,757 g mol−1.

Berechne:

Qth=nF/(3,6M)

Aufgabe 3: Spezifische Energie des aktiven Kathodenmaterials

Verwende die theoretische Kapazität aus Aufgabe 2 und eine mittlere Spannung von 3,40 V.

Berechne die ideale spezifische Energie des aktiven Materials.

Aufgabe 4: C-Rate

Eine Zelle besitzt eine Nennkapazität von 2,00 Ah und wird mit 0,5C betrieben.

Berechne:

  1. den Strom,
  2. die ideale Entladezeit.

Aufgabe 5: Charakteristische Diffusionszeit

Lithium diffundiert über L=10,0 µm mit D=1,00×10−14 m² s−1.

Berechne näherungsweise:

t=L²/D

Aufgabe 6: Ionenleitfähigkeit nach Nernst–Einstein

Ein idealer einwertiger Ionenleiter besitzt n=1,00×1027 mobile Ionen m−3 und D=1,00×10−11 m² s−1 bei 298,15 K.

Berechne:

σ=nq²D/(kBT)

Aufgabe 7: Flächenspezifischer Festelektrolytwiderstand

Ein Festelektrolyt ist 50,0 µm dick und besitzt σ=0,100 S m−1.

Berechne:

RA=L/σ

Gib das Ergebnis in Ω cm² an.

Aufgabe 8: Symmetrische Butler–Volmer-Überspannung

Für αac=0,5, T=298,15 K, i₀=1,00 mA cm−2 und i=10,0 mA cm−2 gilt:

η=2RT/F·asinh[i/(2i₀)]

Berechne η.

Aufgabe 9: Kapazitätserhalt nach vielen Zyklen

Eine Zelle verliert pro Zyklus konstant 0,0500 Prozent ihrer jeweils verbleibenden Kapazität.

Berechne den Kapazitätserhalt nach 500 Zyklen:

R=(1−0,000500)500

Aufgabe 10: Lithium-Metall-Dicke

Eine Lithium-Metall-Anode soll eine Flächenkapazität von 3,00 mAh cm−2 bereitstellen. Verwende:

Qspez,Li=3860 mAh g−1
ρLi=0,534 g cm−3

Berechne die ideale Lithiumdicke.

70. Vollständige Lösungen

Lösung 1
E=−(−350000 J mol−1)/[(1)(96485 C mol−1)]
E=3,627 V

Die positive Spannung entspricht einer spontanen Entladung mit negativer Gibbs-Energie.

Lösung 2
Qth=96485/[3,6(157,757)]
Qth=169,89 mAh g−1

Der praktische Wert ist wegen begrenzter Ausnutzung, Spannungslimits und inaktiver Bestandteile kleiner.

Lösung 3
w≈(169,89 mAh g−1)(3,40 V)
w=577,6 Wh kg−1

Dies ist ein Wert des aktiven Kathodenmaterials, nicht der vollständigen Zelle.

Lösung 4

Strom:

I=(0,5 h−1)(2,00 Ah)=1,00 A

Ideale Zeit:

t=1/(0,5 h−1)=2,00 h

Reale Entladezeiten hängen von Spannungslimit und Ratenkapazität ab.

Lösung 5
t=(10,0×10−6 m)²/(1,00×10−14 m² s−1)
t=1,00×104 s
t=2,78 h

Die Abschätzung zeigt die quadratische Bedeutung der Diffusionslänge.

Lösung 6
σ=(1,00×1027)(1,60218×10−19)²(1,00×10−11)/[(1,38065×10−23)(298,15)]
σ=6,236×10−2 S m−1
σ=0,6236 mS cm−1

Ionenkorrelationen können die reale Leitfähigkeit gegenüber diesem Idealwert verändern.

Lösung 7
RA=(50,0×10−6 m)/(0,100 S m−1)
RA=5,00×10−4 Ω m²

Mit 1 m²=104 cm²:

RA=5,00 Ω cm²

Eine Halbierung der Dicke halbiert diesen ohmischen Beitrag.

Lösung 8
i/(2i₀)=10,0/(2,00)=5,00
η=2(8,31446)(298,15)/96485·asinh(5,00)
η=0,1188 V

Für diese Stromdichte beträgt die reine Aktivierungsüberspannung rund 119 mV.

Lösung 9
R=(0,999500)500
R=0,7788
R=77,88 %

Ein scheinbar kleiner Verlust pro Zyklus summiert sich über viele Zyklen erheblich.

Lösung 10

Flächenmasse:

m/A=(3,00 mAh cm−2)/(3860 mAh g−1)
m/A=7,772×10−4 g cm−2

Dicke:

L=(m/A)/ρ=(7,772×10−4)/(0,534)
L=1,455×10−3 cm
L=14,55 µm

Technische Anoden verwenden oft Lithiumüberschuss, wodurch Masse und Dicke steigen.

71. Zusammenfassung

  • Batteriespannung ist eine Differenz chemischer und elektrochemischer Potentiale.
  • Die reversible Spannung erfüllt E=−ΔG/(nF).
  • Interkalation erhält das Grundgerüst des Wirtsmaterials weitgehend.
  • Festlösungen liefern geneigte, Zweiphasenreaktionen plateauförmige Spannungsverläufe.
  • Grenzflächenenergie und elastische Spannung erzeugen Hysterese und Größenabhängigkeit.
  • Legierungs- und Konversionsanoden bieten hohe Kapazität, aber große Strukturänderungen.
  • Die theoretische Kapazität folgt aus Elektronenzahl und molarer Masse.
  • Praktische Zellenergie enthält die Masse aller inaktiven Komponenten.
  • C-Rate verbindet Nennkapazität, Strom und Zeit.
  • Eine kleine irreversible Ladungsmenge pro Zyklus akkumuliert zu deutlicher Alterung.
  • Festkörperdiffusion skaliert mit L²/D.
  • Die chemische Diffusivität enthält einen thermodynamischen Faktor.
  • Poröse Elektroden koppeln Ionen- und Elektronentransport.
  • Butler–Volmer beschreibt Ladungsübertragung an Grenzflächen.
  • Zellpolarisation besteht aus Aktivierungs-, ohmischen und Konzentrationsanteilen.
  • Flüssige Elektrolyte benetzen Poren gut, sind aber oft brennbar und thermodynamisch begrenzt.
  • Festelektrolyte müssen Ionen leiten und Elektronen blockieren.
  • Polymer-, Oxid-, Sulfid- und Halogenidelektrolyte besitzen unterschiedliche Stärken.
  • Bulk-Leitfähigkeit allein bestimmt nicht den Zellwiderstand.
  • Grenzflächenreaktionen und mechanischer Kontakt können dominieren.
  • SEI und CEI ermöglichen kinetische Stabilität, verbrauchen aber Lithium und erhöhen Widerstand.
  • Lithium-Metall bietet hohe Kapazität, ist aber anfällig für ungleichmäßige Abscheidung.
  • Ein hoher Elastizitätsmodul garantiert keine Dendritenfreiheit.
  • Schichtoxide, LFP und Spinelle verkörpern unterschiedliche Spannungs- und Stabilitätskompromisse.
  • Graphit ist reversibel, Silizium kapazitätsreich und mechanisch anspruchsvoll.
  • LLI, LAM und Impedanzanstieg sind verschiedene Alterungsmodi.
  • Mechanische Risse erzeugen neue Grenzflächen und beschleunigen Nebenreaktionen.
  • Kalender- und Zyklenalterung überlagern sich.
  • Wärme entsteht irreversibel, reversibel und durch Nebenreaktionen.
  • Thermisches Durchgehen ist eine positive Rückkopplung von Temperatur und Reaktionsrate.
  • Sicherheit ist eine Eigenschaft der vollständigen Zellarchitektur.

72. Häufig gestellte Fragen

Warum besitzt ein Zweiphasenmaterial ein Spannungsplateau?

Während der Koexistenz bleibt das chemische Potential des eingelagerten Ions im idealen Gleichgewicht konstant. Daher bleibt auch die Elektrodenspannung nahezu konstant.

Ist die höchste theoretische Kapazität automatisch am besten?

Nein. Spannung, Dichte, Reversibilität, Volumenänderung, Leitfähigkeit, Sicherheit und inaktive Masse bestimmen die praktische Zellleistung.

Warum sind kleine Partikel leistungsfähiger?

Sie verkürzen Diffusionswege und reduzieren Konzentrationsgradienten. Gleichzeitig vergrößern sie die Oberfläche und damit Nebenreaktionen sowie Grenzflächenmasse.

Warum kann ein Elektrolyt außerhalb seines Stabilitätsfensters funktionieren?

Eine elektronisch isolierende, ionisch leitende Zersetzungsschicht kann weitere Reaktion kinetisch hemmen. Das ist das Grundprinzip einer funktionalen SEI.

Ist ein Festkörperakku automatisch nicht brennbar?

Nicht zwingend. Elektroden enthalten gespeicherte chemische Energie, Festelektrolyte können reagieren, und Verbundmaterialien oder Restpolymere können Wärme freisetzen.

Warum steigt der Widerstand an einer Festkörpergrenzfläche?

Chemische Reaktionsschichten, Raumladung, Poren, geringe reale Kontaktfläche, Risse und mechanischer Kontaktverlust können gleichzeitig wirken.

Was unterscheidet Coulomb- und Energie-Wirkungsgrad?

Der Coulomb-Wirkungsgrad vergleicht Ladungsmengen. Der Energie-Wirkungsgrad berücksichtigt zusätzlich die unterschiedliche Lade- und Entladespannung.

Warum ist Lithium-Plating auf Graphit gefährlich?

Metallisches Lithium kann SEI neu bilden, totes Lithium erzeugen und lokal filamentartig wachsen. Es reduziert Kapazität und erhöht Kurzschlussrisiken.

Warum kann eine Zelle bei geringer Rate noch funktionieren, bei hoher Rate aber versagen?

Diffusions-, Elektrolyt- und Kontaktwiderstände erzeugen bei hohem Strom große Polarisation. Die thermodynamisch vorhandene Kapazität wird dann nicht vollständig erreichbar.

Was ist wichtiger: Leitfähigkeit oder Stabilität?

Beides muss gemeinsam optimiert werden. Ein sehr leitfähiger, aber reaktiver Elektrolyt oder ein stabiler, aber hochohmiger Elektrolyt ist für die Zelle unzureichend.

73. Ausblick auf Teil 20

Teil 20 schließt Physikalische Chemie IV ab. Multiphysik, Materialdesign und Skalenkopplung verbinden elektronische Struktur, Defekte, Diffusion, Grenzflächenreaktionen, mechanische Spannung, Wärmeleitung und Bauteilgeometrie. Außerdem werden digitale Materialentwicklung, Unsicherheitsanalyse, Lebensdauerprognose, Nachhaltigkeit und ein Gesamtüberblick über den vierten Kursteil entwickelt.

Literatur und weiterführende Quellen

  1. M. Winter, J. O. Besenhard, M. E. Spahr & P. Novák: Insertion Electrode Materials for Rechargeable Lithium Batteries.
  2. J. Newman & K. E. Thomas-Alyea: Electrochemical Systems. Wiley.
  3. M. Z. Bazant: Thermodynamic Stability of Driven Open Systems and Control of Phase Separation by Electro-autocatalysis.
  4. C. Julien, A. Mauger, A. Vijh & K. Zaghib: Lithium Batteries: Science and Technology. Springer.
  5. N. J. J. de Klerk & M. Wagemaker: Fachliteratur zu Festelektrolyten, Grenzflächen und Ionentransport.

Didaktischer Hinweis: Ideale Nernst-Spannungen, reguläre Lösungen, konstante Diffusionskoeffizienten, Butler–Volmer-Grenzflächen, Nernst–Einstein-Leitung, ebene SEI-Schichten und exponentielle Alterung sind kontrollierte Modelle. Reale Zellen besitzen konzentrierte Elektrolyte, heterogene Poren, bewegte Phasengrenzen, gekoppelte Mechanik, Wärmegradienten, dynamische Grenzflächen und stochastische Fehlerpfade.

Bild- und Nutzungsnachweis: Das Titelbild und die Abbildungen 1–3 wurden eigens für diesen Blogartikel mit selbst erstelltem Programmcode erzeugt. Es wurden keine Grafiken aus Lehrbüchern oder anderen Publikationen übernommen. Die Abbildungen sind für die Verwendung in diesem monetarisierbaren Blog vorgesehen.

Photovoltaik und photoaktive Halbleiter – Ladungserzeugung, Rekombination und Wirkungsgrad




Titelbild: Eigene Illustration von Sonnenlicht, Halbleiter-Solarzelle und zentralen Beziehungen für Photonenergie, Diodenkennlinie, Wirkungsgrad, Diffusionslänge und Quanteneffizienz.

Photovoltaik wandelt Licht unmittelbar in elektrische freie Energie um. Ein Photon erzeugt in einem Halbleiter ein Elektron-Loch-Paar oder zunächst ein gebundenes Exziton. Damit daraus nutzbare Leistung entsteht, müssen die Ladungen getrennt, zu selektiven Kontakten transportiert und schneller extrahiert werden, als sie rekombinieren.

Die Leistung einer Solarzelle ist deshalb keine reine Absorptionsfrage. Sie entsteht aus dem Zusammenspiel von:

  • Bandlücke und Absorptionskoeffizient,
  • Defekten und Ladungsträgerlebensdauer,
  • Dotierung und Raumladungszone,
  • Drift, Diffusion und Kontaktselektivität,
  • Oberflächen- und Grenzflächenpassivierung,
  • optischer Einkopplung und Reflexionsverlusten,
  • Serien- und Parallelwiderständen,
  • thermischer, chemischer und mechanischer Stabilität.

Eine ideale Solarzelle ist gleichzeitig Lichtabsorber, Diode, Ladungstransportstruktur und thermodynamische Wärmekraftmaschine mit einem nichtthermischen Strahlungsfeld als Energiequelle.

Einordnung: Teil 7 behandelte Halbleiter, Dotierung und p-n-Übergänge, Teil 10 optische Eigenschaften und Teil 11 Nanomaterialien. Teil 18 verbindet diese Grundlagen zu photovoltaischen Bauelementen. Teil 19 führt anschließend zu Batterie- und Festelektrolytmaterialien.

Lernziele

Nach diesem Beitrag kannst du:
  • Photonenergie, Bandlücke und Absorptionskante verknüpfen,
  • direkte und indirekte optische Übergänge unterscheiden,
  • Generation mit dem Beer–Lambert-Gesetz beschreiben,
  • Exzitonen und freie Ladungsträger einordnen,
  • p-n-Übergang, Raumladungszone und eingebautes Feld erklären,
  • Drift-, Diffusions- und Kontinuitätsgleichungen verwenden,
  • Quasi-Fermi-Niveaus interpretieren,
  • radiative, SRH-, Auger- und Oberflächenrekombination vergleichen,
  • Jsc, Voc, Füllfaktor und Wirkungsgrad berechnen,
  • EQE, IQE und spektrale Antwort unterscheiden,
  • Shockley–Queisser-Verluste und Tandemprinzip erklären,
  • Silizium-, Dünnschicht-, Perowskit- und organische Solarzellen vergleichen,
  • Degradationspfade und Messfehler erkennen.

1. Energie eines Photons

Eph=hν=hc/λ

Kurzwelliges Licht besitzt eine höhere Photonenergie. Für praktische Rechnungen gilt:

Eph/eV≈1239,84/(λ/nm)

Ein Halbleiter kann ein Photon über einen Band-zu-Band-Übergang absorbieren, wenn dessen Energie mindestens der wirksamen optischen Übergangsenergie entspricht.

2. Bandlücke und Absorptionskante

Die Bandlücke Eg trennt Valenz- und Leitungsband. Die ideale Grenzwellenlänge lautet:

λg=hc/Eg

Photonen mit E<Eg durchlaufen den idealen Kristall weitgehend, sofern keine Defekt-, Exzitonen- oder freien Ladungsträgerübergänge auftreten.

3. Direkte Bandlücke

Bei einer direkten Bandlücke liegen Valenzbandmaximum und Leitungsbandminimum beim gleichen Kristallimpuls. Ein Photon kann den Übergang ohne zusätzliches Phonon ermöglichen.

Direkte Halbleiter besitzen nahe der Bandkante häufig hohe Absorptionskoeffizienten. Dünne Schichten können daher bereits einen großen Lichtanteil absorbieren.

4. Indirekte Bandlücke

Bei einer indirekten Bandlücke unterscheiden sich die Impulse der Bandextrema. Zusätzlich zum Photon wird ein Phonon benötigt, damit Kristallimpuls erhalten bleibt.

Silizium absorbiert deshalb nahe seiner Bandkante schwächer als ein direkter Halbleiter. Optische Wegverlängerung und dickere Absorber kompensieren diesen Nachteil.

5. Absorptionskoeffizient

Für monochromatisches Licht in einem homogenen Material:

I(x)=I0exp(−αx)

α ist der Absorptionskoeffizient. Ohne Reflexion ist der absorbierte Anteil einer Schichtdicke d:

A=1−exp(−αd)

Ein hoher α-Wert erlaubt dünne Absorber, verschärft aber die Bedeutung von Oberflächen und Grenzflächen.

6. Optische Generation

Ist jedes absorbierte Photon zunächst ein Elektron-Loch-Ereignis, lautet die lokale Generationsrate:

G(x)=αΦ0exp(−αx)

Φ0 ist der einfallende Photonenfluss. Reflexion, Interferenz und strukturierte Schichten verändern das einfache exponentielle Profil.

7. Reflexion und Antireflexschichten

Ein Sprung des Brechungsindex erzeugt Fresnel-Reflexion. Antireflexschichten nutzen destruktive Interferenz der reflektierten Teilwellen. Für eine einfache Viertelwellenschicht gilt näherungsweise:

d=λ0/(4n)

Texturierte Oberflächen verlängern zusätzlich den optischen Weg und reduzieren die direkte Rückreflexion.

8. Exzitonen

Elektron und Loch ziehen sich Coulombisch an. Die Bindungsenergie eines wasserstoffähnlichen Wannier-Exzitons ist:

EB≈13,6 eV·(μ/me)/εr²

Ist EB deutlich kleiner als kBT, entstehen bei Raumtemperatur leicht freie Ladungsträger. In organischen Halbleitern sind Exzitonen wegen kleiner Permittivität und lokalisierter Zustände stärker gebunden.

9. Exzitonendissoziation

Gebundene Exzitonen werden an:

  • Donor-Akzeptor-Grenzflächen,
  • Heteroübergängen,
  • starken elektrischen Feldern,
  • energetisch versetzten Kontakten

getrennt. Die Exzitonendiffusionslänge begrenzt, wie weit die Erzeugung von einer trennenden Grenzfläche entfernt stattfinden darf.

10. p-n-Übergang im Gleichgewicht

Werden p- und n-dotierte Bereiche verbunden, diffundieren Elektronen und Löcher über die Grenzfläche. Zurück bleiben ionisierte Dotieratome. Die daraus entstehende Raumladung erzeugt ein eingebautes elektrisches Feld.

11. Eingebaute Spannung

Für einen abrupten, nichtentarteten Homojunction-Übergang:

Vbi=kBT/q·ln(NAND/ni²)

Im thermodynamischen Gleichgewicht ist das Fermi-Niveau räumlich konstant; Leitungs- und Valenzband biegen sich.

12. Verarmungsbreite

W=√{2εs/q·(1/NA+1/ND)(Vbi−V)}

Bei Vorwärtsspannung wird die Raumladungszone schmaler, bei Sperrspannung breiter. Der schwächer dotierte Bereich trägt den größeren Anteil der Breite.

13. Ladungstrennung

Elektron-Loch-Paare, die in oder nahe der Raumladungszone entstehen, werden vom eingebauten Feld getrennt. Elektronen driften zur n-Seite, Löcher zur p-Seite.

Außerhalb der Raumladungszone müssen Minoritätsladungsträger durch Diffusion die Sammelregion erreichen.

14. Driftstrom

Jn,Drift=qμnnE
Jp,Drift=qμppE

μ ist die Beweglichkeit und E das elektrische Feld.

15. Diffusionsstrom

Jn,Diff=qDn∇n
Jp,Diff=−qDp∇p

Die Vorzeichen folgen aus Ladung und Teilchenfluss. Für nichtentartete Träger gilt die Einstein-Beziehung:

D/μ=kBT/q

16. Kontinuitätsgleichungen

∂n/∂t=(1/q)∇·Jn+G−R
∂p/∂t=−(1/q)∇·Jp+G−R

Generation G, Rekombination R und Stromdivergenz bestimmen gemeinsam die lokale Trägerdichte.

17. Diffusionslänge

Ln=√(Dnτn)
Lp=√(Dpτp)

Die Diffusionslänge ist eine zentrale Sammelskala. Ein Absorber darf optisch dick, aber elektrisch nicht zu dick sein.

18. Quasi-Fermi-Niveaus

Unter Beleuchtung sind Elektronen und Löcher nicht mehr durch ein einziges Fermi-Niveau beschreibbar. Es entstehen EFn und EFp.

np=ni²exp[(EFn−EFp)/(kBT)]

Die Aufspaltung der Quasi-Fermi-Niveaus misst die maximal lokal verfügbare elektrochemische Spannung.

19. Selektive Kontakte

Ein Elektronenkontakt soll Elektronen übertragen und Löcher blockieren; ein Lochkontakt entsprechend umgekehrt. Ungeeignete Bandversätze erhöhen Kontaktwiderstand oder Rekombination.

20. Heteroübergänge

Bei einem Heteroübergang treffen verschiedene Halbleiter aufeinander. Bandversätze können Ladungstrennung fördern, aber auch Barrieren bilden. Grenzflächenzustände und chemische Reaktionen müssen passiviert werden.



Abbildung 1. Photonen erzeugen Elektron-Loch-Paare. Eingebautes Feld und Konzentrationsgradienten sammeln die Ladungen; Defekte, Oberflächen und hohe Trägerdichten öffnen konkurrierende Rekombinationskanäle.

21. Strahlende Rekombination

Bei einem direkten Band-zu-Band-Prozess rekombinieren Elektron und Loch unter Photonemission:

Rrad=B(np−ni²)

Strahlende Rekombination ist in Leuchtdioden erwünscht. In Solarzellen ist sie ein Verlustkanal, zugleich aber der unvermeidbare Gegenprozess zur Lichtabsorption. Eine hochwertige Solarzelle nähert sich im Idealfall dem strahlenden Grenzfall.

22. Shockley–Read–Hall-Rekombination

Tiefe Defektniveaus in der Bandlücke fangen nacheinander Elektron und Loch ein:

RSRH=(np−ni²)/[τp0(n+n₁)+τn0(p+p₁)]

Defekte nahe der Bandmitte sind häufig besonders wirksame Rekombinationszentren. Metallverunreinigungen, Leerstellenkomplexe und Grenzflächenzustände können die Lebensdauer stark verkürzen.

23. Auger-Rekombination

Bei der Auger-Rekombination wird die freiwerdende Energie auf einen dritten Ladungsträger übertragen:

RAuger=(Cnn+Cpp)(np−ni²)

Der Prozess wird bei hoher Dotierung oder starker Beleuchtung wichtig. Er begrenzt unter anderem die Spannung sehr hochwertiger Siliziumzellen.

24. Oberflächenrekombination

An einer unpassivierten Oberfläche unterbrechen Bindungen die Periodizität. Zustände in der Bandlücke ermöglichen SRH-ähnliche Rekombination. Ein einfaches Maß ist die Oberflächenrekombinationsgeschwindigkeit S:

Jrec,surf≈qSΔn

Chemische Passivierung reduziert die Zustandsdichte, Feldeffektpassivierung verdrängt eine Trägerart von der Oberfläche.

25. Lebensdauer

Bei kleinen Überschussdichten wird häufig eine effektive Lebensdauer definiert:

1/τeff=1/τrad+1/τSRH+1/τAuger+1/τsurf

Diese additive Form ist eine Näherung, weil einzelne Raten unterschiedlich von Trägerdichte und Ort abhängen.

26. Dunkelkennlinie einer idealen Diode

JD=J₀[exp(qV/(nkBT))−1]

J₀ ist die Sperrsättigungsstromdichte, n der Idealitätsfaktor. Für diffusionsbestimmte Rekombination liegt n häufig nahe 1, für SRH-Rekombination in der Raumladungszone näher bei 2.

27. Beleuchtete Kennlinie

Im idealen Superpositionsmodell:

J(V)=Jph−J₀[exp(qV/(nkBT))−1]

Die Vorzeichenkonvention kann in Diagrammen umgekehrt gewählt werden. Physikalisch verschiebt der Photostrom die Dunkelkennlinie.

28. Kurzschlussstromdichte

Bei V=0 gilt ideal:

Jsc≈Jph

Jsc hängt von Photonenfluss, Absorption und Sammlung ab. Eine hohe Bandlücke reduziert den nutzbaren Spektralbereich; eine niedrige Bandlücke erhöht den Photonenstrom, aber auch thermische und Spannungsverluste.

29. Leerlaufspannung

Bei J=0:

Voc=nkBT/q·ln(Jph/J₀+1)

Voc steigt nur logarithmisch mit Beleuchtungsintensität, sinkt aber stark mit einem erhöhten J₀. Defektpassivierung verbessert deshalb vor allem die Spannung.

30. Quasi-Fermi-Aufspaltung und Voc

Unter ideal selektiven Kontakten entspricht die externe Leerlaufspannung näherungsweise der Quasi-Fermi-Niveau-Aufspaltung zwischen den Kontakten:

qVoc≈EFn−EFp

Kontaktverluste, Bandverbiegung und Grenzflächenrekombination können die extern messbare Spannung reduzieren.

31. Maximum-Power-Point

Die elektrische Leistungsdichte ist:

P(V)=VJ(V)

Am Maximum-Power-Point gelten Vmp und Jmp. Ein Leistungselektroniksystem verfolgt diesen Punkt unter wechselnder Einstrahlung und Temperatur.

32. Füllfaktor

FF=VmpJmp/(VocJsc)

Der Füllfaktor misst die Rechteckigkeit der Kennlinie. Er sinkt durch:

  • Serienwiderstand,
  • Leckpfade und niedrigen Parallelwiderstand,
  • nichtideale Rekombination,
  • Kontaktbarrieren,
  • räumliche Inhomogenität.

33. Wirkungsgrad

η=VocJscFF/Pin

Unter Standardmessbedingungen wird eine definierte spektrale Bestrahlungsstärke verwendet. Temperatur, Spektrum, Einfallswinkel und Stabilisierung der Zellleistung müssen dokumentiert werden.

34. Serienwiderstand

Rs umfasst Widerstände von Absorber, transparenten Elektroden, Metallgittern, Kontakten und Verbindungen. Bei hohem Strom entsteht:

ΔV=JRs

Der Spannungsabfall reduziert vor allem den Füllfaktor.

35. Parallelwiderstand

Rsh beschreibt Leckströme durch Pinholes, Kanten, Defekte oder unerwünschte leitfähige Pfade. Ein kleiner Rsh reduziert Spannung und Füllfaktor, besonders bei niedriger Beleuchtung.

36. Erweitertes Ersatzschaltbild

J=Jph−J₀{exp[q(V+JRs)/(nkBT)]−1}−(V+JRs)/Rsh

Ein einzelnes Diodenmodell ist praktisch, kann aber mehrere Rekombinations- und Transportpfade nur effektiv zusammenfassen.

37. Externe Quanteneffizienz

EQE(λ)=gesammelte Elektronen/einfallende Photonen

EQE enthält Reflexions-, Absorptions- und Sammlungsverluste. Die Kurzschlussstromdichte folgt aus dem Spektrum:

Jsc=q∫EQE(λ)Φ(λ)dλ

38. Interne Quanteneffizienz

IQE(λ)=gesammelte Elektronen/absorbierte Photonen

Mit Reflexion R(λ) und Transmission T(λ) gilt näherungsweise:

IQE=EQE/[1−R−T]

IQE isoliert die elektrische Sammlung, bleibt aber von der Genauigkeit der optischen Messungen abhängig.

39. Spektrale Antwort

Die spektrale Stromantwort besitzt die Einheit A W−1:

SR(λ)=qλEQE(λ)/(hc)

EQE und spektrale Antwort enthalten dieselbe Information in unterschiedlicher Normierung.

40. Temperaturkoeffizienten

Mit steigender Temperatur sinkt die Bandlücke vieler Halbleiter. Jsc kann leicht steigen, Voc sinkt jedoch meist stärker, weil J₀ zunimmt. Der Wirkungsgrad kristalliner Siliziumzellen nimmt deshalb typischerweise mit Temperatur ab.

41. Fundamentale Spektralverluste

Eine Einzelzelle mit Bandlücke Eg verliert Energie auf zwei unvermeidbare Weisen:

  • Photonen mit E<Eg werden nicht absorbiert,
  • bei E>Eg thermalisiert die Überschussenergie auf die Bandkante.

Eine einzelne Bandlücke kann das breite Sonnenspektrum daher nicht vollständig nutzen.

42. Detailliertes Gleichgewicht

Eine ideale absorbierende Zelle muss nach Reziprozität auch emittieren. Im Leerlauf wächst die strahlende Emission, bis Generation und Rekombination ausgeglichen sind. Aus dieser Bilanz folgen J₀, Voc und die fundamentale Wirkungsgradgrenze.

43. Shockley–Queisser-Grenze

Für eine ideale Einzelübergangszelle unter unkonzentrierter Sonneneinstrahlung liegt das Maximum bei einer Bandlücke um etwa 1,3 bis 1,4 eV und einem Wirkungsgrad von ungefähr einem Drittel.

Die Grenze setzt voraus:

  • vollständige Absorption oberhalb Eg,
  • ausschließlich strahlende Rekombination,
  • thermisches Gleichgewicht der Ladungsträger,
  • einen einzigen p-n-Übergang,
  • keine optischen oder ohmschen Zusatzverluste.

44. Spannungseinbuße

Die Differenz zwischen Bandlücke und externer Spannung wird häufig als Spannungseinbuße betrachtet:

ΔEV=Eg−qVoc

Ein Teil ist thermodynamisch unvermeidbar; zusätzliche nichtstrahlende Rekombination vergrößert die Einbuße.

45. Lumineszenz als Qualitätsmaß

Eine hohe externe Lumineszenzeffizienz im Vorwärtsbetrieb zeigt, dass nichtstrahlende Rekombination klein ist und Photonen den Baustein verlassen können. Gute Solarzellen können daher auch effiziente Leuchtdioden im umgekehrten Betrieb sein.



Abbildung 2. Kurzschlussstrom, Leerlaufspannung und Füllfaktor erfassen verschiedene Verlustklassen. Die spektrale Quanteneffizienz verbindet Photonenfluss und messbaren Solarstrom.

46. Kristallines Silizium

Silizium besitzt eine indirekte Bandlücke von rund 1,1 eV. Industrielle Zellen verwenden deshalb vergleichsweise dicke Wafer, Oberflächentexturen und optische Rückreflektoren.

Entscheidende Qualitätsmerkmale sind:

  • hohe Volumenlebensdauer,
  • sehr gute Oberflächenpassivierung,
  • selektive Kontakte,
  • geringer Metall- und Schichtwiderstand,
  • niedrige optische Abschattung.

47. Homojunction- und Heterojunction-Siliziumzellen

Klassische Zellen verwenden dotierte Siliziumbereiche. Heterojunction-Konzepte kombinieren kristallines Silizium mit dünnen amorphen oder nanokristallinen selektiven Schichten. Passivierte Kontaktkonzepte trennen Metallkontakt und kristallines Silizium durch ultradünne dielektrische oder halbleitende Schichten.

48. CdTe- und CIGS-Dünnschichtzellen

CdTe und Cu(In,Ga)Se₂ besitzen direkte Bandlücken und hohe Absorption. Wenige Mikrometer Material können ausreichen. Vorteile sind geringer Materialbedarf und monolithische Modulverschaltung.

Herausforderungen sind Defektchemie, Grenzflächenkontrolle, Elementverfügbarkeit, Recycling und langfristige Kontaktstabilität.

49. Metallhalogenid-Perowskite

Perowskitabsorber kombinieren hohe Absorption, lange effektive Diffusionslängen und abstimmbare Bandlücken. Sie können bei niedrigen Prozesstemperaturen aus Lösung oder Dampf hergestellt werden.

Besondere physikalische Eigenschaften sind:

  • mobile Ionen,
  • weiche polare Gitter,
  • ausgeprägte Defekt- und Grenzflächenchemie,
  • mögliche Strom-Spannungs-Hysterese,
  • Bandlückenabstimmung durch Halogenid- und Kationenzusammensetzung.

50. Organische Solarzellen

Organische Halbleiter erzeugen meist stark gebundene Exzitonen. Eine Donor-Akzeptor-Grenzfläche trennt sie in einen Ladungstransferzustand und anschließend freie Träger.

Eine Bulk-Heterojunction verteilt die Grenzfläche nanoskalig im gesamten Film. Die Morphologie muss gleichzeitig Exzitonendiffusion, Ladungstrennung und durchgängige Transportpfade ermöglichen.

51. Farbstoffsolarzellen

Ein Farbstoff absorbiert Licht und injiziert ein Elektron in ein poröses Oxid. Ein Redoxmediator regeneriert den Farbstoff. Die Funktionen Lichtabsorption, Elektronentransport und Lochtransport sind räumlich auf verschiedene Komponenten verteilt.

52. Quantenpunkt-Solarzellen

Quantenpunkte besitzen größenabhängige Bandlücken. Oberflächenzustände, Liganden und Kopplung zwischen Punkten bestimmen jedoch den Transport. Kurze Liganden verbessern die Leitfähigkeit, können aber Passivierung und Stabilität verschlechtern.

53. Tandemsolarzellen

Tandems stapeln Teilzellen mit unterschiedlichen Bandlücken. Die obere Zelle absorbiert hochenergetische Photonen, die untere nutzt den langwelligen Rest.

Dadurch werden sowohl Transmissions- als auch Thermalisierungsverluste reduziert.

54. Zwei- und Vierterminal-Tandems

  • Zweipolig: Teilzellen sind in Serie; Stromanpassung ist erforderlich, Spannungen addieren sich.
  • Vierpolig: Teilzellen arbeiten elektrisch getrennt; weniger Stromanpassungszwang, aber zusätzliche Kontakte und Optik.

55. Stromanpassung

In einer idealen Serienschaltung fließt durch alle Teilzellen derselbe Strom. Die schwächere Teilzelle begrenzt daher den Tandemstrom:

JTandem≈min(JTop,JBottom)

Spektrum, Schichtdicken und Bandlücken müssen gemeinsam optimiert werden.

56. Rekombinationskontakt im Tandem

Zwischen zwei seriellen Teilzellen müssen Elektronen der einen und Löcher der anderen Teilzelle effizient rekombinieren. Der Zwischenkontakt soll optisch transparent und elektrisch verlustarm sein.

57. Photoelektrochemische Energieumwandlung

Ein Photoelektroden-Halbleiter kann Lichtenergie für elektrochemische Reaktionen bereitstellen, beispielsweise Wasserstoffentwicklung oder Wasseroxidation. Erforderlich sind:

  • geeignete Bandkanten,
  • ausreichende Photospannung,
  • katalytisch schnelle Grenzflächenreaktionen,
  • Beständigkeit im Elektrolyten.

58. Photokorrosion

Photogenerierte Ladungen können nicht nur die gewünschte Reaktion, sondern auch den Halbleiter selbst oxidieren oder reduzieren. Schutzschichten müssen Ladung übertragen, den Elektrolyten aber chemisch vom Absorber trennen.

59. Stabilität und Degradation

Reale Solarzellen altern durch gekoppelte Einflüsse:

  • Feuchte und Sauerstoff,
  • UV- und sichtbares Licht,
  • Temperaturzyklen,
  • elektrisches Feld und Strom,
  • mobile Ionen und Metallmigration,
  • mechanische Spannung und Delamination.

60. Verkapselung

Eine Verkapselung reduziert den Eintritt von Wasser und Sauerstoff, muss aber transparent, elektrisch isolierend, mechanisch beständig und thermisch kompatibel sein. Kantenabdichtung ist häufig besonders kritisch.

61. Lichtinduzierte Veränderungen

Licht kann Defektladungszustände, Ionenverteilung, lokale Zusammensetzung und Grenzflächenreaktionen verändern. Eine anfängliche Leistungssteigerung oder -minderung kann reversibel sein und muss von irreversibler Degradation getrennt werden.

62. Modulverluste

Vom Zell- zum Modulwirkungsgrad entstehen zusätzliche Verluste durch:

  • Zwischenräume und Randflächen,
  • Verschaltung und Kontaktbänder,
  • optische Verkapselung,
  • Zellfehlanpassung,
  • Teilverschattung und Bypassdioden,
  • Temperaturerhöhung im Betrieb.


Abbildung 3. Photovoltaische Technologien verteilen Absorption, Exzitonentrennung und Ladungstransport unterschiedlich. Tandems nutzen mehrere Bandlücken; technische Leistung erfordert außerdem chemisch und mechanisch stabile Schichtstapel.

63. Experimentelle Methoden

MethodePrimäre InformationWichtige Einschränkung
J-V unter SonnensimulatorJsc, Voc, FF und ηSpektralanpassung, Zelltemperatur und Scanprotokoll
EQEspektrale LadungssammlungBiaslicht und elektrische Vorspannung können nötig sein
PhotolumineszenzQuasi-Fermi-Aufspaltung und nichtstrahlende Verlusteoptische Auskopplung und Reabsorption
zeitaufgelöste Lumineszenzeffektive Rekombinationsdynamikmehrere Prozesse und injektionsabhängige Lebensdauer
Kapazität-SpannungRaumladung und scheinbare DotierungFallen, Ionenbewegung und Frequenzdispersion
ImpedanzspektroskopieTransport- und RekombinationszeitkonstantenErsatzschaltbilder sind nicht eindeutig
LBIC/EL-Bildgebungräumliche Defekte und StrominhomogenitätAuflösung und Kontaktgeometrie
StabilitätstestLeistungsentwicklung unter definierter BelastungProtokoll muss reale Belastung repräsentieren

64. Typische Fehler

FehlerKorrektur
jedes absorbierte Photon automatisch als gesammeltes Elektron zählenExzitonentrennung, Rekombination und Kontaktverluste berücksichtigen.
Bandlücke direkt mit Voc gleichsetzenthermodynamische und nichtstrahlende Spannungseinbußen abziehen.
Raumladungsfeld als einzigen Sammelmechanismus ansehenDiffusion in quasineutralen Bereichen einbeziehen.
lange Lebensdauer ohne Beweglichkeit als ausreichend ansehenDiffusionslänge L=√(Dτ) und Kontaktselektivität bewerten.
hohe Lumineszenz grundsätzlich als Solarzellenverlust interpretierenIm strahlenden Grenzfall zeigt Lumineszenz geringe nichtstrahlende Rekombination.
EQE und IQE verwechselnIQE auf absorbierte, EQE auf einfallende Photonen beziehen.
J-V-Scan ohne Stabilisierung als stationären Wirkungsgrad angebenScanrichtung, Geschwindigkeit und stabilisierten MPP dokumentieren.
hohen Serienwiderstand mit Rekombination verwechselnKennlinienform, Suns-Voc und Impedanz ergänzend auswerten.
Shockley–Queisser-Grenze als Grenze aller Photovoltaik ansehenSie gilt für eine ideale Einzelübergangszelle unter definierten Annahmen.
beste Zellleistung direkt auf Module übertragenVerschaltung, Temperatur, Verkapselung und Flächenverluste einbeziehen.
Tandemspannungen addieren, ohne Stromanpassung zu prüfenBei Zweipol-Tandems begrenzt die schwächere Teilzelle den Strom.
kurze beschleunigte Tests direkt in jahrzehntelange Lebensdauer umrechnenMechanismusgleichheit und validierte Beschleunigungsmodelle verlangen.

65. Übungsaufgaben

Aufgabe 1: Photonenergie

Berechne die Energie eines Photons mit λ=550 nm in Joule und Elektronenvolt.

E=hc/λ

Aufgabe 2: Absorbierter Lichtanteil

Ein direkter Halbleiter besitzt bei einer bestimmten Wellenlänge α=1,00×105 cm−1. Die Schicht ist 500 nm dick. Reflexion werde vernachlässigt.

Berechne:

A=1−exp(−αd)

Aufgabe 3: Eingebaute Spannung

Ein abruptes Silizium-p-n-System besitzt bei 300 K:

NA=1,00×1023 m−3
ND=1,00×1022 m−3
ni=1,00×1016 m−3

Berechne Vbi.

Aufgabe 4: Verarmungsbreite

Verwende die Daten aus Aufgabe 3 sowie εr=11,7. Berechne bei V=0:

W=√{2εs/q·(1/NA+1/ND)Vbi}

Bestimme außerdem die Breitenanteile xn und xp mit:

xn=WNA/(NA+ND)
xp=WND/(NA+ND)

Aufgabe 5: Diffusionslänge

Für Minoritätsladungsträger gelten D=2,50×10−3 m² s−1 und τ=1,00 µs.

Berechne:

L=√(Dτ)

Aufgabe 6: Leerlaufspannung

Eine ideale Solarzelle besitzt bei 300 K:

Jsc=35,0 mA cm−2
J₀=1,00×10−12 A cm−2, n=1

Berechne:

Voc=nkBT/q·ln(Jsc/J₀+1)

Aufgabe 7: Idealer Füllfaktor

Verwende Voc aus Aufgabe 6. Für eine ideale Diode kann der Füllfaktor näherungsweise bestimmt werden durch:

voc=qVoc/(nkBT)
FF≈[voc−ln(voc+0,72)]/(voc+1)

Aufgabe 8: Solarzellenwirkungsgrad

Eine Zelle besitzt Jsc=35,0 mA cm−2, Voc=0,650 V und FF=0,800. Die eingestrahlte Leistung beträgt 100 mW cm−2.

Berechne η.

Aufgabe 9: Monochromatischer EQE-Strom

Ein monochromatischer Photonenfluss beträgt 2,00×1017 Photonen cm−2 s−1. Die EQE beträgt 0,850.

Berechne die erzeugte Stromdichte:

J=qΦ·EQE

Aufgabe 10: Serielles Tandem

Eine Zweipol-Tandemzelle besitzt Teilzellspannungen von 0,750 V und 1,10 V. Die Teilzellströme betragen 18,0 und 16,0 mA cm−2. Der Tandemfüllfaktor beträgt 0,820, Pin=100 mW cm−2.

Berechne Strom, Spannung und Wirkungsgrad des ideal serienverschalteten Tandems.

66. Vollständige Lösungen

Lösung 1
E=(6,62607×10−34)(2,99792×108)/(550×10−9)
E=3,608×10−19 J

Umrechnung in Elektronenvolt:

E=(3,608×10−19)/(1,60218×10−19)
E=2,254 eV

Ein Material mit optischer Übergangsenergie über 2,254 eV könnte dieses Photon nicht durch einen einfachen Band-zu-Band-Prozess absorbieren.

Lösung 2

Die Dicke beträgt:

d=500 nm=5,00×10−5 cm

Damit:

αd=(1,00×105)(5,00×10−5)=5,00
A=1−exp(−5,00)=0,9933

Es werden 99,33 Prozent des eintretenden Lichts absorbiert. Reale Verluste durch Reflexion sind hier nicht berücksichtigt.

Lösung 3

Thermische Spannung bei 300 K:

kBT/q=0,025852 V

Dotierungsverhältnis:

NAND/ni²=1045/1032=1013
Vbi=0,025852 ln(1013)
Vbi=0,7738 V

Die Formel setzt nichtentartete Dotierung und vollständige Ionisation voraus.

Lösung 4

Die Permittivität ist:

εs=11,7ε₀=1,036×10−10 F m−1

Einsetzen:

W=√{2(1,036×10−10)/(1,602×10−19)·(10−23+10−22)(0,7738)}
W=3,318×10−7 m=331,8 nm

Aufteilung:

xn=331,8·10/11=301,6 nm
xp=331,8·1/11=30,16 nm

Der schwächer dotierte n-Bereich trägt den größten Anteil der Raumladungszone.

Lösung 5
L=√[(2,50×10−3)(1,00×10−6)]
L=5,00×10−5 m
L=50,0 µm

Ein außerhalb der Raumladungszone erzeugter Träger kann über diese charakteristische Länge diffundieren, bevor er rekombiniert.

Lösung 6
Jsc/J₀=(35,0×10−3)/(1,00×10−12)=3,50×1010
Voc=0,025852 ln(3,50×1010+1)
Voc=0,6277 V

Eine Verringerung von J₀ verbessert Voc logarithmisch. Dafür ist meist eine starke Reduktion nichtstrahlender Rekombination erforderlich.

Lösung 7
voc=0,6277/0,025852=24,279
FF≈[24,279−ln(24,279+0,72)]/(24,279+1)
FF≈0,8331

Dies ist ein idealer Näherungswert ohne Serien- und Parallelwiderstandsverluste.

Lösung 8
η=(0,0350)(0,650)(0,800)/(0,100)
η=0,182
η=18,2 %

Die Einheiten kürzen sich, weil Strom- und Eingangsleistungsdichte auf dieselbe Fläche bezogen sind.

Lösung 9
J=(1,60218×10−19)(2,00×1017)(0,850)
J=2,724×10−2 A cm−2
J=27,24 mA cm−2

Die Rechnung setzt voraus, dass die EQE für diesen gesamten monochromatischen Fluss konstant ist.

Lösung 10

Bei Serienschaltung begrenzt die schwächere Teilzelle den Strom:

J=16,0 mA cm−2

Die Spannungen addieren sich:

V=0,750+1,10=1,85 V

Wirkungsgrad:

η=(0,0160)(1,85)(0,820)/(0,100)
η=0,2427
η=24,27 %

Der zusätzliche Strom der oberen Teilzelle kann in einer idealen Zweipol-Serienstruktur nicht genutzt werden. Eine bessere Stromanpassung würde den Wirkungsgrad erhöhen.

67. Zusammenfassung

  • Photonenergie ist E=hc/λ.
  • Die Bandlücke legt die ideale langwellige Absorptionsgrenze fest.
  • Direkte Halbleiter absorbieren nahe der Bandkante stärker als indirekte.
  • Das Beer–Lambert-Gesetz beschreibt die exponentielle Intensitätsabnahme.
  • Optische Generation hängt von Absorptionskoeffizient und Photonenfluss ab.
  • Antireflexschichten und Texturen reduzieren optische Verluste.
  • Exzitonen sind Coulomb-gebundene Elektron-Loch-Paare.
  • Organische Solarzellen benötigen Donor-Akzeptor-Grenzflächen zur Exzitonentrennung.
  • Der p-n-Übergang erzeugt Raumladung, Bandverbiegung und ein eingebautes Feld.
  • Drift sammelt Träger im Feld, Diffusion entlang von Konzentrationsgradienten.
  • Diffusionslänge verbindet Beweglichkeit, Diffusion und Lebensdauer.
  • Quasi-Fermi-Niveaus beschreiben beleuchtete Nichtgleichgewichtsträger.
  • Selektive Kontakte übertragen eine Trägerart und blockieren die andere.
  • Radiative Rekombination ist der reziproke Prozess zur Absorption.
  • SRH-Rekombination wird durch Defektzustände in der Bandlücke vermittelt.
  • Auger-Rekombination wird bei hoher Trägerdichte wichtig.
  • Oberflächenpassivierung reduziert Grenzflächenrekombination.
  • Die beleuchtete Solarzelle entspricht näherungsweise einer Diode plus Photostromquelle.
  • Jsc misst vor allem spektrale Absorption und Sammlung.
  • Voc reagiert empfindlich auf die Sperrsättigungsstromdichte.
  • Der Füllfaktor misst die Rechteckigkeit der Kennlinie.
  • Serienwiderstand reduziert vor allem FF, Leckpfade reduzieren Spannung und FF.
  • Der Wirkungsgrad ist JscVocFF/Pin.
  • EQE bezieht sich auf einfallende, IQE auf absorbierte Photonen.
  • Spektrale Antwort und EQE sind unterschiedlich normierte Messgrößen.
  • Einzelzellen verlieren langwellige Photonen durch Transmission und Überschussenergie durch Thermalisierung.
  • Die Shockley–Queisser-Grenze folgt aus detailliertem Gleichgewicht.
  • Eine hohe Lumineszenzeffizienz kann geringe nichtstrahlende Rekombination anzeigen.
  • Silizium benötigt wegen seiner indirekten Bandlücke starke Lichtführung.
  • CdTe und CIGS ermöglichen dünne direkte Absorber.
  • Perowskite verbinden starke Absorption mit beweglichen Ionen und anspruchsvoller Stabilität.
  • Organische Zellen benötigen nanoskalig kontrollierte Donor-Akzeptor-Morphologien.
  • Tandems reduzieren Transmissions- und Thermalisierungsverluste.
  • Zweipol-Tandems benötigen Stromanpassung.
  • Photoelektroden müssen Bandkanten, Katalyse und Korrosionsbeständigkeit vereinen.
  • Langzeitstabilität erfordert kontrollierte Chemie, Kontakte, Mechanik und Verkapselung.

68. Häufig gestellte Fragen

Erzeugt jedes Photon in einer Solarzelle Strom?

Nein. Das Photon kann reflektiert, transmittiert oder parasitär absorbiert werden. Erzeugte Ladungsträger können außerdem rekombinieren, bevor sie die Kontakte erreichen.

Warum ist Voc kleiner als Eg/q?

Die Solarzelle benötigt endliche Entropieproduktion und strahlende Emission. Nichtstrahlende Rekombination und Kontaktverluste vergrößern die Spannungseinbuße zusätzlich.

Warum braucht Silizium einen dicken Absorber?

Seine indirekte Bandlücke macht die Absorption nahe der Bandkante phononenabhängig und vergleichsweise schwach. Lichtfangstrukturen verlängern den optischen Weg.

Ist ein starkes elektrisches Feld überall in der Zelle notwendig?

Nein. In quasineutralen Bereichen können Minoritätsladungsträger durch Diffusion gesammelt werden, wenn ihre Diffusionslänge ausreichend groß ist.

Warum kann eine leuchtende Solarzelle hochwertig sein?

Wenn nichtstrahlende Verluste gering sind, nähert sich die Rekombination dem strahlenden Grenzfall. Das verbessert die Quasi-Fermi-Aufspaltung und damit Voc.

Was unterscheidet EQE und Absorption?

Absorption zählt aufgenommene Photonen. EQE zählt Ladungsträger, die tatsächlich im äußeren Stromkreis erscheinen.

Warum sinkt der Wirkungsgrad bei hoher Temperatur?

Die Sperrsättigungsstromdichte steigt stark, wodurch Voc sinkt. Der kleine Gewinn an Photostrom kann diesen Spannungsverlust meist nicht ausgleichen.

Warum sind Tandemzellen effizienter?

Mehrere Bandlücken teilen das Spektrum. Hochenergetische Photonen verlieren weniger Energie durch Thermalisierung, langwellige Photonen können in einer kleineren Bandlücke absorbiert werden.

Warum zeigen Perowskitzellen manchmal Hysterese?

Mobile Ionen, Grenzflächenladungen, Fallenfüllung und zeitabhängige Polarisation können den Zustand während eines Spannungsscans verändern.

Ist der höchste anfängliche Wirkungsgrad automatisch die beste Technologie?

Nein. Energieertrag, Stabilität, Herstellung, Materialverfügbarkeit, Reparierbarkeit, Recycling und Modulzuverlässigkeit sind ebenso entscheidend.

69. Ausblick auf Teil 19

Teil 19 behandelt moderne Batterie- und Festelektrolytmaterialien. Interkalation, chemische Potentiale, Zellspannung, Mehrphasengleichgewichte, Diffusion in Elektroden, poröse Elektroden, Festelektrolyte, Grenzflächenwiderstände, SEI-Bildung, Dendriten, mechanische Spannungen, Degradation und Sicherheitsfragen werden systematisch entwickelt.

Literatur und weiterführende Quellen

  1. P. Würfel & U. Würfel: Physics of Solar Cells. Wiley-VCH.
  2. J. Nelson: The Physics of Solar Cells. Imperial College Press.
  3. M. A. Green: Solar Cells: Operating Principles, Technology and System Applications.
  4. A. Luque & S. Hegedus: Handbook of Photovoltaic Science and Engineering. Wiley.
  5. S. M. Sze & K. K. Ng: Physics of Semiconductor Devices. Wiley.

Didaktischer Hinweis: Beer–Lambert-Absorption, abrupte p-n-Übergänge, konstante Lebensdauern, ideale selektive Kontakte, Superposition von Dunkel- und Photostrom, Ein-Dioden-Modelle und detailliertes Gleichgewicht sind kontrollierte Modelle. Reale Solarzellen besitzen optische Interferenz, inhomogene Defekte, feld- und injektionsabhängige Rekombination, bewegliche Ionen, Kontaktbarrieren und gekoppelte thermische sowie mechanische Alterung.

Bild- und Nutzungsnachweis: Das Titelbild und die Abbildungen 1–3 wurden eigens für diesen Blogartikel mit selbst erstelltem Programmcode erzeugt. Es wurden keine Grafiken aus Lehrbüchern oder anderen Publikationen übernommen. Die Abbildungen sind für die Verwendung in diesem monetarisierbaren Blog vorgesehen.