Planarity and multiple components promote organic photovoltaic efficiency by improving electronic transport

Establishing how the conformation of organic photovoltaic (OPV) polymers affects their electronic and transport properties is critical in order to determine design rules for new OPV materials and in particular to understand the performance enhancements recently reported for ternary blends. We report coupled classical and ab initio molecular dynamics simulations showing that polymer linkage twisting significantly reduces optical absorption efficiency, as well as hole transport rates in donor polymers. We predict that blends with components favoring planar geometries contribute to the enhancement of the overall efficiency of ternary OPVs. Furthermore, our electronic structure calculations for the PTB7-PID2-PC71BM system show that hole transfer rates are enhanced in ternary blends with respect to their binary counterpart. Finally, our results point at thermal disorder in the blend as a key reason responsible for device voltage losses and at the need to carry out electronic structure calculations at finite temperature to reliably compare with experiments.


Introduction
2][3] A relative abundance of semiconducting polymers, the ease of solution-based processing, and the relatively high optical absorption coefficients of OPV materials are all promising properties for grid-scale solar energy harvesting, together with the low production cost. 2,4,5Typical OPV cells consist of a p-conjugated polymeric donor material, commonly poly-3-hexyl-thiophene (P3HT), along with a fullerene acceptor functionalized to enhance solubility, such as phenyl-C61-butyric acid methyl ester (PC61BM) or phenyl-C71-butyric acid methyl ester (PC71BM).
Solar cells constructed from these materials are highly nonuniform by design.They contain a bulk heterojunction (BHJ) in order to maximize interfacial contact between components, which is necessary to promote exciton dissociation and hence harvest electrons and holes; the latter are transported through the acceptor (fullerene) and donor (polymer) regions, respectively.
][17][18][19][20][21] One traditional approach to organic solar cell optimization has involved tuning the HOMO energy level and optical gaps of donor materials, in addition to thin film morphology optimization.This approach has led to developing new donor materials consisting of copolymers made from two different building blocks, 22 which resulted in added flexibility in tuning band positions and exciton binding energies, and thus hole and electron separation. 23Recently, more complex systems have been proposed, consisting of ternary blends of two donor materials (polythieno [3,4-b]-thiophene/benzodithiophene (PTB7), and poly-3-oxothieno [3,4-d]isothiazole-1,1-dioxide/benzodithiophene (PID2)) and one fullerene acceptor (PC71BM), resulting in photoconversion efficiencies which exceeded those of binary solar cells.These improvements were hypothesized to come from changes in charge separation and transport with respect to binary blends and from changes in the blend morphology. 24n order to provide design rules for new OPV materials, as well as understand the performance enhancement reported for ternary blends, it is important to establish how the conformation of the OPV polymers affects the electronic properties of the blend, particularly level alignments, and how the overall morphology of the blend affects hole mobilities in donor regions.
The rational design of OPV materials using computational methods is a challenging task, due to the size and morphology of the blends composing OPV devices and the numerous interdependent properties affecting the photoconversion efficiency.Typical molecular weights of OPV polymers suggest polymer sizes of 10-1000 repeating units, composed of 1000-100 000 atoms, which extend way beyond the regime of accessibility of first principles calculations.6][27] Although useful starting points, isolated oligomers, often represented by their zero temperature optimized geometry, do not account for the inherent disorder present in OPV blends, 28 or for the presence of polydispersity, and of many different stereoisomers resulting from either variation in polymeric linkages or side-chain conformations. 29Changes in polymers (including side-chains, bridging units, and atomic substitutions) may result in surprisingly varied morphologies, electronic properties, and ultimate photoconversion efficiency. 30][33][34][35][36][37][38] In this work, we used classical molecular dynamics (MD) to obtain realistic disordered structures of ternary polymer blends, whose electronic and transport properties were then investigated using first principles methods.In particular, we considered the ternary blend of PTB7, PID2, and PC71BM, recently proposed experimentally 24 and shown in Fig. 1, and we determined the relative energy level alignment of the three components and the transport properties of the composite system, using Marcus theory to describe the hopping transport [39][40][41] of holes.
Our focus on hopping transport builds on the work of Bre ´das et al., 42 Coropceanu, et al., 43 and Ba ¨ssler and Ko ¨hler, 44 who proposed that hole transport occurs primarily via hopping in disordered conjugated polymers, after charge transfer to fullerene acceptors has occurred. 6Measured current versus time curves, obtained by time of flight experiments, are consistent with the proposed picture of transport occurring primarily by hopping. 45][48] Our coupled classical and ab initio molecular dynamics simulations showed that polymer linkage twisting significantly reduces optical absorption efficiency, as well as hole transport rates in donor polymers compared to those of polymer blends with planar, ordered geometries.Hence, we predicted that blends with components favoring planar geometries contribute to the enhancement of the overall efficiency of ternary OPVs.Furthermore, our electronic structure calculations showed that hole transfer rates from PC71BM -PID2 are much larger than those from PC71BM -PTB7, contributing to an improved charge current density in ternary blends with respect to their binary counterpart.
In our analysis of the ternary blend solar cell, we first characterized the equilibrium geometries of PTB7 and PID2 using density functional theory.Computational details are given in Section 2. We started by investigating the twisting of links between subunits in an oligomer picture, followed by the analysis of ground state periodic configurations, as described in Section 3.1.As validation of our methodology and optimized structures, we compared the computed energy levels of the lowest lying configurations to the experimental measurements in Section 3.2.Next, we used molecular dynamics to sample the effects of structural disorder on the properties of PC71BM and PTB7 in Section 3.3.Then, we used classical molecular dynamics simulations to assess the effect of PID2 on the morphology of PTB7 in Section 3.4.Using umbrella sampling 49,50 combined with the weighted histogram analysis method, 51 we find that the presence of PID2 forces PTB7 to adopt a more planar geometry, thereby increasing charge transfer rates.Finally, we computed the charge transport rates at interfaces between donor polymers and PC71BM and these are discussed in Section 3.5.

