Descubriendo la verdadera distribución de vida media de nanocristales de silicio: estudio combinado espectroscópico resuelto en tiempo y frecuencia

Resumen

Los nanocristales de silicio (SiNC) exhiben dinámicas de luminiscencia complejas que surgen de tasas radiativas y no radiativas dependientes del tamaño, ensanchamiento inhomogéneo e interacciones inter-partícula. Estos factores generan perfiles de decaimiento no exponenciales, típicamente modelados por función exponencial estirada (SE). En este trabajo, primero derivamos la función de decaimiento de población para un decaimiento de la forma “exp[−(t/τ)^{β}]”, y comparamos las distribuciones de vida media obtenidas cuando la luminiscencia o la población siguen esta función. Para valores β bien por debajo de 1, ambos enfoques producen vidas medias marcadamente diferentes. Aplicamos ambos modelos SE, una distribución lognormal y un decaimiento bimolecular a datos de luminiscencia de dos conjuntos de SiNC crecidos térmicamente con tamaños medios distintos. Las vidas medias dependen fuertemente de la configuración experimental y el modelo de ajuste elegido, y ninguno de los modelos convencionales captura adecuadamente la dinámica de decaimiento del conjunto. Luego empleamos espectroscopía resuelta en frecuencia (FRS) para extraer la distribución de vida media directamente, revelando un semi-ancho de ~0,5 décadas y una distribución similar a una función lognormal sesgada hacia altas frecuencias. El uso combinado de técnicas resueltas en tiempo y frecuencia proporciona la visión más fiable sobre la dinámica de luminiscencia de nanocristales con amplios espectros de emisión.

Introducción

Las nanopartículas coloidales son integrales en catálisis, biomedicina y optoelectrónica [1–4]. Los nanocristales semiconductores (NC) son especialmente atractivos para emisión de luz, fotovoltaica y fotocatálisis [5,6]. Los nanocristales de silicio (SiNC) han atraído atención por su emisión sintonizable, núcleo de silicio abundante y biocompatible, y potencial para síntesis escalable [7–9]. El conocimiento preciso de sus propiedades optoelectrónicas es esencial para desarrollo de dispositivos, y la espectroscopía resuelta en tiempo es herramienta principal para este propósito.

Las vidas medias de luminiscencia de SiNC se modelan comúnmente por función exponencial estirada (SE), exp[−(λt)^{β}]”, donde 0 < β ≤ 1 es parámetro de dispersión, λ es tasa de decaimiento, y t es tiempo. Esta forma implica distribución asimétrica de tasas de decaimiento que cola hacia vidas medias más largas. Una vez obtenidos β y λ ajustando curva de decaimiento, la distribución subyacente de tasas puede aproximarse [9].

El origen físico del decaimiento SE en silicio y otros puntos cuánticos se debate desde hace dos décadas, con explicaciones propuestas incluyendo tunelado y atrapamiento de portadores en conjuntos densamente empaquetados [11], acoplamiento electrón-fonón dependiente del tamaño [10], y distribución de alturas de barrera no radiativas [13,14]. Una comprensión clara de la distribución de tasas es esencial para elucidar el mecanismo de luminiscencia en SiNC y sistemas relacionados.

Muchos estudios previos asumen a priori el modelo SE, a menudo descuidando distribuciones alternativas y la responsividad del detector sobre el amplio espectro de emisión de SiNC. Esta omisión dificulta la comparabilidad entre estudios y oculta la verdadera dinámica de decaimiento. Además, la espectroscopía resuelta en frecuencia (FRS) no se ha aplicado a SiNC, aunque ofrece ruta libre de modelo a distribuciones de vida media.

El objetivo de este artículo es establecer marco riguroso para medir, modelar e interpretar dinámicas de luminiscencia de SiNC. Buscamos resolver inconsistencias en literatura, proporcionar comparación más fiable de diferentes conjuntos de datos, y profundizar comprensión de mecanismos subyacentes.

Teoría Básica

