Simulación de morfologías de heterounión en volumen con cortes eléctricos 2D

1. Introducción

Las heterouniones en volumen (BHJ) son la microestructura de trabajo para muchos dispositivos electrónicos orgánicos—de forma muy destacada las células solares orgánicas, pero también fotodetectores, sensores y otros sistemas mixtos iónico–electrónicos donde importan las trayectorias percolantes y los cuellos de botella locales. En una BHJ, dos fases se interpenetran a la escala de 10–100 nm, creando una gran interfaz interna y un paisaje de transporte que fundamentalmente no es unidimensional. La contrapartida es la complejidad: una vez que se resuelve la microestructura se heredan inevitablemente límites de percolación, callejones sin salida, constricciones y una fuerte variación espacial en la densidad de portadores y la recombinación. Para la mayor parte del ajuste diario de dispositivos, no empezaría aquí. Normalmente se trabajaría al nivel de dispositivo 1D—rápido, robusto y, a menudo, completamente apropiado—posiblemente incluyendo física adicional como atrapamiento y recombinación. Por ejemplo, el tutorial de inicio rápido de OPV demuestra un flujo de trabajo 1D estándar (incluyendo física de trampas cuando es necesario) que es la herramienta adecuada para la mayoría de los estudios rutinarios y la extracción de parámetros.

En el otro extremo de la escala de abstracción, si su objetivo son cuestiones a nivel de sistema como escalado, pérdidas serie/paralelo o comportamiento de módulo/panel, normalmente pasaría a una descripción de gran área / consciente del circuito en lugar de resolver la morfología a escala nanométrica. El tutorial de gran área PM6:Y6 es un ejemplo de ese enfoque de modelado de más alto nivel.

Este tutorial va deliberadamente en la dirección opuesta: mantiene la microestructura (mediante una morfología BHJ explícita) pero sigue siendo computacionalmente práctico al resolver el problema eléctrico sobre cortes 2D. La idea no es que cada dispositivo deba ajustarse de esta manera; la idea es que disponer de un flujo de trabajo que resuelva la microestructura es científicamente y tecnológicamente valioso. Permite poner a prueba las suposiciones de medio efectivo, ver dónde fallan los modelos “simples” y desarrollar intuición sobre los límites de rendimiento impulsados por la morfología—particularmente cuando dominan la percolación y el atrapamiento/recombinación espacialmente localizados.

2. Cree una simulación de morfología e inspeccione la BHJ

Desde la ventana principal, haga clic en New simulation. Esto abre la biblioteca de tipos de dispositivo mostrada en ??. Haga doble clic en Morphology y, a continuación, guarde el proyecto en su disco local. Guardarlo localmente importa porque las simulaciones resueltas en morfología pueden generar muchos archivos intermedios (mallas, instantáneas, cachés), y el almacenamiento local suele ser mucho más rápido y fiable que las carpetas de red o sincronizadas con la nube.

Tras crear el proyecto, se abre la ventana principal del dispositivo. Debería ver un objeto rojo central que representa una fase de la heterounión en volumen (BHJ), incrustado dentro de la región envolvente de la capa activa, como se muestra en ??. Aunque la vista 3D completa es informativa, la morfología puede ser visualmente densa. Para hacer la estructura interna más clara, haga clic con el botón derecho en el objeto BHJ y seleccione View → Show cut through Y.

Al habilitar la vista en sección se obtiene la representación mostrada en ??, que expone un corte 2D a través de la morfología 3D. Esto facilita mucho ver la conectividad de fase, las interfaces y las trayectorias de percolación dentro de la BHJ.

Un modelo mental útil en esta etapa es el siguiente: la región rectangular envolvente define el volumen de la capa activa, mientras que el objeto rojo incrustado define una fase explícita de la BHJ (por ejemplo, el polímero). La fase complementaria se define implícitamente como el resto del volumen de la capa activa no ocupado por el objeto incrustado. Esta separación se vuelve explícita al examinar las propiedades de los objetos en el editor de objetos en las secciones siguientes.

3. Estructura del dispositivo: capa activa y morfología BHJ incrustada