Computational methods
The characterization of isolated monomers of PTB7 and PID2 was performed within the Gaussian basis set code ORCA 3 52 using the def2-TZVP 53 basis set for density functional calculations (Perdew-Burke-Ernzerhof 54 (PBE), the hybrid PBE0, 55 and the hybrid B3LYP [56][57][58] ) and the Dunning basis cc-pVTZ 59 for resolution of the identity second order Møller-Plesser perturbation theory 60 (MP2) calculations.Unless otherwise noted, all density functional approximations used for Gaussian basis set calculations included the Grimme dispersion correction (-D3) 61 with Becke-Johnson (BJ) damping 62,63 and the geometric This journal is © the Owner Societies 2016 counterpoise correction (gCP) of Kruse, et al. 64 Gaussian basis set calculations used the chain-of-spheres approximation 65,66 and density fitting for Coulomb integrals. 67,68or periodic systems, plane wave density functional theory (DFT) calculations were performed using Quantum Espresso, 69 with Troullier-Martins and Hartwigsen-Goedecker-Hutter pseudopotentials 70,71 and energy cutoffs of 85 and 340 Ry for the wavefunction and charge density, respectively.First principles molecular dynamics simulations were performed with the PBE density functional approximation using Qbox 72 and normconserving pseudopotentials of the HSCV type 73,74 and a cutoff of 85 Ry. van der Waals complexes of PTB7 and PID2 with PC71BM were optimized using the rVV10 van der Waals density functional. 75,76arious conformations of the polymers representing structural disorder were obtained from constant volume and temperature (NVT) classical molecular dynamics (MD) simulations performed with the LAMMPS package 77 and a modified OPLS-AA force field. 78,79A Nose-Hoover 80,81 thermostat was used with a time step of 1 femtosecond and a temperature of 300 K.
To examine the structural effects of PID2 on PTB7, we performed bulk simulations of PTB7 and PID2, which were prepared by initializing systems at low density, then equilibrating at 1 atm and 500 K for 1 ns until condensed.Systems were then equilibrated at 500 K for 1 ns.Sampling was performed over 2.5 ns at 500 K and 1 atm.We find that polymer chains are within the diffusive regime in these timescales, as shown in ESI.† We used PTB7 and PID2 chains with 5 monomer units, allowing for rapid relaxation.Umbrella restraints were applied to the appropriate PTB7 dihedral during the entire simulation.We performed umbrella sampling using 50 evenly space umbrellas with a restraint strength of 200 kcal mol À1 .
Given that the major differences between OPV binary and ternary blends come from variation in the donor material composition, we focused our investigation on hole transport, hypothesized to primarily occur in the donor region of OPVs 8 by hopping.Charge transfer parameters for hopping transport in semiconducting systems are commonly computed using either constrained DFT (CDFT), [82][83][84] excited state wavefunction methods, 85 or model Hamiltonians, [86][87][88] and then subsequently used to estimate transfer rate using, e.g.Marcus theory.
In our calculations of transport parameters, we assumed that holes were localized either on a copolymer unit (consisting of both a thienothiophene and a benzodithiophene building block) or a PC71BM molecule.Indeed, ab initio calculations of hole wavefunctions showed that holes do localize within one repeating unit of copolymers 45 or within four monomer units of the much smaller homopolymer P3HT, 89 where a P3HT tetramer is approximately the same size as one monomer of PTB7 or PID2.Localization of charge carriers in organic polymer semiconductors is further supported by terahertz spectroscopy experiments. 90Previous studies using drift-diffusion models hypothesized that current in these materials is largely determined by the variability in 1-dimensional conductivity along polymer chains; 91 to that end, we calculated charge transport parameters for extended chain configurations.
We computed hole transfer rates hole transfer rates (k AB ) at the interfaces of the donor polymers PTB7 and PID2 with PC71BM using Marcus theory: where A and B denote either a polymer unit or PC71BM, respectively; H AB is the nonadiabatic coupling and DG AB the driving force for the charge transfer process; l AB is the reorganization energy.Other rate equations, similar to eqn (1), are often used in the limit of large (H AB Z l AB /2) nonadiabatic couplings. 92,93Quantum effects from zero-point energies (typically on the order of 0.05-0.075eV for organic semiconductors 93 ) do not substantially modify calculated rates in the case of sufficiently large activation energies (as compared to zero point energies).In our case, activation energies for hole transfer to/from the PC71BM are on the order of 0.2-1.1 eV. 94Hence, we neglected quantum corrections to hole transfer rates in this work.Nonadiabatic couplings and driving forces were calculated using CDFT, following the implementation described by Oberhofer and Blumberger; 95 such implementation was carried out in a developmental version of Quantum Espresso and will be described in a future publication. 96In our implementation of CDFT, differences in integrated charge density (using Hirshfeld partitioning 97 ) were localized in two regions defined by atoms donating or accepting an electron.We note that, as previously suggested in a CDFT study of periodic systems, 98 nonadiabatic couplings H AB between sites are reduced by a factor of two when two equivalent sites B and B 0 (B 0 being the periodic image of B) are equidistant from site A.
In order to understand how conformation changes the electronic properties of OPV polymers, we computed the energy levels and hole transport couplings of 200 thermally disordered configurations of extended polymers, which were generated by classical MD.To determine specific structure-property relationships, we performed a clustering analysis of dihedral angles in these configurations; we identified distinct structural motifs solely from the dihedral angles between polymer subunits, and we computed averaged properties per cluster.First, we used principal component analysis 99 to perform a dimensionality reduction to two principal components.We then identified clusters (i.e.structural motifs) in this reduced space with the affinity propagation method, 100 relying upon a Euclidean distance metric.Our identification of clusters was performed using the scikit-learn python library. 101

