Entschleierung der wahren Lebensdauerverteilung von Silicium-Nanokristallen: Kombinierte zeit- und frequenzaufgelöste spektroskopische Studie

Zusammenfassung

Silicium-Nanokristalle (SiNCs) zeigen komplexe Lumineszenzdynamiken, die aus größenabhängigen radiativen und nicht-radiativen Raten, inhomogener Verbreiterung und Teilchenwechselwirkungen resultieren. Diese Faktoren erzeugen nicht-exponentielle Zerfallsprofile, die typischerweise durch eine gestreckte Exponentialfunktion (SE) modelliert werden. In dieser Arbeit leiten wir zunächst die Populationszerfallsfunktion für einen Zerfall der Form „exp[−(t/τ)^{β}]“ ab und vergleichen die Lebensdauerverteilungen, die sich ergeben, wenn der Lumineszenzzerfall oder der Populationszerfall dieser Funktion folgt. Für β-Werte weit unter 1 liefern die beiden Ansätze deutlich unterschiedliche mittlere Lebensdauern. Wir wenden beide SE-Modelle, eine lognormale Verteilung und einen bimolekularen Zerfall auf Lumineszenzdaten von zwei thermisch gewachsenen SiNC-Ensembles mit unterschiedlichen mittleren Größen an. Die mittleren Lebensdauern hängen stark von der experimentellen Konfiguration und dem gewählten Anpassungsmodell ab, und keines der konventionellen Modelle erfasst die Ensemble-Zerfallsdynamik adäquat. Wir setzen dann frequenzaufgelöste Spektroskopie (FRS) ein, um die Lebensdauerverteilung direkt zu extrahieren, was eine Halbwertsbreite von ~0,5 Dekaden und eine einer hochfrequenz-verzerrten lognormalen Funktion ähnelnde Verteilung offenbart. Die kombinierte Nutzung zeit- und frequenzaufgelöster Techniken liefert die zuverlässigsten Einblicke in die Lumineszenzdynamik von Nanokristallen mit breiten Emissionsspektren.

Einführung

Kolloidale Nanopartikel sind integraler Bestandteil von Katalyse, Biomedizin und Optoelektronik [1–4]. Halbleiter-Nanokristalle (NCs) sind besonders attraktiv für Lichtemission, Photovoltaik und Photokatalyse [5,6]. Silicium-Nanokristalle (SiNCs) haben aufgrund ihrer einstellbaren Emission, des reichlich vorhandenen und biokompatiblen Siliciumkerns und des Potenzials für skalierbare Synthese Aufmerksamkeit erregt [7–9]. Genaues Wissen über ihre optoelektronischen Eigenschaften ist für die Geräteentwicklung essenziell, und zeitaufgelöste Spektroskopie ist ein primäres Werkzeug dafür.

SiNC-Lumineszenzlebensdauern werden häufig durch eine gestreckte Exponentialfunktion (SE) modelliert, exp[−(λt)^{β}], wobei 0 < β ≤ 1 der Dispersionsparameter, λ die Zerfallsrate und t die Zeit ist. Diese Form impliziert eine asymmetrische Verteilung von Zerfallsraten, die zu längeren Lebensdauern hin abfällt. Einmal aus der Anpassung einer Zerfallskurve β und λ erhalten, kann die zugrundeliegende Ratenverteilung approximiert werden [9].

Der physikalische Ursprung des SE-Zerfalls in Silicium und anderen Quantenpunkten wird seit zwei Jahrzehnten diskutiert, mit vorgeschlagenen Erklärungen einschließlich Ladungsträgertunneln und -fang in dicht gepackten Ensembles [11], größenabhängiger Elektron-Phonon-Kopplung [10] und einer Verteilung nicht-radiativer Barrierehöhen [13,14]. Ein klares Verständnis der Ratenverteilung ist somit essenziell für die Aufklärung des Lumineszenzmechanismus in SiNCs und verwandten Systemen.

Viele frühere Studien nehmen das SE-Modell a priori an, oft vernachlässigend alternative Verteilungen und die Detektoreigenschaften über das breite SiNC-Emissionsspektrum. Dieser Überblick erschwert die Vergleichbarkeit zwischen Studien und verschleiert die wahren Zerfallsdynamiken. Darüber hinaus wurde frequenzaufgelöste Spektroskopie (FRS) noch nicht auf SiNCs angewendet, obwohl sie einen modellfreien Weg zu Lebensdauerverteilungen bietet.

Ziel dieser Arbeit ist es, einen rigorosen Rahmen für die Messung, Modellierung und Interpretation von SiNC-Lumineszenzdynamiken zu etablieren. Wir wollen Inkonsistenzen in der Literatur auflösen, einen zuverlässigeren Vergleich verschiedener Datensätze ermöglichen und das Verständnis der zugrundeliegenden Mechanismen vertiefen.

