Coarse-grained simulations of the solution-phase self-assembly of poly(3-hexylthiophene) nanostructures

A Molecular Dynamics Study of Mechanical and Conformational Properties of Conjugated Polymer Thin Films

Understanding and predicting the mechanical and conformational properties of conjugated polymer (CP) thin films are a central focus in flexible electronic device research. Employing molecular dynamics simulations with an architecture-transferable chemistry-specific coarse-grained (CG) model of poly(3-alkylthiophene)s (P3ATs), developed by using an energy renormalization approach, we investigate the mechanical and conformational behavior of P3AT thin films during deformation. The density profiles and measures of local mobility identify a softer interfacial layer for all films, the thickness of which does not depend on M w or side-chain length. Remarkably, Young's modulus measured via nanoindentation is more sensitive to M w than for tensile tests, which we attribute to distinct deformation mechanisms. High-M w thin films show increased toughness, whereas longer side-chain lengths of P3AT resulted in lower Young's modulus. Fractures in low-M w thin films occur through chain pullout due to insufficient chain entanglement and crazing in the plastic region. Importantly, stretching promoted both chain alignment and longer conjugation lengths of P3AT, potentially enhancing its electronic properties. For instance, at room temperature, stretching P3HT thin films to 150% increases the conjugated length of P3HT thin films from 2.7 nm to 4.7 nm, aligning with previous experimental findings and all-atom simulation results. Furthermore, high-M w thin films display elevated friction forces due to the chain accumulation on the indenter, with negligible variations in the friction coefficient across all thin film systems. These findings offer valuable insights that enhance our understanding and guide the rational design of CP thin films in flexible electronics.

Competing single-chain folding and multi-chain aggregation pathways control solution-phase aggregate morphology of organic semiconducting polymers

Understanding the solution-phase behaviour of organic semiconducting polymers is important for systematically improving the performance of devices based on solution-processed thin films of these molecules. Conventional polymer theory predicts that polymer conformations become more compact as solvent quality decreases, but recent experiments have shown the high-performance organic-semiconducting polymer P(NDI2OD-T2) to form extended rod-like aggregates much larger than a single chain in poor solvents, with the formation of these extended aggregates correlated with enhanced electron mobility in films deposited from these solutions. We explain the unexpected formation of extended aggregates using a novel coarse-grained simulation model of P(NDI2OD-T2) that we have developed to study the effect of solvent quality on its solution-phase behaviour. In poor solvents, we find that aggregation through only a few monomers gives effectively inseparable chains, leading to the formation of extended structures of partially overlapping chains via non-equilibrium assembly. This behaviour requires that multi-chain aggregation occurs faster than chain folding, which we show is the case for the chain lengths and concentrations shown experimentally to form rod-like aggregates. This kinetically controlled process introduces a dependence of aggregate structure on concentration, chain length, and chain flexibility, which we show is able to reconcile experimental findings and is generalisable to the solution-phase assembly of other semiflexible polymers.

Coarse-Grained Model Incorporating Short- and Long-Range Effective Potentials for the Fast Simulation of Micelle Formation in Solutions of Ionic Surfactants

A coarse-grained model comprising short- and long-range effective potentials, parametrized with the iterative Boltzmann inversion (IBI) method, is presented for capturing micelle formation in aqueous solutions of ionic surfactants using as a model system sodium dodecyl sulfate (SDS). In the coarse-grained (CG) model, each SDS molecule is represented as a sequence of four beads while each water molecule is modeled as a single bead. The proposed CG scheme involves ten potential energy functions: four of them describe bonded interactions and control the distribution functions of intramolecular degrees of freedom (bond lengths, valence angles, and dihedrals) along an SDS molecule while the other six account for intermolecular interactions between pairs of SDS and water beads and control the radial distribution functions. The nonbonded effective potentials between coarse-grained SDS molecules extend up to about 12 nm and capture structural and morphological features of the micellar solution both at short and long distances. The long-range component of these potentials, in particular, captures correlations between surfactant molecules belonging to different micelles and is essential to describe ordering associated with micelle formation. A new strategy is introduced for determining the effective potentials through IBI by using information (target distribution functions) extracted from independent atomistic simulations of a micellar reference system (a salt-free SDS solution at total surfactant concentration c_T equal to 103 mM, temperature T equal to 300 K, and pressure P equal to 1 atm) obtained through a multiscale approach described in an earlier study. It employs several optimization steps for bonded and nonbonded interactions and a gradual parametrization of the short- and long-range components of the latter, followed by reparametrization of the bonded ones. The proposed CG model can reproduce remarkably accurately the microstructure and morphology of the reference system within only a few hours of computational time. It is therefore very promising for future studies of structural and morphological behavior of various liquid surfactant formulations.

