DOI:
10.1039/C5RA12318F
(Paper)
RSC Adv., 2015,
5, 70939-70948
The prediction of the morphology and PCE of small molecular organic solar cells
Received
25th June 2015
, Accepted 10th August 2015
First published on 11th August 2015
Abstract
The efficiency of bulk heterojunction (BHJ) solar cells depends strongly on the morphology of the electron donors and electron acceptors in the active layer. Here we use Dissipative Particle Dynamics (DPD) simulation to predict the donor–acceptor morphology and graph theory to predict the efficiency of small molecular organic solar cells (SM OSCs). We focus on a recently reported small molecular organic solar cell based on a new molecule donor, DTS(PTTh2)2 and three molecules DR3TBDTT, DR3TBDTT-HD, and DR3TBD2T with a benzo[1,2-b:4,5-b′]dithiophene (BDT) unit. With our theoretical approach, we are able to study the critical factors affecting the morphology and efficiency such as, the chemical structure of the conjugated molecular, fullerene functional group, solvents, additives, blend ratio, and processing conditions (e.g., annealing temperature). Our results are consistent with experimental conclusions and provide useful guidelines for improving the efficiency of small organic solar cells.
1 Introduction
In the circumstance of reduced availability of energy, as a new generation of solar cells, organic solar cells (OSCs) are attracting more and more attention worldwide due to their potential to be cost effective and flexible solar energy conversion devices.1–4 In the past few years, many research efforts have been focused on polymer solar cells and impressive results were achieved.5–10 To date, power conversion efficiencies (PCE) over 10% have been obtained for both solution-processed single-junction polymer solar cells and the tandem devices.2 Small molecule OSCs are of great interest to the research community owing to the distinct advantages of small-molecule semiconductors, such as well-defined structures, and batch-to-batch stability.3 With consistent effort, the PCE of small-molecule OSCs had been improved from 0.03% to 10.08% in the past few years.11–25 However, there is still much work needed to be done to overcome the hurdles of low PCE of small-molecule OSCs.
In order to find ways to overcome these limitations, computational approaches have been utilized for deciphering the roles of different variables that could improve the PCE and lifetime of an OPV cell.26 A number of such studies have focused on some or all of the pertinent steps in solar energy conversion in an OPV device. These steps in chronological order are: (1) the absorption of light and excitons (i.e., electron–hole pair) generation; (2) the diffusion of the excitons to an interface; (3) the separation of the excitons into electrons and holes; (4) the transport of these electrons and holes to their corresponding electrode; and (5) the transfer of electrons and holes from the electrode to an external circuit. It is believed that the third step needs optimization in a bulk heterojunction (BHJ) design.26 Here the design consists of a mixture of an electron donor and an electron acceptor (fullerene based molecules), which tends to form a phase segregated morphology. The phase segregation leads to a marked increase in the interfacial area between the donor and acceptor thus facilitating excitons separation. However, the ultimate goal is to achieve sizes of the donor and electron acceptor domains commensurate with the excitons diffusion length so that a higher number of excitons can reach the donor–acceptor interface. The importance of the morphology of polymer bulk heterojunctions could be found in the reviews, such as, Liu et al.27 and Ruderer and Müller-Buschbaum.28
Recently, three molecules named DR3TBDT-HD, DR3TBDTT and DR3TBD2T with a benzo[1,2-b:4,5-b′]dithiophene (BDT) unit and fullerene derivatives have been used as electron donor and acceptor materials.29,30 The effects of blend ratio,31 additives32,33 have been investigated to improve its morphology and efficiency. However, it is not clear how the chemical structures and the processing conditions affect the morphology and the efficiency. Here we combine DPD simulations and graph theory to predict the morphology and the efficiency of these small molecules organic solar cells. We study the effects of blend ratio of donor/PC71BM, the annealing temperature of system, the additives, and the desolvation process. Our results are consistent with experimental conclusions and provide useful guidelines for improving the efficiency of small organic solar cell.
2 Simulation methods
2.1 DPD simulation
DPD was first proposed by Hoogerbrugge and Koelman34,35 and later improved by Groot et al.36 It is essentially a coarse-grained molecular dynamics (CGMD) method with a stochastic thermostat. In the DPD method, each bead represents a group of atoms and their movements are governed by: |
![image file: c5ra12318f-t1.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t1.gif) | (1) |
Where vi, ri and Fi represent the velocity, position vectors and total force acting on bead, respectively. Due to the significant reduction in the number of degrees of freedom, DPD can be used to simulate complex systems for much longer time. The forces acting on each bead can be divided into three components: conservative forces (FCij), dissipative forces (FDij) and random forces (FRij). |
![image file: c5ra12318f-t2.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t2.gif) | (2) |
All the three forces act only within a cut-off radius rc, which is taken as the unit of length in the simulation. The conservative force represents a soft repulsive interaction in linear form:
|
![image file: c5ra12318f-t3.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t3.gif) | (3) |
Where
rij =
rj −
ri,
rij = |
rij| and
ij =
rij/|
rij|. The parameters
aij, which represent the maximum repulsion between beads
i and
j, can be related to the Flory–Huggins parameter as follows:
|
![image file: c5ra12318f-t4.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t4.gif) | (4) |
Where
ρ represents the bead density of the system and
aii is the repulsion parameter between the same beads. By matching the compressibility of water,
aii and
ρ are related through
|
![image file: c5ra12318f-t5.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t5.gif) | (5) |
The dissipative forces FDij are the drag forces and the random forces FRij represent the thermal noise of the system:
|
![image file: c5ra12318f-t6.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t6.gif) | (6) |
|
![image file: c5ra12318f-t7.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t7.gif) | (7) |
Where
γ is the friction parameter,
σ is the noise amplitude parameter and
vij =
vi −
vj. The following equations should be satisfied to ensure the canonical ensemble equilibrium conditions and the fluctuation–dissipation relation:
|
![image file: c5ra12318f-t8.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t8.gif) | (8) |
We use γ = 4.5, kBT = 1, and thus σ = 3 according to eqn (9). FSi acts between consecutive beads within the same molecule, and has the following form:
|
![image file: c5ra12318f-t9.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t9.gif) | (10) |
Where
C is the spring constant, and
C is set to be 4 in this paper.
37
2.2 Predict the efficiency by using graph theory based on the morphology of the active layer
To improve the efficiency of organic solar cells, it is essential to understand the role of morphology and tailor fabrication process to get desired morphologies. Wodo et al.38 provide a comprehensive suite of physically meaningful morphology descriptors, which reflect the complex nature of the BHJ and the underlying device physics. They formulate and implement a graph-based framework to efficiently construct this comprehensive suite. This approach is motivated by the equivalence between a discretized 2D/3D morphology image. These morphology descriptors are further classified according to the physical processes of the photovoltaic process: excitons diffusion, excitons dissociation into free charge carriers, and charge transport. We obtain the 3D morphology from our DPD simulation firstly. Then we use the following three descriptors38 to characterize the discrete morphology. The descriptors are fabs, fdiss, fout.
fabs Fraction of donor materials. Light is absorbed by the electron donor to generate the excited hole–electron pairs (excitons) in the active layer. We should find the fraction of donor materials.
|
![image file: c5ra12318f-t10.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t10.gif) | (11) |
If there is only donor and acceptor in the morphology, we could use fabs respectively. However, our study involves solvent/additive in addition to donor and acceptor. It is necessary for us to modify descriptor fabs to indicate the fraction of average donor/acceptor density.
|
![image file: c5ra12318f-t11.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t11.gif) | (12) |
(
D: the average density of donor;
A: the average of density of acceptor).
fdiss Fraction of photoactive material within 10 nm distance (excitons diffusion length) to the interface.
|
![image file: c5ra12318f-t12.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t12.gif) | (13) |
By finding the fraction of donor material within this distance to the donor–acceptor interface, we identify the fraction of donor material that can theoretically contribute to this stage. Donor material distributed beyond this distance to the interface has a significantly lower chance of being useful.
fout Fraction of materials with useful domains connected to the electrode.
|
![image file: c5ra12318f-t13.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t13.gif) | (14) |
Only regions with a connection to the electrode constitute paths for charge to travel and consequently can contribute to current generation. If charges are created in islands or “cul-de-sacs” with no direct connection to appropriate electrode, they eventually recombine.
Wodo et al. further define a “performance indicator” as product of these three fractions.
They imply that f is a good marker for Jsc. Jsc is a direct measurement of the device performance. In addition, the magnitude of short circuit current density is a quantity which is strongly affected by the morphology.39,40
The size of the interface and the domain size affect the charge collection efficiency. When the interface area is big, the domain size will be small. There is an optimum interface area of the morphology for the efficiency. Young Min Nam et al.’s experimental results41 show that, when the domain size is about 6 nm, the active layer shows the optimum PCE. Meanwhile, Peter K. Watkins et al.42 use dynamical Monte Carlo model to evaluate the relationship between interface area and efficiency. Watkins et al. evaluated the “excitons dissociation efficiency” and “charge collection efficiency” separately and they conclude that there is an optimum interface area. We use the relationship of interface area and internal quantum efficiency derived from Watkins et al. and the optimum domain size 6 nm based on Young Min Nam et al.’s experimental results. And the descriptor is defined as the following:
|
![image file: c5ra12318f-t14.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t14.gif) | (16) |
A is the interface area.
When we use graph theory to estimate the efficiency, we use the empirical factors (fabs, fdiss, fout, finterface) to represent different factors affecting the efficiency. We don't explicitly consider the process of the generation of the excition, the dissociation of the excition, and the charge transfer. Thus our model doesn't include all the factors affecting the final efficiency and overestimates the final efficiency. For example, Lin et al. explicitly investigate the charge-transfer processes at pentacene-C60 interface by using simulations.43 They found that the positional relationship between donor and acceptor materials plays a critical role for the charge-transfer process. It will be interesting and challenging to simulate the charge transfer process of small molecule OPV explicitly and incorporate the factor into our model to improve the accuracy of the estimated efficiency.
In summary, the power conversion efficiencies based on morphology descriptors is defined as:
|
PCE = f′absfdissfoutfinterface
| (17) |
3 Models and parameters
We model donor and acceptor molecules in a coarse-grained manner to simulate these systems on a length scale that is large enough to examine the blend morphology, while retaining necessary details at molecular level. As shown in Fig. 1, we model a donor using four coarse-grained beads. The acceptor or the solvent is represented by one bead. The additive (PDMS) is represented by three beads.
![image file: c5ra12318f-f1.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-f1.gif) |
| Fig. 1 Chemical structures of donor, acceptor, solvent and additive are represented by coarse-grained beads. (a) Chemical structures of C1 (DR3TBDTT-HD), C2 (DR3TBDTT), C3 (DR3TBDT2T)29 are represented by coarse-grained beads A, B, D and C. (b) The fullerene derivative is represented by one coarse-grained bead Z. (c) The additive PDMS is represented by three coarse-grained beads, G, F, F. (d) The solvent chloroform is represented by one coarse-grained bead S. | |
To obtain the solubility parameters and interaction potential parameters of particles, we construct the amorphous structures of each fragment and perform molecular dynamics with COMPASS force-field44,45 by using the Forcite module of Materials Studio software package to obtain cohesive energy, solubility parameters (δi)46–48 and volume. According to Groot et al.,49 we choose a standard volume of bead Vm, and this DPD bead is equal to Nm (real-space renormalization factor) water molecules (30 Å3). Then we obtain the length scale, time scale and the number of step. The formulas are listed as following:
|
rc = (ρVbead)1/3 = (3 × 100)1/3 = 6.7 Å = 0.67 nm
| (18) |
|
![image file: c5ra12318f-t15.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t15.gif) | (19) |
|
![image file: c5ra12318f-t16.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t16.gif) | (20) |
The Newtonian equation of position and velocity of particles is solved by a modified version of the velocity Verlet algorithm.50 In the simulation, the radius of interaction, the particle mass, and the temperature were chosen as rc = m = kT = 1 and σ = 3.67, while the particle density ρ = 3 (taking into account the computational efficiency, ρ = 3 is a reasonable choice).34,35 The only parameter to be determined is the maximum repulsive force aij, which is chosen according to the linear relation with Flory–Huggins χ parameters.51
|
aij = aii + 3.27χ, aii = 25
| (21) |
The Flory–Huggins χ parameter can be obtained from experiments, free energy of mixing, or solubility parameter.52,53 Here we obtain the χ through the solubility parameter (Table 1). And we obtain the solubility parameters from molecular dynamics (MD) simulations. A high-quality force field COMPASS (condensed-phase optimized molecular potentials for atomistic simulation studies)54–56 is adopted, which is known to accurately reproduce experimental structures, densities, and solubility parameters.57 MD simulations are performed under NPT thermodynamic ensemble. In our MD simulations, we obtain the cohesive energy (δ).58 Then we get the Flory–Huggins χ parameter by the eqn (22).
|
![image file: c5ra12318f-t17.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-t17.gif) | (22) |
Where
V is the arithmetic average of molar volumes of two beads,
δi and
δj are the solubility parameters of beads
i and
j, respectively.
R is perfect gas constant.
T is the temperature of system. Then we obtain the repulsive force
aij according to the
eqn (21). The calculation results are listed in
Table 3.
Table 1 Parameters derived for the coarse-grained beads defined in Fig. 1
|
A |
B |
C1 |
C2 |
C3 |
D |
S |
F |
G |
Z |
Solubility parameter |
17.6 |
14.8 |
15.8 |
15.1 |
16.1 |
9.26 |
15.5 |
11.3 |
10.02 |
20.5 |
Vmol (Å) |
170 |
82 |
923.5 |
579 |
1048.6 |
83.68 |
96.74 |
141.01 |
126.4 |
1003 |
M (g mol−1) |
130 |
84 |
872 |
578 |
956 |
29 |
29 |
89 |
58 |
1030 |
Table 2 The parameters of DPD simulation
Parameter |
Value |
Meaning |
aii |
25 |
Repulsion between identical particle types |
T* |
1.0 |
Reduced temperature |
Vm |
100 Å3 |
The volume of standard bead |
Nm |
3 |
Number of water molecules represented by one bead |
δt |
0.02 |
Time step |
T |
103.55 ps |
Characteristic time scale |
rc |
0.67 nm |
Range of soft repulsive interaction |
N |
25000 step |
The number of step |
CS |
4.0 |
Spring constant |
Table 3 The interaction parameters aij between beads at 346.4 K
|
A |
B |
C |
D |
F |
G |
S |
Z |
A |
25.00 |
|
|
|
|
|
|
|
B |
26.01 |
25.00 |
|
|
|
|
|
|
C |
26.86 |
25.52 |
25.00 |
|
|
|
|
|
D |
34.48 |
37.68 |
47.56 |
25.00 |
|
|
|
|
F |
31.87 |
25.57 |
25.14 |
25.43 |
25.00 |
|
|
|
G |
33.83 |
27.55 |
43.64 |
25.10 |
25.16 |
25.00 |
|
|
S |
25.62 |
25.05 |
25.05 |
28.70 |
27.40 |
28.50 |
25.00 |
|
Z |
28.92 |
37.75 |
43.64 |
58.35 |
45.83 |
27.40 |
35.10 |
25.00 |
In DPD method, the simulation length scale and time scale of molecular motions are much larger than ordinary molecular dynamics (MD) simulations. We choose the particle mass and the temperature and interaction range as units of mass, energy, and length, hence rc = m = kT = 1, and the simulated time is expressed in the natural unit of time.49 The simulation boxes of sizes in our DPD simulation is 32 × 32 × 32 rc3. The parameters of DPD simulations are summarized in Table 2.
4. Result and discussion
4.1 The effect of the blend ratio on the morphology and the efficiency
The blend ratio of donor and acceptor is one of the most important factors affecting the morphology of BHJ solar cells. Jiaoyan Zhou et al.,29 Yanming Sun59 and Aung Ko Ko Kyaw60 illustrated the importance of the ratio of small molecular donor/PC71BM. Here we study the morphology and domain size of small molecular donor/PC71BM for different composition ratio. Fig. 2 shows the equilibrated morphology for different molar percentage of PC71BM: 33.3%, 50%, 66.7%, and 75%. In addition, we studied three different small molecules C1, C2, C3 as the donor, as shown in Fig. 2a–c respectively. The predicted efficiencies are summarized in Fig. 2a–c. The domain sizes are summarized in Fig. 2d–f.
![image file: c5ra12318f-f2.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-f2.gif) |
| Fig. 2 (a) The dependence of Power Conversion Efficiency (PCE) on the molar percentage of PC71BM for C1/PC71BM; (b) the dependence of Power Conversion Efficiency (PCE) on the molar percentage of PC71BM for C2/PC71BM; (c) the dependence of Power Conversion Efficiency (PCE) on the molar percentage of PC71BM for C3/PC71BM; (d) the dependence of the domain size on the molar percentage of PC71BM for C1/PC71BM; (e) the dependence of the domain size on the molar percentage of PC71BM for C2/PC71BM; (f) the dependence of the domain size on the molar percentage of PC71BM for C3/PC71BM; domain size is derived from the volume divided by the interface area. The optimum domain size favorable for efficiency is around 6 nm. | |
Fig. 2a–c show consistent results that, when percentage of PC71BM is 66.7%, the efficiency is the highest. By analyzing the domain sizes as shown in Fig. 2d–f, we can see that when the domain size is around 6 nm, the efficiency is good. In order to rationalize the different factors affecting the efficiency, we list the factors in Table 4.
Table 4 Morphology descriptors depends on the blend ratio of C1/PC71BM
Factor |
Molar percentage of PC71BM |
33.3% |
50% |
66.7% |
75% |
fabs |
Fraction of light absorbing material |
89.6% |
79.5% |
64.8% |
54.8% |
fdiss |
Fraction of photoactive material within d < 10 nm to the interface |
81.8% |
83.0% |
87.5% |
87.4% |
fout |
Fraction of useful domains connected to the electrode |
6.8% |
22.9% |
95.1% |
97.8% |
F |
f = fabs* fdiss* fout |
5.0% |
15.1% |
54.0% |
46.8% |
PCE |
Power conversion efficiency |
0.5% |
2.3% |
10.3% |
9.7% |
As shown in Table 4, among the various factors, fout is the most important factor for different blend ratios. When the percentage of PC71BM is 33.3% or 50%, the islands or ‘cul-de-sacs’ morphology is formed. Such morphology is poor for the fraction of useful domains connected to the electrode and leads to poor efficiency. Thus for 33.3% and 50%, fout is low and the final efficiency is low. Meanwhile, when the percentage of PC71BM is high, such as 75%, fabs is limited, since the donor is responsible for absorbing the light. In conclusion, 66.7% percentage shows the best efficiency.
When Zhou et al.29 made the device of C1, C2, C3 with PC71BM, they chose the blend ratio as (1
:
0.8 w/w). It corresponds to the volume percentage of PC71BM: 66.7%, which is the best blend ratio based on our results. Aung Ko Ko Kyaw et al.60 studied the dependence on the percentage of PCBM. They concluded that, when the weight: weight ratio between DTS(PTTh2)2 and PCBM is 70
:
30 or 60
:
40 shows the best mobility and performance, which is close to our best blend ratio for C1, C2, C3 with PC71BM: (1
:
0.8 w/w).
4.2 The effect of the annealing temperature on the morphology and the efficiency
We then investigated how the annealing temperature affects the morphology, PCE and the domain size of SM OSC (see Fig. 3). In the DPD simulation, we set the reduced temperatures (T*) to the corresponding physical temperature. A rough but simple estimation of the value of this reduced temperature in terms of physical temperature can be given by the equation61,62: T = 133T∗ − 33 to simplify the relationship between T and T*. The reduced temperatures T* = 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, and 1.6 correspond to physical temperature T = −6.4, 20.2, 46.8, 73.4, 100, 126.6, 153.2, and 179.8 °C.
![image file: c5ra12318f-f3.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-f3.gif) |
| Fig. 3 (a) The dependence of Power Conversion Efficiency (PCE) on the temperature of PC71BM for C1/PC71BM; (b) the dependence of Power Conversion Efficiency (PCE) on the temperature for C2/PC71BM; (c) the dependence of Power Conversion Efficiency (PCE) on the temperature for C3/PC71BM; (d) the dependence of the domain size on the temperature for C1/PC71BM; (e) the dependence of the domain size on the temperature for C2/PC71BM; (f) the dependence of the domain size on the temperature for C3/PC71BM; the molar percentage of PC71BM is 66.7%. | |
Fig. 3 shows that the morphology and PCE are the best when the reduced temperatures is 1.0, namely, the physical temperature is 100 °C.63 Similarly, temperature has been found to be an important parameter for morphology and efficiency, experimentally.64,65 Upon changing the heating temperature, the extent of phase separation evolves and finally forms phase-separated network-like morphology, which benefits the efficiency of the device.
4.3 The effect of the additives on the morphology and the efficiency
Kenneth R. Graham et al.63 investigated the effect of macromolecular additive, polydimethylsiloxane (PDMS) on the performance of solution processed molecular bulk heterojunction solar cells. Their results show that PDMS has a strong influence on film morphology, with a significant decrease in film roughness and feature size observed. As a general, the additive has a higher solubility for PC71BM than the donor material.32 This increases donor aggregation and leads to the formation of high ordered D and A phases, thus improving device performance.
Here we investigated how the Degree of Polymerization (DP) of PDMS affects the morphology and PCE of OSC. As shown in Fig. 4, the introduction of PDMS will improve the PCE. For C1, C2, and C3, PCE will increase, when there is PDMS introduced. In addition, our results show that, when DP of PDMS is equal to 20, PCE will be slightly better than other DP values (e.g. DP = 4, or 100). This is consistent with Kenneth R. Graham's66 results. Graham et al. showed that the PCE is enhanced from 1.42% without additive to 2.35% with 410 g mol−1 PDMS, which corresponds to DP around 4.4.
![image file: c5ra12318f-f4.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-f4.gif) |
| Fig. 4 The dependence of the power conversion efficiency (PCE) on the degree of polymerization (DP) of PDMS as additive in the active layer formed by C1, C2, or C3 with PC71BM. We highlight DP = 4, 20, or 100. When DP = 20, it leads to the best PCE. | |
As summarized in Table 5, we can see that the domain size will increase as PCE, when there is PDMS introduced.29 Our results show that, when there is PDMS introduced, f will increases, which corresponds to the short current. Our f values are consistent with the short current from experiments.29 Our approach is useful to illustrate the effect of additives and will be useful for rational design of additives.
Table 5 The effect of the additive on the PCE, domain size, short current, and morphology descriptors
Donor: PC71BM |
C1 |
C2 |
C3 |
Without additive |
With additive |
Without additive |
With additive |
Without additive |
With additive |
PCE (%) |
Exptl.29 |
6.07 |
6.50 |
7.22 |
7.80 |
11.66 |
11.73 |
Modeling |
6.95 |
10.30 |
10.19 |
10.49 |
10.04 |
10.35 |
Domain size (nm) |
Exptl.29 |
15–35 |
20–45 |
10–30 |
15–40 |
8–30 |
15–35 |
Modeling |
5.76 |
6.02 |
6.26 |
6.24 |
5.79 |
5.83 |
J (mA cm−2) |
Exptl.29 |
12.16 |
11.75 |
12.92 |
12.95 |
11.66 |
11.73 |
f (%) |
Modeling |
52.8 |
53.9 |
50.2 |
51.8 |
55.0 |
56.4 |
In addition, we studied the desolvation process based on DPD simulation and graph theory. Fig. 5 shows the change of PCE and morphology during desolvation process. Based on our results, we can see that the aggregation level of PC71BM phase is becoming large with the evaporation of the solvent and percolating network of electron-donor and electron-acceptor materials is formed. And the PCE increases in the desolvation process. Our approach is able to provide the dynamical process during desolvation and will be useful for optimizing the desolvation procedure in the fabrication.
![image file: c5ra12318f-f5.tif](/image/article/2015/RA/c5ra12318f/c5ra12318f-f5.gif) |
| Fig. 5 PCE and morphology of C1/PC71BM evaluated from DPD simulation and graph theory during desolvation process. XS/A represents the percentage of the solvent. | |
5 Conclusions
In this work, we use dissipative particle dynamics simulation and graph theory to predict the morphology and efficiency of small molecular organic solar cell. Our simulations explore the effects of the structural parameters on the morphology. We compare our results with recent experiments and provide a quantitative and qualitative analysis for the experiments studies.
Certainly the Small molecule as the donor is the good choice for organic solar cells. In summary, we combine atomistic simulation with mesoscale simulation to predict the morphology of BHJ solar cells and the optimal parameters for morphology of BHJ solar cells. For the effect of the blend ratio, we find that the molar ratio of PC71BM 66.7% leads to the best performance. In addition, we study the effect of temperature and the additives on the morphology and the performance of small molecular organic solar cell. Our approach is useful to illustrate the mechanism and the important factors for organic solar cell. And it could be used to optimize other important factors to improve the performance of organic solar cell.
Acknowledgements
The work is supported by the National Basic Research Program of China (973 Program, Grant No. 2012CB932400), the National Natural Science Foundation of China (Grant No. 91233115, 21273158, and 91227201), a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). This is also a project supported by the Fund for Innovative Research Teams of Jiangsu Higher Education Institutions, Jiangsu Key Laboratory for Carbon-Based Functional Materials and Devices, Collaborative Innovation Center of Suzhou Nano Science and Technology.
References
- Y.-J. Cheng, S.-H. Yang and C.-S. Hsu, Synthesis of conjugated polymers for organic solar cell applications, Chem. Rev., 2009, 109, 5868–5923 CrossRef CAS PubMed.
- G. Li, R. Zhu and Y. Yang, Polymer solar cells, Nat. Photonics, 2012, 6, 153–161 CrossRef CAS PubMed.
- A. Mishra and P. Bäuerle, Small molecule organic semiconductors on the move: promises for future solar energy technology, Angew. Chem., Int. Ed., 2012, 51, 2020–2067 CrossRef CAS PubMed.
- R. Po, M. Maggini and N. Camaioni, Polymer solar cells: recent approaches and achievements, J. Phys. Chem. C, 2009, 114, 695–706 Search PubMed.
- H. Bronstein, E. Collado-Fregoso, A. Hadipour, Y. W. Soon, Z. Huang, S. D. Dimitrov, R. S. Ashraf, B. P. Rand, S. E. Watkins and P. S. Tuladhar, Thieno [3, 2-b] thiophene-diketopyrrolopyrrole Containing Polymers for Inverted Solar Cells Devices with High Short Circuit Currents, Adv. Funct. Mater., 2013, 23, 5647–5654 CrossRef CAS PubMed.
- J. J. Intemann, K. Yao, Y. X. Li, H. L. Yip, Y. X. Xu, P. W. Liang, C. C. Chueh, F. Z. Ding, X. Yang and X. Li, Highly Efficient Inverted Organic Solar Cells Through Material and Interfacial Engineering of Indacenodithieno [3, 2-b] thiophene-Based Polymers and Devices, Adv. Funct. Mater., 2014, 24, 1465–1473 CrossRef CAS PubMed.
- W. Li, K. H. Hendriks, A. Furlan, W. Roelofs, S. C. Meskers, M. M. Wienk and R. A. Janssen, Effect of the Fibrillar Microstructure on the Efficiency of High Molecular Weight Diketopyrrolopyrrole-Based Polymer Solar Cells, Adv. Mater., 2014, 26, 1565–1570 CrossRef CAS PubMed.
- C. M. Cabanetos, A. El Labban, J. A. Bartelt, J. D. Douglas, W. R. Mateker, J. M. Fréchet, M. D. McGehee and P. M. Beaujuge, Linear side chains in benzo [1, 2-b: 4, 5-b′] dithiophene–thieno [3, 4-c] pyrrole-4, 6-dione polymers direct self-assembly and solar cell performance, J. Am. Chem. Soc., 2013, 135, 4656–4659 CrossRef CAS PubMed.
- K. Li, Z. Li, K. Feng, X. Xu, L. Wang and Q. Peng, Development of large band-gap conjugated copolymers for efficient regular single and tandem organic solar cells, J. Am. Chem. Soc., 2013, 135, 13549–13557 CrossRef CAS PubMed.
- I. Meager, R. S. Ashraf, S. Mollinger, B. C. Schroeder, H. Bronstein, D. Beatrup, M. S. Vezie, T. Kirchartz, A. Salleo and J. Nelson, Photocurrent Enhancement from diketopyrrolopyrrole polymer solar cells through alkyl-chain branching point manipulation, J. Am. Chem. Soc., 2013, 135, 11537–11540 CrossRef CAS PubMed.
- H. Cha, D. S. Chung, S. Y. Bae, M. J. Lee, T. K. An, J. Hwang, K. H. Kim, Y. H. Kim, D. H. Choi and C. E. Park, Complementary Absorbing Star-Shaped Small Molecules for the Preparation of Ternary Cascade Energy Structures in Organic Photovoltaic Cells, Adv. Funct. Mater., 2013, 23, 1556–1565 CrossRef CAS PubMed.
- Y.-H. Chen, L.-Y. Lin, C.-W. Lu, F. Lin, Z.-Y. Huang, H.-W. Lin, P.-H. Wang, Y.-H. Liu, K.-T. Wong and J. Wen, Vacuum-deposited small-molecule organic solar cells with high power conversion efficiencies by judicious molecular design and device optimization, J. Am. Chem. Soc., 2012, 134, 13616–13623 CrossRef CAS PubMed.
- A. Leliège, J. Grolleau, M. Allain, P. Blanchard, D. Demeter, T. Rousseau and J. Roncali, Small D–π–A Systems with o-Phenylene-Bridged Accepting Units as Active Materials for Organic Photovoltaics, Chem.–Eur. J., 2013, 19, 9948–9960 CrossRef PubMed.
- P. Li, H. Tong, J. Ding, Z. Xie and L. Wang, Small molecules based
on 2, 7-carbazole for efficient solution-processed organic solar cells, J. Mater. Chem. A, 2013, 1, 8805–8812 CAS.
- J. Liu, Y. Sun, P. Moonsin, M. Kuik, C. M. Proctor, J. Lin, B. B. Hsu, V. Promarak, A. J. Heeger and T. Q. Nguyen, Tri-Diketopyrrolopyrrole Molecular Donor Materials for High-Performance Solution-Processed Bulk Heterojunction Solar Cells, Adv. Mater., 2013, 25, 5898–5903 CrossRef CAS PubMed.
- X. Liu, Y. Sun, L. A. Perez, W. Wen, M. F. Toney, A. J. Heeger and G. C. Bazan, Narrow-band-gap conjugated chromophores with extended molecular lengths, J. Am. Chem. Soc., 2012, 134, 20609–20612 CrossRef CAS PubMed.
- Y. Liu, C.-C. Chen, Z. Hong, J. Gao, Y. M. Yang, H. Zhou, L. Dou, G. Li and Y. Yang, Solution-processed small-molecule solar cells: breaking the 10% power conversion efficiency, Sci. Rep., 2013, 3, 3356 Search PubMed.
- D. Patra, C.-C. Chiang, W.-A. Chen, K.-H. Wei, M.-C. Wu and C.-W. Chu, Solution-processed benzotrithiophene-based donor molecules for efficient bulk heterojunction solar cells, J. Mater. Chem. A, 2013, 1, 7767–7774 CAS.
- B. Walker, A. B. Tamayo, X. D. Dang, P. Zalar, J. H. Seo, A. Garcia, M. Tantiwiwat and T. Q. Nguyen, Nanoscale Phase Separation and High Photovoltaic Efficiency in Solution-Processed, Small-Molecule Bulk Heterojunction Solar Cells, Adv. Funct. Mater., 2009, 19, 3063–3069 CrossRef CAS PubMed.
- H. Wang, F. Liu, L. Bu, J. Gao, C. Wang, W. Wei and T. P. Russell, The Role of Additive in Diketopyrrolopyrrole-Based Small Molecular Bulk Heterojunction Solar Cells, Adv. Mater., 2013, 25, 6519–6525 CrossRef CAS PubMed.
- K. N. Winzenberg, P. Kemppinen, F. H. Scholes, G. E. Collis, Y. Shu, T. B. Singh, A. Bilic, C. M. Forsyth and S. E. Watkins, Indan-1, 3-dione electron-acceptor small molecules for solution-processable solar cells: a structure–property correlation, Chem. Commun., 2013, 49, 6307–6309 RSC.
- D. Ye, X. Li, L. Yan, W. Zhang, Z. Hu, Y. Liang, J. Fang, W.-Y. Wong and X. Wang, Dithienosilole-bridged small molecules with different alkyl group substituents for organic solar cells exhibiting high open-circuit voltage, J. Mater. Chem. A, 2013, 1, 7622–7629 CAS.
- J. Zhou, X. Wan, Y. Liu, Y. Zuo, Z. Li, G. He, G. Long, W. Ni, C. Li and X. Su, Small molecules based on benzo [1, 2-b: 4, 5-b′] dithiophene unit for high-performance solution-processed organic solar cells, J. Am. Chem. Soc., 2012, 134, 16345–16351 CrossRef CAS PubMed.
- J. Zhou, Y. Zuo, X. Wan, G. Long, Q. Zhang, W. Ni, Y. Liu, Z. Li, G. He and C. Li, Solution-processed and high-performance organic solar cells using small molecules with a benzodithiophene unit, J. Am. Chem. Soc., 2013, 135, 8484–8487 CrossRef CAS PubMed.
- B. Kan, M. M. Li, Q. Zhang, F. Liu, X. J. Wan, Y. C. Wang, W. Ni, G. K. Long, X. Yang, H. R. Feng, Y. Zuo, M. T. Zhang, F. Huang, Y. Cao, T. P. Russell and Y. S. Chen, A Series of Simple Oligomer-like Small Molecules Based on Oligothiophenes for Solution-Processed Solar Cells with High Efficiency, J. Am. Chem. Soc., 2015, 137, 3886–3893 CrossRef CAS PubMed.
- B. G. Sumpter and V. Meunier, Can computational approaches aid in untangling the inherent complexity
of practical organic photovoltaic systems?, J. Polym. Sci., Part B: Polym. Phys., 2012, 50, 1071–1089 CrossRef CAS PubMed.
- F. Liu, Y. Gu, J. W. Jung, W. H. Jo and T. P. Russell, On the morphology of polymer-based photovoltaics, J. Polym. Sci., Part B: Polym. Phys., 2012, 50, 1018–1044 CrossRef CAS PubMed.
- M. A. Ruderer and P. Müller-Buschbaum, Morphology of polymer-based bulk heterojunction films for organic photovoltaics, Soft Matter, 2011, 7, 5482–5493 RSC.
- J. Zhou, Y. Zuo, X. Wan, G. Long, Q. Zhang, W. Ni, Y. Liu, Z. Li, G. He, C. Li, B. Kan, M. Li and Y. Chen, Solution-processed and high-performance organic solar cells using small molecules with a benzodithiophene unit, J. Am. Chem. Soc., 2013, 135, 8484–8487 CrossRef CAS PubMed.
- Y. Liu, C. C. Chen, Z. Hong, J. Gao, Y. M. Yang, H. Zhou, L. Dou, G. Li and Y. Yang, Solution-processed small-molecule solar cells: breaking the 10% power conversion efficiency, Sci. Rep., 2013, 3, 3356 Search PubMed.
- J. Zhou, X. Wan, Y. Liu, Y. Zuo, Z. Li, G. He, G. Long, W. Ni, C. Li, X. Su and Y. Chen, Small molecules based on benzo[1,2-b:4,5-b']dithiophene unit for high-performance solution-processed organic solar cells, J. Am. Chem. Soc., 2012, 134, 16345–16351 CrossRef CAS PubMed.
- J. K. Lee, W. L. Ma, C. J. Brabec, J. Yuen, J. S. Moon, J. Y. Kim, K. Lee, G. C. Bazan and A. J. Heeger, Processing additives for improved efficiency from bulk heterojunction solar cells, J. Am. Chem. Soc., 2008, 130, 3619–3623 CrossRef CAS PubMed.
- J. K. Park, B. Walker and J. H. Seo, Controlling polarity of organic bulk heterojunction field-effect transistors via solvent additives, ACS Appl. Mater. Interfaces, 2013, 5, 4575–4580 CAS.
- P. Hoogerbrugge and J. Koelman, Simulating microscopic hydrodynamic phenomena with dissipative particle dynamics, EPL, 1992, 19, 155 CrossRef.
- J. Koelman and P. Hoogerbrugge, Dynamic simulations of hard-sphere suspensions under steady shear, EPL, 1993, 21, 363 CrossRef CAS.
- R. D. Groot and T. J. Madden, Dynamic simulation of diblock copolymer microphase separation, J. Chem. Phys., 1998, 108, 8713 CrossRef CAS PubMed.
- Y. Li, X. He, X. Cao, Y. Shao, Z. Li and F. Dong, Mesoscopic simulation study on the efficiency of surfactants adsorbed at the liquid/liquid interface, Mol. Simul., 2005, 31, 1027–1033 CrossRef CAS PubMed.
- O. Wodo, S. Tirthapura, S. Chaudhary and B. Ganapathysubramanian, A graph-based formulation for computational characterization of bulk heterojunction morphology, Org. Electron., 2012, 13, 1105–1113 CrossRef CAS PubMed.
- B. Lei, Y. Yao, A. Kumar, Y. Yang and V. Ozolins, Quantifying the relation between the morphology and performance of polymer solar cells using Monte Carlo simulations, J. Appl. Phys., 2008, 104, 024504 CrossRef PubMed.
- G. A. Buxton and N. Clarke, Predicting structure and property relations in polymeric photovoltaic devices, Phys. Rev. B: Condens. Matter Mater. Phys., 2006, 74, 085207 CrossRef.
- Y. Min Nam, J. Huh and W. Ho Jo, Optimization of thickness and morphology of active layer for high performance of bulk-heterojunction organic solar cells, Sol. Energy Mater. Sol. Cells, 2010, 94, 1118–1124 CrossRef PubMed.
- P. K. Watkins, A. B. Walker and G. L. B. Verschoor, Dynamical Monte Carlo Modelling of Organic Solar Cells:
The Dependence of Internal Quantum Efficiency on Morphology, Nano Lett., 2005, 5, 1814–1818 CrossRef CAS PubMed. - B. C. Lin, B. T. Koo, P. Clancy and C. P. Hsu, Theoretical Investigation of Charge-Transfer Processes at Pentacene-C-60 Interface: The Importance of Triplet Charge Separation and Marcus Electron Transfer Theory, J. Phys. Chem. C, 2014, 118, 23605–23613 CAS.
- H. Sun, COMPASS: An ab initio force-field optimized for condensed-phase applications - Overview with details on alkane and benzene compounds, J. Phys. Chem. B, 1998, 102, 7338–7364 CrossRef CAS.
- D. Rigby, H. Sun and B. E. Eichinger, Computer simulations of poly(ethylene oxide): force field, pvt diagram and cyclization behaviour, Polymer. Int., 1997, 44, 311–330 CrossRef CAS.
- J. Xue, T. Hou and Y. Li, Optimal parameters for morphology of bulk heterojunction solar cells from simulations, Appl. Phys. Lett., 2012, 100, 053307 CrossRef PubMed.
- D. R. Kozub, K. Vakhshouri, L. M. Orme, C. Wang, A. Hexemer and E. D. Gomez, Polymer Crystallization of Partially Miscible Polythiophene/Fullerene Mixtures Controls Morphology, Macromolecules, 2011, 44, 5722–5726 CrossRef CAS.
- F. Machui, S. Langner, X. Zhu, S. Abbott and C. J. Brabec, Determination of the P3HT: PCBM solubility parameters via a binary solvent gradient method: Impact of solubility on the photovoltaic performance, Sol. Energy Mater. Sol. Cells, 2012, 100, 138–146 CrossRef CAS PubMed.
- R. Groot and K. Rabone, Mesoscopic simulation of cell membrane damage, morphology change and rupture by nonionic surfactants, Biophys. J., 2001, 81, 725–736 CrossRef CAS.
- N. S. Martys and R. D. Mountain, Velocity Verlet algorithm for dissipative-particle-dynamics-based models of suspensions, Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top., 1999, 59, 3733 CrossRef CAS.
- J. P. Flory, Principles of polymer chemistry, Cornell University Press, 1953 Search PubMed.
- A. Maiti and S. McGrother, Bead–bead interaction parameters in dissipative particle dynamics: relation to bead-size, solubility parameter, and surface tension, J. Chem. Phys., 2004, 120, 1594 CrossRef CAS PubMed.
- Y. Zhao, L.-Y. You, Z.-Y. Lu and C.-C. Sun, Dissipative particle dynamics study on the multicompartment micelles self-assembled from the mixture of diblock copolymer poly (ethyl ethylene)-< i> block</i>-poly (ethylene oxide) and homopolymer poly (propylene oxide) in aqueous solution, Polymer, 2009, 50, 5333–5340 CrossRef CAS PubMed.
- H. Sun, COMPASS: An ab initio force-field optimized for condensed-phase applications overview with details on alkane and benzene compounds, J. Phys. Chem. B, 1998, 102, 7338–7364 CrossRef CAS.
- H. Sun, P. Ren and J. Fried, The COMPASS force field: parameterization and validation for phosphazenes, Comput. Theor. Polym. Sci., 1998, 8, 229–246 CrossRef CAS.
- J. Yang, Y. Ren, A.-M. Tian and H. Sun, COMPASS force field for 14 inorganic molecules, He, Ne, Ar, Kr, Xe, H2, O2, N2, NO, CO, CO2, NO2, CS2, and SO2, in liquid phases, J. Phys. Chem. B, 2000, 104, 4951–4957 CrossRef CAS.
- P. Petrus, M. Lisal and J. K. Brennan, Self-assembly of symmetric diblock copolymers in planar slits with and without nanopatterns: Insight from dissipative particle dynamics simulations, Langmuir, 2009, 26, 3695–3709 CrossRef PubMed.
- Z.-W. Li, Z.-Y. Lu, Z.-Y. Sun, Z.-S. Li and L.-J. An, Calculating the equation of state parameters and predicting the spinodal curve of isotactic polypropylene/poly (ethylene-co-octene) blend by molecular dynamics simulations combined with Sanchez-Lacombe lattice fluid theory, J. Phys. Chem. B, 2007, 111, 5934–5940 CrossRef CAS PubMed.
- Y. Sun, G. C. Welch, W. L. Leong, C. J. Takacs, G. C. Bazan and A. J. Heeger, Solution-processed small-molecule solar cells with 6.7% efficiency, Nat. Mater., 2012, 11, 44–48 CrossRef CAS PubMed.
- A. Ko Ko Kyaw, D. Hwan Wang, H.-R. Tseng, J. Zhang, G. C. Bazan and A. J. Heeger, Electron and hole mobility in solution-processed small molecule-fullerene blend: Dependence on the fullerene content, Appl. Phys. Lett., 2013, 102, 163308 CrossRef PubMed.
- M. Venturoli, B. Smit and M. M. Sperotto, Simulation Studies of Protein-Induced Bilayer Deformations, and Lipid-Induced Protein Tilting, on a Mesoscopic Model for Lipid Bilayers with Embedded Proteins, Biophys. J., 2005, 88, 1778–1798 CrossRef CAS PubMed.
- X. Li, I. V. Pivkin, H. Liang and G. E. Karniadakis, Shape Transformations of Membrane Vesicles from Amphiphilic Triblock Copolymers: A Dissipative Particle Dynamics Simulation Study, Macromolecules, 2009, 42, 3195–3200 CrossRef CAS.
- K. R. Graham, J. Mei, R. Stalder, J. W. Shim, H. Cheun, F. Steffy, F. So, B. Kippelen and J. R. Reynolds, Polydimethylsiloxane as a macromolecular additive for enhanced performance of molecular bulk heterojunction organic solar cells, ACS Appl. Mater. Interfaces, 2011, 3, 1210–1215 CAS.
- W. Ma, C. Yang, X. Gong, K. Lee and A. J. Heeger, Thermally Stable, Efficient Polymer Solar Cells with Nanoscale Control of the Interpenetrating Network Morphology, Adv. Funct. Mater., 2005, 15, 1617–1622 CrossRef CAS PubMed.
- E. Verploegen, R. Mondal, C. J. Bettinger, S. Sok, M. F. Toney and Z. Bao, Effects of Thermal Annealing Upon the Morphology of Polymer–Fullerene Blends, Adv. Funct. Mater., 2010, 20, 3519–3529 CrossRef CAS PubMed.
- K. R. Graham, P. M. Wieruszewski, R. Stalder, M. J. Hartel, J. Mei, F. So and J. R. Reynolds, Improved Performance of Molecular Bulk-Heterojunction Photovoltaic Cells through Predictable Selection of Solvent Additives, Adv. Funct. Mater., 2012, 22, 4801–4813 CrossRef CAS PubMed.
|
This journal is © The Royal Society of Chemistry 2015 |