Thermal eﬀect on the morphology and performance of organic photovoltaics †

The morphology of organic photovoltaics (OPVs) is a significant factor in improving performance, and establishing a method for controlling morphology is necessary. In this study, we propose a device-size simulation model, combining reptation and the dynamic Monte Carlo (DMC) algorithm, to investigate the relationship between the manufacturing process, morphology, and OPV performance. The reptation reproduces morphologies under thermal annealing, and DMC showed morphology-dependence of performance: not only short-circuit current density but also open-circuit voltage had optimal interfacial areas due to competition between exciton dissociation and charge collection. Besides, we performed transient absorption spectroscopy of various BHJ morphologies under realistic conditions, which revealed prompt and delayed dynamics of charge generation—the majority of the charges were from excitons that were generated on interfaces and dissociated within a few picoseconds, and the others from excitons that migrated to interfaces and dissociated on the order of sub-nanoseconds.


Introduction
Organic photovoltaics (OPVs) are promising devices for ubiquitous electricity generation which replace traditional silicon solar cells because of their light weight, flexibility, roll-to-roll printability, 1,2 and short energy payback time. 3 However, the power conversion efficiencies (PCEs) of the OPVs are still low (11% at the highest 4 ) and hence optimization schemes are required. The active layers of OPVs are composed of organic donor and acceptor semiconductors: p-conjugated polymers are frequently used as the donor materials, e.g., poly(3-hexylthiophene) (P3HT), while small molecules such as fullerene derivatives are used as acceptor materials, e.g., [6,6]-phenyl-C 61 -butyric acid methyl ester (PC 61 BM). P3HT:PC 61 BM is one of the most studied pairs of semiconductors for OPVs, the chemical structures of which are shown in Fig. 1.
The morphology of OPVs-the degree of phase separation, miscibility, and crystallinity of the organic materials-is considered as one of the significant factors for the PCE improvement of OPVs because it strongly affects the efficiencies of exciton dissociation and charge collection described above. In the case of OPVs with a smaller interfacial area and a larger domain, charge carriers are collected efficiently because of transport pathways without bottlenecks, while few excitons dissociate before geminate recombination. 5,6 Bulk heterojunction (BHJ) morphology is a breakthrough in OPV development because of its large interfacial area and bicontinuous charge transport pathways. 7 Morphologies are determined not only by the chemical properties of materials-interaction or solubility-but also by fabrication processes. There have been various experimental studies that investigate how to control morphology to improve the PCE: thermal or solvent annealing, [8][9][10][11] and control of the drying conditions. 12,13 For instance, it is known that the annealing temperature has an optimum value for efficiency. 8,9 Transient absorption spectroscopy was also performed to track the efficiency and dynamics of excitons and charge carriers, and influences of morphologies were examined. 14 Meanwhile, various theoretical studies on OPVs have been presented to obtain a detailed understanding of the mechanism of photocurrent generation. 15 Quantum chemical calculations estimate properties such as energy levels and distributions of molecular orbitals, absorbance, [16][17][18][19] and rate of an elementary process. 20 However, such molecular-level investigations alone are not sufficient to provide guiding principles for OPV optimization because photocurrent generation is a complicated macroscopic process. Therefore, several theoretical approaches have been presented to investigate morphologies and efficiencies, which are briefly reviewed in the remainder of this section.
To analyze morphologies, molecular dynamics (MD) calculations are performed, which reveal the molecular packing 21 and molecular orientation of polymers (edge-on and face-on configurations). 22 Although MD simulations can reveal all morphological properties for OPVs in principle, the computational cost of a device-size system is extremely high even if coarse-grained MD techniques are applied. Another approach is simulations using the Cahn-Hilliard equation, which has been established in soft matter physics. This method requires low computational cost and is applicable to device-size simulation. [23][24][25][26][27][28] However, it cannot determine the orientation of polymers because the Cahn-Hilliard equation is a differential equation of concentration of chemical components. Based on the above, we adopt a coarse-grained model called reptation, which has been established to describe polymer melts and solutions. 29,30 Frost et al. applied reptation to OPVs and revealed that the fraction of polymers would influence quantum efficiency through morphologies. 31 Reptation also involves low computational cost, but additionally can reveal morphological properties like orientation and miscibility of polymers, and thermal influences on them. 31 Therefore this method is suitable for studying OPVs, especially the influence of thermal annealing on them.
Several approaches to investigate the PCEs of OPVs have been presented. The simplest one is the equivalent-circuit model, which assumes a solar cell as a circuit of current sources, parallel diodes, and series and parallel resistors. 32,33 Device simulators known as Scharfetter-Gummel models simultaneously solve coupled equations of electric field and carrier dynamics in the OPVs-the Poisson equation, equation of continuity, and drift-diffusion equation. [34][35][36] The last ones are the dynamic Monte Carlo (DMC) algorithm with the first reaction method (FRM) and kinetic Monte Carlo (KMC), which stochastically simulate the dynamics of elementary processes according to the rate. There have been studies which examine morphologydependence of performance-current density and efficiencies of elementary processes such as internal quantum efficiency (IQE)-and dynamics of excitons and charge carriers. 27,28,[37][38][39][40][41][42][43][44][45][46] We adopt DMC because it can treat morphologies explicitly and track carrier dynamics in the morphologies.
In this work, we present a combined model of reptation and DMC, which enables the investigation of the relationship between thermal annealing, morphologies, and OPV performances. Morphologies were generated under various conditions by use of reptation, in which current density-voltage ( J-V) characteristics were simulated and performance parameters such as short-circuit current density ( J SC ), open-circuit voltage (V OC ), and PCE were investigated. We showed the existence of optimal annealing conditions for the PCE because of a trade-off between exciton dissociation and charge recombination.