Development of hybrid coarse-grained atomistic models for rapid assessment of local structuring of polymeric semiconductors

Decades of work in the field of computational study of semiconducting polymers using atomistic models illustrate the challenges of generating equilibrated models for this class of materials. While adopting a coarse-grained model can be helpful, the process of developing a suitable model is particularly non-trivial and time-consuming for semiconducting polymers due to a large number of different interactions with some having an anisotropic nature. This work introduces a procedure for the rapid generation of a hybrid model for semiconducting polymers where atoms of secondary importance (those in the alkyl side chains) are transformed into coarse-grained beads to reduce the computational cost of generating an equilibrated structure. The parameters are determined from easy-to-equilibrate simulations of very short oligomers and the model is constructed to enable a very simple back-mapping procedure to reconstruct geometries with atomistic resolution. The model is illustrated for three related polymers containing DPP (diketopyrrolopyrrole) to evaluate the transferability of the potential across different families of polymers. The accuracy of the model, determined by comparison with the results of fully equilibrated simulations of the same material before and after back-mapping, is fully satisfactory for two out of the three cases considered. We noticed that accuracy can be determined very early in the workflow so that it is easy to assess when the deployment of this method is advantageous. The hybrid representation can be used to evaluate directly the electronic properties of structures sampled by the simulations.

Temperature-Dependent Conformation Behavior of Isolated Poly(3-hexylthiopene) Chains

We use atomistic as well as coarse-grained molecular dynamics simulations to study the conformation of a single poly(3-hexylthiopene) chain as a function of temperature. We find that mainly bundle and toroid structures appear with bundles becoming more abundant for decreasing temperatures. We compare an atomistic and a Martini-based coarse-grained model which we find in very good agreement. We further illustrate how the temperature dependence of P3HT can be connected to that of simple Lennard–Jones model polymers in a vacuum. Upon adding solvent (THF) we observe the occurrence of a prominent swelling of the molecular size at a temperature of about 220 K. This swelling is in close agreement with the interpretation of recent spectroscopic experiments which allows us to explain the experimental observations by an increased frequency of bundle structures.

Strategies for the Development of Conjugated Polymer Molecular Dynamics Force Fields Validated with Neutron and X-ray Scattering

Conjugated polymers (CPs) enable a wide range of lightweight, lower cost, and flexible organic electronic devices, but a thorough understanding of relationships between molecular structure and dynamics and electronic performance is critical for improved device efficiencies and for new technologies. Molecular dynamics (MD) simulations offer in silico insight into this relationship, but their accuracy relies on the approach used to develop the model's parameters or force field (FF). In this Perspective, we first review current FFs for CPs and find that most of the models implement an arduous reparameterization of inter-ring torsion potentials and partial charges of classical FFs. However, there are few FFs outside of simple CP molecules, e.g., polythiophenes, that have been developed over the last two decades. There is also limited reparameterization of other parameters, such as nonbonded Lennard-Jones interactions, which we find to be directly influenced by conjugation in these materials. We further provide a discussion on experimental validation of MD FFs, with emphasis on neutron and X-ray scattering. We define multiple ways in which various scattering methods can be directly compared to results of MD simulations, providing a powerful experimental validation metric of local structure and dynamics at relevant length and time scales to charge transport mechanisms in CPs. Finally, we offer a perspective on the use of neutron scattering with machine learning to enable high-throughput parametrization of accurate and experimentally validated CP FFs enabled not only by the ongoing advancements in computational chemistry, data science, and high-performance computing but also using oligomers as proxies for longer polymer chains during FF development.