Comparamos tres modelos de decaimiento: exponencial estirada (SE), distribución lognormal, y decaimiento bimolecular. Para cualquier modelo, la función de densidad de probabilidad de emisión g(t) se relaciona con la fracción de excitaciones restantes en tiempo t por

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

donde c_t y c₀ son conteos de NC excitados en tiempo t e inicialmente. Para recombinación de primer orden (monomolecular), dc_t/dt = −λc_t, resultando c_t/c₀ = exp(−λt) y g(t) = λ exp(−λt). Incluyendo tasas radiativas (λ_R) y no radiativas (λ_NR), la intensidad medible se vuelve g(t) = λ_R exp[−(λ_R+λ_NR)t].

Cuando la población decae como exponencial estirada, c_t/c₀ = exp[−(λ_SE t)^{β}], la función de probabilidad de emisión es

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

Usando transformadas inversas de Laplace, la distribución de tasas correspondiente H(λ) se ensancha y sesga hacia altas frecuencias al disminuir β. Sin embargo, esta forma no puede descomponerse en contribuciones radiativas y no radiativas separadas, y es estrictamente normalizada solo cuando λ_NR = 0.

Para el modelo de decaimiento de intensidad g(t) = (λ_SE β/Γ(1/β)) exp[−(λ_SE t)^{β}], el decaimiento de población se obtiene vía

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

llevando a vidas medias

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

La distribución lognormal de tasas de decaimiento se define como

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

con tasa media exp(μ), media exp(μ+σ²/2), y vida media más probable exp(μ−σ²). La intensidad resultante es la transformada de Laplace de H(λ):

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

El decaimiento bimolecular sigue dc_t/dt = −λc_t², resultando c_t/c₀ = 1/(λc₀t+1) y ley de potencia I_t/I₀ = A λc₀/(λc₀+1)². Este modelo no tiene vida media finita.

La Figura 1 ilustra diferencias entre las dos formulaciones SE y sus correspondientes distribuciones de tasas.

Espectroscopía Resuelta en Frecuencia en Cuadratura

La FRS en cuadratura (QFRS) emplea haz de bombeo modulado sinusoidalmente a frecuencia angular ω y registra fase y amplitud de respuesta de fotoluminiscencia (PL). El componente de cuadratura Q_PL = Z_PL sin(Δθ_PL) proporciona medida directa de distribución de vida media tras normalizar por amplitud de modulación láser Z_LA.

La señal QFRS para decaimiento exponencial simple es

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

con τ₀ = 1/ω₀. Dado que una distribución realista de decaimiento abarca varias décadas, se requiere deconvolución para extraer la distribución verdadera.

Resultados y Discusión

Caracterización Básica

Imágenes TEM de campo claro estimaron diámetros de partícula por conteo de píxeles. Las distribuciones de tamaño resultantes se describieron bien por funciones lognormales con diámetros medios de 2,9 nm (σ=0,1555) para muestra recocida a 1100 °C y 5,4 nm (σ=0,1917) para 1200 °C. TEM de alta resolución confirmó franjas de red y estimaciones de tamaño. FTIR y XPS verificaron funcionalización exitosa con dodeceno; partículas menores exhibieron mayor oxidación y menor densidad de funcionalización.

Fotoluminiscencia y Espectroscopía Resuelta en Tiempo

Los espectros PL picaron a 660 nm (FWHM = 123 nm) para NC pequeños y 825 nm (FWHM = 198 nm) para grandes. Usando fórmula de banda prohibida dependiente del tamaño E_g = √{E_{g,bulk}² + D/R²} con D = 4,8 eV²/nm², obtuvimos bandas prohibidas predichas de 1,87 eV (pequeños) y 1,37 eV (grandes), consistentes con picos PL medidos. Rendimientos cuánticos absolutos (AQY) fueron 12 % (pequeños) y 56 % (grandes); mediciones independientes dieron 18 % y 48 %, respectivamente, dentro del rango típico de incertidumbre 10 %.

