Predicting the Mechanical Properties of Organic Semiconductors Using Coarse-Grained Molecular Dynamics Simulations

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro-Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density-temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

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.

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.

Temperature-Transferable Coarse-Grained Models for Volumetric Properties of Poly(lactic Acid)

A new coarse-grained (CG) model, for which each monomer is mapped as one bead at its center of mass, was developed for simulating the volumetric properties of the polylactide (PLA) bulk. The three bonded CG potentials are first parametrized against the strain energies of the dimer, trimer, and tetramer, and the nonbonded CG potentials are then optimized to match the melt densities of the decamer. With the derived CG potentials, molecular dynamics (MD) simulations are found to reproduce thermal expansion and glass transition. By rescaling the dihedral and nonbonded potentials with temperature-independent factors, the glass transition temperature (Tg) is also satisfactorily restored with little modifications on the volumetric expansive coefficients at both rubbery and glassy states. Therefore, the finally optimized CG potentials exhibit excellent temperature transferability, as rationalized by the Simha-Boyer relation. Furthermore, it is confirmed that the dihedral torsions and nonbonded interactions play key roles in glass transition. Also, the simulated bulk moduli and conformational properties in a wide temperature range compare well with the referenced data. The proposed multiscale scheme has great potential in simulating thermo-mechanical properties of PLA and other polymers.

Understanding the thermomechanical behavior of graphene-reinforced conjugated polymer nanocomposites via coarse-grained modeling

Graphene-reinforced conjugated polymer (CP) nanocomposites are attractive for flexible and electronic devices, but their mechanical properties have been less explored at a fundamental level. Here, we present a predictive multiscale modeling framework for graphene-reinforced poly(3-alkylthiophene) (P3AT) nanocomposites via atomistically informed coarse-grained molecular dynamics simulations to investigate temperature-dependent thermomechanical properties at a molecular level. Our results reveal reduced graphene dispersion with increasing graphene loading. Nanocomposites with shorter P3AT side chains, lower temperatures, and higher graphene content exhibit stronger mechanical responses, which correlates with polymer dynamics. The elastic modulus increases linearly with the graphene content, which slightly deviates from the "Halpin-Tsai" micromechanical model prediction. Local stiffness analysis shows that graphene possesses the highest stiffness, followed by the P3AT backbone and side chains. Deformation-induced stronger chain alignment of the P3AT backbone compared to graphene may further promote conductive behavior. Our findings provide insights into the dynamical heterogeneity of nanocomposites, paving the way for understanding and predicting their thermomechanical properties.

Enzymatically Polymerized Organic Conductors on Model Lipid Membranes

Seamless integration between biological systems and electrical components is essential for enabling a twinned biochemical-electrical recording and therapy approach to understand and combat neurological disorders. Employing bioelectronic systems made up of conjugated polymers, which have an innate ability to transport both electronic and ionic charges, provides the possibility of such integration. In particular, translating enzymatically polymerized conductive wires, recently demonstrated in plants and simple organism systems, into mammalian models, is of particular interest for the development of next-generation devices that can monitor and modulate neural signals. As a first step toward achieving this goal, enzyme-mediated polymerization of two thiophene-based monomers is demonstrated on a synthetic lipid bilayer supported on a Au surface. Microgravimetric studies of conducting films polymerized in situ provide insights into their interactions with a lipid bilayer model that mimics the cell membrane. Moreover, the resulting electrical and viscoelastic properties of these self-organizing conducting polymers suggest their potential as materials to form the basis for novel approaches to in vivo neural therapeutics.

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.

Technology Roadmap for Flexible Sensors

