Entanglements in P3HT and their influence on thin-film mechanical properties: Insights from molecular dynamics simulations

Impact of doping on the mechanical properties of conjugated polymers

Conjugated polymers exhibit a unique portfolio of electrical and electrochemical behavior, which - paired with the mechanical properties that are typical for macromolecules - make them intriguing candidates for a wide range of application areas from wearable electronics to bioelectronics. However, the degree of oxidation or reduction of the polymer can strongly impact the mechanical response and thus must be considered when designing flexible or stretchable devices. This tutorial review first explores how the chain architecture, processing as well as the resulting nano- and microstructure impact the rheological and mechanical properties. In addition, different methods for the mechanical characterization of thin films and bulk materials such as fibers are summarized. Then, the review discusses how chemical and electrochemical doping alter the mechanical properties in terms of stiffness and ductility. Finally, the mechanical response of (doped) conjugated polymers is discussed in the context of (1) organic photovoltaics, representing thin-film devices with a relatively low charge-carrier density, (2) organic thermoelectrics, where chemical doping is used to realize thin films or bulk materials with a high doping level, and (3) organic electrochemical transistors, where electrochemical doping allows high charge-carrier densities to be reached, albeit accompanied by significant swelling. In the future, chemical and electrochemical doping may not only allow modulation and optimization of the electrical and electrochemical behavior of conjugated polymers, but also facilitate the design of materials with a tunable mechanical response.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

Microstructural and Thermal Transport Properties of Regioregular Poly(3-hexylthiophene-2,5-diyl) Thin Films

Polymeric thin films offer a wide range of exciting properties and applications, with several advantages compared to inorganic counterparts. The thermal conductivity of such thin films ranges typically between 0.1-1 W m-1 K-1. This low thermal conductivity can cause problems with heat dissipation in various applications. Detailed knowledge about thermal transport in polymeric thin films is desired to overcome these shortcomings, especially in light of the multitude of possible microstructures for semi-crystalline thin films. Therefore, poly(3-hexylthiophene-2,5-diyl) (P3HT) is chosen as a model system to analyze the microstructure and optoelectronic properties using X-ray scattering and absorption spectra along with the thermal transport properties using the photoacoustic technique. This combination of analysis methods allows for determining the optoelectronic and thermal transport properties on the same specimen, supplemented by structural information. The effect of different molecular weights and solvents during film preparation is systematically examined. A variation of the optoelectronic properties, mainly regarding molecular weight, is apparent, while no direct influence of the solvent during preparation is discernible. In contrast, the thermal conductivities of all films examined fall within a similar range. Therefore, the microstructural properties in the ordered regions do not significantly affect the resulting thermal properties in the sample space investigated in this work. We conclude that it is mainly the amorphous regions that determine the thermal transport properties, as these represent a bottleneck for thermal transport.

Influence of Asphaltene Modification on Structure of P3HT/Asphaltene Blends: Molecular Dynamics Simulations

Further development and commercialization of bulk heterojunction (BHJ) solar cells require the search for novel low-cost materials. The present study addresses the relations between the asphaltenes' chemical structure and the morphology of the poly(3-hexylthiohene) (P3HT)/asphaltene blends as potential materials for the design of BHJ solar cells. By means of all-atom molecular dynamics simulations, the formation of heterophase morphology is observed for the P3HT-based blends with carboxyl-containing asphaltenes, as well as the aggregation of the asphaltenes into highly ordered stacks. Although the π-π interactions between the polyaromatic cores of the asphaltenes in solutions are sufficient for the molecules to aggregate into ordered stacks, in a blend with a conjugated polymer, additional stabilizing factors are required, such as hydrogen bonding between carboxyl groups. It is found that the asphaltenes' aliphatic side groups may improve significantly the miscibility between the polymer and the asphaltenes, thereby preventing the formation of heterophase morphology. The results also demonstrate that the carboxyl-containing asphaltenes/P3HT ratio should be at least 1:1, as a decrease in concentration of the asphaltenes leads to the folding of the polymer chains, lower ordering in the polymer phase and the destruction of the interpenetrating 3D structure formed by P3HT and the asphaltene phases. Overall, the results of the present study for the first time reveal the aggregation behavior of the asphaltenes of varying chemical structures in P3HT, as well the influence of their presence and concentration on the polymer phase structure and blend morphology, paving the way for future development of BHJ solar cells based on the conjugated polymer/asphaltene blends.