Grundlegende Theorie

Wir vergleichen drei Zerfallsmodelle: die gestreckte Exponentialfunktion (SE), eine lognormale Verteilung und einen bimolekularen Zerfall. Für jedes Modell ist die Emissionswahrscheinlichkeitsdichtefunktion g(t) mit dem Anteil der bei Zeit t verbleibenden Anregungen verknüpft durch

∫₀ᵗ g(t′) dt = 1 − c_t/c₀,

wobei c_t und c₀ die Anzahlen angeregter NCs zum Zeitpunkt t und initial sind. Für eine monomolekulare Rekombination erster Ordnung gilt dc_t/dt = −λc_t, was zu c_t/c₀ = exp(−λt) und g(t) = λ exp(−λt) führt. Einbeziehung sowohl radiativer (λ_R) als auch nicht-radiativer (λ_NR) Raten ergibt die messbare Intensität g(t) = λ_R exp[−(λ_R+λ_NR)t].

Wenn die Population als gestreckte Exponentialfunktion zerfällt, c_t/c₀ = exp[−(λ_SE t)^{β}], ist die Emissionswahrscheinlichkeitsfunktion

g(t) = β λ_SE^{β} t^{β−1} exp[−(λ_SE t)^{β}].

Mittels inverser Laplace-Transformationen erweitert und verzerrt sich die entsprechende Ratenverteilung H(λ) zu hohen Frequenzen hin, je kleiner β wird. Diese Form kann jedoch nicht in separate radiative und nicht-radiative Beiträge zerlegt werden und ist nur strikt normalisiert, wenn λ_NR = 0.

Für das Intensitätszerfallsmodell g(t) = (λ_SE β/Γ(1/β)) exp[−(λ_SE t)^{β}] wird der Populationszerfall über

c_t/c₀ = 1/Γ(1/β) Γ[1/β, (λ_SE t)^{β}],

erhalten, was zu mittleren Lebensdauern führt

⟨τ_SE⟩ = τ_SE Γ(2/β)/Γ(1/β), ⟨t⟩ = τ_SE Γ(3/β)/(2Γ(2/β)).

Die lognormale Zerfallsratenverteilung ist definiert als

H(λ) = (1/λ) (1/(σ√{2π})) exp[−(ln λ − μ)²/(2σ²)],

mit Medianrate exp(μ), Mittelwert exp(μ+σ²/2) und wahrscheinlichster Lebensdauer exp(μ−σ²). Die resultierende Intensität ist die Laplace-Transformation von H(λ):

I_t = A ∫₀^∞ H(λ) exp(−λt) dλ + dc.

Der bimolekulare Zerfall folgt dc_t/dt = −λc_t², was c_t/c₀ = 1/(λc₀t+1) und ein Potenzgesetz für die Intensität I_t/I₀ = A λc₀/(λc₀+1)² liefert. Dieses Modell hat keine endliche mittlere Lebensdauer.

Abbildung 1 illustriert die Unterschiede zwischen den beiden SE-Formulierungen und ihren entsprechenden Ratenverteilungen.

Quadratur-Frequenzaufgelöste Spektroskopie

Quadratur-FRS (QFRS) verwendet einen sinusförmig modulierten Pumpstrahl bei Winkelgeschwindigkeit ω und zeichnet Phase und Amplitude der Photolumineszenz (PL)-Antwort auf. Die Quadraturkomponente Q_PL = Z_PL sin(Δθ_PL) bietet ein direktes Maß für die Lebensdauerverteilung nach Normalisierung durch die Lasermodulationsamplitude Z_LA.

Das QFRS-Signal für einen einfachen exponentiellen Zerfall ist gegeben durch

S_{log10 r} = (ω τ₀)/(1+ω² τ₀²),

mit τ₀ = 1/ω₀. Da eine realistische Zerfallsverteilung mehrere Dekaden umfasst, ist eine Dekonvolution erforderlich, um die wahre Verteilung zu extrahieren.

Ergebnisse und Diskussion

Grundlegende Charakterisierung

Hellfeld-TEM-Bilder wurden verwendet, um Teilchendurchmesser durch Pixelzählung zu schätzen. Die resultierenden Größenverteilungen wurden gut durch lognormale Funktionen mit mittleren Durchmessern von 2,9 nm (σ=0,1555) für die bei 1100 °C geglühte Probe und 5,4 nm (σ=0,1917) für die 1200 °C-Probe beschrieben. Hochauflösendes TEM bestätigte Gitterstreifen und Größenschätzungen. FTIR und XPS verifizierten erfolgreiche Dodecen-Funktionalisierung; die kleineren Teilchen zeigten stärkere Oxidation und geringere Funktionalisierungsdichte.

Photolumineszenz und zeitaufgelöste Spektroskopie

