Dévoilement de la Vraie Distribution des Durées de Vie des Nanocristaux de Silicium : Étude Combinée Spectroscopique Temps-Résolue et Fréquence-Résolue

Résumé

Les nanocristaux de silicium (SiNC) présentent une dynamique de luminescence complexe issue de taux radiatifs et non-radiatifs dépendants de la taille, d'élargissement inhomogène et d'interactions inter-particules. Ces facteurs génèrent des profils de décroissance non exponentiels, typiquement modélisés par une fonction exponentielle étirée (SE). Dans ce travail, nous dérivons d'abord la fonction de décroissance de population pour une décroissance de la forme « exp[−(t/τ)^β] », et comparons les distributions de durées de vie obtenues lorsque la luminescence ou la population suit cette fonction. Pour des valeurs β bien inférieures à 1, les deux approches donnent des durées de vie moyennes nettement différentes. Nous appliquons les deux modèles SE, une distribution log-normale et une décroissance bimoléculaire à des données de luminescence de deux ensembles de SiNC croîts thermiquement avec des tailles moyennes distinctes. Les durées de vie moyennes dépendent fortement de la configuration expérimentale et du modèle d'ajustement choisi, et aucun des modèles conventionnels ne capture adéquatement la dynamique de décroissance d'ensemble. Nous employons ensuite la spectroscopie fréquence-résolue (FRS) pour extraire la distribution de durées de vie directement, révélant une demi-largeur de ~0,5 décade et une distribution ressemblant à une fonction log-normale biaisée vers les hautes fréquences. L'utilisation combinée de techniques temps-résolues et fréquence-résolues fournit l'insight le plus fiable sur la dynamique de luminescence de nanocristaux à larges spectres d'émission.

Introduction

Les nanoparticules colloïdales sont intégrales à la catalyse, la biomédecine et l'optoélectronique. Les nanocristaux semiconducteurs (NC) sont particulièrement attractifs pour l'émission lumineuse, le photovoltaïque et la photocatalyse. Les nanocristaux de silicium (SiNC) ont attiré l'attention grâce à leur émission réglable, leur cœur de silicium abondant et biocompatible, et leur potentiel de synthèse scalable. Une connaissance précise de leurs propriétés optoélectroniques est essentielle pour le développement de dispositifs, et la spectroscopie temps-résolue est un outil principal à cet effet.

Les durées de vie de luminescence des SiNC sont couramment modélisées par une fonction exponentielle étirée (SE), exp[−(λt)^β], où 0 < β ≤ 1 est le paramètre de dispersion, λ le taux de décroissance, et t le temps. Cette forme implique une distribution asymétrique de taux de décroissance s'allongeant vers les durées de vie plus longues. Une fois β et λ obtenus par ajustement d'une courbe de décroissance, la distribution sous-jacente de taux peut être approximée.

L'origine physique de la décroissance SE dans le silicium et autres points quantiques est débattue depuis deux décennies, avec explications proposées incluant tunneling et piégeage de porteurs dans des ensembles densément empilés, couplage électron-phonon dépendant de la taille, et distribution de hauteurs de barrières non-radiatives. Une compréhension claire de la distribution de taux est donc essentielle pour élucider le mécanisme de luminescence dans les SiNC et systèmes apparentés.

De nombreuses études antérieures supposent a priori le modèle SE, négligeant souvent les distributions alternatives et la responsivité du détecteur sur le large spectre d'émission des SiNC. Cet oubli entrave la comparabilité entre études et obscurcit la vraie dynamique de décroissance. De plus, la spectroscopie fréquence-résolue (FRS) n'a pas été appliquée aux SiNC, bien qu'elle offre une voie sans modèle vers les distributions de durées de vie.

Le but de cet article est d'établir un cadre rigoureux pour mesurer, modéliser et interpréter la dynamique de luminescence des SiNC. Nous visons à résoudre les incohérences dans la littérature, fournir une comparaison plus fiable de différents jeux de données, et approfondir la compréhension des mécanismes sous-jacents.

