Monte Carlo Simulation for Investigation of Morphology Dependent

Oct 9, 2018 - Monte Carlo continues time random walk simulations. 1. INTRODUCTION ... studied by continuous time random walk (CTRW) by using the Monte...
0 downloads 0 Views 6MB Size
Article Cite This: J. Phys. Chem. C XXXX, XXX, XXX−XXX

pubs.acs.org/JPCC

Monte Carlo Simulation for Investigation of Morphology Dependent Charge Transport in Bulk-Heterojunction Organic Solar Cells Mahyar Taherpour and Yaser Abdi*

J. Phys. Chem. C Downloaded from pubs.acs.org by REGIS UNIV on 10/19/18. For personal use only.

Nanophysics Research Laboratory, Department of Physics, University of Tehran, North Kargar Avenue, Tehran 1439955961, Iran

ABSTRACT: Morphology, the spatial distribution of traps, interdomain connectivity, and phase separation of the active layer play a critical role in the performance of the bulk heterojunction (BHJ) organic solar cells (OSCs). In this work, we utilize the hopping transport model to simulate the effect of morphological and structural parameters on the diffusion coefficient and efficiency of the polymer-fullerene BHJ solar cells. In BHJ solar cells there are two distinct phases as electron transport material (acceptor) and hole transport material (donor). Here we try to create an almost realistic network containing P3HT polymer chains and PCBM clusters for simulating the charge transport in the active layer. The blend ratio of P3HT:PCBM polymers and alignment of these bicontinuous networks of active layer are considered here as the morphological parameter affecting the cell performance. The dependency of the charge transport on such morphological parameters is obtained in this study by using Monte Carlo continues time random walk simulations.

1. INTRODUCTION Nanostructured solar cells have attracted remarkable interest as an alternative to silicon-based solar cells. Many types of nanostructured solar cells have been developed in recent years based on semiconducting nanoparticles,1 perovskites,2 and polymers.3,4 In all types of the last generation of solar cells, optimization of the morphology is critical for cell performance. Experimental investigations on the efficiency of the fabricated solar cells have shown that the cell performance strongly depends on the morphology of the transport media.5 Electrical performance of the solar cells comes from three successive steps: generation of electron−hole pairs, transferring the electron to electron transporter and hole to hole transporter, and finally charge transport in transport media. After separation of the electron−hole pairs, they need to be conducted to their respective electrodes by a successive transport. Depending on the mobility of the electrons and holes through their transport materials, the charge transport may be limited by electron transport or hole transport. Charge transport in nanostructured media seems to be occurred through the localized states by the hopping mechanism.6,7 Also, the recombination of the charge carriers is believed to be aided by impurities that act as deep traps.8 Morphology of the transport media strongly affects the charge transport, leading to affecting the cell performance.9,10 © XXXX American Chemical Society

Charge transport in the nanostructured solar cells was recently studied by continuous time random walk (CTRW) by using the Monte Carlo (MC) simulations.11 In the RW simulation of charge transport, a random movement of charge carriers in a 3dimensional network of traps is considered.11−14 Arrangement of the traps in the simplest version is carried out on an ordered lattice with the same average distance, although disordered trap configurations can be implemented. CTRW simulation of charge transport is based on diffusion of photocharge carriers through a trap contained media. The precise nature of the localized traps and their spatial location in the transport media are in doubt. Besides, Street et al.15,16 have suggested that the recombination rate through hole-traps in the polymer play an important role in the performance of BHJs. There are two main ideas about the energy distribution of the trap states. Exponentially decreasing tail of trap states below the conduction band is usually given for inorganic media.17 For example, Anta et al. have been utilized a truncated exponential distribution for trap states in dye-sensitized solar cells.18 On the other hand, Gaussian distribution is believed to be observed for localized states in disordered organic materiReceived: May 26, 2018 Revised: August 16, 2018 Published: October 9, 2018 A

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 1. (a) P3HT crystal schematic containing lamellar stacking direction with face-on and edge-on orientations. The chains in the crystalline lamellae are drawn in black. (b) Schematic showing the shift factor and overlap parameter defined in our simulation.

als.19−22 Charge transport in disordered trap contained media as photoanode in solar cells was investigated versus the morphological parameters.23 Ansari et al. have shown by CTRW simulation in a porous TiO2 that the geometry dependency of the diffusion coefficient can be described only as a function of the porosity.24 Also, Abdi et al. have introduced equations describing the morphological dependency of electron transport in nanostructured solar cells.25 The mentioned previous reports on the morphology-dependent electron transport are applicable just for semiconductor media in which the limiting parameter for charge transport is only porosity. On the other hand, there are some reports on the simulation of charge transport in organic solar cells. For example, Koster26 and Woellner et al.27,28 have reported the investigation of the influence of morphology on the charge transport of BHJ solar cells by numerically solving the Pauli master equation. Heiber et al.29 have investigated the effect of the tortuosity of the charge-transport pathways through a bulk heterojunction film on the charge transport using kinetic Monte Carlo simulations. They have used simplified models such as the Ising-based model to generate a porous media for charge transport. Here, we have tried to create a more realistic media to simulate the BHJ active layers. According to the experimental reports, in the organic BHJ solar cells, the limiting cases are a perfect mixture of donor and acceptor.30,31 However, the perfect morphology for organic bulk heterojunction (BHJ) solar cells is unknown because of

experimental limitations. In this work, we investigate the effects of morphology on the charge transport in BHJ solar cells. We show that not only the blend ratio but also the phase separation and the interdomain connectivity can affect the charge transport in this kind of solar cells. There are two main approaches for utilizing in CTRW simulation to explain charge transport. In the first one, which is called multiple trapping (MT) model, transport is affected by trapping and detrapping events.32 In this model, the charge transport occurs by electrons in the conduction band or holes by the valence band. The second approach is hopping transport model in which, the transport can be occurred by the direct transition between the localized trap states.33−36 The hopping events are assumed to be occurred as a consequence of thermal activation and quantum mechanical tunneling.37 Here we utilized hopping transport model in CTRW Monte Carlo simulation, with a Gaussian distribution for trap energies, to describe the morphological dependency of the charge transport in BHJ solar cells.

2. FILM MORPHOLOGY SIMULATION There are some experimental reports on the study of the morphology and structure of organic semiconductors and organic photovoltaics.38−42 Grazing incidence X-ray scattering (GIXS) and scanning tunneling microscope (STM) analyses of P3HT show rectangular crystallites. According to these reports, the lamellar stacking of P3HT has a spacing of ∼1.6 nm and B

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 2. (a) Realistic network including the P3HT clusters and PCBM spherical clusters obtained in the simulation with the parameters of SF = 0.7, Ov = 0.1, and average tortuosity of 1.41. (b) PCBM cluster obtained by random-cluster model and the schematic showing the P3HT cluster with the aromatic (π−π) stacking spacing of ∼0.39 nm.