Results and discussion
In the following, we present our results for structural optimizations and electronic structure calculations.As mentioned in the introduction, we first computed optimized geometries of monomers and periodic conformers using density functional theory, and we compared computed energy levels with experimental electron affinities and ionization potentials.Then, we established the effect of structural disorder on electronic properties using configurations generated by classical molecular dynamics simulations, focusing on hole transport parameters.Finally, we calculated charge-transfer rates for specific interfacial geometries to address the increased charge current density observed experimentally in ternary blends, 24 and we found that that hole transfers from PC71BM -PID2 proceed orders of magnitude faster than from PC71BM -PTB7; these results indicate that ternary blends will have significantly improved charge carrier collection resulting from photoexcitation of PC71BM.Our results are discussed in detail below.

Conformer identification
In order to build realistic model geometries of OPV polymers, we first investigated the total energy of polymer twisting as a function of the dihedral angle between copolymer subunits.3][104] As noted by Bhatta et al., 102 the dependence of the energy of these polymers on structural parameters is well-represented using only one-dimensional dihedral scans since the inter-mode vibrational couplings for dihedral fluctuations are small (0.11 kcal mol À1 ), on the order of the errors expected when using finite basis density functional approximations for conformational energies.In Fig. 2, we show potential energy surfaces as a function of dihedrals (denoted in Fig. 1) for PBE-D3(BJ)-gCP/ def2-TZVP geometries of PTB7 and PID2 constrained at each dihedral value.To remove the effect of side-chains, the 2-ethylhexyl groups were replaced by hydrogens.The MP2, PBE, and B3LYP surfaces present minima in correspondence to approximately aligned sulfurs (hereafter denoted as a cis geometry) or anti-aligned (trans geometry), with an offset of 10-30 degrees for PTB7 and 0-20 degrees for PID2.This preference for nearplanarity is consistent with the conjugated, or quinoidal, nature of these polymers along the backbone.The barriers to twisting the linkage between the subunits are appreciably higher in PID2 (B6 kcal mol À1 in Fig. 2, top) than in PTB7 (B3 kcal mol À1 in Fig. 2, bottom), and these barriers, as well as the relative positions of minima, are within chemical accuracy for the three methods used.In addition to a higher barrier, the energy surface of PID2 presents local minima with a higher energetic separation of 3-4 kcal mol À1 .From these results, we expect the conformations accessed by PID2 at finite temperature to be significantly less diverse than those of PTB7.
From the configurations with minimal energy determined above, we constructed periodic configurations for use in planewave density functional theory (DFT) models in order to determine the energy levels of extended chains.Since each dihedral potential energy surface has two minima, we examined a total of four conformers for extended chains.Each of the four conformers for PTB7 and PID2 were optimized using the PBE density functional approximation with 5 k-points in the direction of periodicity and 20 Å separation between images.Geometry optimizations were performed with hydrogens and methyl groups in place of the full side chains.We report structural parameters (dihedral angles), fundamental gaps, and relative energies in Tables 1 and 2. The lowest energy conformation for PTB7 is cis-cis aligned for the dihedrals, as denoted in Table 1.Similarly, the cis-cis conformer of PID2 is also lowest lying, with the other conformers much higher in energy.We note that our calculations of PTB7 conformer energies are in agreement with those of Bhatta, Perry, and Tsige, 102 while we found different minimum structures from those reported by Khoshkholgh, et al. 103 and Li, 104 and Freeman, 104 possibly due to a very limited configurational sampling in previous works.The variability of our calculated band gaps as a function of the structural motif indicates that models using a single, possibly high energy configuration are inadequate to describe the electronic properties of the blends considered here.Table 1 Dihedral angles a and b (degrees) (see Fig. 1), band gap E g (eV) at the B3LYP level of theory, and relative energies (kcal mol À1 ) for four candidate conformers of PTB7 extended chains, computed using a periodic cell containing two repeating units.Structural parameters and fundamental gaps were determined using geometries with methyl-terminated polymers, and relative energies were obtained with fixed cells of dimension 24 Å with hydrogen-terminated structures to reduce structural changes due to sidechain positions.Geometry relaxations for periodic systems were performed using the PBE density functional approximation This journal is © the Owner Societies 2016 spectroscopy) HOMO and LUMO energy levels for the ternary blend examined by Lu, et al., 24 using our lowest energy configurations for extended polymer chains and the equilibrium geometry of PC71BM.8][109][110][111][112][113][114][115][116][117] Both UPS and CV measurements reported here (corresponding to adiabatic and vertical ionization potentials, respectively) agree on the trends in electronic energy levels, as discussed in literature. 118,119The high variability in UPS threshold ionization potentials and electron affinities may stem from the effects of contaminants and differing surface orientations of materials. 111,120,121The DFT HOMO energy levels are clearly separated for each component, with spacings similar to experiment regardless of the density functional chosen.