Humans rely increasingly on sensors to address grand challenges and to improve quality of life in the era of digitalization and big data. For ubiquitous sensing, flexible sensors are developed to overcome the limitations of conventional rigid counterparts. Despite rapid advancement in bench-side research over the last decade, the market adoption of flexible sensors remains limited. To ease and to expedite their deployment, here, we identify bottlenecks hindering the maturation of flexible sensors and propose promising solutions. We first analyze challenges in achieving satisfactory sensing performance for real-world applications and then summarize issues in compatible sensor-biology interfaces, followed by brief discussions on powering and connecting sensor networks. Issues en route to commercialization and for sustainable growth of the sector are also analyzed, highlighting environmental concerns and emphasizing nontechnical issues such as business, regulatory, and ethical considerations. Additionally, we look at future intelligent flexible sensors. In proposing a comprehensive roadmap, we hope to steer research efforts towards common goals and to guide coordinated development strategies from disparate communities. Through such collaborative efforts, scientific breakthroughs can be made sooner and capitalized for the betterment of humanity.

Fast, Faithful Simulations of Donor-Acceptor Interface Morphology

The local structure of conjugated polymers governs key optoelectronic properties, such as charge conduction and photogeneration at donor-acceptor interfaces. Because conjugated polymers are large, stiff, and relax slowly, all-atom molecular dynamics simulations are computationally expensive. Here, we describe a coarse-graining method that exploits the stiffness of constituent aromatic moieties by representing each moiety as rigidly bonded clusters of atoms wherein virtual sites replace several atoms. This approach significantly reduces the degrees of freedom while faithfully representing the shape and interactions of the moieties, resulting in 10 times faster simulations than all-atom simulations. Simulation of a donor polymer (P3HT) and a non-fullerene acceptor (O-IDTBR) validates the coarse-graining method by comparing structural properties from experiments, such as the density and persistence length. The fast simulation produces equilibrated systems with realistic morphologies. The simulation results of an equimolar mixture of P3HT, with a molecular weight of 1332 g mol^-1, and an O-IDTBR mixture suggest that the interface width must be larger than 7 nm. Also, we investigate the effect of slow cooling on morphologies, particularly the number of close contacts that facilitates carrier transport. Slow cooling increases close contacts, and the effect is more pronounced in crystal-forming P3HT than in O-IDTBR, where bulky side-groups hinder crystal formation.

Light-responsive self-strained organic semiconductor for large flexible OFET sensing array

With the wide application of organic semiconductors (OSCs), researchers are now grappling with a new challenge: design and synthesize OSCs materials with specific functions to satisfy the requirements of high-performance semiconductor devices. Strain engineering is an effective method to improve the semiconductor material’s carrier mobility, which is fundamentally originated from the rearrangement of the atomic packing model of materials under mechanic stress. Here, we design and synthesize a new OSC material named AZO-BTBT-8 based on high-mobility benzo[b]benzo[4,5]thieno[2,3-d]thiophene (BTBT) as the semiconductor backbone. Octane is employed to increase molecular flexibility and solubility, and azobenzene at the other end of the BTBT backbone provides photoisomerization properties and structural balance. Notably, the AZO-BTBT-8 photoisomerization leads to lattice strain in thin-film devices, where exceptional device performance enhancement is realized. On this basis, a large-scale flexible organic field-effect transistor (OFET) device array is fabricated and realizes high-resolution UV imaging with reversible light response.

Strain engineering is effective to improve the carrier mobility of semiconductor materials. Here, the authors demonstrate lattice strain-induced mobility enhancement of an azobenzene compound under photoisomerization and its application in large-scale flexible organic field-effect transistors.

Multiscale modeling of thermomechanical properties of stereoregular polymers

Multiscale coarse-grained (CG) models are expected to play the critical roles in molecular simulations of complex polymers. However, this poses a great challenge for accurately simulating their thermomechanical properties, for which excellent representability and transferability are required for the CG potentials. In this work, virtual sites and elastic network bonds are introduced to improve the structural and volumetric property-based CG models including explicit electrostatic interactions, which is exemplarily applied to the iso- and syndio-tactic poly(methyl methacrylate). A variety of thermomechanical properties of the two stereoregular polymer bulks are reasonably reproduced by the extensive molecular dynamics simulations with the so-parameterized CG potentials. In particular, the attractive nature of electrostatic interactions and tacticity effects on glass transition temperatures (T_g) are well captured. Furthermore, stronger electrostatic interactions lead to higher mass density and bulk modulus, and their effects on Young's modulus, Poisson's ratio, and shear modulus depend upon the chain tacticity. It is also demonstrated that all these elastic constants can be effectively modulated by imposing external electric field. The proposed multiscale scheme can be very valuable to molecular designs of polar polymer materials.