Performance improvement of P3HT nanowire-based organic solar cells by interfacial morphology engineering

Poly (3-hexylthiophene-2,5-diyl) nanowires (nw-P3HT) have been a great interest for organic electronics, including organic field-effect transistors, organic photodetectors, organic photovoltaics, etc due to easy formation in the solution process. Thus both explanations of charge transport dynamics and morphology are crucial for device performance. Here we demonstrated the optoelectronic properties of the P3HT nanowires where the polymer backbones were parallel to the nanowire axis. The nanowires tended to form a bundle due to van der Waals interactions. Nanowire bundles were separated by 1,8-diiodooctane (DIO) additive for photovoltaic fabrication. The bundle separation was visualized by atomic force microscopy. The charge transfer mechanism was evaluated by electrochemical impedance spectroscopy. The electrical analysis showed that short-circuit current density (J _sc) increases to 10.74 mA cm^-2 after the bundle separation. According to impedance analysis, there is a correlation between effective lifetime and DIO ratio. These findings were considered as promising results for improving the transport by forming pathways for charge carriers.

Predicting optical spectra for optoelectronic polymers using coarse-grained models and recurrent neural networks

Coarse-grained modeling of conjugated polymers has become an increasingly popular route to investigate the physics of organic optoelectronic materials. While ultraviolet (UV)-vis spectroscopy remains one of the key experimental methods for the interrogation of these materials, a rigorous bridge between simulated coarse-grained structures and spectroscopy has not been established. Here, we address this challenge by developing a method that can predict spectra of conjugated polymers directly from coarse-grained representations while avoiding repetitive procedures such as ad hoc back-mapping from coarse-grained to atomistic representations followed by spectral computation using quantum chemistry. Our approach is based on a generative deep-learning model: the long-short-term memory recurrent neural network (LSTM-RNN). The latter is suggested by the apparent similarity between natural languages and the mathematical structure of perturbative expansions of, in our case, excited-state energies perturbed by conformational fluctuations. We also use this model to explore the level of sensitivity of spectra to the coarse-grained representation back-mapping protocol. Our approach presents a tool uniquely suited for improving postsimulation analysis protocols, as well as, potentially, for including spectral data as input in the refinement of coarse-grained potentials.

Polaron hopping barriers and rates in semiconducting polymers

Conjugated polymers are potential next-generation materials for organic electronic devices. The ability of these materials to transport charges is a key factor limiting their performance. Charge carriers in conjugated polymers are localized by disorder and polaronic effects. Charge transport in these materials is often described by thermally activated hopping, with a rate given by Marcus theory. The polaron hopping activation energy determines the temperature dependence of the Marcus rate. This energy barrier is dictated by the transition state, in which the charge carrier is equally divided between the initial and final locations. The prefactor for the polaron hopping rate is set by the charge tunneling rate between the initial and final locations. We use a tight-binding polaron model, in which charge carriers are stabilized by both nuclear reorganization and polarization of the surrounding dielectric, to compute the activation energy, charge tunneling rate and overall rate constant for intrachain and interchain charge hopping processes in poly(3-hexylthiophene) (P3HT) crystalline lamellae and amorphous melts. Charge transport in these two environments is limited by interchain hopping processes. Both hopping barriers and rates predicted by the model are in good agreement with experiments on a variety of crystalline and amorphous P3HT materials. Qualitatively, the barriers largely depend on how well the transition state is stabilized by polarization effects, and on the hopping integral between the initial and final locations, both of which penalize hopping over longer distances.

The interplay of interfaces, supramolecular assembly, and electronics in organic semiconductors

