Tutorial de dominio de excitones (Parte A): Ejecutar una simulación 3D de dominio de excitones

1. Introducción

Comprender cómo se generan, transportan y disocian los excitones es central para el funcionamiento de las células solares orgánicas de heterounión volumétrica (BHJ). En estos sistemas, la luz se absorbe predominantemente en la fase donadora, generando pares electrón–hueco ligados conocidos como excitones. Estos excitones deben difundirse a través del material donador y alcanzar una interfaz donador–aceptor antes de decaer. En la interfaz pueden disociarse en electrones y huecos libres, con una probabilidad determinada por la morfología local, la dimensionalidad y las constantes cinéticas relevantes.

El equilibrio entre la difusión de excitones, la disociación interfacial y los procesos de pérdida competidores desempeña por tanto un papel decisivo en la fijación de la eficiencia efectiva de fotogeneración. Sin embargo, históricamente se ha evitado el modelado explícito de excitones no por dificultad conceptual, sino por la incertidumbre de los parámetros. Las longitudes de difusión, los tiempos de vida, las tasas de disociación interfacial, los canales de decaimiento radiativo y no radiativo, y los procesos de aniquilación de excitones estaban a menudo poco restringidos para un sistema material dado. En muchos flujos de trabajo de modelado prácticos, estos efectos se absorbían por tanto en un único factor escalar de eficiencia fotónica, \(\eta_{\mathrm{photon}}\), que representaba las pérdidas geminadas netas sin resolver el transporte y la cinética subyacentes.

Esta situación está cambiando actualmente. Los avances en caracterización experimental están empezando a proporcionar mediciones directas de tiempos de vida de excitones, longitudes de difusión y canales de pérdida en sistemas tipo BHJ. A medida que estos parámetros quedan mejor restringidos, el modelado explícito de dominios de excitones se vuelve cada vez más informativo, no como un medio para predecir eficiencias absolutas del dispositivo, sino como una forma de explorar cómo el tamaño de dominio, la dimensionalidad y la cinética de materiales determinan conjuntamente un rendimiento efectivo de generación de carga. El modelo de dominio de excitones usado en este tutorial fue desarrollado y publicado como parte de este reciente trabajo experimental–de modelado (véase, por ejemplo, Nature Materials 21, 55–61 (2022)).

En este tutorial adoptamos por tanto una geometría de celda unidad idealizada pero totalmente tridimensional: un dominio donador embebido dentro de una matriz aceptora, representado inicialmente como una esfera donadora dentro de una caja aceptora. Aunque la geometría está simplificada, el modelo trata explícitamente la generación de excitones, la difusión, la disociación interfacial y los procesos de pérdida competidores. Esto lo convierte en un marco práctico y físicamente transparente para explorar cómo interactúan parámetros medidos experimentalmente, para comprobar su consistencia interna y para desarrollar intuición sobre cómo la morfología y la cinética controlan conjuntamente la eficiencia efectiva de fotogeneración en tres dimensiones.

2. Ecuaciones gobernantes

En este tutorial, el transporte de excitones se trata explícitamente resolviendo la ecuación tridimensional de difusión de excitones en todo el dominio de simulación. Cuando el modelo de excitones está activado, la absorción óptica alimenta directamente la población de excitones, que luego evoluciona según la siguiente ecuación:

\[ \frac{\partial X}{\partial t} = \nabla \cdot \left( D \nabla X \right) + G_{\mathrm{optical}} - k_{\mathrm{dis}} X - k_{\mathrm{FRET}} X - k_{\mathrm{PL}} X - \alpha X^2. \]

Aquí \(X(\mathbf{r},t)\) es la densidad de excitones (\(\mathrm{m^{-3}}\)) y \(D\) es el coeficiente de difusión de excitones (\(\mathrm{m^2\,s^{-1}}\)). El término fuente \(G_{\mathrm{optical}}\) (\(\mathrm{m^{-3}\,s^{-1}}\)) representa la generación de excitones debida a la absorción óptica. Los términos restantes describen los procesos competidores que eliminan excitones de la población:

