Tutorial de SCLC: Inicio rápido — simular corriente limitada por carga espacial

La corriente limitada por carga espacial (SCLC) es el régimen de transporte donde los portadores inyectados dominan y la corriente está limitada por su movimiento a través de la película, no por la generación. En un dispositivo ideal, libre de trampas, la densidad de corriente sigue la ley de Mott–Gurney: \( J = \frac{9}{8}\,\varepsilon\,\mu\,\frac{V^2}{L^3} \), con constante dieléctrica \( \varepsilon \), movilidad \( \mu \), voltaje \( V \) y espesor \( L \). Las medidas SCLC (a menudo usando diodos solo de huecos o solo de electrones) se utilizan ampliamente para extraer movilidad y evaluar efectos de trampas. En este inicio rápido, configurará una estructura SCLC, ejecutará un barrido JV, localizará la región J ∝ V² y verá cómo las trampas o el espesor desplazan la curva y la movilidad extraída.

En ?? comparamos configuraciones de contacto que controlan qué portadores pueden entrar en el dispositivo. En el panel (a) la estructura estándar tiene contactos selectivos de electrones y huecos, creando un potencial incorporado y permitiendo que tanto electrones como huecos sean inyectados/extraídos. Ajustando la energética de los contactos o añadiendo/seleccionando capas de transporte/bloqueo, se puede imponer inyección de un solo portador: en el panel (b) se forma un dispositivo solo de electrones (SCLC) proporcionando bajas barreras hacia la banda de conducción en ambos contactos mientras se bloquea la banda de valencia (inyección de huecos), y en el panel (c) se forma un dispositivo solo de huecos (SCL) alineando la banda de valencia en ambos contactos mientras se bloquea la banda de conducción (inyección de electrones). En comparación con el dispositivo estándar en (a), los casos de un solo portador (b,c) suprimen la recombinación y fuerzan que la corriente esté gobernada por transporte limitado por carga espacial, lo que es ideal para extraer la movilidad de portadores y los efectos de los contactos.

Paso 1: Crear una nueva simulación

Inicie OghmaNano desde el menú Inicio de Windows. La ventana principal de OghmaNano aparecerá como se muestra en ??.

Paso 2: Comprobar que los contactos están configurados para SCLC

Una vez que haya guardado el ejemplo SCLC, aparecerá la ventana principal (??). Esta muestra la representación 3D del dispositivo. Los dos botones clave que utilizaremos en esta ventana son el botón Run Simulation y el botón Contacts.

Antes de ejecutar la simulación haciendo clic en Run Simulation (o pulsando F9), primero debemos examinar los contactos. En la ventana principal (??), el botón Contacts está resaltado en el recuadro rojo inferior. Al hacer clic en él se abre el editor de contactos. El editor de contactos permite que ambos contactos se configuren con el mismo tipo de portador. Por ejemplo, ambos pueden definirse como contactos de Hole (??) o ambos como contactos de Electron (??). En otras palabras, se selecciona un tipo de portador y se aplica a ambos lados del dispositivo.

Esta configuración es muy diferente de la de una célula solar, donde los dos contactos deben ser distintos: un contacto de electrones y un contacto de huecos. No importa qué lado sea cuál, pero la asimetría es esencial para separar cargas e impulsar la corriente. Sin embargo, en las medidas SCLC la simetría es deliberada. Si ambos contactos se configuran como Electron, el dispositivo mide la movilidad electrónica. Si ambos se configuran como Hole, mide la movilidad de huecos. La lógica es simple: los contactos inyectan el tipo de portador para el que están configurados, y ese portador domina el transporte a través del dispositivo.

Paso 2: Comprobar que la simulación se ejecuta en la oscuridad

SCLC casi siempre se realiza en la oscuridad. Antes de ejecutar el experimento, asegúrese de que la simulación también esté en condiciones de oscuridad. Esto se hace desde la cinta Optical ajustando Light intensity (Suns) a 0.0, como se muestra en ??.

Paso 3: Activar más salida opcional

El paso final antes de ejecutar la simulación SCLC es abrir el JV editor, al que puede acceder desde la cinta Editors → JV editor en la ventana principal (??). Esto abre la ventana de configuración JV (??).

En esta ventana, asegúrese de que Save parameter sweeps esté ajustado a Disk, como se muestra en la parte inferior de ??. Activar esta opción garantiza que parámetros clave del dispositivo, como movilidad, densidades de portadores y otras magnitudes eléctricas, se registren como función del voltaje aplicado. En lugar de producir instantáneas 1D simples, el software integra estos valores a través del dispositivo y los almacena para que posteriormente puedan representarse frente al voltaje. Esto es esencial para el análisis SCLC porque nos permite ver cómo evoluciona la movilidad con el voltaje y determinar los valores extraídos de las curvas SCLC.

De forma predeterminada, esta función suele estar desactivada porque escribir datos de barrido en disco puede ralentizar las simulaciones. En la mayoría de los casos se desea minimizar la salida a disco para mantener las ejecuciones eficientes. Sin embargo, para SCLC merece la pena activarla a pesar del tiempo extra requerido, ya que sin ella no sería posible analizar correctamente las tendencias de movilidad en el dispositivo.

Paso 4: Ejecutar la simulación

Una vez que la simulación SCLC esté completamente preparada — configurando ambos contactos como Hole o Electron, asegurando que la luz esté apagada y activando la opción Save Parameter Sweep — ya está listo para ejecutar. Regrese a la ventana principal e inicie la simulación haciendo clic en el botón Play o pulsando F9.

