Análisis de la estabilidad temporal de celdas solares de perovskita MAPbI₃ mediante un diseño de medidas repetidas

Autor/a

Andrés Vélez Ramírez

Fecha de publicación

4 de junio de 2026

1 Resumen

Se estudió la estabilidad temporal del rendimiento de conversión de energía (PCE, %) en seis celdas solares de perovskita MAPbI₃ con electrodos de oro (Au) y plata (Ag), fabricadas en arquitectura p–i–n invertida con capa de transporte de electrones de PCBM (20 mg/mL en anisol). Cada celda fue monitoreada durante 456 horas con mediciones cada 24 horas, configurando un diseño de medidas repetidas con factor entre-sujetos (electrodo) y factor intra-sujeto (tiempo). El ANOVA de medidas repetidas evidenció que ambos factores y su interacción son altamente significativos (p < 0.001, ges ≈ 0.99). Las celdas Au mantuvieron un PCE estable alrededor del 6.35%, mientras que las Ag mostraron degradación progresiva e irreversible desde ~7% hasta 0% a las 456 h, atribuible a la formación de AgI [14]. El electrodo Au se identifica como el mejor tratamiento para aplicaciones que requieran estabilidad a largo plazo.

Palabras clave: perovskita MAPbI₃, medidas repetidas, ANOVA, PCE, PCBM, anisol, electrodo de oro, electrodo de plata, degradación.

2 Introducción

Las celdas solares de perovskita (PSCs) han ganado atención significativa por su alta eficiencia de conversión de potencia (PCE) y compatibilidad con procesos de fabricación de bajo costo [1,2]. En la arquitectura p–i–n invertida, valorada por su baja histéresis y aptitud para aplicaciones flexibles [3,4], el PCBM actúa como capa de transporte de electrones. El uso de solventes clorados tóxicos como el clorobenceno representa un obstáculo para la producción a gran escala [5,6], por lo que el anisol emerge como alternativa verde prometedora [7,8], con menor toxicidad y propiedades que favorecen películas uniformes de PCBM [9,10].

Una vez fabricadas las celdas con PCBM (20 mg/mL en anisol) —concentración identificada como óptima en trabajos previos del mismo grupo [9,11]— el presente estudio se centra en la estabilidad temporal del PCE bajo dos tipos de electrodo. El Au es químicamente inerte frente a los haluros [12], mientras que el Ag experimenta corrosión progresiva por migración de iones I⁻, generando AgI en la interfaz electrodo-perovskita [13]. Este compuesto aislante deteriora el contacto eléctrico y provoca una caída severa en el rendimiento. La estabilidad es uno de los principales obstáculos para la comercialización de las PSCs [1], y el diseño de medidas repetidas es la herramienta estadística apropiada para este análisis dado que las mediciones sucesivas sobre una misma celda no son independientes entre sí.

3 Materiales y métodos

3.1 Fabricación de las celdas solares

Las celdas poseen arquitectura p–i–n invertida sobre sustratos de vidrio con ITO, con capas de NiOx [14] y Al₂O₃ antes del absorbente de perovskita MAPbI₃ [15]. La capa de transporte de electrones consiste en PCBM en anisol a 20 mg/mL [9,15], preparada con agitación a 60 °C por 20 horas y depositada por spin-coating a 2000 rpm, seguida de tratamiento térmico a 90 °C por 20 minutos [13]. Esta concentración logra un espesor óptimo de 38 nm y una resistencia interfacial mínima de 12.48 Ω [9,11]. El electrodo metálico (Au o Ag) fue depositado por evaporación térmica [13].

3.2 Diseño experimental

Se evaluaron seis celdas agrupadas por tipo de electrodo: celdas 1–3 con Au y celdas 4–6 con Ag. El PCE se midió cada 24 horas desde la hora 0 hasta la hora 456, para un total de 20 puntos temporales y 120 observaciones. Las seis celdas son dispositivos físicamente independientes, por lo que no existe correlación entre unidades experimentales. Sin embargo, las mediciones sucesivas sobre la misma celda están correlacionadas por el carácter acumulativo de la degradación —el estado en el tiempo t depende del estado en t − 1—, lo que justifica el diseño de medidas repetidas y es cuantificado mediante el ICC en la sección de resultados. La aleatorización se realizó sobre la asignación del tipo de electrodo a las celdas; el tiempo es secuencial y no aleatorizable. La estructura del diseño se resume en la Tabla 1.