• Disociación (\(k_{\mathrm{dis}} X\)): conversión de excitones en electrones y huecos libres en interfaces donador–aceptor. Este es el canal productivo que finalmente contribuye a la generación de carga.
• Transferencia tipo Förster (\(k_{\mathrm{FRET}} X\)): transferencia de energía no radiativa de excitones a sitios o dominios vecinos, que redistribuye la energía del excitón sin generar directamente portadores libres.
• Pérdida radiativa (\(k_{\mathrm{PL}} X\)): recombinación de excitones vía emisión de fotones (fotoluminiscencia), que representa una vía intrínseca de pérdida.
• Aniquilación (\(\alpha X^2\)): aniquilación bimolecular excitón–exciton a altas densidades, que se vuelve importante bajo excitación fuerte o en geometrías confinadas.

El coeficiente de difusión \(D\) se parametriza usando una longitud de difusión de excitones \(L\) (m) y un tiempo de vida del excitón \(\tau\) (s) mediante

\[ D = \frac{L^2}{\tau}. \]

En este tutorial ejecutamos una simulación solo de excitones: el solucionador calcula la densidad tridimensional de excitones \(X(\mathbf{r})\) y los términos de reacción asociados, pero no resuelve el transporte de electrones y huecos. El término de disociación \(k_{\mathrm{dis}} X\) se interpreta por tanto como una tasa de generación de carga potencial resuelta espacialmente, que indica dónde y con qué eficiencia los excitones se convertirían en portadores libres en interfaces donador–aceptor.

3. Crear la simulación de dominio de excitones

Abra la ventana New simulation (Figura ??). Puede accederse a ella haciendo clic en New simulation en la ventana principal. Si hace doble clic en Exciton simulations, verá la biblioteca de excitones (Figura ??). Haga doble clic en Exciton domain para abrir el proyecto de ejemplo.

4. Inspeccionar la geometría y los parámetros

Después de abrir el ejemplo, la ventana principal muestra una escena tridimensional simple: una esfera donadora embebida dentro de una caja aceptora (Figura ??). Los rayos verdes indican iluminación incidente desde arriba. Para inspeccionar los parámetros de excitones, haga clic en Electrical parameters en el panel izquierdo (dentro de la pestaña Device structure). Esto abre el editor de parámetros eléctricos. Si se desplaza hacia la parte inferior de esta ventana, encontrará una sección etiquetada Excitons, donde los parámetros específicos de excitones se enumeran por separado para cada objeto de la escena (la capa circundante y la esfera embebida), como se muestra en las Figuras ?? y ??.

La sección de excitones define los parámetros que aparecen en la ecuación de difusión–reacción. Estos incluyen la longitud de dispersión del excitón \(L\), que fija el coeficiente de difusión mediante \(D = L^{2}/\tau\), y el tiempo de vida del excitón \(\tau\), que controla la escala temporal global del transporte en ausencia de canales de pérdida adicionales. La recombinación radiativa se describe mediante la tasa de fotoluminiscencia \(k_{\mathrm{PL}}\), mientras que los procesos de transferencia de energía no radiativa se capturan mediante la tasa tipo Förster \(k_{\mathrm{FRET}}\). A densidades de excitones más altas, la aniquilación bimolecular se incluye mediante el coeficiente \(\alpha\), correspondiente al término \(-\alpha X^{2}\) en la ecuación gobernante. Finalmente, la tasa de disociación \(k_{\mathrm{dis}}\) especifica con qué eficiencia los excitones se convierten en estados de transferencia de carga en interfaces donador–aceptor y es el parámetro clave que controla la eficiencia efectiva de generación de carga en este modelo.

La característica conceptual clave de este ejemplo es que la esfera (donador) y la capa (matriz aceptora) usan física de excitones diferente. En las capturas de pantalla anteriores, la esfera donadora incluye un \(k_{\mathrm{PL}}\) no nulo (un canal de pérdida). La región circundante está configurada con un \(k_{\mathrm{dis}}\) muy grande, de modo que una vez que los excitones alcanzan (o salen hacia) esa región se disocian rápidamente. Esto hace que la frontera donador/aceptor se comporte como un sumidero eficiente, produciendo el perfil característico de “alto en el centro, agotado cerca de la frontera”.

5. Ejecutar la simulación