Organic semiconductors, which include a diverse range of carbon-based small molecules and polymers with interesting optoelectronic properties, offer many advantages over conventional inorganic semiconductors such as silicon and are growing in importance in electronic applications. Although these materials are now the basis of a lucrative industry in electronic displays, many promising applications such as photovoltaics remain largely untapped. One major impediment to more rapid development and widespread adoption of organic semiconductor technologies is that device performance is not easily predicted from the chemical structure of the constituent molecules. Fundamentally, this is because organic semiconductor molecules, unlike inorganic materials, interact by weak non-covalent forces, resulting in significant structural disorder that can strongly impact electronic properties. Nevertheless, directional forces between generally anisotropic organic-semiconductor molecules, combined with translational symmetry breaking at interfaces, can be exploited to control supramolecular order and consequent electronic properties in these materials. This review surveys recent advances in understanding of supramolecular assembly at organic-semiconductor interfaces and its impact on device properties in a number of applications, including transistors, light-emitting diodes, and photovoltaics. Recent progress and challenges in computer simulations of supramolecular assembly and orientational anisotropy at these interfaces is also addressed.

Solving the Trivial Crossing Problem While Preserving the Nodal Symmetry of the Wave Function

In an adiabatic mixed quantum-classical simulation, the avoided crossing of weakly coupled eigenstates can lead to unphysical discontinuities in wave function dynamics, otherwise known as the trivial crossing problem. A standard solution to the trivial crossing problem eliminates spatial discontinuities in wave function dynamics by imposing changes to the eigenstate of the wave function. In this paper, we show that this solution has the side effect of introducing transient discontinuities in the nodal symmetry of the wave function. We present an alternative solution to the trivial crossing problem that preserves both the spatial and nodal structure of the adiabatic wave function. By considering a model of exciton dynamics on conjugated polymer systems, we show that failure to preserve wave function symmetry yields exciton dynamics that depends unphysically on polymer system size. We demonstrate that our symmetry-preserving solution to the trivial crossing problem yields more realistic dynamics and can thus improve the accuracy of simulations of larger systems that are prone to the trivial crossing problem.

Theoretical and experimental studies of an oseltamivir-triazole-based thermoresponsive organogel

Low-molecular weight organic gelators have been of significant interest in recent years because of their interesting properties and potential applications in sensing technology, biomedicine and drug delivery. Herein, the synthesis, characterization and gelation properties of new oseltamivir conjugates are reported. The oseltamivir-triazole conjugate 1 was synthesized via a click-reaction in a 75% yield. The key features of this conjugate include the presence of amide, flexible ester linkages and a triazole scaffold linking a hydrophobic alkyl chain. The conjugate 1, possessing a long alkyl chain, showed gelation properties in various apolar organic solvents. This gelation behavior was not observed in the case of the deesterified conjugate 2; this indicated the necessity of the alkyl chain for gelation. The gelator 1 showed thermoreversible gelation properties in a range of linear alkane solvents (from n-pentane to n-dodecane). A scanning electron microscopic study suggests that the gelator 1 exists as cross-linked structures, which are self-aggregated in the range of submicrometers, as supported by extensive ¹H-NMR studies. The rheological parameters supported the occurrence of a soft gelation process, and the gel formed in n-decane was found to be stiffer than that formed in n-hexane. Computational studies suggested that the gelation behavior was indeed due to micelle formation and dependent on the lipophilicity of solvents.

Assessment of molecular dynamics simulations for amorphous poly(3-hexylthiophene) using neutron and X-ray scattering experiments

The molecular morphology and dynamics of conjugated polymers in the bulk solid state play a significant role in determining macroscopic charge transport properties. To understand this relationship, molecular dynamics (MD) simulations and quantum mechanical calculations are used to evaluate local electronic properties. In this work, we investigate the importance of system and simulation parameters, such as force fields and equilibration methods, when simulating amorphous poly(3-hexylthiophene) (P3HT), a model semiconducting polymer. An assessment of MD simulations for five different published P3HT force fields is made by comparing results to experimental wide-angle X-ray scattering (WAXS) and to a broad range of quasi-elastic neutron scattering (QENS) data. Moreover, an in silico analysis of force field parameters reveals that atomic partial charges and torsion potentials along the backbone and side chains have the greatest impact on structure and dynamics related to charge transport mechanisms in P3HT.