Methodology
In this article, we present a model of OPVs by combining the two simulation techniques. Morphologies and OPV characteristics were simulated by reptation and DMC, respectively. Various parameters used in these techniques were set to model OPVs of indium tin oxide (ITO)/P3HT:PC 61 BM/aluminum, and the simulations were performed at a device scale with a system size of 150 3 nm 3 to perform realistic analyses. The simulation details of reptation and DMC are described below.

Reptation
Reptation, which was introduced by de Gennes 29 and refined by Doi and Edwards, 30 is a coarsening model of entangled polymer systems, such as melts or concentrated solutions. The polymer is enclosed in a space called a ''tube'' and can move easily along the tube, but not perpendicularly to it like a slithering snake.
In this work, polymers were modeled as self-avoiding chains of a certain number of segments connected by bonds, and the excluded volume effect was considered. The system was composed of 150 3 sites, each of which was occupied by either a polymer segment or a small molecule, and the lattice constant was set to 1 nm. Periodic boundary conditions were applied for the x, y, and z directions.
The reptation was simulated stochastically with the Metropolis Monte Carlo (MMC) algorithm. 31 The scheme is shown in Fig. 2. In each MMC step, an end of a polymer tries to move to one of the nearest neighboring empty lattices at random, followed by the rest of the polymer. The acceptance probability p is determined by the temperature (T) and energy change (DE) due to the trial reptation. The total energy of the system, E, is given by the sum of interactions between the nearest neighboring sites according to the materials occupying them, while DE is calculated by the interaction change of only the two sites, where the polymer enters and exits, because the other segments of the polymer are unchanged.
Let n p and n m be the coordination numbers of polymer segments and small molecules for a site, respectively, and e pp , e pm , and e mm be polymer-polymer, polymer-small molecule, and small molecule-small molecule interactions, respectively. Interaction between polymer segments connected by a bond was assumed to be constant and ignored in the calculation of the energy because the length of the polymers was fixed in the simulation. The DE due to reptation is calculated as: DE en/ex = (n en/ex p e pp + n en/ex m e pm ) À (n en/ex p e pm + n en/ex m e mm ).
Here, the upper indices ''en'' and ''ex'' represent those for the site where the polymer enters or exits, respectively. If the coordination numbers of the above site are the same, i.e., n en p + n en m = n ex p + n ex m , the DE of reptation is simply calculated as with De e pp À 2e pm + e mm ; negative or positive De correspond to mutually attractive or mutually repulsive polymers, respectively. The simulation scheme is as follows: (i) configure polymers randomly in the system ( Fig. 5(a)). (ii) Randomly and independently choose a polymer end and direction to reptate that are possible-rechoose them if the pointed site is already occupied by another polymer. (iii) Calculate the energy difference DE of reptation. (iv) Accept the reptation with probability p calculated as To investigate the dependence of morphologies on annealing temperature and other conditions, reptation was simulated for various bDe, number of polymer segments n seg , and volumetric ratios. The number of polymer segments, n seg , was set to 25, 50, 75, and 100, and the volume fraction of polymers to 33% or 50%; the number of polymers was adjusted to achieve those volume fractions in the simulation volume of 150 3 nm 3 . The dimensionless parameter bDe was set to nonpositive values to simulate mutually attractive or noninteracting polymers or infinite temperature.