Energy levels of ternary blend systems
The dielectric constant e N for thin films of PTB7 and PC71BM is B4; 122 hence, PBE0 (25% exact exchange) or B3LYP (20% exact exchange) are expected to yield good trends in ionization potentials and electron affinities since the optimal percent of exact exchange is inversely proportional to e N in condensed phase materials. 123

Thermal configurational disorder
As a first step toward assessing the effect of thermal disorder on polymer energy levels, we performed a first principles molecular dynamics simulation of benzodithiophene.This simulation allowed us to isolate the influence of the thermal disorder of side-chains on the electronic properties of organic polymers, and our results provided qualitative validation of classical MD simulations of the full PTB7 polymer.Fig. 4 presents the PBE HOMO and LUMO energy levels for periodic benzodithiophene configurations with varying side-chain lengths.Geometries were obtained from a first principles NVT molecular dynamics trajectory at 300 K using the Bussi-Donadio-Parrinello thermostat, 124 with 80 snapshots taken every 1.2 femtoseconds over a trajectory of 400 steps, sampled with a timestep of 0.242 femtoseconds.Increasing side-chain length shifts the average position of the HOMO energy level up and broadens the energy level positions.The presence of highly branched side-chains induces significantly less broadening, attributable to a loss of flexibility and steric hindrances.From these calculations, we concluded that the different thermal configurations of side chains induce significant band broadening, and hence they need to be taken into account in calculations of electronic energy levels.
To estimate the energy level broadening of PC71BM, we performed first principles molecular dynamics simulations of PC71BM, whose B3LYP HOMOÀ2 through LUMO+2 energy levels (calculated using B3LYP/def2-TZVP) are presented in Fig. 5. Eighty configurations were obtained in the same manner as for the benzodithiophene configurations above.The distribution of HOMO energy levels (À5.86 AE 0.06 eV) displays two peaks and a shoulder, with a dispersion of ca.0.25 eV.The distribution of LUMO energy levels (À3.47 AE 0.10 eV) displays four peaks, with a dispersion of ca.0.45 eV.The apparent width of these distributions exceeds two standard deviations, indicating that different atomic configurations accessible due to thermal disorder exhibit different electronic properties.The energy onset of the LUMO (ca.À3.7 eV) is 0.2 eV lower than the mean energy level for these configurations, which suggests that disorder of PC71BM geometries may result in a voltage loss of 0.2 V for the entire device.
Our results indicated that side chain dynamics must be included to describe PTB7, as well as changes in torsions and Table 2 Dihedral angles a and b (degrees) (see Fig. 1), band gap E g (eV) at the B3LYP level of theory, and relative energies (kcal mol À1 ) for four candidate conformers of PID2 extended chains, computed using a periodic cell containing two repeating units.Structural parameters and fundamental gaps were determined using geometries with methyl-terminated polymers, and relative energies were obtained with fixed cells of dimension 24 Å with hydrogen-terminated structures to reduce structural changes due to sidechain positions.Geometry relaxations for periodic systems were performed using the PBE density functional approximation   out-of-plane bending angles, whose timescales extend beyond those sampled by first principles molecular dynamics.Hence, we performed classical NVT molecular dynamics calculations on periodic cells with two polymer repeating units of PTB7, and we considered 200 snapshots taken every 0.5 nanoseconds after 100 nanoseconds of equilibration.Fig. 6 shows the broadening of the HOMO and LUMO energy levels at the PBE level of theory for those 200 configurations.We note the large spread for the HOMO and LUMO energy levels, with band onsets of 0.2-0.3eV from the average values of À4.95 eV (HOMO) and À3.65 eV (LUMO), possibly accounting for 0.2 V in operating voltage losses for these materials.The voltage losses due to thermal disorder will be cumulative, adding up to 0.4 V for the PTB7/ PC71BM binary mixture from the band widths of the PTB7 HOMO and the PC71BM LUMO.Since the theoretical maximum voltage for that binary mixture is 1.45 V, and the experimental open-circuit voltage is 0.72 V, 24 , voltage losses of 1.45 À 0.72 V = 0.73 V must be explained.We attribute part of these losses (0.40 V) to thermal broadening of energy levels.The remainder of the voltage losses (B0.33 V) can be reasonably attributed to the energy required to separate charges (bound by large excitonic binding energies, ca.0.3 eV, in organic materials 125 ).
In order to establish correlations between electronic and transport properties and specific structural motifs, we performed a cluster analysis (principal component analysis, 99 followed by clustering with the affinity propagation method 100 ) of the inter-unit dihedral angles of these configurations.From the 200 configurations, we identified 11 clusters.We present our results for the electronic properties (PBE HOMO, LUMO, band gap, and H AB nonadiabatic coupling for hole transfer) of each cluster within Table 3.
The calculated electronic properties for clustered PTB7 configurations in Table 3 are statistically distinguishable, as determined by Kruskal-Wallis 126 and one-way 127 analysis of variance p-values (p o 1 Â 10 À12 for all calculated values).The clusters with the highest hole nonadiabatic couplings have average dihedral angles approximately 30-40 degrees out of plane, with sequential dihedrals typically either both aligned (cis) or antialigned (trans), corresponding to the ground state models examined in Section 5.1.By comparison, clusters with a lower degree of planarity (e.g.B and G) exhibit much lower hole transport couplings and higher band gaps.Cluster I includes all-trans model geometries, and cluster F all-cis geometries; however, several motifs with low band gap and high hole transfer couplings have a cis-cis-trans-trans character (e.g.clusters D and K).These results suggest that the local pairwise alignment of dihedrals  Table 3 Dihedral angles a and b (degrees) (see Fig. 1), variation from planarity g (degrees), calculated PBE HOMO and LUMO energy levels and electronic band gap E g (eV), transfer integrals H AB (meV) (see eqn (1)), and frequency of values for different clusters generated from 300 K configurations of PTB7, as well as the average values over all the 200 configurations (overall) and overall planarity are key descriptors to optimize opto-electronic properties of semiconducting polymers.As a further simplification, we reduced the number of structural parameters to g ¼ 1 4 . the variation of dihedrals from 0 or 180 deg.(the planar configurations).Fig. 7 shows the dependence of the HOMO-LUMO gap and the alongchain hole transfer nonadiabatic couplings on g.As the chains deform, conjugation is broken, raising the HOMO-LUMO gap linearly as a function of g, and thus shifting the absorption threshold towards values less optimal for solar energy absorption.This breaking of conjugation also interrupts the charge transport along the chain, reducing the hole transfer integrals by a factor of four.Since rates depend quadratically on the transfer integrals (within Marcus theory), such a reduction translates into a decrease of rates by a factor of 16.These results highlight the key role of microphase structural changes in determining optical absorption and charge transport in OPV blends.In particular, structural aggregation motifs (induced by solvent, additives, or ternary components) are expected to change the planarity of polymer chains and thus hole conductivities and the efficiency of light harvesting.Our conclusion from this analysis is that the planarity of the backbone is highly correlated with optimal electronic gaps and efficient hole transport.