Cuando la ejecución finalice, abra la pestaña Output (??). Aquí verá la salida estándar de la simulación, incluido jv.csv y el directorio sweep. Hacer doble clic en jv.csv abre la gráfica JV (??), que muestra una curva SCLC típica. Al pulsar L seguido de Shift+L se cambian ambos ejes a escala logarítmica, facilitando la identificación del régimen SCLC en los datos.

Paso 5: Extraer la movilidad SCLC por el camino difícil

Para determinar la movilidad del portador de carga a partir de una medida SCLC, utilizamos la ley de Mott–Gurney, que relaciona la densidad de corriente medida con el voltaje aplicado bajo condiciones libres de trampas. El paso clave es identificar el régimen SCLC en la curva JV, medida en la oscuridad (??). El régimen SCLC ocurre una vez que la inyección de portadores desde los contactos es suficientemente eficiente como para que la corriente ya no esté limitada por portadores generados térmicamente, sino por la acumulación de carga espacial en el dispositivo. En este régimen la densidad de corriente aumenta cuadráticamente con el voltaje, siguiendo la relación \( J \propto V^2 \).

En una gráfica log–log de densidad de corriente frente a voltaje, el régimen SCLC puede reconocerse por una región rectilínea con pendiente de aproximadamente 2.0. A voltajes más bajos la pendiente es más cercana a 1.0, reflejando conducción óhmica dominada por portadores en equilibrio. A voltajes más altos la corriente puede desviarse nuevamente del comportamiento cuadrático si el llenado de trampas, la resistencia en serie o los efectos de alto campo se vuelven importantes. Por tanto, la movilidad debe extraerse específicamente de la región de voltaje intermedio donde la pendiente es cercana a 2, ya que esto corresponde a la condición SCLC libre de trampas asumida en la ley de Mott–Gurney.

Para una corriente limitada por carga espacial libre de trampas, la densidad de corriente viene dada por

$$ J = \frac{9}{8}\,\varepsilon \mu \frac{V^2}{L^3}. $$

Reordenando para la movilidad:

$$ \mu = \frac{8}{9} \cdot \frac{J L^3}{\varepsilon V^2}. $$

Sustituyendo valores representativos:

• \( L = 100~\text{nm} = 1.0\times10^{-7}~\text{m} \;\;\Rightarrow\;\; L^3 = 1.0\times10^{-21}~\text{m}^3 \)
• \( \varepsilon_r = 3.0 \), \( \varepsilon_0 = 8.85\times10^{-12}~\text{F·m}^{-1} \), por lo que \( \varepsilon = 2.655\times10^{-11}~\text{F·m}^{-1} \)
• \( V = 1.0~\text{V} \)
• \( J \approx 1.0\times10^{3}~\text{A·m}^{-2} \)

Ahora:

$$ \mu = \frac{8}{9} \cdot \frac{ (1.0\times10^{3})(1.0\times10^{-21}) } { (2.655\times10^{-11})(1.0^2) } = 3.35\times10^{-8}~\text{m}^2\text{V}^{-1}\text{s}^{-1}. $$

En unidades cgs, esto corresponde a

Este ejemplo resuelto muestra cómo la movilidad puede extraerse directamente de la curva JV en el régimen SCLC, usando la relación de Mott–Gurney.

Paso 6: Representar la movilidad real dentro del dispositivo

Una de las verdaderas fortalezas de la simulación es que no está restringido a las predicciones de modelos analíticos. En su lugar, puede mirar directamente dentro del dispositivo simulado y examinar cómo evolucionan las magnitudes físicas a medida que se aplica el voltaje. En este caso, podemos explorar los resultados almacenados en la carpeta sweep, mostrada en ??. Abrir esta carpeta muestra una lista de archivos (??) que contienen parámetros del dispositivo guardados como función del voltaje. Estos incluyen tasas de generación, densidades de portadores, tasas de recombinación y muchas otras magnitudes. Para el análisis de movilidad, los archivos clave son mun.csv y mup.csv, que informan de las movilidades electrónicas y de huecos simuladas. El archivo que elija depende de si configuró el dispositivo anteriormente con contactos de electrones o de huecos. En este ejemplo, nos interesa el transporte electrónico, por lo que examinamos la salida mun.csv, mostrada en ??. Esta gráfica revela cómo cambia la movilidad con el voltaje aplicado.

La tendencia creciente aparece en dispositivos que contienen trampas. A medida que crece la polarización, se llenan más trampas y se liberan portadores adicionales a estados libres, lo que conduce a un aumento aparente de la movilidad. Este comportamiento va más allá de la simple descripción de Mott–Gurney y resalta el valor de la simulación: no solo reproduce la curva corriente–voltaje, sino que también nos permite ver los procesos físicos subyacentes que configuran el rendimiento del dispositivo.

Comparando el resultado analítico de la Sección 5 con el resultado numérico de la Sección 6, obtenemos una movilidad analítica de aproximadamente \(3.35\times10^{-8}\ \mathrm{m^2\,V^{-1}\,s^{-1}}\) frente a un valor simulado de aproximadamente \(2\times10^{-8}\ \mathrm{m^2\,V^{-1}\,s^{-1}}\) para el caso SCLC. La discrepancia es relativamente pequeña aquí, pero resalta un resultado común al contrastar enfoques analíticos y numéricos: los modelos analíticos simplificados (por ejemplo, suposiciones ideales y libres de trampas de Mott–Gurney, inyección perfectamente óhmica, campos uniformes, sin resistencia en serie ni dependencia de campo) pueden diferir de simulaciones completas del dispositivo que resuelven trampas, variaciones espaciales y contactos no ideales. Como resultado, los valores analíticos y numéricos no tienen por qué coincidir exactamente — incluso en SCLC — aunque sus tendencias deben ser consistentes.