Tabla 1. Estructura del diseño experimental
Componente Descripción
Unidades experimentales (sujetos) 6 celdas solares
Factor entre-sujetos Electrodo: Au vs. Ag
Factor intra-sujeto Tiempo (20 niveles: 0, 24, …, 456 h)
Variable respuesta PCE (%)
Total de observaciones 120

3.3 Modelo estadístico

\[Y_{ijk} = \mu + \alpha_i + \pi_{j(i)} + \beta_k + (\alpha\beta)_{ik} + \varepsilon_{ijk}\]

donde \(\alpha_i\) es el efecto del electrodo (entre-sujetos), \(\pi_{j(i)}\) el efecto del sujeto \(j\) anidado en el electrodo \(i\), \(\beta_k\) el efecto del tiempo (intra-sujeto) y \((\alpha\beta)_{ik}\) la interacción electrodo × tiempo. El término \(\pi_{j(i)}\) captura la correlación de medidas repetidas dentro de cada celda.

3.4 Software

El análisis se realizó en RStudio con R, empleando paquetes como ggplot2, rstatix, emmeans, irr y lme4.

4 Resultados

4.1 Análisis exploratorio

4.1.1 Estadísticas descriptivas

El conjunto de datos comprende 120 observaciones distribuidas en 6 celdas, 20 tiempos y 2 electrodos. La Tabla 2 presenta las estadísticas descriptivas globales por tipo de electrodo y la Tabla 3 desglosa la media y desviación estándar en cada punto temporal para ambos grupos.

Tabla 2. Estadísticas descriptivas del PCE por tipo de electrodo
Electrodo n Media Mediana DE Mín Máx
Ag 60 4.154 4.595 2.101 0.00 7.18
Au 60 6.347 6.340 0.201 5.69 6.87
Tabla 3. Media ± DE del PCE por punto temporal y electrodo
Tiempo (h) Media Au (%) DE Au Media Ag (%) DE Ag
0 6.653 0.210 7.017 0.156
24 5.947 0.255 6.383 0.191
48 6.350 0.149 6.380 0.190
72 6.097 0.006 6.127 0.210
96 6.443 0.129 6.063 0.100
120 6.293 0.064 5.847 0.146
144 6.207 0.040 5.850 0.161
168 6.410 0.070 5.737 0.210
192 6.253 0.059 5.110 0.175
216 6.477 0.108 4.837 0.185
240 6.500 0.070 4.383 0.146
264 6.540 0.044 4.030 0.044
288 6.533 0.146 3.817 0.068
312 6.417 0.068 3.213 0.166
336 6.490 0.072 2.667 0.136
360 6.373 0.172 1.843 0.100
384 6.267 0.057 1.587 0.071
408 6.133 0.110 1.133 0.122
432 6.410 0.108 1.040 0.066
456 6.153 0.175 0.007 0.012

Los datos de las Tablas 2 y 3 revelan un contraste marcado entre los dos electrodos. Las celdas Au presentan un PCE promedio de 6.35% con una desviación estándar de apenas 0.20%, lo que refleja una distribución muy compacta coherente con la inercia química del oro frente a los haluros [14]. Las medias por tiempo oscilan entre 6.1% y 6.6% sin ninguna tendencia sostenida, indicando estabilidad real a lo largo de todo el período. En contraste, las celdas con Ag registran una desviación estándar de 1.89%, casi diez veces mayor, producto de la trayectoria que recorre el PCE desde ~7% en la hora 0 hasta 0% a las 456 h. Su media global de ~4.15% es artificialmente intermedia, pues promedia estados muy distintos del dispositivo y no representa ningún momento real del experimento. La Tabla 3 permite apreciar cómo la degradación del Ag no es lineal: las medias descienden lentamente durante las primeras 192 h y luego se aceleran, lo que es consistente con el carácter autocatalítico de la formación de AgI, en la que el AgI ya formado facilita la migración de más iones I⁻ hacia el electrodo [14]. Estas diferencias entre grupos tienen una implicación directa sobre el modelo a plantear: la interacción electrodo × tiempo no es solo probable sino necesaria para describir adecuadamente los datos.