Thermal Fluctuations Lead to Cumulative Disorder and Enhance Charge Transport in Conjugated Polymers

All conjugated polymers examined to date exhibit significant cumulative lattice disorder, although the origin of this disorder remains unclear. Using atomistic molecular dynamics (MD) simulations, the detailed structures for single crystals of a commonly studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) are obtained. It is shown that thermal fluctuations of thiophene rings lead to cumulative disorder of the lattice with an effective paracrystallinity of about 0.05 in the π-π stacking direction. The thermal-fluctuation-induced lattice disorder can in turn limit the apparent coherence length that can be observed in diffraction experiments. Calculating mobilities from simulated crystal structures demonstrates that thermal-fluctuation-induced lattice disorder even enhances charge transport in P3HT. The mean inter-chain charge transfer integral is enhanced with increasing cumulative lattice disorder, which in turn leads to pathways for fast charge transport through crystals.

Generic Model for Lamellar Self-Assembly in Conjugated Polymers: Linking Mesoscopic Morphology and Charge Transport in P3HT

We develop a generic coarse-grained model of soluble conjugated polymers, capable of describing their self-assembly into a lamellar mesophase. Polymer chains are described by a hindered-rotation model, where interaction centers represent entire repeat units, including side chains. We introduce soft anisotropic nonbonded interactions to mimic the potential of mean force between atomistic repeat units. The functional form of this potential reflects the symmetry of the molecular order in a lamellar mesophase. The model can generate both nematic and lamellar (sanidic smectic) molecular arrangements. We parametrize this model for a soluble conjugated polymer poly(3-hexylthiophene) (P3HT) and demonstrate that the simulated lamellar mesophase matches morphologies of low molecular weight P3HT, experimentally observed at elevated temperatures. A qualitative charge-transport model allows us to link local chain conformations and mesoscale order to charge transport. In particular, it shows how coarsening of lamellar domains and chain extension increase the charge carrier mobility. By modeling large systems and long chains, we can capture transport between lamellar layers, which is due to rare, but thermodynamically allowed, backbone bridges between neighboring layers.

Brownian orientational lath model (BOLD): A computational model relating the self-assembly in a fluid of lath like particles with its rheology and gelation

By means of a computational model, we study the relation between two complementary views of gelation, rheological tests against the characterization of a network of consecutive particles. The model we propose consists of slender, plane, colloidal sized particles, which we name laths, which self-assemble into long ordered aggregates of several particles; called whiskers in the literature. Within a whisker, the interaction potential is a minimum when: the planes of two consecutive laths are aligned, thus favoring their alignment; when the center of three consecutive laths lie in a straight line, thus favoring stacking; and when the center of two consecutive laths are located at a certain distance, which mimics excluded volume. A threshold value of the potential gives a condition for sticking free laths into whiskers, and for the breaking of whiskers. The simplicity of the model allows the simulation to reach large enough times, of the order of minutes, needed to simulate numerical rheology tests. We are able to characterize the whisker formation, as well as to simulate the gel transition, by means of an oscillatory shear numerical experiment. We conclude that according to the usual rheological definition a gel transition occurs at about 250K, even though there is no branching and less than 10% of whiskers are long enough as to percolate the system.

Polaron formation mechanisms in conjugated polymers

In semiconducting polymers, interactions with conformational degrees of freedom can localize charge carriers, and strongly affect charge transport. Polarons can form when charges induce deformations of the surrounding medium, including local vibrational modes or dielectric polarization. These deformations then interact attractively with the charge, tending to localize it. First we investigate vibrational polaron formation in poly(3-hexylthiophene) [P3HT], with a tight-binding model for charges hopping between adjacent rings, coupled to ring distortions. We use density functional theory (DFT) calculations to determine coupling constants, including the "spring constant" for ring distortions and the coupling to the charge carrier. On single chains, we find only broad, weakly bound polarons by this mechanism. In 2d crystalline layers of P3HT, even weak transverse hopping between chains destabilizes this polaron. Then, we consider polarons stabilized by dielectric polarization, described semiclassically with a polarizable continuum interacting with the carrier wavefunction. In contrast to vibrational polarons, we find dielectrically stabilized polarons in P3HT are narrower, more strongly bound, and stable in 2d crystalline layers.