Dynamic Monte Carlo
2.2.1 Basic algorithm. DMC simulations were performed to investigate morphology dependencies of photovoltaic characteristics. The system size was unchanged at 150 3 nm 3 , while the lattice constant a was increased to 3 nm, 37,39,44 i.e., the number of sites was decreased to 50 3 to examine three-dimensional devicescale morphologies. We investigated some of the morphologies generated by the reptation (BHJ), and a cube of 3 3 sites was coarsened to one donor or acceptor site depending on the fraction of polymeric donors and molecular acceptors. In addition, some ordered morphologies were investigated for comparisonstriped morphologies with widths of 75 nm and 15 nm, and checkered morphologies with periods of 75 nm Â 75 nm, 15 nm Â 75 nm, and 15 nm Â 15 nm (shown in Fig. 3). Periodic boundary conditions were applied for the x and y directions, and an anode (ITO) and a cathode (aluminum) were set on the top and bottom in the z-direction, respectively. Temperature T was fixed at 300 K.
In DMC simulations with the FRM, waiting times of all possible events, calculated from their rate, are stored in a queue; the event holding the minimum waiting time is executed in each DMC step. 37,47 This work considers three particles, i.e., excitons, holes, and electrons, which are allowed single occupancy of a site. A donor site would be occupied by an exciton or a hole, while an acceptor site by an electron. The simulation scheme is summarized as follows: (i) initialize the system: the first events under photo and dark conditions are exciton generation and charge (hole or electron) injection, respectively. (ii) Insert waiting times of all possible events associated with the particles into the queue; waiting time t of each process is calculated as where W is the rate of the process and X is a random number uniformly distributed in (0,1). (iii) Pop the event that has the shortest waiting time from the queue and execute it. (iv) Advance the simulation time and reduce waiting times in the queue by the expired time. Insert waiting times of newly possible events into the queue, and remove impossible events from the queue, for example events where single occupancy becomes unsatisfied. The next subsection explains the details of the stochastic modeling of elementary processes.

Elementary processes.
Elementary processes considered in the DMC simulations are (i) generation, (ii) hopping, (iii) dissociation, and (iv) recombination of excitons; (v) injection, (vi) hopping, (vii) collection, (viii) leak, and (ix) recombination of free carriers, i.e., holes and electrons. These processes are schematically shown in Fig. 4(a). In particular, events related to carrier transition between the active layer (P3HT:PC 61 BM) and the electrodes (ITO and aluminum) are summarized in Fig. 4(b). The parameters used in this simulation are listed in Table 1.