Todas las muestras exhibieron perfiles de decaimiento no exponenciales. Ajustes usando Eqs. 5, 9, 16 y 17 (SE, lognormal, bimolecular) se realizaron mediante minimización de mínimos cuadrados. Residuos oscilaron para todos los modelos, indicando captura incompleta de forma de decaimiento; el SE simple (Eq. 9) y lognormal (Eq. 16) proveyeron menores sumas de cuadrados. Las vidas medias variaron fuertemente con modelo elegido (Tabla 2). El método Higashi-Kastner, que localiza pico de I_t·t versus t, produjo constantes de decaimiento similares a (1/β)^{1/β}·τ_SE derivado de Eq. 9. El modelo bimolecular describió pobremente los datos, consistente con NC aislados sin altas densidades de excitones.

Condiciones de excitación se estimaron usando secciones transversales de absorción (~10⁻¹⁴ cm²) e irradiancia de 4500 W m⁻² a 352 nm. Excitaciones promedio por NC fueron < 1 para NC grandes y ≈ 0,2 para pequeños, indicando efectos de sobre-excitación despreciables, confirmado por mediciones de vida media dependientes de potencia mostrando < 2 % variación.

Para sondear dependencia de longitud de onda, decaimientos PL se midieron en longitudes de onda de emisión discretas (banda 3 nm) usando monocromador. El parámetro de dispersión β se acercó a unidad para longitudes de onda mayores, y vidas medias aumentaron con longitud de onda (Fig. 5, Tabla 3). Partículas menores exhibieron consistentemente vidas medias más cortas que mayores a misma longitud de onda, reflejando menor AQY y mayores contribuciones no radiativas superficiales.

Espectroscopía Resuelta en Frecuencia

FRS se validó usando circuito RC conocido (decaimiento mono-exponencial, τ = 78,9 µs) y microesfera dopada con quelato Eu (τ ≈ 670 µs). Los espectros QFRS coincidieron con función de respuesta esperada (Eq. 18) con frecuencias pico de 12,7 kHz y 1570 Hz, respectivamente.

Datos FRS de SiNC fueron solo ligeramente más anchos que respuesta intrínseca, necesitando deconvolución via algoritmo Richardson-Lucy para recuperar distribución de vida media. Las distribuciones deconvolucionadas (Fig. 6) abarcaron ~0,5 décadas, con NC pequeños mostrando forma simétrica en escala log y pico a 19,9 kHz (50,3 µs), y grandes pico a 6,28 kHz (159,2 µs). Ajustes de modelos SE y lognormal produjeron formas de distribución y frecuencias pico diferentes, todas menores que picos QFRS, reflejando sensibilidad de TRS de estado estacionario a duración de pulso y responsividad de detector.

La responsividad del detector afectó significativamente resultados TRS; mediciones APD y PMT difirieron en factor ~2 en tiempos medios de decaimiento. APD pico a 600 nm, PMT a 850 nm, llevando a ajustes de “mejor” modelo divergentes. Corregir responsividad en FRS es factible si se conoce dependencia espectral de distribución de tasas de decaimiento; TRS carece de método de corrección directo.

Conclusiones

Derivamos el decaimiento de población correspondiente al decaimiento de luminiscencia SE comúnmente usado y proporcionamos expresiones para vidas medias. Dos muestras SiNC crecidas térmicamente, funcionalizadas con dodeceno (diámetros medios 2,9 nm y 5,4 nm) fueron sintetizadas y caracterizadas. Datos TRS ajustados con modelos SE, lognormal y bimolecular no capturaron completamente forma de decaimiento; vidas medias dependieron de responsividad de detector y elección de modelo. Espectroscopía resuelta en frecuencia, tras deconvolución, produjo distribución de vida media directamente, revelando dispersión ~0,5 décadas y sesgo alta frecuencia para NC grandes. FRS es por tanto método preferido para sondear dinámicas de luminiscencia SiNC, aunque deconvolución cuidadosa y calibración de detector son esenciales. Estudios futuros deben incorporar correcciones de responsividad de detector y evitar default a modelos SE sin validar distribución subyacente.