Morphology of a self-doped conducting oligomer for green energy applications

A recently synthesized self-doped conducting oligomer, salt of bis[3,4-ethylenedioxythiophene]3thiophene butyric acid, ETE-S, is a novel promising material for green energy applications. Recently, it has been demonstrated that it can polymerize in vivo, in plant systems, leading to a formation of long-range conducting wires, charge storage and supercapacitive behaviour of living plants. Here we investigate the morphology of ETE-S combining the experimental characterisation using Grazing Incidence Wide Angle X-ray Scattering (GIWAXS) and atomistic molecular dynamics (MD) simulations. The GIWAXS measurements reveal a formation of small crystallites consisting of π-π stacked oligomers (with the staking distance 3.5 Å) that are further organized in h00 lamellae. These experimental results are confirmed by MD calculations, where we calculated the X-ray diffraction pattern and the radial distribution function for the distance between ETE-S chains. Our MD simulations also demonstrate the formation of the percolative paths for charge carriers that extend throughout the whole structure, despite the fact that the oligomers are short (6-9 rings) and crystallites are thin along the π-π stacking direction, consisting of only two or three π-π stacked oligomers. The existence of the percolative paths explains the previously observed high conductivity in in vivo polymerized ETE-S. We also explored the geometrical conformation of ETE-S oligomers and the bending of their aliphatic chains as a function of the oligomer lengths.

Simulating charge transport in organic semiconductors and devices: a review

Charge transport simulation can be a valuable tool to better understand, optimise and design organic transistors (OTFTs), photovoltaics (OPVs), and light-emitting diodes (OLEDs). This review presents an overview of common charge transport and device models; namely drift-diffusion, master equation, mesoscale kinetic Monte Carlo and quantum chemical Monte Carlo, and a discussion of the relative merits of each. This is followed by a review of the application of these models as applied to charge transport in organic semiconductors and devices, highlighting in particular the insights made possible by modelling. The review concludes with an outlook for charge transport modelling in organic electronics.

Controlled growth mechanism of poly (3-hexylthiophene) nanowires

Synthesis of 1D-polymer nanowires by a self-assembly method using marginal solvents is an attractive technique. While the formation mechanism is poorly understood, this method is essential in order to control the growth of nanowires. Here we visualized the time-dependent assembly of poly (3-hexyl-thiophene-2,5-diyl) (P3HT) nanowires by atomic force microscopy and scanning tunneling microscopy. The assembly of P3HT nanowires was carried out at room temperature by mixing cyclohexanone (CHN), as a poor solvent, with polymer solution in 1,2-dichlorobenzene (DCB). Both π-π stacking and planarization, obtained at the mix volume ratio of P3HT (in DCB):CHN (10:7), were considered during the investigation. We find that the length of nanowires was determined by the ordering of polymers in the polymer repetition direction. Additionally, our density functional theory calculations revealed that the presence of DCB and CHN molecules that stabilize the structural distortions due to tail group of polymers was essential for the core-wire formation.

Tight binding model of conformational disorder effects on the optical absorption spectrum of polythiophenes

Semiconducting polymers are soft materials with many conformational degrees of freedom. The limited understanding of how conformational disorder affects their optoelectronic properties is a key source of difficulties that limits their widespread usage in electronic devices. We develop a coarse-grained approach based on the tight binding approximation to model the electronic degrees of freedom of polythiophene chains, taking into account conformational degrees of freedom. Particularly important is dihedral disorder, which disrupts extended electronic states. Our tight binding model is parameterized using density functional theory (DFT) calculations of the one-dimensional band structures for chains with imposed periodic variations in dihedral angles. The model predicts valence and conduction bands for these chain conformations that compare well to DFT results. As an initial application of our model, we compute the optical absorption spectrum of poly(3-hexylthiophene) chains in solution. We observe a broadening of the absorption edge resulting from dihedral disorder, just shy of the experimental broadening. We conclude that the effects of molecular disorder on the optoelectronic properties of conjugated polymer single chains can be mostly accounted for by torsional disorder alone.