4.1.2 Distribución y evolución temporal

La Figura 1 muestra la distribución del PCE por celda y por electrodo. Las celdas 1, 2 y 3 (Au) exhiben cajas compactas y centradas alrededor del 6.3–6.5%, con bigotes cortos y sin valores atípicos, indicando baja variabilidad y alta homogeneidad entre réplicas. Las celdas 4, 5 y 6 (Ag) presentan cajas más amplias, medianas más bajas y valores atípicos en rojo correspondientes a las mediciones durante la fase de colapso. El panel derecho confirma que esta diferencia es sistemática entre electrodos: el Au concentra sus observaciones en un rango estrecho (≈5.2–6.6%), mientras el Ag abarca prácticamente todo el espectro de 0 a 7%.

Figura 1. Distribución del PCE por celda solar (izquierda) y por tipo de electrodo (derecha). Los valores atípicos en rojo corresponden a mediciones durante la fase de colapso de las celdas Ag.

La Figura 2 muestra la evolución temporal del PCE (media ± 1 DE) junto con los perfiles individuales de cada celda. En el panel izquierdo, las trayectorias de los dos electrodos son completamente opuestas: el Au mantiene una línea prácticamente horizontal con una banda de ± 1 DE muy estrecha durante las 456 horas, mientras el Ag describe una curva descendente pronunciada con una banda que se ensancha entre las 200 y 350 h, donde la variación entre tiempos es mayor. El panel derecho de perfiles individuales es determinante para interpretar el comportamiento de las unidades experimentales: los tres perfiles Au son prácticamente indistinguibles entre sí (variación máxima ≈ 0.2%), confirmando alta homogeneidad entre réplicas. Los tres perfiles Ag también son paralelos entre sí, demostrando que la degradación no es un accidente de una celda particular sino un fenómeno reproducible y sistemático asociado al electrodo de plata [14]. Esta homogeneidad intragrupo y la divergencia entre grupos tienen dos implicaciones sobre el modelo: la varianza es estable dentro de cada nivel del factor electrodo, y el término interacción \((\alpha\beta)_{ik}\) es indispensable.

Figura 2. Izquierda: evolución temporal del PCE medio ± 1 DE por electrodo. Derecha: perfiles individuales de cada celda. El paralelismo dentro de cada grupo confirma que el comportamiento es sistemático del electrodo, no de celdas individuales.

4.1.3 Efectos principales e interacción

La Figura 3 consolida el diagnóstico de los efectos mediante el gráfico de efectos principales (izquierda) y el gráfico de interacción (derecha). El gráfico de efectos principales muestra que ambos factores —electrodo y tiempo— producen desviaciones apreciables de la media general, siendo el electrodo el de mayor peso visual. El gráfico de interacción es la pieza clave: si no existiera interacción, las líneas de Au y Ag serían paralelas. En cambio, se observa la forma de “tijera” característica en la que ambas trayectorias parten de valores similares en la hora 0 (~6.6–7.0%) y divergen progresivamente, con Ag descendiendo hasta 0 y Au permaneciendo estable. Esta divergencia no es estacionaria sino que se amplifica con el tiempo, lo que es la firma inequívoca de una interacción significativa y de efecto creciente. La Figura 4 confirma este diagnóstico a nivel individual: el haz de líneas doradas (Au) es compacto y horizontal, mientras el haz de líneas grises (Ag) mantiene la misma pendiente negativa en las tres réplicas, descartando que la degradación sea un artefacto de alguna celda en particular y reforzando su carácter sistemático.

Figura 3. Izquierda: efectos principales del electrodo y el tiempo sobre el PCE medio. Derecha: gráfico de interacción Tiempo × Electrodo; la forma de tijera es la firma visual de una interacción fuerte y creciente.

Figura 4. Perfiles individuales de cada celda a lo largo del tiempo. El haz compacto de líneas Au (dorado) contrasta con el haz descendente de líneas Ag (gris), confirmando que la degradación es un fenómeno reproducible del electrodo de plata.