Templating of PTB7 planarity by PID2
To assess the effect of PID2 on the morphology of PTB7, we used the weighted histogram analysis method combined with umbrella sampling to determine the free energy of rotation for the middle dihedral in a PTB7 5-mer, which is the b dihedral (see Fig. 1).We determined the free energy of this rotation in two different solvents: a solvent of PTB7, and a solvent of PID2.The solvent consisted of 25 PTB7 or PID2 5-mers, depending on the sample-a system large enough to prevent self-interaction of the PTB7 5-mer of interest.We sampled the free energy of the entire dihedral rotation using 50 uniformly spaced umbrellas, each with a restraint strength of 200 kcal mol À1 .Sampling was performed at 500 K, well above the glass transition temperature of PTB7.Samples were prepared by first initializing molecules on a low-density lattice, then evolving the system with a barostat set to 1 atmosphere at 500 K until volume reached equilibrium (5 ns).A temperature above the glass transition temperature was used to ensure that the system was equilibrated.
The dihedral restraint for each umbrella was active during the entire simulation.Free energies of rotation for PTB7 in a solvent of PTB7 and PID2 solvent at 500 K are shown in Fig. 8.The free energies of rotation reproduce the distinguishing characteristics of the dihedral potential energy surfaces shown in Fig. 2, in that planar configurations are strongly favored.Comparing the free energies of rotation in a solvent of PTB7 to that in a solvent of PID2, the planar configurations centered around 0 and 180 degrees are stabilized by PID2 solvent.While the energy maxima at 90 and 270 degrees are of nearly identical magnitude in a solvent of PTB7, PID2 imparts an asymmetry, with the energy maxima centered around 90 degrees 0.7 kcal mol À1 lower than that centered around 270 degrees.We calculated the average deviation from planarity to be 36.4and 34.6 degrees for PTB7 in a solvent of PTB7 and PID2, respectively, showing that PID2 imparts a mild degree of planarity onto PTB7.Interestingly, we performed the same experiment at 400 K and calculated the average deviation from planarity to be 32.3 and 29.2 degrees.These free energies are shown in the ESI.† Planarity increases significantly as temperature drops, demonstrating that planarity is the enthalpically favorable state.Note also that PTB7 at 500 K is on average 1.8 degrees more planar in a  solvent of PID2 than in PTB7.At 400 K, this difference increases to 3.1 degrees.Our results show that PID2 leads to a higher degree of planarity of PTB7 in bulk heterojunctions, increasing charge transfer rates in donor regions and thus device efficiency.

