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.

Structure-Property Relationships of Granular Hybrid Hydrogels Formed through Polyelectrolyte Complexation.

Hybrid hydrogels are hydrogels that exhibit heterogeneity in the network architecture by means of chemical composition and/or microstructure. The different types of interactions, together with structural heterogeneity, which can be created on different length scales, determine the mechanical properties of the final material to a large extent. In this work, the microstructure-mechanical property relationships for a hybrid hydrogel that contains both electrostatic and covalent interactions are investigated. The hybrid hydrogel is composed of a microphase-separated polyelectrolyte complex network (PEC) made of poly(4-styrenesulfonate) and poly(diallyldimethylammonium chloride) within a soft and elastic polyacrylamide hydrogel network. The system exhibits a granular structure, which is attributed to the liquid-liquid phase separation into complex coacervate droplets induced by the polymerization and the subsequent crowding effect of the polyacrylamide chains. The coacervate droplets are further hardened into PEC granules upon desalting the hydrogel. The structure formation is confirmed by a combination of electron microscopic imaging and molecular dynamics simulations. The interpenetration of both networks is shown to enhance the toughness of the resulting hydrogels due to the dissipative behavior of the PEC through the rupture of electrostatic interactions. Upon cyclic loading-unloading, the hydrogels show recovery of up to 80% of their original dissipative behavior in less than 300 s of rest with limited plasticity. The granular architecture and the tough and self-recoverable properties of the designed hybrid networks make them good candidates for applications, such as shape-memory materials, actuators, biological tissue mimics, and elastic substrates for soft sensors.

Bidispersity Improves the Toughness and Impact Resistance of Star-Polymer Thin Films.

Branched polymer architectures are used to tune the mechanical properties of impact-resistant thin films through parameters, such as chain length and grafting density. While chain dispersity affects molecular properties, such as interpenetration and entanglements, structure-property relationships accounting for dispersity are challenging to obtain experimentally and are often neglected in computational models. We employ molecular dynamics simulations to model the high-rate tensile elongation and nanoballistic impact of thin films composed of bidisperse star polymers with varying arm lengths. We find that, at fixed molecular weight, high dispersity can significantly enhance the toughness and impact resistance of the films without decreasing their elastic modulus. Bidisperse stars with fewer longer arms are less entangled, but stretch and interpenetrate for longer times during crazing, leading to increased toughness. These findings highlight controlled dispersity as a design strategy to improve the mechanical properties of polymer composites across Pareto fronts.

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.

The effect of nanoparticle softness on the interfacial dynamics of a model polymer nanocomposite.

The introduction of soft organic nanoparticles (NPs) into polymer melts has recently expanded the material design space for polymer nanocomposites, compared to traditional nanocomposites that utilize rigid NPs, such as silica, metallic NPs, and other inorganic NPs. Despite advances in the fabrication and characterization of this new class of materials, the effect of NP stiffness on the polymer structure and dynamics has not been systematically investigated. Here, we use molecular dynamics to investigate the segmental dynamics of the polymer interfacial region of isolated NPs of variable stiffness in a polymer matrix. When the NP-polymer interactions are stronger than the polymer-polymer interactions, we find that the slowing of segmental dynamics in the interfacial region is more pronounced for stiff NPs. In contrast, when the NP-polymer interaction strength is smaller than the matrix interaction, the NP stiffness has relatively little impact on the changes in the polymer interfacial dynamics. We also find that the segmental relaxation time τα of segments in the NP interfacial region changes from values lower than to higher than the bulk material when the NP-polymer interaction strength is increased beyond a "critical" strength, reminiscent of a binding-unbinding transition. Both the NP stiffness and the polymer-surface interaction strength can thus greatly influence the relative segmental relaxation and interfacial mobility in comparison to the bulk material.

Systematic Coarse-graining of Epoxy Resins with Machine Learning-Informed Energy Renormalization.