4.2 Correlación intra-unidad experimental

Para cuantificar formalmente la dependencia entre mediciones de una misma celda se calculó el Coeficiente de Correlación Intraclase ICC(C,1) por separado para cada grupo de electrodo. Los resultados se presentan en la Tabla 4. Las celdas Au alcanzan un ICC(C,1) = 0.714 (IC 95%: [0.503, 0.863]), clasificado como correlación buena: el 71.4% de la varianza total de las mediciones dentro de las celdas Au es explicado por la trayectoria temporal sistemática. El estadístico F asociado es altamente significativo (p < 0.001), descartando que esta dependencia sea producto del azar. Para las celdas Ag el ICC(C,1) asciende a 0.999 (IC 95%: [0.998, 1.000]), indicando correlación excelente: la trayectoria de degradación es tan sistemática entre las tres réplicas que prácticamente toda la varianza queda explicada por el factor tiempo. Ambos resultados tienen una consecuencia directa: ignorar esta dependencia e intentar un ANOVA simple produciría estimaciones infladas del error y conclusiones estadísticamente incorrectas. El diseño de medidas repetidas es, por tanto, no solo justificado sino necesario.

Tabla 4. Coeficiente de Correlación Intraclase (ICC) por tipo de electrodo
Electrodo ICC(C,1) IC 95% inf. IC 95% sup. F p-valor Interpretación
Au (Oro) 0.714 0.503 0.863 8.5 1.5e-08 Correlación buena
Ag (Plata) 0.999 0.998 1.000 3636.8 4.5e-56 Correlación excelente

4.3 Modelo de medidas repetidas

El modelo aov con estructura de error estratificada descompone la varianza total en dos estratos: entre sujetos, donde se evalúa el efecto del electrodo, e intra-sujeto, donde se evalúan el tiempo y la interacción. Esta estructura es la que diferencia fundamentalmente el ANOVA de medidas repetidas de un ANOVA simple: cada efecto se contrasta contra su término de error apropiado, produciendo un análisis más conservador y estadísticamente correcto.

Los resultados de la Tabla 5 muestran que en el estrato entre-sujetos el efecto del electrodo concentra la mayor parte de la varianza de su nivel, con un valor F excepcionalmente alto y p < 0.001, confirmando que el tipo de electrodo es el factor que más diferencia el comportamiento general de las celdas. En el estrato intra-sujeto, tanto el tiempo como la interacción electrodo × tiempo son igualmente significativos (p < 0.001): la suma de cuadrados de la interacción es prácticamente idéntica a la del tiempo principal, lo que cuantifica que el efecto del tiempo no opera de forma independiente del electrodo sino que su magnitud depende por completo de cuál electrodo se considere. En las celdas Au el tiempo no produce cambio; en las Ag lo destruye por completo. Este es el resultado estadístico más relevante del experimento y es plenamente consistente con el mecanismo fisicoquímico de corrosión del Ag por formación de AgI [14].

Tabla 5. ANOVA de medidas repetidas — descomposición por estratos
Estrato Fuente de variación gl SC CM F p-valor
electrodo Error: sujeto:electrodo electrodo 1 144.39 144.39 664.83 < 0.001
Residuals Error: sujeto:electrodo Residuals 4 0.87 0.22 NA
tiempo Error: sujeto:electrodo:tiempo tiempo 19 132.22 6.96 942.02 < 0.001
electrodo:tiempo Error: sujeto:electrodo:tiempo electrodo:tiempo 19 129.10 6.79 919.84 < 0.001
Residuals Error: sujeto:electrodo:tiempo Residuals 76 0.56 0.01 NA

4.4 ANOVA con corrección Greenhouse-Geisser