Charge transfer at interfaces
Differences in hole transfer properties at interfaces between components in organic solar cell systems have been proposed as one reason for improved efficiency of the PTB7, PID2, and PC71BM system. 24Upon photo-excitation of PC71BM, electron hole pairs may be harvested if holes can successfully be extracted from the PC71BM microphases.To this end, we examined the hole transfer rates from PC71BM -PTB7 and -PID2 using constrained DFT for interfacial structures containing two repeating units of a polymer and one PC71BM molecule.Fig. 9 shows isosurfaces of the difference between the charged and neutral densities for a positively ionized interfacial structure of the PTB7/ PC71BM interface.Our constrained DFT implementation provided us with direct access to the nonadiabatic couplings (H AB ), as well as with differences in energy between diabatic states (DG AB ).We approximated our reorganization energy (l AB = l int.AB + l ext.AB ) using the internal reorganization energy (l int.
AB ) of the isolated, periodic polymer and the isolated PC71BM, with a 0.11 eV estimate of the external reorganization energy (l ext.AB ) in these systems (corresponding to the response of the bulk environment), which was taken from a literature study on similar materials. 128Our total reorganization energies (the sum of l int. AB and l ext.AB ) for these systems are of the order 0.2-0.3eV.The resulting rates for hole transfer at the PTB7/ PC71BM interface are presented in Table 4.In addition to Marcus theory, we also present results using a different rate equation 129 proposed by Troisi 93 for reactions with large nonadiabatic coupling (H AB Z l AB /2).A comparison of the magnitudes of our calculated nonadiabatic couplings and reorganization energies, showing H AB Z l AB /2, suggests that the use of the Troisi rate equation may be appropriate.However, we computed H AB using semi-local functionals and based on calculations reported in the literature for small molecules, 130,131 we may expect our values to be overestimated compared to hybrid functional results.Hence, which rate equation is applicable may depend upon the level of theory used to compute H AB .
Therefore, we calculated hole transfer rates using both rate equations (Marcus and Troisi) and examined trends.The different rate equations yielded substantially different results for alongchain hole transfer rates; however, the qualitative trends for hole transfer rates are similar.The calculated hole transfer rates suggest extremely efficient hole transfer along the polymer backbone (k B 10 14 -10 15 second À1 within Marcus theory), while transfer of holes from PC71BM -PTB7 is calculated to occur 6-8 orders of magnitude more slowly; this suggests that few charge carrier collection events will occur resulting from photoexcitation of the PC71BM in binary PTB7/PC71BM solar cells.
The hole transfer dynamics at the PID2/PC71BM interface are presented in Table 5.Here, hole transfer rates along the backbone of PID2 are calculated to occur at similar rates to within PTB7 (k B 10 13 -10 15 s À1 within Marcus theory), while the hole transfer from the PC71BM can occur at comparatively fast rates (k B 10 13 s À1 ) to the nearest PID2 subunit.As suggested by Lu, et al., 24 a lower driving force, DG AB , of the reaction is the primary difference for this hole extraction process, indicating this reaction is within the inverted Marcus regime. 132Thus the difference in HOMO energy level between PTB7 and PID2 is responsible for different hole transport rates at interfaces with PC71BM.Hence, our results predict that only a small portion of PID2 at interfaces with PC71BM would be necessary for significant increases in charge collection following PC71BM excitation.