Die PL-Spektren peaked bei 660 nm (FWHM = 123 nm) für kleine NCs und 825 nm (FWHM = 198 nm) für große NCs. Mit der größenabhängigen Bandlückenformel E_g = √{E_{g,bulk}² + D/R²} mit D = 4,8 eV²/nm² erhielten wir vorhergesagte Bandlücken von 1,87 eV (klein) und 1,37 eV (groß), konsistent mit den gemessenen PL-Peaks. Absolute Quantenausbeuten (AQY) betrugen 12 % (klein) und 56 % (groß); unabhängige Messungen ergaben 18 % und 48 %, jeweils innerhalb der typischen 10 % Unsicherheit.

Alle Proben zeigten nicht-exponentielle Zerfallsprofile. Anpassungen mittels Gleichungen 5, 9, 16 und 17 (SE, lognormal, bimolekular) wurden über Kleinste-Quadrate-Minimierung durchgeführt. Residuen oszillierten für alle Modelle, was auf unvollständige Erfassung der Zerfallsform hinweist; das einfache SE (Gl. 9) und lognormal (Gl. 16) lieferten die geringsten Summenquadrate. Mittlere Lebensdauern variierten stark mit dem gewählten Modell (Tabelle 2). Die Higashi-Kastner-Methode, die das Maximum von I_t·t versus t lokalisiert, ergab Zerfallskonstanten ähnlich zu (1/β)^{1/β}·τ_SE aus Gl. 9. Das bimolekulare Modell beschrieb die Daten schlecht, konsistent mit isolierten NCs, die keine hohen Excitondichten erfahren.

Anregungsbedingungen wurden mit Absorptionsquerschnitten (~10⁻¹⁴ cm²) und einer Bestrahlungsstärke von 4500 W m⁻² bei 352 nm geschätzt. Die durchschnittlichen Anregungen pro NC lagen < 1 für große NCs und ≈ 0,2 für kleine NCs, was vernachlässigbare Überanregungseffekte anzeigt, bestätigt durch leistungsabhängige Lebensdauermessungen mit < 2 % Variation.

Um Wellenlängenabhängigkeit zu untersuchen, wurden PL-Zerfälle an diskreten Emissionswellenlängen (3 nm Bandbreite) mit einem Monochromator gemessen. Der Dispersionsparameter β näherte sich für längere Wellenlängen der Eins an, und Lebensdauern nahmen mit der Wellenlänge zu (Abb. 5, Tabelle 3). Kleinere Teilchen zeigten konsistent kürzere Lebensdauern als größere bei gleicher Wellenlänge, was niedrigere AQY und höhere nicht-radiative Oberflächenbeiträge reflektiert.

Frequenzaufgelöste Spektroskopie

FRS wurde mit einem bekannten RC-Kreis (mono-exponentieller Zerfall, τ = 78,9 µs) und einem Eu-Chelat-dotierten Mikrosphäre (τ ≈ 670 µs) validiert. Die QFRS-Spektren stimmten mit der erwarteten Antwortfunktion (Gl. 18) mit Spitzenfrequenzen von 12,7 kHz und 1570 Hz überein.

SiNC-FRS-Daten waren nur geringfügig breiter als die intrinsische Antwort, was eine Dekonvolution via Richardson-Lucy-Algorithmus zur Extraktion der Lebensdauerverteilung erforderte. Die dekonvolvierten Verteilungen (Abb. 6) erstreckten sich über ~0,5 Dekaden, wobei die kleinen NCs eine symmetrische Log-Skalen-Form mit Peak bei 19,9 kHz (50,3 µs) und die großen NCs einen Peak bei 6,28 kHz (159,2 µs) zeigten. SE- und lognormale Modellanpassungen ergaben unterschiedliche Verteilungsformen und Spitzenfrequenzen, alle niedriger als die QFRS-Peaks, was die Empfindlichkeit von stationärer TRS gegenüber Pulsdauer und Detektoreigenschaften widerspiegelt.

Detektoreigenschaften beeinflussten TRS-Ergebnisse signifikant; APD- und PMT-Messungen unterschieden sich um einen Faktor ~2 in mittleren Zerfallszeiten. Die APD peaked bei 600 nm, die PMT bei 850 nm, was zu divergenten „besten“ Modellanpassungen führte. Eine Korrektur der Eigenschaften in FRS ist machbar, wenn die spektrale Abhängigkeit der Zerfallsratenverteilung bekannt ist; TRS bietet keine einfache Korrekturmethode.

Fazit