Abra la pestaña Device structure y haga clic en Layer editor. El dispositivo consta de tres capas: un contacto superior, una única capa activa y un contacto inferior (Figura ??). Esta pila deliberadamente simple facilita entender cómo entra la morfología en el problema eléctrico, sin introducir interfaces adicionales ni capas de transporte. El punto clave es que la capa activa en sí misma es un volumen rectangular simple definido en el editor de capas. En este ejemplo está asignada al material molecular Y6, un aceptor, y-por sí sola-es simplemente una caja. Toda la complejidad estructural asociada a la heterounión en volumen se introduce por separado, incrustando un segundo objeto dentro de este volumen de capa activa. Ese objeto incrustado es una malla CAD compleja importada de la base de datos de formas, que representa la fase polimérica PM6. La malla polimérica define un volumen tridimensional cerrado dentro de la capa activa. La fase complementaria (la región rica en aceptor) se define entonces implícitamente como la parte de la capa activa que no está ocupada por la malla polimérica. En otras palabras, la BHJ se representa mediante exclusión de volumen en lugar de definir explícitamente dos mallas entrelazadas.

Esta separación se hace clara al inspeccionar los objetos en la vista 3D. Hay dos objetos relevantes: la caja de la capa activa envolvente (normalmente mostrada como una estructura de alambre azul), y la morfología del polímero incrustada (mostrada en rojo). Para inspeccionar sus propiedades, haga clic con el botón derecho en la región envolvente azul y seleccione Edit object, como se muestra en ??.

En el editor de objetos, la distinción es explícita. Para el objeto envolvente Active, la Object shape es una caja simple (con las mallas complejas deshabilitadas). Para el objeto Polymer, la Object shape apunta a una entrada de Shape DB, que es la malla triangulada que define la morfología PM6. Cambiar esta referencia intercambia la morfología del polímero mientras deja inalterada la capa activa circundante.

Esta construcción—caja simple de capa activa más volumen polimérico incrustado—es central para la forma en que se manejan las morfologías BHJ en OghmaNano. Hace que la geometría sea fácil de razonar, garantiza que la fase polimérica forme una región cerrada bien definida y le permite estudiar cómo la percolación, la conectividad y las constricciones en la red polimérica influyen en el transporte de carga una vez que se ejecuta el solucionador eléctrico.

4. Parámetros eléctricos: movilidades y bandas de trampas SRH

Abra el editor de Electrical parameters desde la ventana principal. Se proporcionan paneles de parámetros separados para las dos fases de la capa activa: la región envolvente Active (aceptor / matriz) y la morfología incrustada Polymer, mostradas en ?? y ??. Aquí es donde las dos fases interpenetrantes se hacen eléctricamente distintas.

Ambas fases utilizan el mismo formalismo básico de transporte: drift–diffusion con estadísticas de portadores Maxwell–Boltzmann y tensores de movilidad simétricos. La diferencia clave es el contraste de movilidad intencionado. En la fase polimérica, la movilidad de huecos es alta mientras que la movilidad electrónica está suprimida; en la fase activa complementaria ocurre lo contrario. Esto impone transporte selectivo: los huecos percolan a través de la red polimérica, mientras que los electrones percolan a través de la matriz.

Esta simple asimetría es suficiente para generar un comportamiento BHJ físicamente significativo. Incluso antes de ejecutar la simulación, define lo que debería esperarse ver en las salidas espaciales: trayectorias continuas de percolación, callejones sin salida donde los portadores quedan atrapados geométricamente y acumulación local de carga en constricciones o regiones de la morfología mal conectadas. Además del transporte de portadores libres, se habilita en ambas fases el atrapamiento Shockley–Read–Hall (SRH) no en equilibrio. Como se muestra en los paneles de parámetros, la densidad de estados se representa mediante distribuciones exponenciales de cola tanto para electrones como para huecos, con densidades de trampas y pendientes energéticas idénticas en cada fase. El solucionador sigue cinco bandas de trampas discretas en cada punto de la malla.

Cada banda de trampas intercambia portadores dinámicamente con las poblaciones libres mediante términos explícitos de captura y emisión (las constantes de acoplamiento libre-a-atrapado visibles en el editor). Como resultado, el sistema que se resuelve incluye no solo densidades de electrones y huecos libres, sino también múltiples poblaciones de portadores atrapados en cada localización espacial. Esto incrementa sustancialmente la rigidez numérica y el coste computacional, pero es esencial para capturar transporte dispersivo, recombinación asistida por trampas y comportamiento transitorio realista en semiconductores orgánicos desordenados. Los valores numéricos utilizados aquí están pensados como un ejemplo estable e ilustrativo más que como un modelo de material calibrado. Una vez comprendido el flujo de trabajo, estos parámetros pueden refinarse, validarse frente al experimento o simplificarse si el objetivo es aislar los efectos de la morfología y la percolación más que la física de trampas.