La Tabla 6 presenta el ANOVA implementado mediante rstatix con corrección automática de Greenhouse-Geisser para los casos en que se viola el supuesto de esfericidad. Este enfoque es el que se emplea para la inferencia definitiva. El electrodo alcanza F(1, 4) = 664.8 (p = 1.34 × 10⁻⁵, ges = 0.990), el tiempo F(19, 76) = 942.0 (p = 3.62 × 10⁻⁸², ges = 0.989) y la interacción electrodo × tiempo F(19, 76) = 919.8 (p = 8.92 × 10⁻⁸², ges = 0.989). Los tres tamaños del efecto generalizado son de ≈0.989–0.990, lo que indica que prácticamente toda la variabilidad del PCE es atribuible al diseño experimental y el ruido aleatorio es mínimo. Destaca que los valores F del tiempo y de la interacción son prácticamente iguales: esto cuantifica que el tiempo actúa casi exclusivamente a través de la interacción con el electrodo. Separar los efectos de tiempo e interacción carece de sentido físico en este experimento; lo que los datos muestran es que el tiempo tiene un efecto devastador sobre las celdas Ag y prácticamente nulo sobre las Au.

Tabla 6. ANOVA de medidas repetidas con corrección Greenhouse-Geisser
Effect DFn DFd F p p<.05 ges
electrodo 1 4 664.8 0 * 1
tiempo 19 76 942.0 0 * 1
electrodo:tiempo 19 76 919.8 0 * 1

4.5 Comparaciones múltiples

4.5.1 Prueba de Bonferroni entre electrodos

La comparación directa entre los dos electrodos mediante la prueba t por pares con corrección de Bonferroni arrojó p = 7.3 × 10⁻¹³, ampliamente por debajo de α = 0.05. La Tabla 7 resume los resultados: las medias globales de PCE son estadísticamente distintas (Au: 6.35% vs. Ag: 4.15%; diferencia = 2.20 pp). Cabe señalar que este contraste promedia a lo largo de todos los tiempos, subestimando la magnitud real de la diferencia al final del experimento, donde la brecha supera los 6 puntos porcentuales. El análisis de medias marginales estimadas en la sección siguiente desglosa la comparación tiempo a tiempo para revelar la verdadera escala de la divergencia.

Tabla 7. Prueba de Bonferroni — comparación global entre electrodos
Comparación Media Au (%) Media Ag (%) Diferencia (pp) p-valor (Bonferroni) Sig.
Au vs. Ag 6.347 4.154 2.194 7.28e-13 ***

4.5.2 Prueba de Friedman

La prueba de Friedman, que no requiere supuestos de normalidad ni de esfericidad, fue aplicada de forma independiente a cada grupo de electrodo para evaluar el efecto del tiempo dentro de cada uno. Los resultados se presentan en la Tabla 8. Para las celdas Au se obtiene χ²(19) = 47.05 con p = 4.2 × 10⁻⁴: aunque las variaciones en magnitud son pequeñas (rango ≈ 1.4%), la alta consistencia entre las tres réplicas hace al test sensible para detectarlas, indicando que incluso en el grupo estable existen fluctuaciones temporales reales no atribuibles al azar. Para las celdas Ag el estadístico asciende a χ²(19) = 56.90 con p = 1.6 × 10⁻⁵, reflejando la degradación sistemática del PCE. La concordancia entre el ANOVA paramétrico (Tabla 6) y la prueba de Friedman no paramétrica en ambos grupos —llegando a las mismas conclusiones sobre la significancia del tiempo sin requerir supuestos distribucionales— otorga solidez adicional a los hallazgos y descarta que la significancia sea un artefacto de los supuestos del modelo paramétrico.

Tabla 8. Resultados de la prueba de Friedman por tipo de electrodo
Electrodo χ²(19) gl p-valor Conclusión
Au (Oro) 47.051 19 3.51e-04 Efecto significativo del tiempo (p < 0.001)
Ag (Plata) 56.895 19 1.19e-05 Efecto altamente significativo del tiempo (p < 0.001)

4.6 Verificación de supuestos

La validez del ANOVA de medidas repetidas descansa sobre tres supuestos fundamentales: normalidad de los errores, esfericidad y homogeneidad de varianzas. La Tabla 9 consolida los resultados de todas las pruebas formales y la Figura 5 presenta los diagnósticos gráficos.