Excitons.
Excitons are bound electron-hole pairs generated by an absorbed photon and have no charge.
(i) In this work AM 1.5 irradiation (power of 100 W cm À2 ) was assumed and the exciton generation rate W gen ex was set to 8 Â 10 À10 ps À1 nm À2 , which yields a current density of 12.8 mA cm À2 when the internal quantum efficiency (IQE) reaches 100%. To save computational resources we insert an event of exciton generation into the queue at a time, without considering where to generate an exciton; at execution of exciton generation an unoccupied donor site is randomly chosen. The probability p is weighted by depth according to the Lambert-Beer law, where a is the absorption coefficient of the active layer, set to 3 Â 10 4 cm À1 . 49 (ii) Once generated, an exciton will hop to a nearby donor site; in this simulation, hopping to all possible sites is considered and the rate W hop ex is calculated by the Förster energy-transfer theory, 52 where t life is the exciton lifetime, DE is the energy change due to the hopping, and R 0 and R are the Förster radius and hopping distance, respectively. Because an exciton is neutral (always DE = 0), the rate is simply calculated as The lifetime, t life , and Förster radius, R 0 , were set to 330 ps, 48 and 3 nm, respectively, so that the diffusion length of excitons was in the order of 10 nm. A cutoff was introduced to reduce the computational cost, i.e., exciton hopping was considered only for sites within the cutoff distance R cut = 15 nm. (iii) When an exciton reaches the interface, it is able to dissociate into a hole and an electron. In the simulation all possible donor-acceptor nearest neighboring pairs were considered, and after exciton dissociation a hole was created on the donor site that the dissociated exciton occupied, and an electron on the acceptor site across the interface. The exciton dissociation rate, W dis ex , was set to 0.15 ps À1 , which was estimated by quantum chemical calculations. 20 (iv) Otherwise, an exciton decayed to the ground state and geminate recombination occurs. The exciton recombination rate is the reciprocal of the exciton lifetime: W rec ex = 1/t life .

