Polycyanurates via Molecular Dynamics: In Situ Crosslinking from Di(Cyanate Ester) Resins and Model Validation through Comparison to Experiment

Optimal Molecular Dynamics System Size for Increased Precision and Efficiency for Epoxy Materials

Molecular dynamics (MD) simulation is an important tool for predicting thermo-mechanical properties of polymer resins at the nanometer length scale, which is particularly important for efficient computationally driven design of advanced composite materials and structures. Because of the statistical nature of modeling amorphous materials on the nanometer length scale, multiple MD models (replicates) are typically built and simulated for statistical sampling of predicted properties. Larger replicates generally provide higher precision in the predictions but result in higher simulation times. Unfortunately, there is insufficient information in the literature to establish guidelines between MD model size and the resulting precision in predicted thermo-mechanical properties. The objective of this study was to determine the optimal MD model size of epoxy resin to balance efficiency and precision. The results show that an MD model size of 15,000 atoms provides for the fastest simulations without sacrificing precision in the prediction of mass density, elastic properties, strength, and thermal properties of epoxy. The results of this study are important for efficient computational process modeling and integrated computational materials engineering (ICME) for the design of next-generation composite materials for demanding applications.

Efficient exploration of compositional space for high-performance copolymers via Bayesian optimization

The traditional approach employed in copolymer compositional design, which relies on trial-and-error, faces low-efficiency and high-cost obstacles when attempting to simultaneously improve multiple conflicting properties. For example, designing co-cured polycyanurates that exhibit both moisture and thermal resistance, along with high modulus, is a long-term challenge because of the intrinsic trade-offs between these properties. In this work, to surmount these barriers, we developed a Bayesian optimization (BO)-guided method to expedite the discovery of co-cured polycyanurates exhibiting low water uptake, coupled with higher glass transition temperature and Young's modulus. By virtue of the knowledge of molecular simulations, benchmarking studies were carried out to develop an effective BO-guided method. Propelled by the developed method, several copolymers with improved comprehensive properties were obtained experimentally in a few iterations. This work provides guidance for efficiently designing other high-performance copolymers.

Directly Using Paraffin as the Toughening Agent of Epoxy Composites: An Experimental and Molecular Dynamics Simulation Study

It is still a challenge in studying the toughening mechanism by well combining the experimental and atomistic molecular dynamics (MD) simulation study. This article directly introduced eicosane (C20, model compound of paraffin) into the epoxy matrix (DGEBA) by using a special epoxy resin with alkyl side chains (D12) as a compatibilizer, which was synthesized through thiol-ene click chemistry. The toughening mechanism of the ternary DGEBA/D12/C20 (EPDA-X) systems was systematically investigated by experimental and MD simulation methods. Though C20 can be well dispersed in the curing mixture, the huge polarity difference between C20 and DGEBA can be the driving force for C20 to stay away from DGEBA, demonstrating the self-assembly effect of C20 around the alkyl side chains of D12 because of the good compatibility of D12 and C20. The soft alkyl chains of D12 and C20 as well as the self-assembly effect of C20 around the D12 molecules can simultaneously improve the strength, modulus, and toughness of the EPDA-2.5 system. This article not only provides a brand new toughening strategy by directly using nonfunctional alkyl derivatives as the toughening agent of epoxy composites with superior mechanical properties but also provides a systematic MD simulation method to evaluate whether there is the interaction or not and the strength of interaction between different molecular chains so as to provide a theoretical basis for the cause of the microphase separation structure and related toughening mechanism in cross-linking networks on the atomic and molecular levels.

Conformational Heterogeneity and Interchain Percolation Revealed in an Amorphous Conjugated Polymer

Conjugated polymers are employed in a variety of application areas due to their bright fluorescence and strong biocompatibility. However, understanding the structure of amorphous conjugated polymers on the nanoscale is extremely challenging compared to their related crystalline phases. Using a bespoke classical force field, we study amorphous poly(9,9-di-n-octylfluorene-alt-benzothiadiazole) (F8BT) with molecular dynamics simulations to investigate the role that its nanoscale structure plays in controlling its emergent (and all-important) optical properties. Notably, we show that a giant percolating cluster exists within amorphous F8BT, which has ramifications in understanding the nature of interchain species that drive the quantum yield reduction and bathochromic shift observed in conjugated polymer-based devices and nanostructures. We also show that distinct conformations can be unravelled from within the disordered structure of amorphous F8BT using a two-stage machine learning protocol, highlighting a link between molecular conformation and ring stacking propensity. This work provides predictive understanding by which to enhance the optical properties of next-generation conjugated polymer-based devices and materials by rational, simulation-led design principles.

A Molecular Dynamics Study of Monomer Melt Properties of Cyanate Ester Monomer Melt Properties

The objective of this work was to computationally predict the melting temperature and melt properties of thermosetting monomers used in aerospace applications. In this study, we applied an existing voids method by Solca. to examine four cyanate ester monomers with a wide range of melting temperatures. Voids were introduced into some simulations by removal of molecules from lattice positions to lower the free-energy barrier to melting to directly simulate the transition from a stable crystal to amorphous solid and capture the melting temperature. We validated model predictions by comparing melting temperature against previously reported literature values. Additionally, the torsion and orientational order parameters were used to examine the monomers’ freedom of motion to investigate structure–property relationships. Ultimately, the voids method provided reasonable estimates of melting temperature while the torsion and order parameter analysis provided insight into sources of the differing melt properties between the thermosetting monomers. As a whole, the results shed light on how freedom of molecular motions in the monomer melt state may affect melting temperature and can be utilized to inspire the development of thermosetting monomers with optimal monomer melt properties for demanding applications.