Théorie de Base

Nous comparons trois modèles de décroissance : l'exponentielle étirée (SE), une distribution log-normale, et une décroissance bimoléculaire. Pour tout modèle, la fonction de densité de probabilité d'émission g(t) est liée à la fraction d'excitations restantes au temps t par

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

où c_t et c₀ sont les comptes de NC excités au temps t et initialement. Pour une recombinaison monomoléculaire (premier ordre), dc_t/dt = −λc_t, menant à c_t/c₀ = exp(−λt) et g(t) = λ exp(−λt). En incluant taux radiatif (λ_R) et non-radiatif (λ_NR), l'intensité mesurable devient g(t) = λ_R exp[−(λ_R+λ_NR)t].

Lorsque la population décroît comme une exponentielle étirée, c_t/c₀ = exp[−(λ_SE t)^β], la fonction de probabilité d'émission est

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

En utilisant les transformations de Laplace inverses, la distribution de taux correspondante H(λ) s'élargit et se biaise vers les hautes fréquences quand β diminue. Cependant, cette forme ne peut être décomposée en contributions radiatives et non-radiatives séparées, et n'est strictement normalisée que quand λ_NR = 0.

Pour le modèle de décroissance d'intensité g(t) = (λ_SE β/Γ(1/β)) exp[−(λ_SE t)^β], la décroissance de population est obtenue via

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

menant aux durées de vie moyennes

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

La distribution log-normale de taux de décroissance est définie comme

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

avec taux médian exp(μ), moyenne exp(μ+σ²/2), et durée de vie la plus probable exp(μ−σ²). L'intensité résultante est la transformée de Laplace de H(λ) :

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

La décroissance bimoléculaire suit dc_t/dt = −λc_t², donnant c_t/c₀ = 1/(λc₀t+1) et une intensité en loi de puissance I_t/I₀ = A λc₀/(λc₀+1)². Ce modèle n'a pas de durée de vie moyenne finie.

La Figure 1 illustre les différences entre les deux formulations SE et leurs distributions de taux correspondantes.

Spectroscopie Fréquence-Résolue en Quadrature

La FRS en quadrature (QFRS) utilise un faisceau pompe modulé sinusoïdalement à la fréquence angulaire ω et enregistre la phase et l'amplitude de la réponse de photoluminescence (PL). La composante de quadrature Q_PL = Z_PL sin(Δθ_PL) fournit une mesure directe de la distribution de durées de vie après normalisation par l'amplitude de modulation du laser Z_LA.

Le signal QFRS pour une décroissance mono-exponentielle est donné par

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

avec τ₀ = 1/ω₀. Puisqu'une distribution de décroissance réaliste s'étend sur plusieurs décennies, une déconvolution est requise pour extraire la vraie distribution.

Résultats et Discussion

Caractérisation de Base

Des images MET en champ clair ont servi à estimer les diamètres des particules par comptage de pixels. Les distributions de tailles résultantes étaient bien décrites par des fonctions log-normales avec diamètres moyens de 2,9 nm (σ=0,1555) pour l'échantillon recuit à 1100 °C et 5,4 nm (σ=0,1917) pour l'échantillon à 1200 °C. La MET haute résolution a confirmé les franges de réseau et les estimations de taille. FTIR et XPS ont vérifié la fonctionnalisation dodécène réussie ; les particules plus petites présentaient une oxydation plus grande et une densité de fonctionnalisation plus faible.

Photoluminescence et Spectroscopie Temps-Résolue

Les spectres PL culminent à 660 nm (LHM = 123 nm) pour les petits NC et 825 nm (LHM = 198 nm) pour les grands NC. En utilisant la formule de bande interdépendante de la taille E_g = √{E_{g,bulk}² + D/R²} avec D = 4,8 eV²/nm², nous avons obtenu des bandes interdites prédites de 1,87 eV (petits) et 1,37 eV (grands), cohérentes avec les pics PL mesurés. Les rendements quantiques absolus (AQY) étaient de 12 % (petits) et 56 % (grands) ; des mesures indépendantes ont donné 18 % et 48 %, respectivement, dans la plage d'incertitude typique de 10 %.