Thermomechanical Behavior of Poly(3-hexylthiophene) Thin Films on the Water Surface

The thermomechanical behavior of a conjugated polymer (CP) in a thin film state has rarely been studied despite the importance of understanding the polymer morphologies and optimizing the thermal processes of organic semiconductors. Moreover, the seamless integration of multilayers without mechanical failures in CP-based electronic devices is crucial for determining their operational stability. Large differences in the coefficients of thermal expansion (CTEs) between the multilayers can cause serious degradation of devices under thermal stress. In this study, we measure the intrinsic thermomechanical properties of poly(3-hexylthiophene) (P3HT) thin films in a pseudo-freestanding state on the water surface. The as-cast P3HT thin films exhibited a large thermal shrinkage (-1001 ppm K^-1) during heating on the water surface. Morphological analyses revealed that the thermal shrinkage of the polymer films was caused by the rearrangement of the polymer chain networks accompanied by crystallization, thus indicating that preheating the polymer films is essential for estimating their intrinsic CTE values. Moreover, the rigidity of the substrate significantly influences the thermomechanical behavior of the polymer films. The polymer films that were preheated on the glass substrate showed nonlinear thermal expansion due to the substrate constraint inhibiting sufficient relaxation of the polymer chains. In comparison, a linear expansion behavior is observed after preheating the films on the water surface, exhibiting a consistent CTE value (185 ppm K^-1) regardless of the number of thermal strain measurements. Thus, this work provides a direct method for measuring in-plane CTE values and an in-depth understanding of the thermomechanical behaviors of CP thin films to design thermomechanically reliable organic semiconductors.

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.

Transfer Learned Designer Polymers For Organic Solar Cells

Organic photovoltaic (OPV) materials have been examined extensively over the past two decades for solar cell applications because of the potential for device flexibility, low-temperature solution processability, and negligible environmental impact. However, discovery of new candidate OPV materials, especially polymer-based electron donors, that demonstrate notable power conversion efficiencies (PCEs), is nontrivial and time-intensive exercise given the extensive set of possible chemistries. Recent progress in machine learning accelerated materials discovery has facilitated to address this challenge, with molecular line representations, such as Simplified Molecular-Input Line-Entry Systems (SMILES), gaining popularity as molecular fingerprints describing the donor chemical structures. Here, we employ a transfer learning based recurrent neural (LSTM) model, which harnesses the SMILES molecular fingerprints as an input to generate novel designer chemistries for OPV devices. The generative model, perfected on a small focused OPV data set, predicts new polymer repeat units with potentially high PCE. Calculations of the similarity coefficient between the known and the generated polymers corroborate the accuracy of the model predictability as a function of the underlying chemical specificity. The data-enabled framework is sufficiently generic for use in accelerated machine learned materials discovery for various chemistries and applications, mining the hitherto available experimental and computational data.

Elasto-morphology of P3HT:PCBM bulk heterojunction organic solar cells

Predicting the mechanical properties of organic semiconductors is important when using these materials in flexible electronics applications. For instance, knowledge of the mechanical and thermal stability of thin film organic solar cells (OSCs) is critical for the roll-to-roll production of photovoltaic devices and their use under various operating conditions. Here, we examine the thermal and elasto-mechanical properties of the conjugated donor polymer poly-(3-hexylthiophene) (P3HT) and the interpenetrating mixtures of P3HT and phenyl-C₆₁-butyric acid methyl (PCBM) ester bulk heterojunction (BHJ) active layers under the application of unidirectional tensile deformation using coarse-grained molecular dynamics (CGMD) simulations. The predictions are validated against previous experimental reports as well as with earlier modeling results derived using different intermolecular force fields. Our results reveal that PCBM molecules behave as anti-plasticizers when mixed with P3HT and tend to increase the tensile modulus and glass transition temperature, while decreasing the crack-onset strain relative to pure P3HT. The variations in the mechanical properties with the composition of the BHJ active layer suggest that, in the presence of small oligomers as additives in the BHJ, the P3HT:PCBM mixture resists the anti-plasticizing effect of PCBM molecules due to the low tensile modulus of the short polymer chains.