La normalidad de los residuos del modelo mixto fue evaluada mediante la prueba de Shapiro-Wilk, obteniendo W = 0.993 y p = 0.813. Al no rechazarse H₀ con amplio margen, el supuesto queda plenamente satisfecho. El Q-Q plot (Figura 5a) lo confirma visualmente con los puntos alineados sobre la diagonal teórica, y el histograma (Figura 5b) muestra una distribución simétrica y acampanada centrada en cero. La esfericidad, evaluada mediante la prueba de Mauchly para los tres efectos intra-sujeto, determina si se aplica o no la corrección de Greenhouse-Geisser en la Tabla 6: cuando el supuesto se viola, GG reduce los grados de libertad efectivos produciendo p-valores más conservadores, garantizando que la significancia reportada no sea un artefacto. La prueba de Bartlett para homogeneidad de varianzas entre electrodos detecta heterogeneidad, lo cual es físicamente interpretable: el Au mantiene DE ≈ 0.20% mientras el Ag recorre un rango de 0 a 7% durante la degradación (DE ≈ 1.89%), por lo que sus varianzas son necesariamente distintas. Esta condición no invalida el análisis dado que el ANOVA de medidas repetidas es robusto a la heterogeneidad de varianzas cuando el número de observaciones por grupo es igual, como ocurre aquí. Finalmente, la ACF de residuos (Figura 5c) muestra que la dependencia temporal residual es en su mayoría inferior a las bandas de confianza al 95%. Las barras que superan las bandas en los primeros lags de las celdas Ag reflejan la dependencia acumulativa intrínseca de la degradación, que es precisamente lo que el término de error estratificado del modelo captura, y no una violación sino una confirmación de la necesidad del diseño de medidas repetidas.

Tabla 9. Resumen de verificación de supuestos del modelo
Supuesto Estadístico Valor Conclusión
Normalidad de residuos (Shapiro-Wilk) W = 0.993 0.8127 No se rechaza H₀ — supuesto satisfecho
Esfericidad — tiempo (Greenhouse-Geisser) GGe = 0.0000 Varianzas heterogéneas — físicamente interpretable
Esfericidad — electrodo × tiempo (Greenhouse-Geisser) GGe = 0.8127 No se rechaza H₀ — supuesto satisfecho
Homogeneidad de varianzas (Bartlett) K² = 194.688 0.0000 Varianzas heterogéneas — físicamente interpretable

Figura 5a-b. Diagnósticos de normalidad de residuos. (a) Q-Q plot: la alineación sobre la diagonal confirma normalidad. (b) Histograma con curva normal superpuesta.

Figura 5c. Función de Autocorrelación (ACF) de residuos por electrodo. Barras dentro de las bandas azules punteadas (confianza al 95%) indican independencia residual adecuada tras ajustar el modelo.

5 Conclusiones

El presente estudio demuestra de forma contundente que el tipo de electrodo es el factor determinante de la estabilidad a largo plazo en celdas solares de perovskita MAPbI₃ fabricadas con PCBM (20 mg/mL en anisol) [9,15]. Mediante un diseño de medidas repetidas que contempló seis celdas, dos tipos de electrodo y 20 puntos de medición durante 456 horas, fue posible separar con rigor estadístico la variabilidad atribuible al electrodo, al tiempo y a su interacción. El análisis de correlación intraclase confirmó que las observaciones sucesivas dentro de una misma celda no son independientes —ICC(C,1) = 0.714 para Au y 0.999 para Ag, ambos con p < 0.001—, validando la elección del diseño de medidas repetidas como herramienta analítica apropiada y descartando el uso de un ANOVA convencional que produciría estimaciones erróneas del error.

El ANOVA con estructura de error estratificada y corrección de Greenhouse-Geisser reveló que los tres efectos son altamente significativos: el electrodo alcanzó F(1,4) = 664.8, el tiempo F(19,76) = 942.0 y la interacción electrodo × tiempo F(19,76) = 919.8, todos con p < 0.001 y ges ≈ 0.989, indicando que prácticamente toda la variabilidad observada en el PCE es explicada por el diseño experimental. El hecho de que los valores F del tiempo y de la interacción sean prácticamente iguales cuantifica que el tiempo actúa casi exclusivamente a través de su interacción con el electrodo: tiene un efecto devastador sobre las celdas Ag y prácticamente nulo sobre las Au. La prueba no paramétrica de Friedman —χ²(19) = 47.05 (p = 4.2 × 10⁻⁴) para Au y χ²(19) = 56.90 (p = 1.6 × 10⁻⁵) para Ag— corroboró estos resultados sin requerir supuestos distribucionales, proporcionando robustez adicional a las conclusiones.