the aromatic (π−π) stacking has a spacing of ∼0.39 nm. Verploegen et al.43 have investigated the effects of the thermal annealing on the orientation of P3HT rectangular crystallites relative to the substrate. Sirringhaus et al.44 showed that, in the P3HT films with high regioregularity and low molecular weight, the preferential orientation of ordered domains is with the (100)-axis normal to the substrate and that, in the films with low regioregularity and high molecular weight the preferential orientation of ordered domains is with the (010)-axis normal to the substrate. Also, scanning tunneling microscope images of P3HT clearly show an ordered polymer domain with an average size of 20 nm.45−47 Besides this, it has been shown that after annealing, this polymer chain is extended in the lateral side to create a P3HT crystal along (100) and (010).45,47−50 In order to simulate P3HT as donor species, we have used a cubic box containing flat sheets as P3HT face-on and edge-on crystallites51 with spacing of 0.4 nm between the sheets, as shown in Figure 1a. In order to generate the network for the simulation of the charge transport, the prepared clusters are randomly distributed at the surface of simulation box with the size of 150 × 150 × 150 nm at z = 0 with the overlap parameter less than a desired value. Then, the same clusters are placed randomly on top of the initial clusters in the manner that the overlap parameter (OV) and shift factor (SF) for all of the neighboring clusters be less than the specific values. Shift factor and overlap parameter are defined as SF = shift factor =

Sy Sx = Lx Ly

OV = overlap parameter = 1 −

of the clusters in the vertical direction can be achieved by a minimum value of shift factor. This process is repeated so that the structure grows up in the z-direction. In the literature, such alignment is quantified as a well-known parameter called tortuosity.29,52 Tortuosity quantitatively indicates how convoluted a transport pathway is relative to the shortest straight path. In order to obtain the average tortuosity for the networks with different shift factors, we calculated the shortest pathways through each phase from each corresponding point on the bottom of the simulation box to the top. As calculated pathway was then divided by the box height (150 nm) to achieve the tortuosity for that specific pathway. This calculation was repeated to get an average tortuosity for a specific network. In order to create PCBM clusters as acceptor species, we used a random-cluster model.53 The procedure is started by putting a sphere (as PCBM) with radius r = 1.5 nm at the void spaces between the P3HT network. Then, spheres with the same radius are being brought randomly around the initial sphere on the condition that they have to be connected to at least one other sphere. This process is repeated enough to fill the free spaces between the P3HT layered clusters. It must be mentioned that conflicted clusters with P3HT clusters were removed in this process. A realistic media produced by this algorithm is shown in (Figure 2). The values of the constants which have been chosen from the previous experimental and simulation reports are listed in Table 1. The proper references are included in the table.

(1)

Oy Ox =1− Lx Ly

Table 1. Constants Used in the Morphology Simulation (2)

parameters

where Sx (Sy) is the distance between the edges of the upper cube and initial one along the x (y) direction and Ox (Oy) is the distance between the edges of the cubes located at the same height along the x (y) direction (see Figure 1b). Decreasing the shift factor in the produced network leads to have an ordered continuous region and a maximum alignment

simulation box size cluster size (P3HT) shift factor (P3HT) overlap parameter (P3HT) sphere size (PCBM) C

value

references

Lx, Ly, Lz SF Ov

150 × 150 × 150 nm 20,12,10 nm 0.1 to 0.9 0.1 to 0.25

ref 24 ref 50 − −

r

1.5 nm

ref 54

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 3. (a) Realistic 3-dimensional network of traps produced by the r3 model with Ntv = 0.3 nm−3 for the network with average tortuosity of 1.15. The picture has been visualized by the VMD software. (b) Schematic showing the localized states distributed in the volume of P3HT clusters (r3 model) and illustration of the random walk simulation employed in this work. Some of these traps are occupied by charge carriers (black circles). The possible jump to the neighboring states is selected between the unoccupied traps (open circles) located within the cutoff radius, rcut, with minimum hopping time.

Figure 4. (a) Band structure of the P3HT:PCBM heterojunction in steady-state illumination under open-circuit conditions (V = Voc). (b) Schematic showing the Gaussian energy distribution of the localized states with disorder parameter of σ around the transport level,20 E0. Arrows in the schematic show the realization of the hopping through the localized states.

3. MONTE CARLO RANDOM WALK SIMULATION

g (E , E 0 ) =

3.A. Transport Model Implementation. To simulate the charge transport through the trap states in the P3HT polymer, static traps are distributed throughout the clusters of P3HT (as shown in Figure 3). These static traps were placed on the P3HT sheets based on the well-known r3 models23 (volume diffusion), in which the traps are distributed through the material with a volume density of Ntv. We assume that the energy distribution of the localized states has an exponential form such as36,55

É ÄÅ ÅÅ (E − E )2 ÑÑÑ Ntv 0 Ñ Å ÑÑ Å expÅÅ− 2 ÑÑ ÅÅÅ 2πσn σ 2 ÑÑÖ n Ç

(3)

where Ntv is a total density per unit volume, σn is the half width of the distribution, and E0 is the transport level (Figure 4,parts a and b). The Gaussian peak of occupied states is reported to be shifted down by σ2/kT with respect to the Gaussian peak of localized state.56 Charge transport through the generated traps is simulated here by single-carrier random walk simulation. Anta et al. have shown that the single-carrier simulation can reproduce the results of multicarrier simulation.13 The random walk simulation based on hopping transport is started by jumping a hole (or an electron) to a neighboring D

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C trap with the minimum hopping time.57 Collection efficiency in the BHJ solar cells is affected by both electron transport in acceptor species and hole transport in donor species. However, in the case that the mobilities of electron and hole in their media are different, electron or hole transport can limit the transport parameters. For example, the hole mobility in annealed P3HT:PCBM (50 wt %) is lower than the mobility of electrons.58−60 Consequently, investigation of the hole transport in such BHJ solar cells is more important than the electron transport. According to this discussion, we consider the hole transport in P3HT in this work. In order to calculate the release time (tij) of each hopping step from trap i to the neighboring trap j in a cutoff radius of rcut, we have used the Miller−Abrahams equation:36,61−65 ÄÅ É ÅÅ 2rij Ej − Ei + |Ej − Ei| ÑÑÑÑ 1 Å + tij = −ln(P) expÅÅÅ ÑÑÑ ÅÅÇ αl ÑÑÖ 2kBT υ0 (4)

