Simulación de célula solar orgánica de gran área (PM6:Y6): De redes de contactos 3D a dispositivos iluminados
1. Introducción
En este tutorial continuamos la serie sobre simulación de dispositivos de gran área, centrándonos en dispositivos que han ido más allá de pequeñas celdas de laboratorio y han pasado a geometrías a escala de aplicación. Se sitúa entre dos flujos de trabajo existentes en el manual: diseño de contactos de gran área (sin dispositivos activos), y módulos de perovskita 3D complejos. Este ejemplo conecta ambos, introduciendo diodos y absorción de luz en una estructura de contacto previamente diseñada, pero sin llegar a un módulo interconectado completo.
Aquí simulamos una única célula solar orgánica de gran área basada en una arquitectura tipo PM6:Y6 con contactos metálicos impresos hexagonales. El objetivo es comprender cómo se comporta un dispositivo que funciona bien a pequeña área cuando se escala, donde el flujo lateral de corriente, la caída de tensión y la resistencia de contacto comienzan a dominar. Aunque este tutorial utiliza materiales OPV, el flujo de trabajo es general: sustituyendo los datos ópticos n/k y ajustando los valores de los componentes, el mismo modelo puede representar su dispositivo de investigación.
Para aplicar este flujo de trabajo a su propio dispositivo, simplemente proporciona dos conjuntos de entradas: resistividades de capa (normalmente tomadas de la literatura o de estructuras de prueba sencillas), y un pequeño número de parámetros del diodo, como el factor de idealidad, la fotocorriente y la corriente de saturación. Una vez definidos, el modelo permite escalar virtualmente un dispositivo de laboratorio e identificar qué mecanismos de pérdida limitarán su eficiencia al pasar a geometrías reales de gran área.
Aquí no utilizamos el solucionador completo de deriva-difusión 3D, aunque OghmaNano sí lo proporciona. Para dispositivos de gran área, la deriva-difusión 3D completamente resuelta se vuelve rápidamente impracticable desde el punto de vista computacional y a menudo es innecesaria cuando las limitaciones dominantes provienen del transporte lateral y la resistencia de contacto, en lugar del movimiento de portadores a través de cada capa. En su lugar, OghmaNano construye un modelo de circuito 3D dependiente de la geometría, donde el transporte lateral y los contactos se representan mediante elementos resistivos y el comportamiento fotovoltaico se captura utilizando diodos integrados.
Mientras que el problema eléctrico se resuelve usando una formulación de circuito, la óptica se trata rigurosamente. La fotogeneración se calcula utilizando un solucionador óptico de matriz de transferencia evaluado en cortes unidimensionales a través de la geometría tridimensional, manteniendo los efectos de interferencia y absorción de película delgada mientras el solucionador de circuitos gestiona la recolección de corriente en gran área.
2. Crear la simulación
En la ventana principal, haga clic en New simulation. Verá una lista de categorías de simulación (véase ??). Después haga doble clic en Large area PM6:Y6 solar cell (véase ??).
Después de crear la simulación, la ventana principal se abre y muestra la geometría del dispositivo en tres dimensiones (véase ??). El contacto superior es una malla metálica hexagonal (tipo panal) etiquetada como Ag. Esta capa de plata actúa como una rejilla de recolección de corriente altamente conductora. Como está implementada como una malla, la mayor parte de la superficie permanece abierta, permitiendo que la luz atraviese el contacto. Debajo de la malla de plata se encuentra PEDOT:PSS, que sirve como una capa lateral de dispersión de corriente. Como en el tutorial de contactos de gran área, la carga se recoge primero lateralmente en la capa polimérica y luego se transfiere a la rejilla de plata de baja resistencia para una extracción eficiente.
La capa fotoactiva es PM6:Y6. Eléctricamente, esta capa se representa mediante un elemento diodo con generación óptica habilitada, mientras que componentes resistivos adicionales garantizan que el solucionador tenga en cuenta su conductividad finita. En la parte inferior de la pila hay un contacto continuo de Ag que cubre toda el área del dispositivo y completa la conexión eléctrica.
La barra verde visible en la ventana principal puede parecer flotante, pero es simplemente un indicador visual de la condición de contorno de extracción. Marca la región donde el solucionador aplica el contacto eléctrico utilizado para extraer corriente de la red de circuito.
3. Inspeccionar la pila de capas epitaxiales y los contactos
Esta simulación (y el ejemplo anterior de contacto de gran área) se define utilizando una estructura epitaxial en capas. El dispositivo se construye como una pila de capas de película delgada usando el editor de capas, lo que garantiza que cada capa se asiente limpiamente sobre la siguiente, sin huecos y con un tratamiento consistente del espesor. Este enfoque es muy adecuado para dispositivos planares de gran área, donde la pila vertical es fija y el comportamiento lateral se introduce a través del modelo de circuito. Se diferencia de simulaciones posteriores de tipo módulo (véase el tutorial de módulo interconectado), donde dedos, busbars y diseños más complejos se definen explícitamente usando geometría 3D completa.
Para inspeccionar la estructura, abra el Layer editor para ver la pila de capas (véase ??), y luego abra el panel Contacts en el editor de estructura del dispositivo (véase ??). Geométricamente, la configuración es simple pero importante: el contacto inferior cubre toda el área del dispositivo, mientras que el contacto superior es mucho más pequeño y define la región donde se extrae la corriente, en consonancia con el marcador verde de extracción mostrado en la vista 3D. Comprender esta geometría de contacto es esencial, ya que establece los caminos laterales de corriente que dominan el comportamiento del dispositivo de gran área.
4. Inspección de la malla
Cambie a la pestaña Circuit diagram y haga clic en el icono refresh de la parte inferior izquierda (el botón tipo reciclaje). Esta acción construye la representación en circuito del dispositivo directamente a partir de la geometría definida y de la pila de capas. La malla de circuito resultante se muestra en ??. Cada enlace de la malla corresponde a un elemento de circuito, como una resistencia o un diodo, y el color indica la capa o región de material de la que procede.
Al rotar la vista, puede identificar claramente las regiones de contacto superior e inferior. En ?? y ??, los puntos azules indican los nodos donde se extrae corriente de la malla de circuito. Estos nodos definen las condiciones de contorno del contacto y representan las ubicaciones en las que el solucionador de circuito recoge corriente del dispositivo.
5. Parámetros eléctricos y generación óptica
Abra Electrical parameters desde el panel de estructura del dispositivo. Este editor controla las propiedades de resistencias y diodos utilizadas por el modelo de circuito de gran área. Los ajustes de ejemplo se muestran en ?? y ??. Para la capa PM6:Y6 verá parámetros específicos de un elemento diodo, incluyendo el factor de idealidad n y la corriente en polarización inversa (de saturación) I0. Si sustituye estos valores por valores de sus propias curvas experimentales verá cómo se escalará su dispositivo. También puede habilitar la generación óptica de portadores de carga. Cuando la generación óptica está habilitada, el solucionador óptico de matriz de transferencia proporciona un término de fotogeneración que entra en el modelo eléctrico como una fotocorriente Iph en la ecuación del diodo iluminado:
I(V) = I0 [ exp( qV / (n k T) ) − 1 ] − Iph
En este tutorial, el conmutador controla efectivamente el término Iph. Con la generación óptica deshabilitada, el circuito se comporta como una red de diodos en oscuridad. Con ella habilitada, la iluminación inyecta corriente en el circuito a través del elemento diodo, con Iph determinada por el flujo de fotones absorbidos procedente de la óptica por matriz de transferencia.
💡 Detalle de implementación: la capa activa puede discretizarse en múltiples elementos resistivos a través de su espesor para capturar el transporte en volumen y la resistencia en serie. Sin embargo, el modelo de circuito siempre contiene una única capa de diodo, normalmente colocada en la parte inferior de la región activa. Esto refleja el hecho de que el diodo representa la respuesta JV global y la fotogeneración del dispositivo, mientras que los elementos resistivos tienen en cuenta la caída de tensión y la dispersión de corriente dentro del material.
6. Ejecución de la simulación
Para iniciar la simulación, haga clic en Run simulation (el triángulo azul de la barra de herramientas principal), o pulse F9. El solucionador comenzará inmediatamente a ejecutarse, y el terminal de salida se actualizará en tiempo real (véase ??).
Las líneas coloreadas que aparecen al principio de la salida corresponden al cálculo óptico. En esta etapa, el solucionador de matriz de transferencia está evaluando la fotogeneración a través del dispositivo, corte a corte, para construir el perfil de generación resuelto espacialmente. Una vez completado esto, el solucionador pasa al problema eléctrico y comienza a resolver la red de circuito.
Vale la pena prestar atención a la salida del terminal mientras se ejecuta la simulación. Esta información no es ruido diagnóstico: proporciona una visión útil de si la simulación se está comportando físicamente. En la instantánea mostrada, podemos ver que aproximadamente −9.48 × 102 A m−2 está fluyendo fuera de un contacto, mientras que −9.48 × 101 A m−2 está fluyendo fuera del otro. Al mismo tiempo, el residual reportado del solucionador F = 1.7 × 10−14 indica que la continuidad de corriente se ha satisfecho con una precisión extremadamente alta, con el primer punto JV resuelto en aproximadamente 17 ms.
La diferencia de densidad de corriente entre los dos contactos es totalmente esperable. El contacto inferior cubre toda el área del dispositivo, mientras que el contacto superior ocupa solo una pequeña región. Como resultado, la misma corriente total se extrae a través de un área mucho menor en el contacto superior, lo que conduce a una densidad de corriente significativamente mayor allí. Leer la salida del terminal de esta manera proporciona una comprobación rápida y eficaz de que la geometría y las condiciones de contorno de la simulación son físicamente consistentes. Si el error residual aumentara hasta orden unidad, o si el tiempo de ejecución se volviera inusualmente largo, normalmente indicaría un problema como una malla rota o un contacto desconectado. Una vez completada la simulación, abra la pestaña Output (véase ??). Las salidas disponibles se describen con más detalle en el tutorial anterior y se exploran más a fondo en el siguiente tutorial. Aquí, nos centramos solo en las curvas JV resueltas por contacto.
Cada contacto eléctrico tiene su propio archivo JV, jv_contact0.csv y jv_contact1.csv, que registran la corriente que sale de ese contacto específico. Aunque las densidades de corriente difieren, la corriente total extraída es la misma una vez que se tienen en cuenta las áreas de contacto. Abrir jv_contact0.csv produce la curva JV mostrada a continuación.
7. Edición de parámetros de material
Las propiedades de los materiales pueden modificarse directamente a través del editor de objetos. Haciendo clic con el botón derecho sobre cualquier capa en la vista 3D y seleccionando Edit object, se abre el editor de objetos (véanse ?? y ??).
Desde aquí, puede cambiar el material óptico asignado al objeto, modificar su forma geométrica (por ejemplo al utilizar contactos con patrón como mallas tipo panal), o ajustar atributos visuales como el color. En esta simulación de estilo epitaxial, la posición absoluta del objeto está en gran medida definida por la pila de capas, por lo que parámetros como traslación y rotación suelen tener poco efecto. Lo que más importa son la asignación del material óptico y la geometría del objeto.
💡 Pruebe esto: haga clic con el botón derecho sobre la capa PM6:Y6 y sustituya su material óptico por una perovskita. Luego vuelva a ejecutar la simulación y observe cómo cambia la curva JV. Si lo desea, también puede actualizar los parámetros del diodo — como el factor de idealidad, la corriente de saturación y la fotocorriente — usando valores reportados en la literatura sobre perovskitas.