5. Mallado

Antes de ejecutar la simulación, conviene comprobar la malla eléctrica para asegurarse de que es razonable. En la ventana principal, abra la cinta Electrical (Figura ??) y haga clic en Electrical mesh. Esto abre el editor de malla mostrado en ??.

En este ejemplo la malla es de 40 × 40 en el plano xy, lo que significa que el problema eléctrico se resuelve en dos dimensiones. La dirección z se omite intencionadamente. Las simulaciones eléctricas completamente tridimensionales son posibles en OghmaNano, pero requieren una verdadera morfología 3D y conllevan un aumento sustancial del coste computacional y de la complejidad. La construcción e interpretación de tales modelos queda fuera del alcance de este tutorial.

En su lugar, este cálculo debe entenderse como una sección transversal 2D a través del dispositivo. El corte pasa por el centro de la morfología y conserva la microestructura lateral completa de la BHJ dentro de ese plano. Este enfoque captura percolación, cuellos de botella y conectividad de fase manteniéndose rápido, interpretable y bien adaptado a una exploración a escala de tutorial.

6. Ejecución de la simulación

Haga clic en el botón azul Run simulation en la ventana principal para iniciar el solucionador. Mientras la ejecución está en curso, el terminal imprime el punto de polarización y la densidad de corriente en cada paso (??). En este ejemplo la polarización aplicada al top contact se incrementa desde 0.00 V hasta aproximadamente 0.54 V para generar un barrido JV. A baja polarización la densidad de corriente informada es negativa bajo iluminación, y a medida que la polarización se aproxima al régimen de circuito abierto la corriente neta cruza por cero y cambia de signo, de acuerdo con el comportamiento de una célula solar.

Cuando finalice la ejecución, abra la pestaña Output (??) para inspeccionar los archivos generados. El resumen del terminal informa del tiempo total de ejecución, el número de ecuaciones resueltas y estadísticas básicas de E/S de disco. En simulaciones resueltas en morfología conviene vigilar la salida a archivo: escribir un gran número de instantáneas puede dominar el tiempo de ejecución incluso cuando el solucionador en sí es rápido.

Para una comprobación de consistencia geométrica, haga doble clic en device.csv para ver la representación triangulada del dispositivo utilizada internamente por el solucionador (??). Esto confirma que la morfología BHJ se ha importado y discretizado como se esperaba.

Para ver el rendimiento eléctrico, abra jv_contact_0.csv o jv_contact_1.csv para representar la curva JV (??). Para esta simulación 2D, deben utilizarse las salidas JV por contacto; jv.csv puede ignorarse.

7.0 Examen de instantáneas de la simulación

Tras completarse la ejecución, las salidas más informativas se encuentran en el directorio snapshots/. Los archivos de instantáneas guardan variables internas de estado del sistema drift–diffusion–SRH, escritas en función de la polarización aplicada. Exponen la solución espacialmente resuelta de las ecuaciones de transporte en cada paso de voltaje, en lugar de resúmenes posprocesados. Cuando se abre una instantánea, están disponibles dos controles: un deslizador de polarización, que recorre el barrido JV, y un deslizador de posición del corte, que selecciona la sección 2D a través de la morfología. A lo largo de este tutorial el problema eléctrico se resuelve sobre un corte 2D, y todos los datos de instantáneas corresponden a ese corte.

Las Figuras ??– ?? muestran instantáneas 2D representativas a una polarización fija. Estas cantidades deben interpretarse conjuntamente. El paisaje de la banda de conducción (LUMO), Ec.csv (??), define el entorno energético para los electrones. En una simulación resuelta en morfología este paisaje varía lateralmente debido tanto a las energéticas dependientes de la fase como a la respuesta electrostática a la acumulación local de carga. La densidad de electrones libres, Q_nfree.csv (??), resalta regiones donde los electrones pueden moverse eficientemente a través de la red aceptor/matriz. Las fuertes variaciones espaciales indican trayectorias de percolación, constricciones y cuellos de botella inducidos por la morfología.