Conclusions
In this paper, we addressed the question of increased efficiency of ternary OPV blends with respect to binary ones; in particular we focused on a promising system recently studied experimentally: Fig. 9 Isosurfaces of the difference between charged and neutral densities (r = r neutral À r charged ) corresponding to constrained DFT diabatic states with a hole localized on two PTB7 subunits (1-red, 2-blue) and PC71BM (green).the ternary blend of PTB7, PID2, and PC71BM, and we investigated the reasons underlying its measured absorption and transport efficiency.We used a combination of classical and ab initio MD and electronic structure calculations, and we found that sampling of conformational disorder is necessary for predictive computational design of the electronic and transport properties of OPV materials our results showed that fundamental band gaps and hole transfer nonadiabatic couplings are highly correlated with backbone dihedrals of PTB7, where planar configurations were found to exhibit superior hole transport rates and absorption onsets.
Our simulations of different solvation environments found that the presence of PID2 resulted in more planar PTB7 geometries, which improves charge transport and overall device efficiency.We also found that structural disorder is a likely candidate for voltage losses of up to 0.4 V in the ternary blends investigated here, due to the dispersion of the HOMO energy level of PTB7 and the LUMO energy level of PC71BM as a function of temperature.These fundamental voltage losses are intrinsic to disordered organic polymer semiconductors and not only to the systems considered here.We found that improved current in ternary blends can be attributed to the superior hole transfer rates from PC71BM -PID2 relative to those from PC71BM -PTB7, in agreement with the experimental results and operational picture of Lu, et al. 24 In summary, our computational study showed that the higher planarity of PID2, coupled with the improved hole extraction from PC71BM when PID2 is present, explains the nonadditive improvements seen in ternary solar cell systems. 24As a direction forward in design of OPV materials, we recommend targeting planar structures to minimize voltage and current losses.These structures may be attainable through additives and ternary components or through computationally-aided design of more disorder-resistant (e.g.rigid or planar) OPV materials.which has been modified to include the reaction driving force.We employ a frequency o = 1000 cm À1 /h , corresponding to an optical phonon frequency which is representative for most organic semiconductors.Since o enters primarily as a prefactor, it serves to rescale rates, and errors here will not change trends.