Haga clic en el triángulo azul Run simulation para iniciar el solucionador. La salida del terminal aparecerá en la ventana principal, como se muestra en la Figura ??. Las líneas iniciales informan información general sobre la configuración de la simulación y la geometría. El diagnóstico clave aparece una vez que el solucionador de excitones empieza a ejecutarse y tiene la forma Exciton f() = … , steps = …. Esta línea informa el progreso del solucionador tridimensional de excitones. La cantidad f() es un residuo que mide cuán lejos está el campo actual de densidad de excitones de la autoconsistencia; a medida que el solucionador itera, este valor debería disminuir. El contador steps asociado indica el número de iteración. En el ejemplo mostrado, el residuo desciende desde un valor inicial hasta 2.97 × 104 en el paso 2, tras lo cual el solucionador informa que la ecuación ha sido resuelta y termina.

En un portátil moderno típico, este ejemplo debería completarse en unos 5–10 segundos. Si el tiempo de ejecución se prolonga a minutos, esto suele indicar que la malla es excesivamente fina o que la geometría o los parámetros se han modificado de una forma que incrementa mucho el número de puntos de malla. En ese caso, conviene revisar la configuración de la malla o repetir los pasos previos de configuración antes de continuar.

6. Graficar la salida

Una vez terminada la ejecución, cambie a la pestaña Output (Figura ??). Los elementos más relevantes para este tutorial son el directorio exciton_output/, que contiene los resultados de excitones resueltos espacialmente, y el archivo exciton_sim_info.json, que resume estadísticas globales de generación y pérdidas y se utilizará más adelante en la Parte C. Al hacer doble clic en exciton_output/ se muestra el contenido de la Figura ??.

Abra G.csv para visualizar la tasa de generación de excitones producida por el modelo óptico (Figura ??). En la ventana de gráfico 3D puede rotar la escena con el ratón y usar los deslizadores Z e Y en la parte inferior para cortar el volumen. Esta suele ser la forma más rápida de verificar que la geometría, la dirección de iluminación y la malla se comportan como se pretende.

A continuación, abra exciton.csv (Figura ??), que muestra la densidad de excitones en estado estacionario. Debería observar la firma característica de difusión hacia un sumidero interfacial: una alta densidad de excitones en el interior de la esfera donadora, con agotamiento hacia la frontera donador–aceptor, donde la disociación es fuerte. Use los deslizadores de corte para confirmar que este agotamiento sigue la geometría de la interfaz y no es un artefacto del ángulo de visualización.

Finalmente, abra Gn.csv (Figura ??), que muestra la tasa de generación de electrones resuelta espacialmente que surge de la disociación de excitones. En este ejemplo debería ver un anillo (o capa) distintivo de generación localizado cerca de la interfaz donador–aceptor. Use los deslizadores Z e Y para mover el plano de corte a través del volumen y explorar cómo evoluciona este perfil de disociación con la posición.

7. ¿Qué son los archivos en exciton_output/?

El directorio exciton_output/ contiene dos tipos distintos de archivos. El primer grupo consiste en campos calculados producidos por el solucionador de excitones, como la densidad de excitones y las tasas de generación derivadas de la disociación. El segundo grupo contiene mapas de parámetros: copias tridimensionales de los parámetros de excitones definidos en la GUI. Estos se escriben principalmente como ayuda de consistencia y depuración, permitiéndole verificar que los parámetros de material previstos se han aplicado correctamente en el espacio. Las Tablas 1 y 2 resumen estas dos categorías. Todas las cantidades se informan en unidades SI. (Algunos archivos pueden no escribirse si los submodelos correspondientes están desactivados.)

Los archivos de mapas de parámetros de la Tabla 2 no son salidas del solucionador en el sentido habitual. Simplemente reflejan los parámetros de excitones definidos en la sección Electrical → Excitons de la GUI, escritos sobre la misma malla tridimensional que los campos solución. Su propósito es proporcionar una forma transparente de comprobar que los parámetros se asignan a los objetos y regiones correctos del dominio, particularmente cuando se trabaja con geometrías complejas.

8.1 Interpretación de exciton_sim_info.json (eficiencia de disociación)