Holes and electrons.
Holes and electrons have charge and hence their transport is affected by the difference of site energy, DE, before and after the event. The site energy is a summation of the following three factors: 1. On-site energies of charge carriers (holes/electrons), which are a summation of the orbital energy (HOMO of the acceptor/LUMO of the donor) and Gaussian disorder of density of states (DOS) caused by residual solvents and structural disorder. The HOMO of the acceptor, LUMO of the donor,  and Gaussian width s were set to À5.0 eV, À3.7 eV, and 0.062 eV, 39 respectively. 2. The Coulomb interaction, U, with other charge carriers and image charges, including that of itself. The Coulomb interaction of carrier i is calculated as where e is permittivity, q i is the charge of carrier i, and r ij is the distance between carriers i and j; the minimum image convention was applied, namely where n x,y is the integer in {À1, 0, + 1}, and L x,y is the length of the simulation box along x and y (L x = L y = 150 nm). The relative permittivity e r was set to 3.5 as a typical value for OPVs. 3. The electrostatic built-in potential, V built , which is induced by the applied voltage V and the difference in the work functions of the electrodes (DW = 0.5 V): with elementary charge e.
(v) A charge carrier was created in the active layer by either exciton dissociation or injection from the electrode; a hole (an electron) is injected from the anode (cathode) to one of the donor (acceptor) sites at the surface. These injected carriers contribute to the dark current, which is enhanced by the applied voltage, V. The charge injection rate W inj h/e is calculated by the Miller-Abraham equation 53 as where m h/e is the hole/electron mobility of pristine donor/ acceptor materials, V T is the thermal voltage (V T k B T/e), and a is the lattice constant. Hole and electron mobilities, m h and m e , were set to 0.1 cm 2 V À1 s À1 50 and 0.025 cm 2 V À1 s À1 , 51 respectively. The prefactor is derived from the Einstein relationship. 37,54 In calculations of the energy difference DE, the Gaussian disorder of an electrode was not considered. In addition to the above three factors, DE included the charge injection barrier of +0.4 eV, which is derived from the difference between the LUMO of the acceptor (À3.7 eV) and the Fermi level of the cathode (À4.1 eV). Although an injected carrier could be collected from the electrode, it was forced to hop at once before extraction if possible.
(vi) Holes and electrons inside the active layer are allowed to hop to one of the nearest neighboring donor and acceptor sites, respectively. The charge hopping rate W hop h/e is calculated using the Marcus formula as where the reorganization energy l, equivalent to twice the polaron-binding energy, was set to 0.11 eV. 20 (vii) When a charge carrier reaches the electrode surface, it is able to be collected and then the current flows. The charge collection rate is also calculated by the Marcus formula (eqn (14)), but the negative of the charge injection barrier is considered for DE. Collection of holes from the cathode and electrons from the anode contribute to the photocurrent.
(viii) Leakage of charge carriers, i.e., extraction of holes/ electrons from the anode/cathode, is simulated in the same way as above, but it contributes to the dark current.
(ix) Charge recombinations are considered when a hole and an electron are located on adjacent sites. The rate W rec ch was set to 1.93 ns À1 . 20 2.2.3 OPV characterization. DMC simulations were performed for both photocurrent and dark current. No exciton is generated in the latter case while excitons are generated in the former case according to the procedure explained previously. Performance parameters such as J SC , V OC , IQE, and PCE were calculated for the morphologies. Charge injection (holes/electrons from the anode/cathode) and leakage (holes/electrons extracted to the cathode/anode) contribute to the dark current. Therefore, the net number of collected holes at the anode n net coll The efficiencies of exciton dissociation Z diss ex and charge collection Z coll ch were calculated as respectively, where n gen ex and n dis ex are the number of generated and dissociated excitons, respectively. The IQE is the product of Z diss ex and Z coll ch . The net current density J is calculated as the number of net collected charges divided by the incident area (150 2 nm 2 ) and simulation time. The J-V curves were obtained by repeating the calculation of current density under each applied voltage. Open-circuit voltage V OC , PCE, and fill factor (FF) were obtained by linear interpolation of the J-V curve.
Notably, this study also analyzed transient absorption spectroscopy. In the simulation, a certain number of excitons were created at the start under the conditions that no charge injection, leak, nor collection and no applied voltage were imposed. The initial number of excitons was set to 100 and 1000: the latter corresponds to a density of 3 Â 10 17 cm À3 , which is consistent with the experimental value; 14 the former was for comparison, and note that we can simulate dilute conditions that experimental observation is difficult. There have been studies that discuss the dynamics of carriers-the effect of mobility on charge separation, 45 and the kinetic model of charge separation and recombination. 46 In this work we examined the dependence of transient absorption spectroscopy on BHJ morphology in consideration of annealing temperature and analyzed spacetime dynamics.
Five simulation runs were performed independently for each condition-0.5 ms and 10 ns for a simulation of current density and transient absorption spectroscopy, respectively-and the results were averaged.

Reptation
The generated morphologies with a 1 : 1 volume ratio and 75 polymer segments are shown in Fig. 5, whose phase separation is consistent with an experimental report 55 and theoretical studies using the Cahn-Hilliard equation. [23][24][25][26][27][28] See Fig. S1 and S2 in the ESI † for the other morphologies-those of different volume ratios and number of segments. Generally, the domain size grew as bDe increased. The temperature increased under the negative constant De, while the polymers dispersed when the temperature was too high. These morphologies are characterized numerically: interfacial area per volume and mean distribution of domain size, shown in Fig. 6(a) and (b). Here, the domain size is defined as the length of continuous donor/ acceptor domains along x, y, or z. Both characteristics depended strongly on bDe but only slightly on the polymer length. The interfacial area and domain size indicate the minimum and maximum values, respectively, around bDe = À1 because of the degree of diffusion. At low temperature (bDe = À4.0, shown in Fig. 5(i)), the polymers are hindered in diffusing and trapped in a local minimum because of the attractive interaction. As the temperature rises, the polymer can diffuse and relax, therefore donor domains are formed; however, free diffusion breaks the domains when the temperature is too high. The polymer length also affects diffusion and therefore shows a similar dependence. See Section S1.2 in the ESI † for the results of a polymer/small molecule volumetric ratio of 1 : 2. The interfacial area and domain size were smaller and their dependence on bDe became milder than those of a 1 : 1 ratio, due to the lower polymer fraction and weaker entanglement effects.
The domain-size distributions of the polymers and small molecules are shown in Fig. 6(c) and (d), respectively. Essentially, the distributions of the domain size for both polymers and small molecules exhibit exponential decay and there is no significant difference between polymers and small molecules. However, the morphology of bDe = À0.5 has significantly different characters: the domains of small molecules in that morphology are larger than those of polymers, and exhibit exceptionally oscillating properties. This was because polymer aggregation precluded acceptors and formed large acceptor domains but polymer domains were broken due to diffusion.
The results in this subsection propose that the morphologies of OPVs strongly depend on the temperature by using the physical modeling of blends of polymers and small molecules.