To simulate the recombination process, i.e., charge transfer from the donor polymer to acceptor species, we have used the Marcus theory of charge transfer. The recombination rate Krt , using the Marcus theory can be written as ij (ΔG 0 + λ)2 yz zz K tr ∝ exp( −β ·r ) expjjj− zz j 4 λ k T B k {

(5)

The first distance depended exponential term represents the electronic coupling between the wave functions of donor and acceptor species. The exponent factor, β, is a measure of the electrical coupling between the donor and acceptor. The second exponential term represents the thermal activation over an activation barrier. ΔG0 in this term is free energy change resulting from the electron-transfer reaction, and the positive parameter of λ is called the reorganization energy, which depends on polarization and relaxation in the surrounding environment. Calculating the probability of charge-transfer at the interface, recombination rate, and carrier lifetime is difficult problem in BHJ organic solar cells.73−75 Recently, a simple and general approach was developed based on the assumption of the Marcus model to formulate the recombination rate in such solar cells.67,76−83 Using this approach, we have calculated the recombination rate Krt at each step of the simulation by the following equation:

Here Ei (Ej) is the energy of the starting (target) trap, P is a random number distributed uniformly between 0 and 1, υ0 is the attempt to jump frequency, rij is the distance between the starting and target traps, and αl is the localization radius. Simulation of the hole transport is started by the distribution of the trap states on the produced network. At the beginning of the simulation, a hole randomly selects an initial trap state of i and the hopping time for the possible jump to an empty site j among neighbor traps in the cutoff radius is calculated by the above equation. After that, the hole is jumped to a neighboring trap with the minimum hopping time (Figure 3b). Such a procedure is repeated in each step of the simulation, and the calculated hopping times are added to the total time of the simulation. 3.B. Recombination Model Implementation. Mechanism of the recombination in the disordered nanostructured materials is not well-defined yet. But, the most common mechanism for the recombination of the charge carriers in such media is the recombination via intraband localized states. This trap-assisted recombination, which is referred to the nongeminate recombination of the electron−hole pairs in the presence of a density of localized states is assumed to be the dominant mechanism for recombination in the bulk heterojunction organic solar cells.9,50,66−70 Nongeminate recombination defines the loss process of two charges which originate from different generation processes. There are some reports describing the mechanism of recombination in P3HT: PCBM solar cells. For example, MacKenzie et al.71 show that the dominant mechanism for recombination in such solar cells is the nongeminate recombination of the free carriers with trapped carriers because the number of free carriers is smaller than the trapped carriers. Also Nalwa et al.72 show that the trap-assisted recombination dominates at low light intensities and that, independent of the light intensity, the trap-assisted recombination is the dominant mechanism for the fast grown active layers in which the trap density is very high. In this model of recombination, the rate of recombination is proportional to the depth of trap states. Charge carriers in the shallow traps with small ionization energies can be released with high probabilities into delocalized states. On the other hand, charges located in deep traps, those having much higher ionization energies, are not released easily. According to this assumption, the traps can act as recombination centers in the disordered materials.

2 ji (E − En + λ) zyz zz K tr = K 0i expjjj− i j z 4λkBT k {

(6)

where Ki0 is the time constant of tunneling and assumed to be a constant with units of s−1. Ei is the energy of ith trap state in the HOMO distribution of donor polymer, and En is the most probable energy for the electrons state in LUMO distribution of acceptor polymer. For realization of the recombination in our simulation, we have calculated the recombination probability as Γ = tijK tr

(7)

We have assumed that this probability is proportional to the recombination rate of the holes and the average time that hole spends in localized states. This probability is calculated for each hopping event and as obtained probability is compared with a produced random number. When the random number is equal to or smaller than the calculated probability, the recombination occurs. If it occurs, the total time of the simulation and the distance traveled by the holes are recorded and saved. Finally, carrier lifetime and diffusion length can be computed by statistical averaging over large enough simulations. For a normal diffusion, after a long enough time t, the diffusion coefficient can be calculated by the following expression:55,84,85 Dj =

rt − 1 r0 )2 ⟩ ⟨(1 6t

(8)

Dj is the jump-diffusion coefficient and the ⟨⟩ symbol denotes a statistical average. In order to overcome the finite size effects, the periodic boundary condition is applied to the nonrecombined holes when reaching to the boundary of the simulation box. We have averaged the results over up to 1500 realizations to reduce the statistical error. In order to calculate the collection efficiency, after 1000 times realizations, the number of carriers reaching to the conducting electrode is E

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C counted. The percentage of the collected charge is assumed to be collection efficiency (ηcoll) of the solar cell. The values of the constants have been chosen from some theoretical and experimental reports. They are listed in Table 2. The proper references are included in the table.

VD =

value ‑3

density of traps

Ntv

0.3 nm

attempt to jump frequency

υ0

2 THZ

system temperature traps characteristic temperature most probable energy in acceptor LUMO distribution quasi Fermi level center of the DOS donor HOMO energy constant rate of recombination from trap states localization radius reorganization energy

T T0 En

300 K 600 K −4.1 eV

Ef 0 E0 Ep K0

−4.75 eV −5 eV −5.4 eV 3.6 × 102 s‑1

αl λ

0.5 nm 0.7 eV

cutoff radius half-width of the distribution

Rcut σn

2.5 nm 0.08 eV

(9)

where Ntot v is the total number of traps that are distributed homogeneously in the volume of donor species. In addition, the shift factor has a direct impact on the tortuosity. In order to achieve the different tortuosities, the networks with different shift factors were created and the average tortuosity of the networks were calculated using the previously mentioned algorithm. As obtained tortuosities for the different shift factors are shown in Figure 5b. Also, the domain size which is defined as the total volume of the active layer over the interfacial area is a good indication for phase separation. Parts a and b of Figure 6 show the domain size and the normalized interface area achieved for different average tortuosity, and Figure 6c shows the normalized interface area at constant blend ratio 50%. 4.B. Effect of Morphology on Charge Transport. At the first step of the simulation, we have calculated the diffusion coefficient as a function of the blend ratio. We expect that the diffusion coefficient decreases by decreasing the blend ratio. Figure 7a shows the diffusion coefficient versus the blend ratio. Results of this figure are in agreement with the results of the previous works24,100−103 reporting that the diffusion coefficient is proportional to|P − PC|δ, where PC is the critical porosity obtained to be around 78% and the conductivity exponent of δ is around 1. Ansari et al. have defined the porosity as the volume ratio of the voids in the transport media over the total volume, where (1 − P) is a measure of the volume percentage of donor polymer. According to this definition, the blend ratio can be a good substitution for (1 − P) in the BHJ solar cells. Also, they show that the diffusion coefficient of the photoanodes in nanostructured solar cells only depends on the porosity and the effects of other morphological parameters can be interpreted in the context of porosity.24 Experimental studies of the impact of blending composition on device performance suggested that an optimum P3HT blending ratio was around 50−60% for the P3HT:PCBM system.104−107 So, in order to investigate the effects of other morphological parameters on the diffusion coefficient, we have prepared some simulation network with the same blend ratio of 50% and the different

Table 2. Constants Used in Simulation parameters

Nvtot Ntv

references refs 6 and 86−90 refs 91 and 92 room temp refs93 and 94 ref 74 − − ref 74 refs 76 and 95 ref 96 refs 74 and 76 ref 97 refs 21, 98, and 99

4. RESULTS AND DISCUSSION 4.A. Film Characterization. As mentioned previously, the blend ratio (volume percentage of P3HT) and the interfacial area (total surface area of P3HT polymer) are important morphological parameters affecting the cell performance. In order to achieve a wide range of the blend ratio, we have changed the overlap parameter and shift factor. As shown in Figure 5a, different blend ratios ranging from 25 to 60% can be obtained by changing the shift factor from 0.1 to 0.9. Also, in a constant shift factor, an increase in the overlap parameter leads to an increase in the P3HT blend ratio. Practically blend ratio can be obtained by computing the volume of donor over the total volume of the active layer. Also, the volume of donor species, VD, can be obtained by

Figure 5. (a) Blend ratio versus the shift factor for different overlap parameters. (b) Shift factor versus the average tortuosity obtained by morphological simulations. F

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 6. Domain size (a) and normalized interface area (b) with overlap parameter of OV = 0.1 versus the average tortuosity. (c) Normalized interface area at constant blend ratio (50%) versus the average tortuosity (different morphologies with constant blend ratio obtained by changing the overlap parameter).

Figure 7. (a) Diffusion coefficient versus the blend ratio. The red solid line shows the fitted curves corresponding to |B − Bc|δ. The critical blend ratio, Bc, is ∼22% which is in good agreement with the critical porosity of ∼78% obtained by Ansari et al. (b) Diffusion coefficient and its components versus the average tortuosity at the constant blend ratio of 50% (SF parameter was changed from 0.1 to 0.6).

component of the diffusion coefficient, Dz, decreases by increasing the average tortuosity, but the X and Y components of the diffusion coefficient remain almost constant. Also, the total diffusion coefficient decreases by decreasing the average

average tortuosity. Results of diffusion coefficient versus the average tortuosity are presented in Figure 7b. X, Y, and Z components of the diffusion coefficient are also shown in this figure separately. As we can see, at a constant blend ratio, the Z G

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 8. (a) Diffusion coefficient versus the interfacial area and (b) domain size obtained for the constant blend ratio of 50%. Domain size is defined as the volume divided by the interface area.

Figure 9. (a) Visualization of the actual path of a successive collection of a charge carrier. Each point is the indication of 100 jumps in the diffusion process. (b) Normalized collection efficiency as a function of average tortuosity at the constant blend ratio of 50% for the P3HT network with faceon (circles) and edge-on (squares) orientations.

tortuosity, as a direct consequence of decreasing the Dz. At the much smaller tortuosities, Dz and D begin to grow dramatically. Such increment in the diffusion coefficient is a direct consequence of the directional transport which can be occurred in columnar structures. Javadi and Abdi108 have shown in their previous work that the electron diffusion coefficient in a disordered material which contained localized states in a columnar structure is significantly higher than the diffusion coefficient in the conventional isotropic structure. They believe that the high diffusion coefficient in columnar structures is a consequence of confinement in x and y directions, which leads to vanish the diffusion coefficient in the confined directions and to have a directional transport in the unconfined direction. The diffusion coefficient versus the interfacial area and domain size are provided in Figure 8 to demonstrate the effects of the phase separation on the hole transport. The simulations are carried out at a constant blend ratio. According to the definition of the interfacial area and the domain size, a higher phase separation occurs at higher interfaces and lower domain

sizes. Figure 8b shows that improving the phase separation due to decreasing the domain size leads to a higher diffusion coefficient. 4.C. Effect of Morphology on Collection Efficiency. As an important electrical parameter showing the performance of a solar cell, the experimental works mostly report the efficiency of the fabricated cells. We calculated the collection efficiency of the BHJ solar cells using the calculation of the percentage of carriers reaching to the conducting electrode. The Collection efficiency of the solar cell strongly depends on the thickness of its active layer. Our simulations for calculating the collection efficiency are carried out for the constant thickness of 150 nm. A charge carrier generated in this 150 nm thickness can be collected if it reaches to the bottom electrode at z = 0. After 1000 times simulations, the percentage of carriers reaching to the bottom electrode is calculated to obtain collection efficiency. A successive collection of a carrier is visualized in Figure 9a. Figure 9b shows the collection efficiency of the cells as a function of the average tortuosity. As shown in this figure, decreasing the average tortuosity leads to a significant H

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

5. CONCLUSION Charge transport in BHJ solar cells was investigated by continuous time random walk Monte Carlo simulation based on a hopping model. Morphology, spatial distribution of traps, interdomain connectivity, and phase separation of active layer were considered in the simulations. Diffusion coefficient as the main parameter in charge transport and collection efficiency as an output characteristic were calculated. The blend ratio of the donor and acceptor species as a mixing parameter was considered in the simulations. It can be concluded from the simulation results that, by decreasing the blend ratio below the 50%, the diffusion coefficient decreases. Such a decrement, which is in agreement with experimental results, is a direct consequence of decreasing the accessible space for charge carriers in P3HT:PCBM. Tortuosity of the active layer which quantitatively indicates the convolution of the transport pathway was considered as an affecting parameter in charge transport. Tortuosity is a measure of ordering the bicontinuous structure of active layer in BHJ organic solar cells. Results showed that an increment in the tortuosity leads to a decrement in the diffusion coefficient and the collection efficiency of the cell. Accordingly, it can be concluded that the directional deposition of the polymers can improve the cell performance. Results showed that by decreasing the interdomain connectivity, the collection efficiency decreases dramatically which is in agreement with the experimental results. Orientation of the P3HT clusters relative to the substrate was also considered and it was obtained that the perpendicular P3HT lamellae orientation relative to the substrate has a lower collection efficiency comparing to the parallel P3HT lamellae orientation relative to the substrate.

increment in the collection efficiency. Such increment in the collection efficiency is a direct consequence of the ordering the bicontinuous structure of active layer and realizing the directional charge transport. Figure 9b also shows the effect of the P3HT clusters orientation on the collection efficiency. As shown in this figure, the perpendicular P3HT lamellae orientation relative to the substrate (edge-on) has lower collection efficiency comparing the parallel P3HT lamellae orientation relative to the substrate (face-on). The results of Figure 9b is in agreement with the previously reported experiments. For example, Crossland et al.109 have measured anisotropic charge-transport in P3HT films using field-effect transistor channels placed within well-aligned domains of quasi-parallel lamellae. They have found that the charge mobility is three times higher in the direction parallel to the ordered π-stacks comparing the direction perpendicular to them. The same results have been reported by Wu et al.110 Their results exhibit an improvement in the performance of the field effect transistor in which the direction of most of the (π−π) overlaps was consistent with that of the charge transport from source to drain. It is worth mentioning that the collection efficiency shown in Figure 9b is not the total efficiency of the cell. The measured efficiency in the experiments is affected by not only the charge collection but also the light absorption and successive injection. On the other word, the total efficiency of a solar cell is a multiplication of the collection, injection and absorption efficiencies. In order to investigate the impact of spatial distribution of the traps and the interdomain connectivity on the performance of the BHJ solar cells, we have carried out a set of simulations, and the results are provided in Figure 10. In this regard, we



AUTHOR INFORMATION

Corresponding Author

*(Y.A.) E-mail: [email protected]. ORCID

Yaser Abdi: 0000-0002-7583-7687 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We would like to thank the Iran National Science Foundation (INSF) for partial financial support.



REFERENCES

(1) Archer, M. D. Nanostructured and Photoelectrochemical Systems for Solar Photon Conversion; Imperial College Press: 2008; Vol. 3. (2) Snaith, H. J. Perovskites: The Emergence of a New Era for LowCost, High-Efficiency Solar Cells. J. Phys. Chem. Lett. 2013, 4, 3623− 3630. (3) Huang, H.; Huang, J. Organic and Hybrid Solar Cells; Springer: 2014. (4) Tress, W. Organic Solar Cells. Springer Ser. Mater. Sci. 2014, 208, 67−214. (5) Liu, Y.; Zhao, J.; Li, Z.; Mu, C.; Ma, W.; Hu, H.; Jiang, K.; Lin, H.; Ade, H.; Yan, H. Aggregation and Morphology Control Enables Multiple Cases of High-Efficiency Polymer Solar Cells. Nat. Commun. 2014, 5, 5293. (6) Arkhipov, V.; Heremans, P.; Emelianova, E.; Adriaenssens, G.; Bässler, H. Weak-Field Carrier Hopping in Disordered Organic Semiconductors: The Effects of Deep Traps and Partly Filled Densityof-States Distribution. J. Phys.: Condens. Matter 2002, 14, 9899.

Figure 10. Collection efficiency as a function of the percentage of the removed clusters for the network of face-on P3HT clusters with tortuosity of 1.01 and 1.32.

have removed some of the clusters randomly, and the collection efficiency of the cells was calculated as a function of the percentage of the removed clusters. Figure 10 shows that by decreasing the interdomain connectivity (increasing the percentage of the removed clusters), the collection efficiency decreases dramatically. Low efficiencies of BHJ solar cells reported in the experiments may come from such poor connectivity in the organic polymers. I

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C (7) Arkhipov, V.; Emelianova, E.; Adriaenssens, G.; Bässler, H. Equilibrium Carrier Mobility in Disordered Organic Semiconductors. J. Non-Cryst. Solids 2002, 299-302, 1047−1051. (8) Forrest, S. R. The Path to Ubiquitous and Low-Cost Organic Electronic Appliances on Plastic. Nature 2004, 428, 911. (9) Howard, I. A.; Mauer, R.; Meister, M.; Laquai, F. Effect of Morphology on Ultrafast Free Carrier Generation in Polythiophene: Fullerene Organic Solar Cells. J. Am. Chem. Soc. 2010, 132, 14866− 14876. (10) Watkins, P. K.; Walker, A. B.; Verschoor, G. L. Dynamical Monte Carlo Modelling of Organic Solar Cells: The Dependence of Internal Quantum Efficiency on Morphology. Nano Lett. 2005, 5, 1814−1818. (11) Nelson, J. Continuous-Time Random-Walk Model of Electron Transport in Nanocrystalline Tio 2 Electrodes. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 15374. (12) Nelson, J.; Chandler, R. E. Random Walk Models of Charge Transfer and Transport in Dye Sensitized Systems. Coord. Chem. Rev. 2004, 248, 1181−1194. (13) Anta, J. A.; Mora-Seró, I.; Dittrich, T.; Bisquert, J. Interpretation of Diffusion Coefficients in Nanostructured Materials from Random Walk Numerical Simulation. Phys. Chem. Chem. Phys. 2008, 10, 4478−4485. (14) Anta, J. A. Random Walk Numerical Simulation for Solar Cell Applications. Energy Environ. Sci. 2009, 2, 387−392. (15) Street, R.; Schoendorf, M.; Roy, A.; Lee, J. Interface State Recombination in Organic Solar Cells. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 81, 205307. (16) Beiley, Z. M.; Hoke, E. T.; Noriega, R.; Dacuña, J.; Burkhard, G. F.; Bartelt, J. A.; Salleo, A.; Toney, M. F.; McGehee, M. D. Morphology-Dependent Trap Formation in High Performance Polymer Bulk Heterojunction Solar Cells. Adv. Energy Mater. 2011, 1, 954−962. (17) Anta, J. A.; Nelson, J.; Quirke, N. Charge Transport Model for Disordered Materials: Application to Sensitized Tio 2. Phys. Rev. B: Condens. Matter Mater. Phys. 2002, 65, 125324. (18) Anta, J. A.; Casanueva, F.; Oskam, G. A Numerical Model for Charge Transport and Recombination in Dye-Sensitized Solar Cells. J. Phys. Chem. B 2006, 110, 5372−5378. (19) Baranovski, S. Charge Transport in Disordered Solids with Applications in Electronics; John Wiley & Sons: 2006; Vol. 17. (20) Nenashev, A.; Oelerich, J.; Baranovskii, S. Theoretical Tools for the Description of Charge Transport in Disordered Organic Semiconductors. J. Phys.: Condens. Matter 2015, 27, 093201. (21) Oelerich, J.; Huemmer, D.; Baranovskii, S. How to Find out the Density of States in Disordered Organic Semiconductors. Phys. Rev. Lett. 2012, 108, 226403. (22) Baranovskii, S. Theoretical Description of Charge Transport in Disordered Organic Semiconductors. Phys. Status Solidi B 2014, 251, 487−525. (23) Anta, J. A.; Morales-Flórez, V. Combined Effect of Energetic and Spatial Disorder on the Trap-Limited Electron Diffusion Coefficient of Metal-Oxide Nanostructures. J. Phys. Chem. C 2008, 112, 10287−10293. (24) Ansari-Rad, M.; Abdi, Y.; Arzi, E. Monte Carlo Random Walk Simulation of Electron Transport in Dye-Sensitized Nanocrystalline Solar Cells: Influence of Morphology and Trap Distribution. J. Phys. Chem. C 2012, 116, 3212−3218. (25) Abdi, N.; Abdi, Y.; Alemipour, Z.; NedaaeeOskoee, E. Chemical Diffusion Coefficient in Dye Sensitized Solar Cells as a Function of Porosity and Surface Roughness. Sol. Energy 2016, 135, 506−511. (26) Koster, L. Charge Carrier Mobility in Disordered Organic Blends for Photovoltaics. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 81, 205318. (27) Woellner, C. F.; Freire, J. A. Impact of the Intermixed Phase and the Channel Network on the Carrier Mobility of Nanostructured Solar Cells. J. Chem. Phys. 2016, 144, 084119.

(28) Woellner, C. F.; Li, Z.; Freire, J. A.; Lu, G.; Nguyen, T.-Q. Charge Carrier Mobility in a Two-Phase Disordered Organic System in the Low-Carrier Concentration Regime. Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 125311. (29) Heiber, M. C.; Kister, K.; Baumann, A.; Dyakonov, V.; Deibel, C.; Nguyen, T.-Q. Impact of Tortuosity on Charge-Carrier Transport in Organic Bulk Heterojunction Blends. Phys. Rev. Appl. 2017, 8, 054043. (30) Kim, Y.; Cook, S.; Tuladhar, S. M.; Choulis, S. A.; Nelson, J.; Durrant, J. R.; Bradley, D. D.; Giles, M.; McCulloch, I.; Ha, C.-S.; et al. A Strong Regioregularity Effect in Self-Organizing Conjugated Polymer Films and High-Efficiency Polythiophene: Fullerene Solar Cells. Nat. Mater. 2006, 5, 197−203. (31) Brabec, C. J.; Heeney, M.; McCulloch, I.; Nelson, J. Influence of Blend Microstructure on Bulk Heterojunction Organic Photovoltaic Performance. Chem. Soc. Rev. 2011, 40, 1185−1199. (32) Bisquert, J. Chemical Diffusion Coefficient of Electrons in Nanostructured Semiconductor Electrodes and Dye-Sensitized Solar Cells. J. Phys. Chem. B 2004, 108, 2323−2332. (33) Monroe, D. Hopping in Exponential Band Tails. Phys. Rev. Lett. 1985, 54, 146. (34) Bisquert, J. Hopping Transport of Electrons in Dye-Sensitized Solar Cells. J. Phys. Chem. C 2007, 111, 17163−17168. (35) Tiedje, T.; Rose, A. A Physical Interpretation of Dispersive Transport in Disordered Semiconductors. Solid State Commun. 1981, 37, 49−52. (36) Baranovskii, S.; Zvyagin, I.; Cordes, H.; Yamasaki, S.; Thomas, P. Electronic Transport in Disordered Organic and Inorganic Semiconductors. J. Non-Cryst. Solids 2002, 299-302, 416−419. (37) Jansson, F.; Baranovskii, S.; Gebhard, F.; Ö sterbacka, R. Effective Temperature for Hopping Transport in a Gaussian Density of States. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 77, 195211. (38) Rivnay, J.; Mannsfeld, S. C.; Miller, C. E.; Salleo, A.; Toney, M. F. Quantitative Determination of Organic Semiconductor Microstructure from the Molecular to Device Scale. Chem. Rev. 2012, 112, 5488−5519. (39) Chen, D.; Nakahara, A.; Wei, D.; Nordlund, D.; Russell, T. P. P3ht/Pcbm Bulk Heterojunction Organic Photovoltaics: Correlating Efficiency and Morphology. Nano Lett. 2011, 11, 561−567. (40) Liu, F.; Gu, Y.; Shen, X.; Ferdous, S.; Wang, H.-W.; Russell, T. P. Characterization of the Morphology of Solution-Processed Bulk Heterojunction Organic Photovoltaics. Prog. Polym. Sci. 2013, 38, 1990−2052. (41) Moulé, A. J.; Meerholz, K. Morphology Control in SolutionProcessed Bulk-Heterojunction Solar Cell Mixtures. Adv. Funct. Mater. 2009, 19, 3028−3036. (42) Yoshikawa, S.; Saeki, A.; Saito, M.; Osaka, I.; Seki, S. On the Role of Local Charge Carrier Mobility in the Charge Separation Mechanism of Organic Photovoltaics. Phys. Chem. Chem. Phys. 2015, 17, 17778−17784. (43) Verploegen, E.; Mondal, R.; Bettinger, C. J.; Sok, S.; Toney, M. F.; Bao, Z. Effects of Thermal Annealing Upon the Morphology of Polymer−Fullerene Blends. Adv. Funct. Mater. 2010, 20, 3519−3529. (44) Sirringhaus, H.; Brown, P. J.; Friend, R. H.; Nielsen, M. M.; Bechgaard, K.; Langeveld-Voss, B. M. W.; Spiering, A. J. H.; Janssen, R. A. J.; Meijer, E. W.; Herwig, P.; et al. Two-Dimensional Charge Transport in Self-Organized, High-Mobility Conjugated Polymers. Nature 1999, 401, 685. (45) Dennler, G.; Scharber, M. C.; Brabec, C. J. Polymer-Fullerene Bulk-Heterojunction Solar Cells. Adv. Mater. 2009, 21, 1323−1338. (46) Liao, H.-C.; Tsao, C.-S.; Shao, Y.-T.; Chang, S.-Y.; Huang, Y.C.; Chuang, C.-M.; Lin, T.-H.; Chen, C.-Y.; Su, C.-J.; Jeng, U.-S.; et al. Bi-Hierarchical Nanostructures of Donor−Acceptor Copolymer and Fullerene for High Efficient Bulk Heterojunction Solar Cells. Energy Environ. Sci. 2013, 6, 1938−1948. (47) Liao, H.-C.; Tsao, C.-S.; Lin, T.-H.; Chuang, C.-M.; Chen, C.Y.; Jeng, U.-S.; Su, C.-H.; Chen, Y.-F.; Su, W.-F. Quantitative Nanoorganized Structural Evolution for a High Efficiency Bulk J

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C Heterojunction Polymer Solar Cell. J. Am. Chem. Soc. 2011, 133, 13064−13073. (48) Huang, Y.; Kramer, E. J.; Heeger, A. J.; Bazan, G. C. Bulk Heterojunction Solar Cells: Morphology and Performance Relationships. Chem. Rev. 2014, 114, 7006−7043. (49) Erb, T.; Zhokhavets, U.; Gobsch, G.; Raleva, S.; Stühn, B.; Schilinsky, P.; Waldauf, C.; Brabec, C. J. Correlation between Structural and Optical Properties of Composite Polymer/Fullerene Films for Organic Solar Cells. Adv. Funct. Mater. 2005, 15, 1193− 1196. (50) Huang, F.; Yip, H.-L.; Cao, Y. Polymer Photovoltaics: Materials, Physics, and Device Engineering; Royal Society of Chemistry: 2015. (51) Osaka, I.; Kakara, T.; Takemura, N.; Koganezawa, T.; Takimiya, K. Naphthodithiophene−Naphthobisthiadiazole Copolymers for Solar Cells: Alkylation Drives the Polymer Backbone Flat and Promotes Efficiency. J. Am. Chem. Soc. 2013, 135, 8834−8837. (52) Dijkstra, E. W. A Note on Two Problems in Connexion with Graphs. Numerische mathematik 1959, 1, 269−271. (53) Morales-Flórez, V.; Pinero, M.; De La Rosa-Fox, N.; Esquivias, L.; Anta, J. A.; Primera, J. The Cluster Model: A HierarchicallyOrdered Assemblage of Random-Packing Spheres for Modelling Microstructure of Porous Materials. J. Non-Cryst. Solids 2008, 354, 193−198. (54) Berriman, G. A.; Holmes, N. P.; Holdsworth, J. L.; Zhou, X.; Belcher, W. J.; Dastoor, P. C. A New Model for Pcbm Phase Segregation in P3ht: Pcbm Blends. Org. Electron. 2016, 30, 12−17. (55) Bässler, H. Charge Transport in Disordered Organic Photoconductors a Monte Carlo Simulation Study. Phys. Status Solidi B 1993, 175, 15−56. (56) Deibel, C.; Dyakonov, V. Polymer−Fullerene Bulk Heterojunction Solar Cells. Rep. Prog. Phys. 2010, 73, 096401. (57) Kohary, K.; Cordes, H.; Baranovskii, S.; Thomas, P.; Yamasaki, S.; Hensel, F.; Wendorff, J.-H. One-Dimensional Hopping Transport in Disordered Organic Solids. Ii. Monte Carlo Simulations. Phys. Rev. B: Condens. Matter Mater. Phys. 2001, 63, 094202. (58) Laquai, F.; Andrienko, D.; Mauer, R.; Blom, P. W. Charge Carrier Transport and Photogeneration in P3ht: Pcbm Photovoltaic Blends. Macromol. Rapid Commun. 2015, 36, 1001−1025. (59) Bartelt, J. A.; Lam, D.; Burke, T. M.; Sweetnam, S. M.; McGehee, M. D. Charge-Carrier Mobility Requirements for Bulk Heterojunction Solar Cells with High Fill Factor and External Quantum Efficiency> 90%. Adv. Energy Mater. 2015, 5, 1500577. (60) Foster, S.; Deledalle, F.; Mitani, A.; Kimura, T.; Kim, K. B.; Okachi, T.; Kirchartz, T.; Oguma, J.; Miyake, K.; Durrant, J. R.; et al. Electron Collection as a Limit to Polymer: Pcbm Solar Cell Efficiency: Effect of Blend Microstructure on Carrier Mobility and Device Performance in Ptb7: Pcbm. Adv. Energy Mater. 2014, 4, 1400311. (61) Tessler, N.; Preezant, Y.; Rappaport, N.; Roichman, Y. Charge Transport in Disordered Organic Materials and Its Relevance to ThinFilm Devices: A Tutorial Review. Adv. Mater. 2009, 21, 2741−2761. (62) Mott, N. F.; Twose, W. The Theory of Impurity Conduction. Sir Nevill Mott−65 Years in Physics 1995, 12, 259−316. (63) Miller, A.; Abrahams, E. Impurity Conduction at Low Concentrations. Phys. Rev. 1960, 120, 745. (64) Walker, A.; Kambili, A.; Martin, S. Electrical Transport Modelling in Organic Electroluminescent Devices. J. Phys.: Condens. Matter 2002, 14, 9825. (65) Baranovskii, S.; Rubel, O.; Thomas, P. Theoretical Description of Hopping Transport in Disordered Materials. Thin Solid Films 2005, 487, 2−7. (66) Mandoc, M.; Kooistra, F.; Hummelen, J.; De Boer, B.; Blom, P. Effect of Traps on the Performance of Bulk Heterojunction Organic Solar Cells. Appl. Phys. Lett. 2007, 91, 263505. (67) Ansari-Rad, M.; Garcia-Belmonte, G.; Bisquert, J. Consistent Formulation of the Crossover from Density to Velocity Dependent Recombination in Organic Solar Cells. Appl. Phys. Lett. 2015, 107, 073301.

(68) Wetzelaer, G. J. A.; Kuik, M.; Blom, P. W. Identifying the Nature of Charge Recombination in Organic Solar Cells from ChargeTransfer State Electroluminescence. Adv. Energy Mater. 2012, 2, 1232−1237. (69) Wagenpfahl, A. Mobility Dependent Recombination Models for Organic Solar Cells. J. Phys.: Condens. Matter 2017, 29, 373001. (70) Sun, L.; Sun, J.-X.; Xiong, C.-H.; Shi, X.-H. Trap-Assisted Recombination in Disordered Organic Semiconductors Extended by Considering Density Dependent Mobility. Sol. Energy 2016, 135, 308−316. (71) MacKenzie, R. C.; Kirchartz, T.; Dibb, G. F.; Nelson, J. Modeling Nongeminate Recombination in P3ht: Pcbm Solar Cells. J. Phys. Chem. C 2011, 115, 9806−9813. (72) Nalwa, K. S.; Kodali, H. K.; Ganapathysubramanian, B.; Chaudhary, S. Dependence of Recombination Mechanisms and Strength on Processing Conditions in Polymer Solar Cells. Appl. Phys. Lett. 2011, 99, 263301. (73) Arkhipov, V.; Heremans, P.; Bässler, H. Why Is Exciton Dissociation So Efficient at the Interface between a Conjugated Polymer and an Electron Acceptor? Appl. Phys. Lett. 2003, 82, 4605− 4607. (74) Ripolles-Sanchis, T.; Raga, S. R.; Guerrero, A.; Welker, M.; Turbiez, M.; Bisquert, J.; Garcia-Belmonte, G. Molecular Electronic Coupling Controls Charge Recombination Kinetics in Organic Solar Cells of Low Bandgap Diketopyrrolopyrrole, Carbazole, and Thiophene Polymers. J. Phys. Chem. C 2013, 117, 8719−8726. (75) Brückner, C.; Würthner, F.; Meerholz, K.; Engels, B. Atomistic Approach to Simulate Processes Relevant for the Efficiencies of Organic Solar Cells as a Function of Molecular Properties. Ii. Kinetic Aspects. J. Phys. Chem. C 2017, 121, 26−51. (76) Guerrero, A.; Marchesi, L. F.; Boix, P. P.; Bisquert, J.; GarciaBelmonte, G. Recombination in Organic Bulk Heterojunction Solar Cells: Small Dependence of Interfacial Charge Transfer Kinetics on Fullerene Affinity. J. Phys. Chem. Lett. 2012, 3, 1386−1392. (77) Kuik, M.; Koster, L.; Wetzelaer, G.; Blom, P. Trap-Assisted Recombination in Disordered Organic Semiconductors. Phys. Rev. Lett. 2011, 107, 256805. (78) Kirchartz, T.; Pieters, B. E.; Kirkpatrick, J.; Rau, U.; Nelson, J. Recombination Via Tail States in Polythiophene: Fullerene Solar Cells. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 83, 115209. (79) MacKenzie, R. C.; Shuttle, C. G.; Chabinyc, M. L.; Nelson, J. Extracting Microscopic Device Parameters from Transient Photocurrent Measurements of P3ht: Pcbm Solar Cells. Adv. Energy Mater. 2012, 2, 662−669. (80) Durrant, J. R.; Haque, S. A.; Palomares, E. Photochemical Energy Conversion: From Molecular Dyads to Solar Cells. Chem. Commun. 2006, 3279−3289. (81) Tachibana, Y.; Haque, S. A.; Mercer, I. P.; Moser, J. E.; Klug, D. R.; Durrant, J. R. Modulation of the Rate of Electron Injection in DyeSensitized Nanocrystalline Tio2 Films by Externally Applied Bias. J. Phys. Chem. B 2001, 105, 7424−7431. (82) Scully, S. R. Exciton Transport, Charge Extraction, and Loss Mechanisms in Organic Photovoltaics; Stanford University: 2009. (83) Garcia-Belmonte, G.; Boix, P. P.; Bisquert, J.; Sessolo, M.; Bolink, H. J. Simultaneous Determination of Carrier Lifetime and Electron Density-of-States in P3ht: Pcbm Organic Solar Cells under Illumination by Impedance Spectroscopy. Sol. Energy Mater. Sol. Cells 2010, 94, 366−375. (84) Jansson, F.; Nenashev, A.; Baranovskii, S.; Gebhard, F.; Ö sterbacka, R. Effect of Electric Field on Diffusion in Disordered Materials. Ann. Phys. 2009, 18, 856−862. (85) Mohan, S. R.; Singh, M. P.; Joshi, M.; Kukreja, L. Monte Carlo Simulation of Carrier Diffusion in Organic Thin Films with Morphological Inhomogeneity. J. Phys. Chem. C 2013, 117, 24663− 24672. (86) Martens, H.; Blom, P.; Schoo, H. Comparative Study of Hole Transport in Poly (P-Phenylene Vinylene) Derivatives. Phys. Rev. B: Condens. Matter Mater. Phys. 2000, 61, 7489. K

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C (87) Tanase, C.; Blom, P.; De Leeuw, D.; Meijer, E. Charge Carrier Density Dependence of the Hole Mobility in Poly (P-Phenylene Vinylene). Phys. Status Solidi A 2004, 201, 1236−1245. (88) Arkhipov, V.; Heremans, P.; Emelianova, E.; Baessler, H. Effect of Doping on the Density-of-States Distribution and Carrier Hopping in Disordered Organic Semiconductors. Phys. Rev. B: Condens. Matter Mater. Phys. 2005, 71, 045214. (89) Craciun, N.; Wildeman, J.; Blom, P. Universal Arrhenius Temperature Activated Charge Transport in Diodes from Disordered Organic Semiconductors. Phys. Rev. Lett. 2008, 100, 056601. (90) Germs, W. C.; Van der Holst, J.; Van Mensfoort, S.; Bobbert, P.; Coehoorn, R. Modeling of the Transient Mobility in Disordered Organic Semiconductors with a Gaussian Density of States. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, 165210. (91) Bässler, H. Localized States and Electronic Transport in Single Component Organic Solids with Diagonal Disorder. Phys. Status Solidi B 1981, 107, 9−54. (92) Brédas, J.-L.; Beljonne, D.; Coropceanu, V.; Cornil, J. ChargeTransfer and Energy-Transfer Processes in Π-Conjugated Oligomers and Polymers: A Molecular Picture. Chem. Rev. 2004, 104, 4971− 5004. (93) Guha, S.; Rice, J.; Yau, Y.; Martin, C. M.; Chandrasekhar, M.; Chandrasekhar, H. R.; Guentner, R.; De Freitas, P. S.; Scherf, U. Temperature-Dependent Photoluminescence of Organic Semiconductors with Varying Backbone Conformation. Phys. Rev. B: Condens. Matter Mater. Phys. 2003, 67, 125204. (94) Bisquert, J.; Fabregat-Santiago, F.; Mora-Sero, I.; GarciaBelmonte, G.; Barea, E. M.; Palomares, E. A Review of Recent Results on Electrochemical Determination of the Density of Electronic States of Nanostructured Metal-Oxide Semiconductors and Organic Hole Conductors. Inorg. Chim. Acta 2008, 361, 684−698. (95) Abdi, N.; Abdi, Y.; Alemipour, Z. Effects of Morphology and Thickness of Al2o3 Scaffold on Charge Transport in Perovskite-Based Solar Cells. Sol. Energy 2017, 153, 379−382. (96) Rubel, O.; Baranovskii, S.; Thomas, P.; Yamasaki, S. Concentration Dependence of the Hopping Mobility in Disordered Organic Solids. Phys. Rev. B: Condens. Matter Mater. Phys. 2004, 69, 014206. (97) Gonzalez-Vazquez, J.; Anta, J. A.; Bisquert, J. Random Walk Numerical Simulation for Hopping Transport at Finite Carrier Concentrations: Diffusion Coefficient and Transport Energy Concept. Phys. Chem. Chem. Phys. 2009, 11, 10359−10367. (98) Van der Auweraer, M.; De Schryver, F. C.; Borsenberger, P. M.; Bässler, H. Disorder in Charge Transport in Doped Polymers. Adv. Mater. 1994, 6, 199−213. (99) Mihailetchi, V. D.; van Duren, J. K.; Blom, P. W.; Hummelen, J. C.; Janssen, R. A.; Kroon, J. M.; Rispens, M. T.; Verhees, W. J. H.; Wienk, M. M. Electron Transport in a Methanofullerene. Adv. Funct. Mater. 2003, 13, 43−46. (100) Benkstein, K. D.; Kopidakis, N.; Van de Lagemaat, J.; Frank, A. Influence of the Percolation Network Geometry on Electron Transport in Dye-Sensitized Titanium Dioxide Solar Cells. J. Phys. Chem. B 2003, 107, 7759−7767. (101) Ben-Avraham, D.; Havlin, S. Diffusion and Reactions in Fractals and Disordered Systems; Cambridge University Press: 2000. (102) Ofir, A.; Dor, S.; Grinis, L.; Zaban, A.; Dittrich, T.; Bisquert, J. Porosity Dependence of Electron Percolation in Nanoporous Ti O 2 Layers. J. Chem. Phys. 2008, 128, 064703. (103) Dittrich, T.; Ofir, A.; Tirosh, S.; Grinis, L.; Zaban, A. Influence of the Porosity on Diffusion and Lifetime in Porous Ti O 2 Layers. Appl. Phys. Lett. 2006, 88, 182110. (104) Chiu, M.-Y.; Jeng, U.-S.; Su, M.-S.; Wei, K.-H. Morphologies of Self-Organizing Regioregular Conjugated Polymer/Fullerene Aggregates in Thin Film Solar Cells. Macromolecules 2010, 43, 428−432. (105) Van Bavel, S. S.; Bärenklau, M.; de With, G.; Hoppe, H.; Loos, J. P3ht/Pcbm Bulk Heterojunction Solar Cells: Impact of Blend Composition and 3d Morphology on Device Performance. Adv. Funct. Mater. 2010, 20, 1458−1463.

(106) Zhao, J.; Swinnen, A.; Van Assche, G.; Manca, J.; Vanderzande, D.; Mele, B. V. Phase Diagram of P3ht/Pcbm Blends and Its Implication for the Stability of Morphology. J. Phys. Chem. B 2009, 113, 1587−1591. (107) Müller, C.; Ferenczi, T. A.; Campoy-Quiles, M.; Frost, J. M.; Bradley, D. D.; Smith, P.; Stingelin-Stutzmann, N.; Nelson, J. Binary Organic Photovoltaic Blends: A Simple Rationale for Optimum Compositions. Adv. Mater. 2008, 20, 3510−3515. (108) Javadi, M.; Abdi, Y. Monte Carlo Random Walk Simulation of Electron Transport in Confined Porous Tio2 as a Promising Candidate for Photo-Electrode of Nano-Crystalline Solar Cells. J. Appl. Phys. 2015, 118, 064304. (109) Crossland, E. J.; Tremel, K.; Fischer, F.; Rahimi, K.; Reiter, G.; Steiner, U.; Ludwigs, S. Anisotropic Charge Transport in Spherulitic Poly (3-Hexylthiophene) Films. Adv. Mater. 2012, 24, 839−844. (110) Wu, F.-C.; Lu, C.-C.; Ruan, J.; Tang, F.-C.; Cheng, H.-L.; Chou, W.-Y. Enhanced and Anisotropic Charge Transport in Polymer-Based Thin-Film Transistors by Guiding Polymer Growth. Cryst. Growth Des. 2017, 17, 629−636.

L

DOI: 10.1021/acs.jpcc.8b04962 J. Phys. Chem. C XXXX, XXX, XXX−XXX