Investigation on mechanical behaviors of Cu-Ni binary alloy nanopillars: a molecular dynamics study

Cu-Ni binary alloy has become the attention of scientific world for its potentials in nanodevices. It is indispensable to investigate on the mechanical properties of this material due to lack of previous work done regarding this binary alloy. Molecular dynamics (MD) studies were performed on nanopillar (NP) structures comprised of Cu-Ni binary alloy having an FCC unit cell with Cu atoms selectively replaced by Ni atoms. This selective replacement resulted in a better stress behavior than the randomly replaced alloy structure when both tension and compression load were applied. The effect of crystal orientation, NP dimensions, temperature, and strain rate on the stress-strain curve of Cu-Ni binary alloy NPs was thoroughly investigated under tensile loading. This investigation reveals significant influence of crystal orientation on ultimate strength and flow stress region. Among four different crystal orientations, <111> orientation shows maximum strength behavior under tensile loading, while <110> shows highest strength under compressive load. However, in both cases, i.e. tension and compression, the poorest stress behavior was observed for <001> orientation. Under tensile load, <111>-oriented binary alloy fails due to the formation of Shockley partials followed by formation of complex dislocation network. On the other hand, <110>-oriented binary alloy fails due to the formation of Lomer-Cottrell (LC) lock from the Shockley partials. Total dislocation length is calculated, and its effect on the stress-strain behavior of the Cu-Ni binary alloy is discussed. Highest Young's modulus and yield stress are observed on <111>-oriented binary alloy among other orientations, and these values for <111>-oriented NP was found to decrease with the increment of temperature. If the temperature is increased, yield stress and Young's modulus decrease. The effect of cross section width was also investigated in this study, and it was found that yield stress decreases with the increment of cross section width due to the effect of surface atom fraction. Increasing the strain rate causes the initiation of amorphous structure, resulting in superplastic behavior of the <111>-oriented Cu-Ni binary alloy NP.

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.

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.

Measurement of Cohesion and Adhesion of Semiconducting Polymers by Scratch Testing: Effect of Side-Chain Length and Degree of Polymerization

Most advantages of organic electronic materials are enabled by mechanical deformability, as flexible (and stretchable) devices made from these materials must be able to withstand roll-to-roll printing and survive mechanical insults from the external environment. Cohesion and adhesion are two properties that dictate the mechanical reliability of a flexible organic electronic device. In this paper, progressive-load scratch tests are used for the first time to correlate the cohesive and adhesive behavior of poly(3-alkylthiophenes) (P3ATs) with respect to two molecular parameters: length of the alkyl side chain and molecular weight. In contrast to metrological techniques based on buckling or pull testing of pseudofreestanding films, scratch tests reveal information about both the cohesive and adhesive properties of thin polymeric films from a single procedure. Our data show a decrease in cohesion and adhesion, that is, a decrease in overall mechanical robustness, with increasing length of the side chain. This behavior is likely due to increases in free volume and concomitant decreases in the glass transition temperature. In contrast, we observe increases in both the cohesion and adhesion with increasing molecular weight. This behavior is attributed to an increased density of entanglements with high molecular weight, which manifests as increased extensibility. These observations are consistent with the results of molecular dynamics simulations. Interestingly, the normal (applied) forces associated with cohesive and adhesive failure are directly proportional to the average degree of polymerization, as opposed to simply the molecular weight, as the length of the alkyl side chain increases the molecular weight without increasing the degree of polymerization.

Coarse-Grained Simulations Using a Multipolar Force Field Model