Dynamic Monte Carlo
For BHJ morphologies, we chose those with a 1 : 1 volumetric ratio and 75-segment polymers, n seg = 75, to investigate wide ranges of interfacial area and domain size (see Fig. 5 and 6). It is expected that other morphology groups-different n seg and/or volumetric fraction-show a similar but milder temperature dependence.  For example, morphologies with a 1 : 2 volume ratio have smaller donor domains, which probably enhances the exciton dissociation efficiency, but lower donor fraction yields fewer excitons and thus their effect on current density is likely canceled.
3.2.1 Performance parameters. Current density-voltage ( J-V) characteristics are shown in Fig. 7. Since the sample standard deviations were small (ranging from 0.03 mA cm À2 to 0.3 mA cm À2 , see Section S2.1 in the ESI †), the error bars were omitted for clarity. Both photo-and dark currents depended on the morphologies. Interfacial areas of the coarsened morphologies used in these DMC simulations are listed in Table 2, and interfacial area dependencies of performance parameters-short-circuit current density J SC , open-circuit voltage V OC , PCE, and FF-are shown in Fig. 8. The J SC , V OC , and PCE exhibit maximum values related to the interfacial area, while the FF decreases monotonically as the interfacial area increases. Different BHJ morphologies show different dark currents (Fig. 7(c)), while the ordered ones show similar tendencies (Fig. 7(d)). This is because all of ordered morphologies have suitable charge transport paths and charge recombination merely occurs. Photocurrent, shown in Fig. 7(a) and (b), is determined by exciton dissociation and charge collection, whose details are described in the next paragraph. This study shows the morphology dependence of V OC , which arises from the morphology dependence of the photo-and dark currents.
Efficiencies of exciton dissociation and charge collection, and the IQE under short-circuit conditions are shown in Fig. 8(e), (f), and (g), respectively. Their sample standard deviations were smaller than 0.02. In the morphologies with large interfacial areas, the exciton dissociation is enhanced, but the charge collection is discouraged by a longer transport path. Therefore, the IQE Fig. 7 Current densities for different morphologies and conditions: photocurrent densities of (a) the BHJ morphologies generated by the reptation and (b) the ordered morphologies; dark current densities of (c) the BHJ morphologies generated by the reptation and (d) the ordered morphologies.  exhibits a maximum value because of the trade-off. These tendencies are consistent with previous studies. 27,37,39,40,44 In addition, we theoretically demonstrated that the opencircuit voltage V OC and the logarithm of the IQE have a linear correlation (Fig. 8(h)), which has been reported by an experimental work. 18 The extrapolated V OC to 100% IQE is a little less than 0.5 V, which equals the difference of electrode work functions. Therefore, we revealed that improving the IQE enhances not only J SC but also V OC .
The dependencies of performance parameters of BHJ morphologies on bDe, corresponding to finite temperature, are shown in Fig. 9. The morphologies generated under higher bDe tend to give higher J SC and FF with lower V OC . It is shown that bDe had an optimal value for the PCE-the morphology with bDe of À2 showed the highest PCE of 1.4%. Remarkably, we revealed that morphology only, without parameter change, determines open-circuit voltage. These effects of thermal annealing are consistent with some experimental results. 56,57 3.2.2 Transient absorption spectroscopy. The time dependence of the number and peak normalized number of generated charge carriers with 1000 initial excitons is shown in Fig. 10(a) and (b) corresponding to the dynamics of charge generation and recombination, respectively. The simulated number of excitons is shown in Fig. S5 (ESI †). The standard error of the number of excitons and charge carriers was 10 at the maximum. The tendency was independent of the number of created excitons at the start (see Fig. S6 for the results of 100 initial excitons, ESI †). The ratio of dissociated to generated excitons varied in a range from 40% to 90%, consistent with the exciton dissociation efficiencies shown in Fig. 8(e). The time when the number of charge carriers takes the maximum also depended on morphologies: those with lower bDe tended to show a sudden increase. The morphologies generated with higher temperature tended to have larger interfacial areas, and created more charge carriers, whereas they had more bimolecular recombinations, which is due to their large interfacial area. These tendencies are the same as the experimental observations. 14 The space-time dynamics of charge generation-diffusion length and lifetime of the dissociated excitons-is shown in Fig. 11; see Fig. S7 (ESI †) for the results of other BHJ morphologies. Note that the lifetime of excitons equals the simulation time. The simulation revealed prompt and delayed dynamics of the charge generation: the majority of the charges (40-95%, varying for morphologies) were from excitons that were generated on interfaces and promptly (within a few picoseconds) dissociated, corresponding to the sudden increase as shown in Fig. 10; the others were from excitons that migrated to interfaces and dissociated later (up to several hundred nanoseconds), corresponding to the slow increase. The tendency depends on the morphologies because in the case of morphologies of lower temperature, those have higher   interfacial area and smaller domains, excitons are more likely generated on a site on an interface and otherwise could reach interfaces easily, resulting in prompt increase of carriers; morphologies of higher temperature show fewer prompt and more delayed exciton dissociation, therefore the number of charges takes the maximum value at a few 100 ps.