Aunque los gráficos tridimensionales son invaluables para visualizar dónde se generan, pierden y disocian los excitones, por sí solos no proporcionan resúmenes numéricos concisos. El archivo exciton_sim_info.json cumple esta función informando tasas de generación y pérdidas integradas espacialmente para cada objeto de la simulación (la capa circundante, la esfera embebida y cualquier forma adicional que pueda añadir).

Para ver esta información, abra la pestaña principal Output y haga doble clic en exciton_sim_info.json. El archivo contiene una estructura JSON anidada que enumera, para cada objeto, la tasa total de generación de excitones y las tasas totales asociadas con cada canal competidor de pérdida o disociación. Estas cantidades le permiten cuantificar cómo los cambios en parámetros materiales, geometría o morfología se traducen en cambios en la eficiencia de generación de carga.

La clave de nivel superior segment0 simplemente marca el inicio del informe. Dentro de ella, el campo name identifica el objeto envolvente (aquí, la layer), mientras que el bloque anidado shape0 corresponde al objeto embebido (la sphere). Cada bloque informa la contribución total de ese objeto a la generación, la disociación y los procesos de pérdida.

La Tabla 2 resume los campos que encontrará en este archivo y cómo interpretarlos.

En conjunto, los gráficos tridimensionales y los resúmenes numéricos en exciton_sim_info.json proporcionan una visión complementaria de la simulación: los gráficos muestran dónde ocurren los procesos, mientras que el archivo JSON cuantifica cuánto contribuye cada proceso. Esto lo convierte en una herramienta útil para analizar sistemáticamente cómo afectan los cambios en parámetros o geometría a la disociación de excitones y a la generación efectiva de carga.

9. Integración final

En conjunto, los resultados tridimensionales cuentan una historia física coherente. La absorción óptica genera excitones principalmente en el interior del dominio donador. Estos excitones luego se difunden a través del donador, se agotan cerca de la interfaz donador–aceptor y se disocian eficientemente una vez que alcanzan la región aceptora circundante. Los gráficos espaciales muestran dónde ocurren estos procesos; el informe resumen exciton_sim_info.json cuantifica cuánto contribuye cada proceso.

En este ejemplo, el informe resumen separa limpiamente la física subyacente por región espacial, permitiendo interpretar la generación, las pérdidas y la disociación de una manera transparente.

• Esfera (donador): La mayor parte de la absorción óptica — y por tanto de la generación de excitones — ocurre dentro de la esfera donadora, lo que se refleja en un valor grande de tot_G en el bloque shape0/sphere. Dentro del donador, el canal dominante de pérdida competidora en este ejemplo es el decaimiento radiativo, cuantificado por tot_Rk_pl. La razón entre tot_Rk_pl y tot_G proporciona una medida directa de cuánto suprimen las pérdidas radiativas la población de excitones antes de que pueda alcanzar la interfaz donador–aceptor.
• Capa (matriz aceptora): La capa circundante absorbe solo débilmente en este ejemplo, de modo que su contribución a la generación de excitones (tot_G) es pequeña. Su papel principal es en cambio la conversión de excitones: la capa está configurada con un fuerte canal de disociación, haciendo de tot_Rk_dis el término dominante en el bloque layer. Esto corresponde directamente a la disociación interfacial observada en los gráficos espaciales, donde Gn.csv muestra una región de generación de electrones en forma de anillo o capa concentrada en la frontera donador–aceptor.

La métrica final dis_eff informa la eficiencia global de disociación del sistema (aquí aproximadamente \(89\%\)). Físicamente, esto indica que la mayoría de los excitones generados en el donador se difunden hasta la interfaz y se disocian, con solo una fracción relativamente pequeña perdida por decaimiento radiativo. Otros canales, como la transferencia Förster o la aniquilación excitón–exciton, están inactivos en este ejemplo particular.

En otras palabras, la simulación opera en un régimen donde la disociación interfacial supera al decaimiento de excitones, dando lugar a un alto rendimiento efectivo de conversión de fotones en portadores libres. Esta combinación de visualización tridimensional e informe cuantitativo proporciona una imagen compacta pero físicamente completa de cómo la geometría y los parámetros cinéticos determinan conjuntamente la eficiencia de generación de carga en un modelo de dominio de excitones.