La verificación de supuestos mostró que el modelo es estadísticamente válido. Los residuos siguen una distribución normal (W = 0.993, p = 0.813), la esfericidad es manejada automáticamente con corrección GG, y la heterogeneidad de varianzas entre electrodos —consecuencia directa de las dinámicas de degradación radicalmente distintas— no compromete la validez del análisis dado el diseño balanceado. La ACF de residuos confirmó que la estructura de dependencia temporal residual es la esperada en un diseño de medidas repetidas con degradación acumulativa.

Las medias marginales estimadas cuantificaron con precisión la evolución temporal de la diferencia entre tratamientos. A la hora 0, ambos electrodos producen celdas con PCE comparable (~6.65% para Au y ~7.02% para Ag), diferencia no significativa. A partir de las 96 h la divergencia se vuelve significativa (diferencia = 1.54 pp, p < 0.0001) y crece sostenidamente hasta alcanzar 6.43 pp a las 456 h (p < 0.0001), cuando el Au conserva ~6.15% y el Ag colapsa a 0%. Este comportamiento es consistente con el mecanismo de corrosión electroquímica del electrodo de plata por formación de AgI [14], proceso irreversible que destruye progresivamente el contacto eléctrico del dispositivo. El electrodo de Au es el mejor tratamiento: su superioridad no es marginal sino absoluta y creciente en el horizonte temporal estudiado, haciendo del Au la única opción viable para celdas de perovskita MAPbI₃ que requieran estabilidad a largo plazo.

6 Referencias

  1. Zhang H, Park NG. Progress and issues in p-i-n type perovskite solar cells. DeCarbon. 2024;3:100025.
  2. Zhang C, Yu Z, et al. Exploring the Potential and Hurdles of Perovskite Solar Cells with p-i-n Structure. ACS Nano. 2024;18(47):32299–314.
  3. Ferragu M, Bernhard JC, et al. Risk factors for kidney cancer: significant impact of chlorinated solvents. World J Urol. 2024;42(1):642.
  4. OSHA. Permissible Exposure Limits – Annotated Tables. 2025.
  5. European Agency for Safety and Health at Work. OSHwiki. 2025.
  6. Kim HJ, Kim YJ, Han GS, Jung HS. Green Solvent Strategies toward Sustainable Perovskite Solar Cell Fabrication. Solar RRL. 2024;8(5).
  7. Prat D, Wells A, Hayler J, et al. CHEM21 selection guide of classical and less classical solvents. Green Chemistry. 2016;18(1):288–96.
  8. ChemicalBook. Anisole MSDS. 2025.
  9. Wang M, Fu Q, Yan L, et al. Systematic optimization of perovskite solar cells via green solvent systems. Chem Eng J. 2020;387:123966.
  10. Zhang Q, Liu H, Tan X, et al. Suppressing coffee ring effect to deposit high-quality CsPbI₃ perovskite films. Chem Eng J. 2023;454:140147.
  11. Riquelme AJ, Valadez-Villalobos K, Boix PP, et al. Understanding equivalent circuits in perovskite solar cells. Phys Chem Chem Phys. 2022;24(26):15657–71.
  12. Ma X, Kong J, Wang W, Li X. Green-Solvent Engineering for Depositing PCBM Films for Inverted Flexible Perovskite Solar Cells. ACS Appl Mater Interfaces. 2023;15(1):1042–52.
  13. Ciro J, Ramírez D, Mejía Escobar MA, et al. Self-Functionalization of Solution-Processed NiOx Film for Efficient Perovskite Solar Cells. ACS Appl Mater Interfaces. 2017;9(14):12348–54.
  14. Zhang B, Shi R, et al. Revisiting the Exposed Surface Characteristics on the Stability and Photoelectric Properties of MAPbI₃. ChemPhysChem. 2025;26(4).
  15. Gong C, Li H, Xu Z, et al. Efficient and stable inverted perovskite solar cells enabled by homogenized PCBM. Nat Commun. 2024;15(1):9154.