The Effect of Molecular Structure and Environment on the Miscibility and Diffusivity in Polythiophene-Methanofullerene Bulk Heterojunctions: Theory and Modeling with the RISM Approach

Although better means to model the properties of bulk heterojunction molecular blends are much needed in the field of organic optoelectronics, only a small subset of methods based on molecular dynamics- and Monte Carlo-based approaches have been hitherto employed to guide or replace empirical characterization and testing. Here, we present the first use of the integral equation theory of molecular liquids in modelling the structural properties of blends of phenyl-C61-butyric acid methyl ester (PCBM) with poly(3-hexylthiophene) (P3HT) and a carboxylated poly(3-butylthiophene) (P3BT), respectively. For this, we use the Reference Interaction Site Model (RISM) with the Universal Force Field (UFF) to compute the microscopic structure of blends and obtain insight into the miscibility of its components. Input parameters for RISM, such as optimized molecular geometries and charge distribution of interaction sites, are derived by the Density Functional Theory (DFT) methods. We also run Molecular Dynamics (MD) simulation to compare the diffusivity of the PCBM in binary blends with P3HT and P3BT, respectively. A remarkably good agreement with available experimental data and results of alternative modelling/simulation is observed for PCBM in the P3HT system. We interpret this as a step in the validation of the use of our approach for organic photovoltaics and support of its results for new systems that do not have reference data for comparison or calibration. In particular, for the less-studied P3BT, our results show that expectations about its performance in binary blends with PCBM may be overestimated, as it does not demonstrate the required level of miscibility and short-range structural organization. In addition, the simulated mobility of PCBM in P3BT is somewhat higher than what is expected for polymer blends and falls into a range typical for fluids. The significance of our predictive multi-scale modelling lies in the insights it offers into nanoscale morphology and charge transport behaviour in multi-component organic semiconductor blends.

High aspect ratio conjugated polymer nanowires for high performance field-effect transistors and phototransistors

We synthesized a highly crystalline DPP-based polymer, DPPBTSPE, which contained 1,2-bis(5-(thiophen-2-yl)selenophen-2-yl)ethene as a planar and rigid electron donating group. High- and low-molecular weight (MW) DPPBTSPE fractions were collected by Soxhlet extraction and were employed to investigate their unique charge transport properties in macroscopic films and single crystalline polymer nanowire (SC-PNW), respectively. The low-MW polymer could provide well-isolated and high aspect ratio SC-PNWs, in which the direction of π-π stacking was perpendicular to the wire growing axis. The field effect transistors made of SC-PNWs exhibited remarkably high carrier mobility of 24 cm(2) V(-1) s(-1). In addition, phototransistors (PTs) made of SC-PNW showed very high performance in terms of photoresponsivity (R) and photoswitching ratio (P). The average R of the SC PNW-based PTs were in the range of 160-170 A W(-1) and the maximum R was measured at 1920 A W(-1), which is almost three orders higher than that of thin film-based PT device.

Molecular-Scale Understanding of Cohesion and Fracture in P3HT:Fullerene Blends

Quantifying cohesion and understanding fracture phenomena in thin-film electronic devices are necessary for improved materials design and processing criteria. For organic photovoltaics (OPVs), the cohesion of the photoactive layer portends its mechanical flexibility, reliability, and lifetime. Here, the molecular mechanism for the initiation of cohesive failure in bulk heterojunction (BHJ) OPV active layers derived from the semiconducting polymer poly(3-hexylthiophene) [P3HT] and two monosubstituted fullerenes is examined experimentally and through molecular-dynamics simulations. The results detail how, under identical conditions, cohesion significantly changes due to minor variations in the fullerene adduct functionality, an important materials consideration that needs to be taken into account across fields where soluble fullerene derivatives are used.