Predicting Organic Solar Cell Performance and Stability from Fast, Morphology-aware Current-Voltage Modeling
Pith reviewed 2026-05-20 04:03 UTC · model grok-4.3
The pith
A fast method predicts current-voltage curves of organic solar cells from their complex multi-phase morphologies using morphology-aware descriptors.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors establish a unique method capable of calculating the current-voltage (JV) curve of complex heterojunction morphologies containing up to five phases with a very low computation time using morphology-aware descriptors of light absorption, exciton dissociation, non-geminate recombination and free charge carrier mobilities. This is validated against Monte Carlo and 3D drift-diffusion simulations and applied to P3HT:PCBM and PM6:Y6 systems to shed light on physical compromises for optimizing device performance and lifetime, with the finding that the morphology-performance relationship depends on the materials system studied.
What carries the argument
morphology-aware descriptors of light absorption, exciton dissociation, non-geminate recombination and free charge carrier mobilities that enable fast JV curve calculation for multi-phase morphologies
If this is right
- The method supports analysis of morphologies with up to five phases: donor amorphous, donor crystalline, acceptor amorphous, acceptor crystalline, and mixed amorphous.
- Validation shows agreement with more computationally intensive Monte Carlo and 3D drift-diffusion simulations.
- Application to P3HT:PCBM and PM6:Y6 reveals specific physical compromises between performance and stability.
- The relationship between morphology and device metrics varies depending on the materials system under consideration.
Where Pith is reading between the lines
- Designers could use this to rapidly iterate on morphology to balance efficiency and longevity in new solar cell formulations.
- Similar descriptor-based approaches might accelerate modeling in related fields like perovskite or organic light-emitting devices with phase separation.
- Further validation on experimental devices with controlled multi-phase structures would strengthen the predictive power.
Load-bearing premise
The morphology-aware descriptors accurately capture the combined effects of light absorption, exciton dissociation, recombination, and carrier mobilities across all phases without needing full microscopic simulation for each new morphology.
What would settle it
A significant mismatch between the method's predicted JV curves and direct experimental measurements on a solar cell with a verified five-phase morphology would indicate the descriptors fail to capture the necessary physics.
read the original abstract
Understanding the relationship between morphology and performance in organic solar cells is essential for developing devices that are both high performing and resilient to aging. This work introduces a unique method capable of calculating the current-voltage (JV) curve of complex heterojunction morphologies containing up to five phases (donor amorphous, donor crystalline, acceptor amorphous, acceptor crystalline, mixed amorphous) with a very low computation time using morphology-aware descriptors of light absorption, exciton dissociation, non-geminate recombination and free charge carrier mobilities. The method is validated against Monte Carlo and 3D drift-diffusion simulations and applied to P3HT:PCBM and PM6:Y6 systems, shedding light on the physical compromises encountered to optimize device performance and lifetime. Finally, we show that the morphology-performance relationship is dependent on the materials system studied.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a method to compute current-voltage curves for organic solar cell heterojunctions containing up to five phases (donor amorphous, donor crystalline, acceptor amorphous, acceptor crystalline, mixed amorphous) at low computational cost. It employs morphology-aware descriptors for light absorption, exciton dissociation, non-geminate recombination, and free charge carrier mobilities. The approach is validated against Monte Carlo and 3D drift-diffusion simulations on P3HT:PCBM and PM6:Y6 systems, then applied to examine morphology-performance and stability trade-offs, with the conclusion that these relationships are materials-system dependent.
Significance. If the descriptors prove transferable and the validation holds with quantitative accuracy, the method could enable rapid screening of complex multi-phase morphologies for both efficiency and lifetime in organic photovoltaics, a longstanding bottleneck in the field. The explicit demonstration of system dependence strengthens the case for morphology-aware modeling over generic effective-medium approximations.
major comments (3)
- [Abstract] Abstract: validation is stated against Monte Carlo and 3D drift-diffusion simulations, yet no quantitative error metrics (RMSE on JV curves, predicted PCE deviation, or similar) are supplied. Without these numbers it is impossible to judge whether the low-computation-time advantage is achieved at acceptable accuracy loss for the central predictive claim.
- [Method description (paragraph on method)] Method description (paragraph on method): the claim that the small set of morphology-aware descriptors captures the combined effects of absorption, dissociation, recombination and mobilities across all five phases without requiring full microscopic simulation for each new morphology rests on unshown transferability. No test on deliberately varied fifth-phase fractions, new interface geometries, or out-of-sample morphologies is described, which is load-bearing for the assertion that the method generalizes beyond the two validation systems.
- [Abstract] Abstract: descriptors for mobilities and recombination are presented as morphology-aware, but the text does not clarify whether their extraction involves any fitting or calibration to the same JV data used for validation. If any key parameter is adjusted to the target curves, the reported 'prediction' becomes a fitted interpolation by construction and undermines the claimed advantage over direct simulation.
minor comments (2)
- [Abstract] Abstract: the statement that the method 'sheds light on the physical compromises encountered to optimize device performance and lifetime' would benefit from a one-sentence indication of how stability/aging is incorporated into the descriptor framework.
- Ensure that any timing comparisons with Monte Carlo or drift-diffusion solvers are reported with explicit wall-clock values and hardware specifications so readers can reproduce the claimed computational savings.
Simulated Author's Rebuttal
We thank the referee for their constructive comments on our manuscript. We address each major point below and will revise the manuscript to improve clarity on validation metrics, transferability, and parameter derivation.
read point-by-point responses
-
Referee: [Abstract] Abstract: validation is stated against Monte Carlo and 3D drift-diffusion simulations, yet no quantitative error metrics (RMSE on JV curves, predicted PCE deviation, or similar) are supplied. Without these numbers it is impossible to judge whether the low-computation-time advantage is achieved at acceptable accuracy loss for the central predictive claim.
Authors: We agree that quantitative error metrics are necessary to evaluate the accuracy of our approach. In the revised manuscript, we will add RMSE values for the JV curves and deviations in predicted PCE for comparisons against both Monte Carlo and 3D drift-diffusion simulations on the P3HT:PCBM and PM6:Y6 systems. revision: yes
-
Referee: [Method description (paragraph on method)] Method description (paragraph on method): the claim that the small set of morphology-aware descriptors captures the combined effects of absorption, dissociation, recombination and mobilities across all five phases without requiring full microscopic simulation for each new morphology rests on unshown transferability. No test on deliberately varied fifth-phase fractions, new interface geometries, or out-of-sample morphologies is described, which is load-bearing for the assertion that the method generalizes beyond the two validation systems.
Authors: We acknowledge that explicit tests on varied fifth-phase fractions, new interface geometries, and out-of-sample morphologies would strengthen the transferability claim. The descriptors are derived from general physical principles applicable to multi-phase systems, but the current manuscript focuses on the two validation cases. We will add a dedicated discussion on expected transferability and the design of the descriptors in the revised version; new computational tests can be included if they do not require prohibitive additional resources. revision: partial
-
Referee: [Abstract] Abstract: descriptors for mobilities and recombination are presented as morphology-aware, but the text does not clarify whether their extraction involves any fitting or calibration to the same JV data used for validation. If any key parameter is adjusted to the target curves, the reported 'prediction' becomes a fitted interpolation by construction and undermines the claimed advantage over direct simulation.
Authors: We thank the referee for this clarification request. The descriptors for mobilities and recombination are extracted directly from morphological features and material properties using independent physical relations, with no fitting or calibration performed against the JV curves used for validation. The comparisons with Monte Carlo and drift-diffusion simulations therefore represent genuine out-of-sample predictions. We will add explicit statements clarifying this procedure in the method section and abstract of the revised manuscript. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper presents a method that computes JV curves for up to five-phase morphologies via precomputed morphology-aware descriptors for absorption, dissociation, recombination and mobilities, then validates the outputs against independent Monte Carlo and 3D drift-diffusion simulations on P3HT:PCBM and PM6:Y6. No equations or sections in the supplied text show a parameter fitted to the target JV data and then relabeled as a prediction, nor any self-definitional loop where an output is used to define its own input. The validation step supplies external benchmarks rather than reducing the result to the descriptors by construction, so the derivation chain remains self-contained.
Axiom & Free-Parameter Ledger
free parameters (1)
- mobility and recombination parameters
axioms (1)
- domain assumption The chosen morphology descriptors sufficiently represent the physics of charge generation and transport in multi-phase blends.
Reference graph
Works this paper leans on
-
[1]
Validation of morphology-aware electronic descriptors To validate our simplified approach for evaluating morphology -dependent electronic properties , the values for exciton dissociation, bimolecular recombination, and mobilities have been calculated using Equations 6, 7, 13 and 14 above and compared with previously published 15 simulation results. We ide...
-
[2]
[4] [5] [6] [7] [9] [11] [12] [14] [15] [37] [42] where the simula tion results can be compared to the present approach. [4] [6] [7] [13] They discuss different morphological descriptors and types of morphology, but all these studies treat only two -phase morphologies (either amorphous demixed donor - acceptor BHJ or amorphous -crystalline single -materia...
-
[3]
Results 4.1. Morphology-dependent JV curves and their evolution upon thermal loading To illustrate the interest and the power of the proposed approach, we apply the method to a set of various organic photoactive layer BHJ morphologies, for which we identify the phase distribution, calculate the optical and electronic properties, and finally the JV curve o...
-
[4]
Conclusions With this work, we introduced a new method to study the morphology-performance relationship of OSC. The approach relies on a topological analysis of the BHJ morphology and on the basic physics of the optoelectronic processes involved in power conversion , namely light absorption, free charge carrier generation, transport and recombination. Com...
-
[5]
Conflict of interests There are no conflicts to declare
-
[6]
Data availability The simulation data used for this article is publicly accessible on Zenodo (see DOI:10.5281/zenodo.20272551)
-
[7]
Acknowledgements The authors acknowledge financial support by the German Research Foundation (DFG, project HA 4382/14-1, project BR 4031/21-2 and within the Collaborative Research Center ‘ChemPrint’ , project 538767711, CRC 1719), the European Commission (H2020 Program, Project 101008701/ EMERGE), the Helmholtz Association (SolarTAP Innovation Platform) ,...
-
[8]
O. Wodo, S. Tirthapura, S. Chaudhary, B. Ganapathysubramanian, Organic Electronics 2012, 13, 1105
work page 2012
-
[9]
X. Xu, Y . Ji, C. Du, T. Hou, Y . Li, RSC Advances 2015, 5, 70939
work page 2015
-
[10]
P . K. Watkins, A. B. Walker, G. L. B. Verschoor, Nano Lett. 2005, 5, 1814
work page 2005
-
[12]
S. Khodakarimi, M. H. Hekmatshoar, F. Abbasi, Applied Physics A 2016, 122
work page 2016
- [15]
-
[16]
E. Kawashima, M. Fujii, K. Yamashita, Phys. Chem. Chem. Phys. 2016, 18, 26456
work page 2016
- [17]
-
[18]
R. A. Marsh, C. Groves, N. C. Greenham, Journal of Applied Physics 2007, 101, 083509
work page 2007
- [19]
-
[20]
Y . Geng, Y . Han, L. Xiong, H. Li, Organic Electronics 2024, 127, 107006
work page 2024
-
[21]
E. Van, M. Jones, E. Jankowski, O. Wodo, Mol. Syst. Des. Eng. 2018, 3, 853
work page 2018
-
[22]
M. C. Heiber, K. Kister, A. Baumann, V. Dyakonov, C. Deibel, T.-Q. Nguyen, Phys. Rev. Applied 2017, 8, 054043
work page 2017
-
[23]
L. Ye, H. Hu, M. Ghasemi, T. Wang, B. A. Collins, J.-H. Kim, K. Jiang, J. H. Carpenter, H. Li, Z. Li, T. McAfee, J. Zhao, X. Chen, J. L. Y . Lai, T. Ma, J.-L. Bredas, H. Yan, H. Ade, Nature Mater 2018, 17, 253
work page 2018
- [24]
-
[25]
L. Ye, S. Li, X. Liu, S. Zhang, M. Ghasemi, Y . Xiong, J. Hou, H. Ade, Joule 2019, 3, 443. 26
work page 2019
-
[26]
J. A. Bartelt, Z. M. Beiley, E. T. Hoke, W. R. Mateker, J. D. Douglas, B. A. Collins, J. R. Tumbleston, K. R. Graham, A. Amassian, H. Ade, J. M. J. Fréchet, M. F. Toney, M. D. McGehee, Advanced Energy Materials 2013, 3, 364
work page 2013
-
[27]
Z. Wang, Z. Peng, Z. Xiao, D. Seyitliyev, K. Gundogdu, L. Ding, H. Ade, Advanced Materials 2020, 32, 2005386
work page 2020
-
[28]
M. B. Isichenko, Rev. Mod. Phys. 1992, 64, 961
work page 1992
-
[29]
C. H. Y . Ho, S. H. Cheung, H.-W. Li, K. L. Chiu, Y . Cheng, H. Yin, M. H. Chan, F. So, S.-W. Tsang, S. K. So, Advanced Energy Materials 2017, 7, 1602360
work page 2017
-
[30]
K. Vakhshouri, D. R. Kozub, C. Wang, A. Salleo, E. D. Gomez, Physical Review Letters 2012, 108
work page 2012
-
[31]
R. Wang, Y . Jiang, W. Gruber, Y . He, M. Wu, P . Weitz, K. Zhang, L. Lüer, K. Forberich, T. Unruh, E. Spiecker, C. Deibel, N. Li, C. J. Brabec, Advanced Materials Interfaces 2022, 9, 2200342
work page 2022
-
[32]
J. Liu, Y . Yin, K. Wang, P . Wei, H. Lu, C. Song, Q. Liang, W. Huang, iScience 2022, 25
work page 2022
-
[33]
S. M. Hosseini, N. Tokmoldin, Y . W. Lee, Y . Zou, H. Y . Woo, D. Neher, S. Shoaee, Solar RRL 2020, 4, 2000498
work page 2020
-
[34]
L. Lüer, I. M. Peters, A. S. Smith, E. Dorschky, B. M. Eskofier, F. Liers, J. Franke, M. Sjarov, M. Brossog, D. M. Guldi, A. Maier, C. J. Brabec, Joule 2024, 8, 295
work page 2024
-
[35]
T. Osterrieder, F. Schmitt, L. Lüer, J. Wagner, T. Heumüller, J. Hauch, C. J. Brabec, Energy Environ. Sci. 2023, 16, 3984
work page 2023
- [36]
-
[37]
X. Rodríguez-Martínez, E. Pascual-San-José, M. Campoy-Quiles, Energy Environ. Sci. 2021, 14, 3301
work page 2021
-
[39]
H. K. Kodali, B. Ganapathysubramanian, Solar Energy Materials and Solar Cells 2013, 111, 66
work page 2013
-
[40]
M. C. Heiber, A. Dhinojwala, Physical Review Applied 2014, 2
work page 2014
-
[41]
Y . Min Nam, J. Huh, W. Ho Jo, Solar Energy Materials and Solar Cells 2010, 94, 1118
work page 2010
-
[42]
U. Neupane, B. Bahrami, M. Biesecker, M. F. Baroughi, Q. Qiao, Nano Energy 2017, 35, 128
work page 2017
-
[43]
G. A. Nemnes, S. Iftimie, A. Palici, A. Nicolaev, T. L. Mitran, A. Radu, S. Antohe, Applied Surface Science 2017, 424, 264
work page 2017
-
[44]
O. Wodo, J. Zola, B. S. Sarath Pokuri, P . Du, B. Ganapathysubramanian, Materials Discovery 2015, 1, 21
work page 2015
-
[45]
F. Kaka, R. K. Singh, P . C. Ramamurthy, A. Choudhury, AIP Advances 2020, 10, 065304
work page 2020
-
[46]
F. Kaka, S. Khanna, P . C. Ramamurthy, A. Choudhury, Journal of Applied Physics 2020, 128, 145501
work page 2020
-
[47]
I. Maqsood, L. D. Cundy, M. Biesecker, J.-H. Kimn, E. Darlington, E. P . Hettwer, S. Schill, V. Bommisetty, In 2014 IEEE 40th Photovoltaic Specialist Conference (PVSC), 2014, pp. 1758–1761
work page 2014
-
[48]
K. Bergermann, C. Deibel, R. Herzog, R. C. I. MacKenzie, J.-F. P . Null, M. Stoll, CiCP 2023, 34, 1
work page 2023
- [49]
-
[50]
R. S. Gebhardt, P . Du, O. Wodo, B. Ganapathysubramanian, Computational Materials Science 2017, 129, 220
work page 2017
-
[51]
H. Li, G. Sini, J. Sit, A. J. Moulé, J.-L. Bredas, Energy Environ. Sci. 2020, 13, 601
work page 2020
-
[52]
R. G. E. Kimber, E. N. Wright, S. E. J. O’Kane, A. B. Walker, J. C. Blakesley, Phys. Rev. B 2012, 86, 235206
work page 2012
-
[53]
O. J. J. Ronsin, J. Harting, ACS Appl. Mater. Interfaces 2022, 14, 49785−49800
work page 2022
-
[54]
O. J. J. Ronsin, J. Harting, Advanced Theory and Simulations 2022, 2200286
work page 2022
-
[55]
F. C. Spano, The Journal of Chemical Physics 2005, 122, 234701
work page 2005
-
[56]
F. C. Spano, Acc. Chem. Res. 2010, 43, 429
work page 2010
-
[57]
R. Wang, L. Lüer, S. Langner, T. Heumueller, K. Forberich, H. Zhang, J. Hauch, N. Li, C. J. Brabec, ChemSusChem 2021, 14, 3590
work page 2021
-
[58]
X. Du, L. Lüer, T. Heumueller, J. Wagner, C. Berger, T. Osterrieder, J. Wortmann, S. Langner, U. Vongsaysy, M. Bertrand, N. Li, T. Stubhan, J. Hauch, C. J. Brabec, Joule 2021, 5, 495. 27
work page 2021
-
[59]
L. A. A. Pettersson, L. S. Roman, O. Inganäs, Journal of Applied Physics 1999, 86, 487
work page 1999
-
[60]
R. Kerremans, C. Kaiser, W. Li, N. Zarrabi, P . Meredith, A. Armin, Advanced Optical Materials 2020, 8, 2000319
work page 2020
-
[61]
C. F. Bohren, Eur. J. Phys. 2010, 31, 573
work page 2010
-
[62]
M. Campoy-Quiles, C. Müller, M. Garriga, E. Wang, O. Inganäs, M. I. Alonso, Thin Solid Films 2014, 571, 371
work page 2014
-
[63]
T. Sai, M. Saba, E. R. Dufresne, U. Steiner, B. D. Wilts, Faraday Discuss. 2020, 223, 136
work page 2020
-
[64]
O. Wodo, S. Tirthapura, S. Chaudhary, B. Ganapathysubramanian, Journal of Applied Physics 2012, 112, 064316
work page 2012
-
[65]
O. Wodo, J. D. Roehling, A. J. Moulé, B. Ganapathysubramanian, Energy Environ. Sci. 2013, 6, 3060
work page 2013
-
[66]
M. Rubinstein, R. H. Colby, Polymer Physics, Oxford University Press, Oxford, New York, 2003
work page 2003
- [67]
-
[68]
M. Koopmans, V. M. L. Corre, L. J. A. Koster, Journal of Open Source Software 2022, 7, 3727
work page 2022
-
[69]
J. A. Nelson, The Physics Of Solar Cells, World Scientific Publishing Company, 2003
work page 2003
-
[70]
G. F. Burkhard, E. T. Hoke, M. D. McGehee, Advanced Materials 2010, 22, 3293
work page 2010
-
[71]
M. L. Schubert, H. Metni, J. D. Fischbach, B. Zerulla, M. Krstić, U. W. Paetzold, S. Orooji, O. J. J. Ronsin, Y . Ameslon, J. Harting, T. Kirchartz, S. Ravishankar, C. Dreessen, E. Kim, C. Sprau, M. Hussein, A. Colsmann, K. Forberich, K. Jäger, P . Friederich, C. Rockstuhl, Towards a fully differentiable digital twin for solar cells, arXiv, 2025
work page 2025
- [73]
- [74]
- [76]
-
[77]
C. Wang, C. Wöpke, T. Seiler, J. Faisst, M. List, M. Kuhn, B. Joseph, A. Ehm, D. R. T. Zahn, Y . Vaynzof, E. M. Herzig, R. C. I. Mackenzie, U. Würfel, M. Saladina, C. Deibel, Rethinking Charge Transport and Recombination in Donor-diluted Organic Solar Cells, arXiv, 2026. 1
work page 2026
-
[78]
Validation of morphology-aware electronic descriptors ....................................................................... 2 1.1. Comparison to published Monte Carlo, master equation and 2D-3D drift diffusion simulation results. 2 1.1.1. Exciton dissociation for layered morphologies of binary blends......................................... 2 1.1.2. Exc...
-
[79]
Drift-diffusion simulation of a real PM6-Y6 solar cell; determination of 𝜼𝒅𝒊, 𝝁𝒉𝒊, 𝝁𝒆𝒊, 𝒌𝒓𝒘𝒄 and 𝑵𝒕𝒆𝒉𝒘𝒄 ...................................................................................................................................................... 10
-
[80]
Detailed data to the structure-property investigations for morphologies M1 to M6 ......................... 13 3.1. Morphology and electronic properties of the as cast films ....................................................... 13 3.2. Green equation and simulated fill factor for a simplified drift-diffusion modelling setup ........ 14 3.3. Evolution of ...
-
[81]
Performance evolution under thermal loading for morphologies M1 to M6 calculated with various model parameters ....................................................................................................................................... 18 4.1. PM6:Y6 with different model parameters ..................................................................
-
[82]
References ........................................................................................................................................... 21 2
-
[83]
LG” for lateral gradient) or crystalline (“LGcr
Validation of morphology-aware electronic descriptors 1.1. Comparison to published Monte Carlo, master equation and 2D- 3D drift diffusion simulation results. 1.1.1. Exciton dissociation for layered morphologies of binary blends Yang and coworkers performed Monte -Carlo simulations to calculate the exciton dissociation efficiency for different type of pla...
-
[84]
Drift-diffusion simulation of a real PM6-Y6 solar cell; determination of 𝜼𝒅𝒊, 𝝁𝒉𝒊, 𝝁𝒆𝒊, 𝒌𝒓 𝒘𝒄 and 𝑵𝒕𝒆𝒉 𝒘𝒄 The objective here is to get reasonable reference values for 𝜂𝑑𝑖, 𝜇ℎ𝑖, 𝜇𝑒𝑖, 𝑘𝑟 𝑤𝑐 and 𝑁𝑡𝑒ℎ 𝑤𝑐, which reflect the optoelectronic properties of the considered donor-acceptor blend. For this, we choose to fit the JV curve of real PM6-Y6 solar cells in or...
-
[85]
Detailed data to the structure-property investigations for morphologies M1 to M6 3.1. Morphology and electronic properties of the as cast films The following table show morphological parameters describing the as -cast films and the ir corresponding performance results. We introduce in this section the following additional parameters 𝜙𝐶𝑟,𝑑, 𝜙𝐶𝑟,𝑎, 𝜙𝐴𝑚,𝑑, 𝜙...
-
[86]
Performance evolution under thermal loading for morphologies M1 to M6 calculated with various model parameters 4.1. PM6:Y6 with different model parameters In the main text, it is assumed that bimolecular recombination is non -Langevin. A strong penalty has also been used for crystal to amorphous phase hopping. The figure below shows the sensitivity of the...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.