Material Design and Device Fabrication Strategies for Stretchable Organic Solar Cells

Recent advances in the power conversion efficiency (PCE) of organic solar cells (OSCs) have greatly enhanced their commercial viability. Considering the technical standards (e.g., mechanical robustness) required for wearable electronics, which are promising application platforms for OSCs, the development of fully stretchable OSCs (f-SOSCs) should be accelerated. Here, a comprehensive overview of f-SOSCs, which are aimed to reliably operate under various forms of mechanical stress, including bending and multidirectional stretching, is provided. First, the mechanical requirements of f-SOSCs, in terms of tensile and cohesion/adhesion properties, are summarized along with the experimental methods to evaluate those properties. Second, essential studies to make each layer of f-SOSCs stretchable and efficient are discussed, emphasizing strategies to simultaneously enhance the photovoltaic and mechanical properties of the active layer, ranging from material design to fabrication control. Key improvements to the other components/layers (i.e., substrate, electrodes, and interlayers) are also covered. Lastly, considering that f-SOSC research is in its infancy, the current challenges and future prospects are explored.

Molecular Dynamics Study of the Thermodynamics of Integer Charge Transfer vs Charge-Transfer Complex Formation in Doped Conjugated Polymers

Molecular dopants such as 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ) can interact with conjugated polymers such as poly(3-hexylthiophene-2,5-diyl) (P3HT) in two different ways: they can undergo integer charge transfer (ICT) or they can form a partial-charge-transfer complex (CTC). Both are seen experimentally, but the CTC has been challenging to characterize, making it difficult to answer questions such as the following. Which polymorph is more stable? Do they have similar barriers for formation? Is there a thermodynamic route to convert one to the other? Here, we study the structure and the thermodynamics of bulk F4TCNQ-doped P3HT with all-atom molecular dynamics simulations, using thermodynamic integration to calculate the relative free energies. We find that the ICT and CTC polymorphs have similar thermodynamic stabilities. The barrier to create the ICT polymorph, however, is lower than that to make the CTC polymorph, because the ICT polymorph has a small critical nucleus, but the critical nucleus for the CTC polymorph is larger than what we can simulate. Moreover, simulated thermal annealing shows that the activation barrier for converting the CTC polymorph to the ICT polymorph is relatively modest. Overall, the simulations explain both the observed structures and the thermodynamics of F4TCNQ-doped P3HT and offer guidelines for targeting the production of a desired polymorph for different applications.

Molecular Weight Tuning of Organic Semiconductors for Curved Organic-Inorganic Hybrid X-Ray Detectors

Curved X-ray detectors have the potential to revolutionize diverse sectors due to benefits such as reduced image distortion and vignetting compared to their planar counterparts. While the use of inorganic semiconductors for curved detectors are restricted by their brittle nature, organic-inorganic hybrid semiconductors which incorporated bismuth oxide nanoparticles in an organic bulk heterojunction consisting of poly(3-hexylthiophene-2,5-diyl) (P3HT) and [6,6]-phenyl C71 butyric acid methyl ester (PC₇₀ BM) are considered to be more promising in this regard. However, the influence of the P3HT molecular weight on the mechanical stability of curved, thick X-ray detectors remains less well understood. Herein, high P3HT molecular weights (>40 kDa) are identified to allow increased intermolecular bonding and chain entanglements, resulting in X-ray detectors that can be curved to a radius as low as 1.3 mm with low deviation in X-ray response under 100 repeated bending cycles while maintaining an industry-standard dark current of <1 pA mm^-2 and a sensitivity of ≈ 0.17 μC Gy^-1 cm^-2 . This study identifies a crucial missing link in the development of curved detectors, namely the importance of the molecular weight of the polymer semiconductors used.