Fig. 1
Fig. 1 Structures of compounds included in the ternary blend examined in this paper.Two different dihedral angles a and b (S-C-C-S) between subunits are labeled for use in distinguishing configurations.

Fig. 4
Fig. 4 Density of states (DOS) of calculated PBE HOMO and LUMO energy levels (eV) for 80 benzodithiophene configurations with increasing R group complexity.Plotted densities of states have been vertically shifted as a visual aide.A Gaussian broadening of 0.02 eV was applied to computed energy levels.

Fig. 5
Fig. 5 Density of states (DOS) of calculated B3LYP HOMOÀ2, HOMOÀ1, HOMO, LUMO, LUMO+1, and LUMO+2 energy levels (eV) for 80 PC71BM configurations, with each energy level marked by color.A Gaussian broadening of 0.02 eV was applied to computed energy levels.

Fig. 6
Fig. 6 Density of states (DOS) of calculated PBE HOMO and LUMO energy levels for PTB7 extended chains.A Gaussian broadening of 0.02 eV was used in smoothing the results.

Fig. 7 4 Pi
Fig. 7 Electronic band gap obtained at the PBE level of theory (top) and nonadiabatic coupling (bottom) versus the variation of dihedrals from planarity g ¼ 1 4 P i

Fig. 8
Fig.8Free energies (kcal mol À1 ) of rotation for a central PTB7 b dihedral angle in degrees (see Fig.1) in two different solvents (PTB7 and PID2) at 500 K.
This journal is © the Owner Societies 2016