This paper presents a coarse-grained molecular simulation for fullerenes based on a multipolar expansion method developed previously. The method is enabled by the construction of transferable united atoms potentials that approximate the full atomistic intermolecular interactions, as obtained from ab initio electronic structure calculations supplemented by empirical force fields and experimental data, or any combination of the above. The resultant series contains controllable moment tensors that allow to estimate the errors, and approaches the all-atom intermolecular potential as the expansion order increases. We can compute the united atoms potentials very efficiently with a few interaction moment tensors, in order to implement a parallel algorithm on molecular interactions. Our simulations describe the mechanism for the condensation of fullerenes, and they produce excellent agreement with benchmark fully atomistic molecular dynamics simulations.

Atomistic investigations on the mechanical properties and fracture mechanisms of indium phosphide nanowires

The mechanical properties of indium phosphide (InP) nanowires are an emerging issue due to the promising applications of these nanowires in nanoelectromechanical and microelectromechanical devices. In this study, molecular dynamics simulations of zincblende (ZB) and wurtzite (WZ) crystal structured InP nanowires (NWs) are presented under uniaxial tension at varying sizes and temperatures. It is observed that the tensile strengths of both types of NWs show inverse relationships with temperature, but are independent of the size of the nanowires. Moreover, applied load causes brittle fracture by nucleating cleavage on ZB and WZ NWs. When the tensile load is applied along the [001] direction, the direction of the cleavage planes of ZB NWs changes with temperature. It is found that the {111} planes are the cleavage planes at lower temperatures; on the other hand, the {110} cleavage planes are activated at elevated temperatures. In the case of WZ NWs, fracture of the material is observed to occur by cleaving along the (0001) plane irrespective of temperature when the tensile load is applied along the [0001] direction. Furthermore, the WZ NWs of InP show considerably higher strength than their ZB counterparts. Finally, the impact of strain rate on the failure behavior of InP NWs is also studied, and higher fracture strengths and strains at higher strain rates are found. With increasing strain rate, the number of cleavages also increases in the NWs. This paper also provides in-depth understanding of the failure behavior of InP NWs, which will aid the design of efficient InP NWs-based devices.

Influence of amide-containing side chains on the mechanical properties of diketopyrrolopyrrole-based polymers

An efficient strategy to modify the mechanical properties of conjugated polymers has been developed through the incorporation of amide moieties.

Stretchable and Soft Electronics using Liquid Metals

The use of liquid metals based on gallium for soft and stretchable electronics is discussed. This emerging class of electronics is motivated, in part, by the new opportunities that arise from devices that have mechanical properties similar to those encountered in the human experience, such as skin, tissue, textiles, and clothing. These types of electronics (e.g., wearable or implantable electronics, sensors for soft robotics, e-skin) must operate during deformation. Liquid metals are compelling materials for these applications because, in principle, they are infinitely deformable while retaining metallic conductivity. Liquid metals have been used for stretchable wires and interconnects, reconfigurable antennas, soft sensors, self-healing circuits, and conformal electrodes. In contrast to Hg, liquid metals based on gallium have low toxicity and essentially no vapor pressure and are therefore considered safe to handle. Whereas most liquids bead up to minimize surface energy, the presence of a surface oxide on these metals makes it possible to pattern them into useful shapes using a variety of techniques, including fluidic injection and 3D printing. In addition to forming excellent conductors, these metals can be used actively to form memory devices, sensors, and diodes that are completely built from soft materials. The properties of these materials, their applications within soft and stretchable electronics, and future opportunities and challenges are considered.

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.

Stretchable Organic Semiconductor Devices

Stretchable electronics are essential for the development of intensely packed collapsible and portable electronics, wearable electronics, epidermal and bioimplanted electronics, 3D surface compliable devices, bionics, prosthesis, and robotics. However, most stretchable devices are currently based on inorganic electronics, whose high cost of fabrication and limited processing area make it difficult to produce inexpensive, large-area devices. Therefore, organic stretchable electronics are highly attractive due to many advantages over their inorganic counterparts, such as their light weight, flexibility, low cost and large-area solution-processing, the reproducible semiconductor resources, and the easy tuning of their properties via molecular tailoring. Among them, stretchable organic semiconductor devices have become a hot and fast-growing research field, in which great advances have been made in recent years. These fantastic advances are summarized here, focusing on stretchable organic field-effect transistors, light-emitting devices, solar cells, and memory devices.