Tous les échantillons ont présenté des profils de décroissance non exponentiels. Des ajustements utilisant les équations 5, 9, 16 et 17 (SE, log-normale, bimoléculaire) ont été effectués par minimisation des moindres carrés. Les résidus oscillaient pour tous les modèles, indiquant une capture incomplète de la forme de décroissance ; le SE simple (Éq. 9) et la log-normale (Éq. 16) ont fourni les plus faibles sommes des carrés. Les durées de vie moyennes variaient fortement selon le modèle choisi (Tableau 2). La méthode de Higashi-Kastner, qui localise le pic de I_t·t versus t, a donné des constantes de décroissance similaires à (1/β)^{1/β}·τ_SE dérivée de l'Éq. 9. Le modèle bimoléculaire a mal décrit les données, cohérent avec des NC isolés ne subissant pas de densités d'excitons élevées.

Les conditions d'excitation ont été estimées en utilisant des sections efficaces d'absorption (~10⁻¹⁴ cm²) et une irradiance de 4500 W m⁻² à 352 nm. Les excitations moyennes par NC étaient < 1 pour les grands NC et ≈ 0,2 pour les petits NC, indiquant que les effets de sur-excitation étaient négligeables, confirmé par des mesures de durée de vie en fonction de la puissance montrant < 2 % de variation.

Pour sonder la dépendance en longueur d'onde, les décroissances PL ont été mesurées à longueurs d'onde d'émission discrètes (bande passante 3 nm) via un monochromateur. Le paramètre de dispersion β tendait vers l'unité pour les longueurs d'onde plus longues, et les durées de vie augmentaient avec la longueur d'onde (Fig. 5, Tableau 3). Les particules plus courtes présentaient systématiquement des durées de vie plus courtes que les plus grandes à la même longueur d'onde, reflétant un AQY plus faible et des contributions non-radiatives de surface plus élevées.

Spectroscopie Fréquence-Résolue

La FRS a été validée utilisant un circuit RC connu (décroissance mono-exponentielle, τ = 78,9 µs) et une microsphère dopée chélate d'europium (τ ≈ 670 µs). Les spectres QFRS correspondaient à la fonction de réponse attendue (Éq. 18) avec des fréquences de pic de 12,7 kHz et 1570 Hz, respectivement.

Les données FRS des SiNC n'étaient que légèrement plus larges que la réponse intrinsèque, nécessitant une déconvolution via l'algorithme de Richardson-Lucy pour récupérer la distribution de durées de vie. Les distributions déconvolvées (Fig. 6) s'étendaient sur ~0,5 décade, avec les petits NC montrant une forme symétrique à l'échelle log et un pic à 19,9 kHz (50,3 µs), et les grands NC culminant à 6,28 kHz (159,2 µs). Les ajustements modèles SE et log-normale donnaient des formes de distribution et fréquences de pic différentes, toutes inférieures aux pics QFRS, reflétant la sensibilité de la TRS d'état stationnaire à la durée d'impulsion et à la responsivité du détecteur.

La responsivité du détecteur affectait significativement les résultats TRS ; les mesures APD et PMT différaient d'un facteur ~2 en temps de décroissance moyens. L'APD culminait à 600 nm, le PMT à 850 nm, menant à des ajustements « meilleurs » divergents. Corriger la responsivité en FRS est faisable si la dépendance spectrale de la distribution de taux de décroissance est connue ; la TRS manque de méthode de correction directe.

Conclusions