Concluding remarks
We revealed a relationship between the annealing temperature, which is one of the controllable parameters in the manufacturing process, and the performance of OPVs via morphologies of donor and acceptor domains in the active layers. To this end, we proposed a model that combines reptation and DMC. The reptation and DMC were performed to investigate the physical mechanism of formation of the morphologies and the electrochemical profile of photo-and dark currents, and transient absorption spectroscopy, respectively. Then, the reptation with a 1 nm-lattice showed temperature dependencies of morphology. The interfacial area exhibited a minimum value related to the annealing temperature because of the trade-off between thermal diffusion and attractive interaction. The DMC simulations, in which the generated BHJ morphologies were coarsened to a 3 nm-lattice, revealed morphology dependencies of the photoand dark currents, and OPV performance: short-circuit current density J SC , open-circuit voltage V OC , and PCE exhibit maximum values related to the interfacial area because of the trade-off between exciton dissociation and charge collection. This work theoretically demonstrated the linear correlation between opencircuit voltage V OC and the logarithm of IQE, which agrees with experiment. Besides, we performed transient absorption spectroscopy of various BHJ morphologies, which revealed two different dynamics of charge generation-prompt (a few picoseconds) generation from excitons that were generated on interfaces and delayed (order of sub-nanoseconds) generation from excitons which migrated to interfaces; the former is dominant but its ratio depends on the domain size of morphologies.
Coarsening of lattices in the DMC simulations implicitly assumes that electron-hole separation immediately reaches 3 nm, which likely leads to overestimation of the charge separation efficiency and thus shall be verified in future. It is, actually, numerically observed that charge recombination is suppressed in the case of the coarsened morphologies. A mechanism, for example the charge hopping rate increases locally on an interface due to factors like a hot process, 58 is necessary to compensate it.
The assumption and mechanism require further study, which is beyond the scope of this work.