La densidad de electrones atrapados, Q_ntrap.csv (??), revela dónde se acumulan los portadores en estados trampa. Comparar las densidades libres y atrapadas suele ser la forma más rápida de identificar regiones de extracción deficiente o fuerte atrapamiento local vinculado a la morfología BHJ. La salida relacionada con la movilidad, mu_p_ft.csv (??), muestra una movilidad efectiva promediada sobre portadores libres y atrapados. Las variaciones espaciales reflejan tanto el contraste intrínseco de fase como el atrapamiento local, y a menudo se correlacionan fuertemente con la conectividad de la red subyacente.

Aunque las vistas 2D suelen ser la forma más clara de analizar la variación espacial, a veces es más fácil interpretar los mismos datos en espacio real. Al habilitar el modo 3D en la ventana de instantáneas, el corte 2D se proyecta de nuevo sobre la geometría del dispositivo. Las Figuras ??– ?? muestran ejemplos de esta proyección. Las movilidades promedio de electrones y huecos (mu_n_ft.csv, mu_p_ft.csv), el cuasi-nivel de Fermi de electrones libres (Fn.csv) y la densidad de corriente electrónica en la dirección de transporte (Jn.csv) son particularmente útiles para visualizar trayectorias de percolación y flujo de corriente a través de la BHJ.

Estas salidas de instantáneas proporcionan el principal valor científico de las simulaciones resueltas en morfología: permiten vincular directamente observables a escala de dispositivo con energéticas locales, atrapamiento y conectividad dentro de la heterounión en volumen.

8. Tareas para probar

La forma más rápida de desarrollar intuición es perturbar el modelo y observar qué cambia. Las tareas siguientes están diseñadas para ser rápidas: cambie la malla de morfología, vuelva a ejecutar y compare la curva JV (y, si tiene tiempo, un par de instantáneas a la misma polarización).

Tarea 1 — cambiar la malla de morfología

En la vista principal del dispositivo, haga clic con el botón derecho en el objeto BHJ rojo (la malla polimérica incrustada) y seleccione Mesh editor. En el editor, la opción actual debe ser Shape from database. Cambie la forma seleccionada a morphology/1, morphology/2 y morphology/3 (Figuras ??, ??, ??). Vuelva a ejecutar la simulación después de cada cambio y compare la curva JV.

Lo que debe buscar es cómo la conectividad y los cuellos de botella alteran la extracción y la recombinación: algunas morfologías producen trayectorias de percolación más limpias (corriente más alta), mientras que otras fuerzan a los portadores a pasar por cuellos estrechos y callejones sin salida (más acumulación, más atrapamiento y, normalmente, peor rendimiento JV).

Como experimento deliberadamente "incorrecto pero instructivo", también encontrará una forma llamada teapot en la base de datos. Cambie la malla BHJ a teapot, vuelva a ejecutar y genere una curva JV para un dispositivo con una inclusión con forma de tetera. No es una BHJ físicamente sensata, pero demuestra que el flujo de trabajo está realmente impulsado por la geometría: cualquier malla CAD cerrada puede insertarse y analizarse eléctricamente.

Tarea 2 — cambiar la intensidad de atrapamiento (bandas de trampas SRH)

En Electrical parameters, ajuste la configuración de non-equilibrium SRH traps y vuelva a ejecutar. Pruebe un cambio cada vez:

• Aumente/disminuya la densidad de trampas en 1-2 órdenes de magnitud y compare Q_nfree frente a Q_ntrap.

El objetivo de esta tarea es ver cómo el transporte dispersivo/limitado por trampas compite con la morfología: el atrapamiento tiende a amplificar las “malas” regiones (callejones sin salida, constricciones) porque la carga no puede vaciarse con tanta rapidez.

Tarea 3 — cambiar las constantes de velocidad de recombinación

Aún en Electrical parameters, modifique la recombinación y vuelva a ejecutar. Dos perturbaciones simples:

• Cambie la constante de velocidad de recombinación free-to-free (n_free a p_free) arriba/abajo en 1–2 órdenes de magnitud. Observe cómo responden la curva JV y la distribución espacial de carga.
• Si su proyecto expone parámetros de captura/acoplamiento SRH (constantes de acoplamiento libre↔trampa), escálelos arriba/abajo conjuntamente y compare el equilibrio de Q_nfree y Q_ntrap a la misma polarización.

Un hábito útil para las Tareas 2–3 es mantener fija la morfología, cambiar un único parámetro y comprobar siempre las mismas dos cosas: (i) cómo se desplaza la curva JV, y (ii) si las trayectorias de corriente (por ejemplo Jn y los campos de cuasi-nivel de Fermi) se vuelven más suaves o más “fracturadas” por la microestructura.