Wir leiteten den Populationszerfall ab, der der häufig verwendeten SE-Lumineszenzzerfallsfunktion entspricht, und gaben Ausdrücke für mittlere Lebensdauern an. Zwei thermisch gewachsene, dodecen-funktionalisierte SiNC-Proben (mittlere Durchmesser 2,9 nm und 5,4 nm) wurden synthetisiert und charakterisiert. TRS-Daten, angepasst mit SE, lognormalen und bimolekularen Modellen, erfassten die Zerfallsform nicht vollständig; mittlere Lebensdauern hingen von Detektoreigenschaften und Modellwahl ab. Frequenzaufgelöste Spektroskopie lieferte nach Dekonvolution die Lebensdauerverteilung direkt, offenbarte eine ~0,5-Dekaden-Spreizung und eine Hochfrequenz-Verzerrung für große NCs. FRS ist daher die bevorzugte Methode zur Untersuchung von SiNC-Lumineszenzdynamiken, wenngleich sorgfältige Dekonvolution und Detektorkalibrierung essenziell sind. Zukünftige Studien sollten Detektoreigenschaftskorrekturen einbeziehen und vermeiden, standardmäßig SE-Modelle ohne Validierung der zugrundeliegenden Verteilung zu verwenden.

Methoden

SiNCs wurden durch Glühen von 4 g Hydridsilsesquioxan (HSQ) bei 1100 °C oder 1200 °C für 1 h in 5 % H₂/95 % Ar hergestellt. Die resultierenden SiNC/SiO₂-Komposite wurden gemahlen, HF-geätzt zur Entfernung der Siliciumdioxid-Matrix und mit 1-Dodecen (≈ 4,6 mmol) und AIBN (20 mg) unter Ar funktionalisiert. Nach 24 h bei 70 °C wurde das Reaktionsgemisch mit Methanol/Ethanol gefällt, zentrifugiert und in Toluol redispersiert. Die endgültigen dodecyl-terminierten SiNCs wurden in 5 ml trockenem Toluol (≈ 0,5 mg ml⁻¹) gelagert.

TEM-Proben wurden auf ultradünnen (≈ 3 nm) kohlenstoffbeschichteten Kupfergittern präpariert und mit einem JEOL JEM-2010 (Hellfeld) und JEOL JEM-ARM200CF (HRTEM) abgebildet. FTIR-Spektren wurden an einem Nicolet 8700 aufgenommen, XPS an einem SPECS-System (Mg Kα-Quelle).

PL-Anregung erfolgte mit einem 352-nm Ar⁺-Laser, gepulst (50 % Tastverhältnis, 50–250 Hz) über einen akusto-optischen Modulator (AOM). Der Strahl wurde durch einen Strahlteiler geteilt; der Hauptteil beleuchtete die Probenküvette, und die emittierte PL wurde von einer 0,22 NA-Faser gesammelt, durch ein 450-nm Langpassfilter gefiltert und zum Detektor geleitet. Das PL-Spektrum wurde an einem Ocean Optics Miniaturspektrometer aufgezeichnet, kalibriert mit einer HL-3+-CAL-Quelle. Quantenausbeute wurde mit einer Integriersphäre (405-nm Anregung) an einer verdünnten Lösung (A ≈ 0,15) gemessen.

Zeitaufgelöste PL wurde entweder mit einer Thorlabs 120A2 Lawinenphotodiode (APD) oder einer Hamamatsu h7422-50 Photomultiplierröhre (PMT) aufgezeichnet, verbunden mit einem Moku:Lab-Oszilloskop oder einem Becker-Hickl PMS400 Multiscaler. Detektoreigenschaften wurden gemessen und während der Analyse angewendet. QFRS-Experimente nutzten einen sinusförmig modulierten AOM; das Referenzsignal stammte von einer Thorlabs PDA10A Photodiode. Die PL-Antwort wurde von der APD gesammelt und mit einem Lock-in-Verstärker verarbeitet, um In-Phase- und Quadraturkomponenten zu extrahieren.

Für die Detektion schneller Zerfallskomponenten erregte ein 405-nm Pikosekunden-Diodenlaser (Alphalas GmbH) die NCs, während ein Becker-Hickl HPM-100-50 PMT (100 ps Ansprechzeit) und ein SPC-130 Pulszähler die PL aufzeichneten. Kein Nanosekunden-Zerfall wurde beobachtet.

Abkürzungen

APD

Lawinenphotodiode (Avalanche photodiode)

AQY

Absolute Quantenausbeute (Absolute quantum yield)

FRS

Frequenzaufgelöste Spektroskopie (Frequency-resolved spectroscopy)

LN

Lognormal

NCs

Nanokristalle (Nanocrystals)

PL

Photolumineszenz (Photoluminescence)

PMT

Photomultiplierröhre (Photomultiplier tube)

QFRS

Quadratur-frequenzaufgelöste Spektroskopie (Quadrature frequency-resolved spectroscopy)

SE

Gestreckte Exponentialfunktion (Stretched exponential)

SiNCs

Silicium-Nanokristalle (Silicon nanocrystals)

TRS

Zeitaufgelöste Spektroskopie (Time-resolved spectroscopy)