Nous avons dérivé la décroissance de population correspondant à la décroissance de luminescence SE couramment utilisée et fourni des expressions pour les durées de vie moyennes. Deux échantillons de SiNC croîts thermiquement et fonctionnalisés dodécène (diamètres moyens 2,9 nm et 5,4 nm) ont été synthétisés et caractérisés. Les données TRS ajustées avec modèles SE, log-normale et bimoléculaire n'ont pas pleinement capturé la forme de décroissance ; les durées de vie moyennes dépendaient de la responsivité du détecteur et du choix du modèle. La spectroscopie fréquence-résolue, après déconvolution, a fourni la distribution de durées de vie directement, révélant un étalement de ~0,5 décade et un biais vers les hautes fréquences pour les grands NC. La FRS est donc la méthode préférée pour sonder la dynamique de luminescence des SiNC, bien que déconvolution soignée et calibration du détecteur soient essentielles. Les futures études devraient incorporer des corrections de responsivité du détecteur et éviter de recourir par défaut aux modèles SE sans valider la distribution sous-jacente.

Méthodes

Les SiNC ont été préparés par recuit de 4 g de silsesquioxane d'hydrogène (HSQ) à 1100 °C ou 1200 °C pendant 1 h sous 5 % H₂/95 % Ar. Les composites SiNC/SiO₂ résultants ont été broyés, attaqués par HF pour retirer la matrice silice, et fonctionnalisés avec du 1-dodécène (≈ 4,6 mmol) et AIBN (20 mg) sous Ar. Après 24 h à 70 °C, le mélange réactionnel a été précipité avec méthanol/éthanol, centrifugé, et redispersé dans le toluène. Les SiNC finaux terminés dodécyle étaient stockés dans 5 mL de toluène sec (≈ 0,5 mg mL⁻¹).

Les échantillons MET ont été préparés sur grilles cuivre revêtues de carbone ultramince (≈ 3 nm) et imageries avec un JEOL JEM-2010 (champ clair) et JEOL JEM-ARM200CF (HRTEM). Les spectres FTIR ont été enregistrés sur un Nicolet 8700, et le XPS réalisé sur un système SPECS (source Mg Kα).

L'excitation PL a utilisé un laser Ar⁺ à 352 nm pulsé (cycle de travail 50 %, 50–250 Hz) via un modulateur acousto-optique (AOM). Le faisceau était divisé par un séparateur de faisceau ; la partie principale illuminait la cuvette d'échantillon, et la PL émise était collectée par une fibre 0,22 NA, filtrée par un filtre passe-haut 450 nm, et dirigée vers le détecteur. Le spectre PL a été enregistré sur un spectromètre miniature Ocean Optics, calibré avec une source HL-3+-CAL. Le rendement quantique a été mesuré avec une sphère intégrante (excitation 405 nm) sur une solution diluée (A ≈ 0,15).

La PL temps-résolue a été enregistrée avec soit une photodiode avalanche Thorlabs 120A2 (APD) soit un tube photomultiplicateur Hamamatsu h7422-50 (PMT), interfacés à un oscilloscope Moku:Lab ou un multiscaler Becker-Hickl PMS400. Les fonctions de réponse des détecteurs ont été mesurées et appliquées lors de l'analyse. Les expériences QFRS ont employé un AOM modulé sinusoïdalement ; le signal de référence provenait d'une photodiode Thorlabs PDA10A. La réponse PL a été collectée par l'APD et traitée avec un amplificateur lock-in pour extraire les composantes en phase et en quadrature.

Pour la détection de composantes de décroissance rapide, un laser diode picoseconde 405 nm (Alphalas GmbH) a excité les NC tandis qu'un PMT Becker-Hickl HPM-100-50 (réponse 100 ps) et un compteur d'impulsions SPC-130 enregistraient la PL. Aucune décroissance nanoseconde n'a été observée.

Abréviations

APD

Photodiode avalanche

AQY

Rendement quantique absolu

FRS

Spectroscopie fréquence-résolue

LN

Log-normale

NCs

Nanocristaux

PL

Photoluminescence

PMT

Tube photomultiplicateur

QFRS

Spectroscopie fréquence-résolue en quadrature

SE

Exponentielle étirée

SiNCs

Nanocristaux de silicium

TRS

Spectroscopie temps-résolue