CONI-Net: Machine Learning of Separable Intermolecular Force Fields

Noncovalent interactions (NCIs) play an essential role in soft matter and biomolecular simulations. The ab initio method symmetry-adapted perturbation theory allows a precise quantitative analysis of NCIs and offers an inherent energy decomposition, enabling a deeper understanding of the nature of intermolecular interactions. However, this method is limited to small systems, for instance, dimers of molecules. Here, we present a scale-bridging approach to systematically derive an intermolecular force field from ab initio data while preserving the energy decomposition of the underlying method. We apply the model in molecular dynamics simulations of several solvents and compare two predicted thermodynamic observables-mass density and enthalpy of vaporization-to experiments and established force fields. For a data set limited to hydrocarbons, we investigate the extrapolation capabilities to molecules absent from the training set. Overall, despite the affordable moderate quality of the reference ab initio data, we find promising results. With the straightforward data set generation procedure and the lack of target data in the fitting process, we have developed a method that enables the rapid development of predictive force fields with an extra dimension of insights into the balance of NCIs.

Directly Probing the Fracture Behavior of Ultrathin Polymeric Films

Understanding fracture mechanics of ultrathin polymeric films is crucial for modern technologies, including semiconductor and coating industries. However, up to now, the fracture behavior of sub-100 nm polymeric thin films is rarely explored due to challenges in handling samples and limited testing methods available. In this work, we report a new testing methodology that can not only visualize the evolution of the local stress distribution through wrinkling patterns and crack propagation during the deformation of ultrathin films but also directly measure their fracture energies. Using ultrathin polystyrene films as a model system, we both experimentally and computationally investigate the effect of the film thickness and molecular weight on their fracture behavior, both of which show a ductile-to-brittle transition. Furthermore, we demonstrate the broad applicability of this testing method in semicrystalline semiconducting polymers. We anticipate our methodology described here could provide new ways of studying the fracture behavior of ultrathin films under confinement.

Molecular Origin of Electrical Conductivity in Gold-Polythiophene Hybrid Particle Films

Hybrid electronic materials combine inorganic metals and semiconductors with π-conjugated polymers. The orientation of the polymer molecules in relation to the inorganic components is crucial for electrical material properties and device performance, but little is known of the configuration of π-conjugated polymers that bind to inorganic surfaces. Highly curved surfaces are common when using nanoscale components, for example, metal nanocrystal cores covered with conductive polymers. It is important to understand their effect on molecular arrangement. Here, we compare the molecular structures and electrical conductivities of well-defined nanoscale gold spheres and rods with shells of the covalently bound polythiophene PTEBS (poly[2-(3-thienyl)-ethyloxy-4-butylsulfonate]). We prepared aqueous sinter-free inks from the particles and printed them. The particles formed highly conductive films immediately after drying. Films with spherical metal cores consistently had 40% lower conductivities than films based on nanorods. Raman and X-ray photoelectron spectroscopy revealed differences in the gold-sulfur bonds of PTEBS on rods and spheres. The fractions of bond sulfur groups implied differences in the alignment of PTEBS with the surface. More polymer molecules were bound in an edge-on configuration on spheres than on rods, where almost all polymers aligned "face-on" with the metal surface. This leads to different interface resistances: gold-polythiophene-gold interfaces between rods with π-π-tacked face-on PTEBS apparently foster electron transport along the surface-normal direction, while edge-on PTEBS does not. Molecular confinement thus increases the conductivity of hybrid inks based on highly curved nanostructures.

Hole-Transport Layer Molecular Weight and Doping Effects on Perovskite Solar Cell Efficiency and Mechanical Behavior