A persistent challenge in predictive molecular modeling of thermoset polymers is to capture the effects of chemical composition and degree of crosslinking (DC) on dynamical and mechanical properties with high computational efficiency. We established a new coarse-graining (CG) approach that combines the energy renormalization method with Gaussian process surrogate models of the molecular dynamics simulations. This allows a machine-learning informed functional calibration of DC-dependent CG force field parameters. Taking versatile epoxy resins consisting of Bisphenol A diglycidyl ether combined with curing agent of either 4,4-Diaminodicyclohexylmethane or polyoxypropylene diamines, we demonstrated excellent agreement between all-atom and CG predictions for density, Debye-Waller factor, Young's modulus and yield stress at any DC. We further introduce a surrogate model enabled simplification of the functional forms of 14 non-bonded calibration parameters by quantifying the uncertainty of a candidate set of high-dimensional/flexible calibration functions. The framework established provides an efficient methodology for chemistry-specific, large-scale investigations of the dynamics and mechanics of epoxy resins.

Predictive relation for the α-relaxation time of a coarse-grained polymer melt under steady shear.

We examine the influence of steady shear on structural relaxation in a simulated coarse-grained unentangled polymer melt over a wide range of temperature and shear rates. Shear is found to progressively suppress the α-relaxation process observed in the intermediate scattering function, leading ultimately to a purely inertially dominated ß-relaxation at high shear rates, a trend similar to increasing temperature. On the basis of a scaling argument emphasizing dynamic heterogeneity in cooled liquids and its alteration under material deformation, we deduce and validate a parameter-free scaling relation for both the structural relaxation time τα from the intermediate scattering function and the "stretching exponent" ß quantifying the extent of dynamic heterogeneity over the entire range of temperatures and shear rates that we can simulate.

Influence of polymer topology on crystallization in thin films.

We investigate how varying molecular topology of polymers influences crystallization in thin polymer films. In particular, we simulate linear and star polymers of fixed mass having a progressively increasing number of arms (f ≤ 16) in a system where the linear polymer exhibits crystallization in a thin film geometry, but no apparent crystallization in the corresponding bulk material. The degree of crystallization of the polymer film at long times decreases progressively with increasing f, and no crystallization is observed beyond f = 8. Crystallization for smaller values of f develops as a sigmoidally shaped wavefront initiating from the supporting crystalline interface. We suggest that large shape fluctuations and the competition of length scales of star polymers with high f lead to inhibited crystallization.

Mutual information does not detect growing correlations in the propensity of a model molecular liquid.

The dynamical spatial correlations detected by the mutual information (MI) in the isoconfigurational particle displacements of a monodisperse molecular viscous liquid are studied via molecular-dynamics simulations by changing considerably both the molecular mobility and the degree of dynamical heterogeneity. Different from atomic liquids, the MI correlation length does not grow on approaching the glass transition by considering the liquid both in full detail as a collection of monomers and as a coarse-grained ensemble of molecular centers of mass. In the detailed picture, it is found that: (i) the MI correlations between monomers are largely due to inter-molecular correlations, (ii) the MI length scale is numerically identical, within the errors, to the correlation length scale of the displacement direction, as drawn by conventional correlation functions. The time evolution of the MI spatial correlations complies with the scaling between the fast vibrational dynamics and the long-time relaxation. Our findings suggest that the characteristics of the MI length scale are markedly system-dependent and not obviously related to dynamical heterogeneity.

Molecular layers in thin supported films exhibit the same scaling as the bulk between slow relaxation and vibrational dynamics.

We perform molecular-dynamics simulations of a supported molecular thin film. By varying thickness and temperature, we observe anisotropic mobility as well as strong gradients of both the vibrational motion and the structural relaxation through film layers with monomer-size thickness. We show that the gradients of the fast and the slow dynamics across the layers (except the adherent layer to the substrate) comply, without any adjustment, with the same scaling between the structural relaxation time and the Debye-Waller factor originally observed in the bulk [Larini et al., Nat. Phys., 2008, 4, 42]. The scaling is not observed if the average dynamics of the film is inspected. Our results suggest that the solidification process of each layer may be tracked by knowing solely the vibrational properties of the layer and the bulk.