The effect of tuning molecular weight ( M_n) in poly[bis(4-phenyl)(2,4,6-trimethylphenyl)amine] (PTAA) to increase both mechanical properties of the film and electrical properties of perovskite solar cells is reported. Perovskite solar cell devices are fabricated to investigate the effect of M_n on power conversion efficiency. Moisture stability for various M_n is also studied in PTAA films exposed to mechanical loads in humid environments. Furthermore, cohesion and tensile tests are employed to determine the mechanical properties of PTAA, where higher M_n leads to more robust films. To elucidate the effect of M_n on the debonding kinetics, a viscoelastic fracture kinetic model is proposed as a function of M_n, and the debonding mechanism is found to be dependent on M_n. Finally, the effect of small-molecule-based dopants on the mechanical stability of PTAA is investigated.

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.

Enhancement of the Mechanical Properties of Polyimide Film by Microwave Irradiation

Polyimide films have conventionally been prepared by thermal imidization of poly(amic acid)s (PAAs). Here we report that the improvement of tensile strength while increasing (or maintaining) film flexibility of polyimide films was accomplished by simple microwave (MW) irradiation of the PAAs. This improvement in mechanical properties can be attributed to the increase in molecular weight of the polyimides by MW irradiation. Our results show that the mechanical properties of polyimide films can be improved by MW irradiation, which is a green approach that requires relatively low MW power, very short irradiation time, and no incorporation of any additional inorganic substance.

Comparison of Methods for Determining the Mechanical Properties of Semiconducting Polymer Films for Stretchable Electronics

This paper describes a comparison of two characterization techniques for determining the mechanical properties of thin-film organic semiconductors for applications in soft electronics. In the first method, the film is supported by water (film-on-water, FOW), and a stress-strain curve is obtained using a direct tensile test. In the second method, the film is supported by an elastomer (film-on-elastomer, FOE), and is subjected to three tests to reconstruct the key features of the stress-strain curve: the buckling test (tensile modulus), the onset of buckling (yield point), and the crack-onset strain (strain at fracture). The specimens used for the comparison are four poly(3-hexylthiophene) (P3HT) samples of increasing molecular weight (M_n = 15, 40, 63, and 80 kDa). The methods produced qualitatively similar results for mechanical properties including the tensile modulus, the yield point, and the strain at fracture. The agreement was not quantitative because of differences in mode of loading (tension vs compression), strain rate, and processing between the two methods. Experimental results are corroborated by coarse-grained molecular dynamics simulations, which lead to the conclusion that in low molecular weight samples (M_n = 15 kDa), fracture occurs by chain pullout. Conversely, in high molecular weight samples (M_n > 25 kDa), entanglements concentrate the stress to few chains; this concentration is consistent with chain scission as the dominant mode of fracture. Our results provide a basis for comparing mechanical properties that have been measured by these two techniques, and provide mechanistic insight into fracture modes in this class of materials.

Bulk Heterojunction Morphologies with Atomistic Resolution from Coarse-Grain Solvent Evaporation Simulations

Control over the morphology of the active layer of bulk heterojunction (BHJ) organic solar cells is paramount to achieve high-efficiency devices. However, no method currently available can predict morphologies for a novel donor–acceptor blend. An approach which allows reaching relevant length scales, retaining chemical specificity, and mimicking experimental fabrication conditions, and which is suited for high-throughput schemes has been proven challenging to find. Here, we propose a method to generate atom-resolved morphologies of BHJs which conforms to these requirements. Coarse-grain (CG) molecular dynamics simulations are employed to simulate the large-scale morphological organization during solution-processing. The use of CG models which retain chemical specificity translates into a direct path to the rational design of donor and acceptor compounds which differ only slightly in chemical nature. Finally, the direct retrieval of fully atomistic detail is possible through backmapping, opening the way for improved quantum mechanical calculations addressing the charge separation mechanism. The method is illustrated for the poly(3-hexyl-thiophene) (P3HT)–phenyl-C61-butyric acid methyl ester (PCBM) mixture, and found to predict morphologies in agreement with experimental data. The effect of drying rate, P3HT molecular weight, and thermal annealing are investigated extensively, resulting in trends mimicking experimental findings. The proposed methodology can help reduce the parameter space which has to be explored before obtaining optimal morphologies not only for BHJ solar cells but also for any other solution-processed soft matter 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.