Impact of the Force Field on the Calculation of Density and Surface Tension of Epoxy-Resins

The molecular simulation of interfacial systems is a matter of debate because of the choice of many input parameters that can affect significantly the performance of the force field of reproducing the surface tension and the coexisting densities. After developing a robust methodology for the calculation of the surface tension on a Lennard-Jones fluid, we apply it with different force fields to calculate the density and surface tension of pure constituents of epoxy resins. By using the model that best reproduces the experimental density and surface tension, we investigate the impact of composition in mass fraction on uncured epoxy resins and the effects of degree of cross-linking on cured resins.

Fast-Processable Non-Flammable Phthalonitrile-Modified Novolac/Carbon and Glass Fiber Composites

Phthalonitrile resins (PN) are known for their incredible heat resistance and at the same time poor processability. Common curing cycle of the PN includes dozens hours of heating at temperatures up to 375 °C. This work was aimed at reducing processing time of phthalonitrile resin, and with this purpose, a novolac oligomer with hydroxyl groups fully substituted by phthalonitrile moieties was synthesized with a quantitative yield. Formation of the reaction byproducts was investigated depending on the synthesis conditions. The product was characterized by ¹H NMR and FT-IR. Curing of the resins with the addition of different amounts of novolac phenolic as curing agent (25, 50 and 75 wt.%) was studied by rheological and DSC experiments. Based on these data, a curing program was developed for the further thermosets' investigation: hot-pressing at 220 °C and 1.7 MPa for 20 min. TGA showed the highest thermal stability of the resin with 25 wt.% of novolac (T_5% = 430 °C). The post-curing program was developed by the use of DMA with different heating rates and holding for various times at 280 or 300 °C (heating rate 0.5 °C/min). Carbon and glass fiber plastic laminates were fabricated via hot-pressing of prepregs with T_g's above 300 °C. Microcracks were formed in the CFRP, but void-free GFRP were fabricated and demonstrated superior mechanical properties (ILSS up to 86 MPa; compressive strength up to 620 MPa; flexural strength up to 946 MPa). Finally, flammability tests showed that the composite was extinguished in less than 5 s after the flame source was removed, so the material can be classified as V-0 according to the UL94 ratings. For the first time, fast-curing phthalonitrile prepregs were presented. The hot-pressing cycle of 20 min with 150 min free-standing post-curing yielded composites with the unique properties. The combination of mechanical properties, scale-up suitable fast-processing and inflammability makes the presented materials prospective for applications in the electric vehicle industries, fast train construction and the aerospace industry.

Molecular Simulations of Thermomechanical Properties of Epoxy-Amine Resins

All-atom molecular dynamics (MD) simulations were performed with the CHARMM force field to characterize various epoxy resins, such as aliphatic and bisphenol-based resins. A multistep cross-linking algorithm was established, and key properties such as density, glass temperature, and elastic modulus were calculated. A quantitative comparison was made and was proven to be in good agreement with experimental data, with average absolute deviations between experiments and molecular simulation comprised between 2% and 12%. Additional findings on structure-property relationships were highlighted such as the effect of the cross-linking rate and oligomerization of the resin.

Arbitrarily-shaped microgels composed of chemically unmodified biopolymers

Biohydrogels, composed of naturally occurring biopolymers are typically preferred over their synthetic analogues in bioapplications thanks to their biocompatibility, bioactivity, mechanical or degradation properties. Shaping biohydrogels on the single-cell length scales (micrometers) is a key ability needed to create bioequivalent artificial cell/tissue constructs and cannot be achieved with current methods. This work introduces a method for photolithographic synthesis of arbitrarily shaped microgels composed purely of a biopolymer of choice. The biopolymer is mixed with a sacrificial photocrosslinkable polymer, and the mixture is photocrosslinked in a lithographic process, yielding anisotropic microgels with the biopolymer entrapped in the network. Subsequent ionic or covalent biopolymer crosslinking followed by template cleavage yields a microgel composed purely of a biopolymer with the 3D shape dictated by the photocrosslinking process. Method feasibility is demonstrated with two model polysaccharide biopolymers (alginate, chitosan) using suitable crosslinking methods. Next, alginate microgels were used as microtaggants on a pharmaceutical oral solid dose formulation to prevent its counterfeiting. Since the alginate is approved as an additive in the food and pharmaceutical industries, the presented tagging system can be implemented in practical use much easier than systems comprising synthetic polymers.

MARTINI-based simulation method for step-growth polymerization and its analysis by size exclusion characterization: a case study of cross-linked polyurethane

Simulation studies of step-growth polymerization, e.g., polymerization of polyurethane systems, hold great promise due to having complete control over the reaction conditions and being able to perform an in-depth analysis of network structures. In this work, we developed a (completely automated) simulation method based on a coarse-grained (CG) methodology, i.e., the MARTINI model, to study the cross-linking reaction of a diol, a tri-isocyanate molecule and one-hydroxyl functional molecule to form a polyurethane network without and with dangling chains. This method is capable of simulating the cross-linking reactions not only up to very high conversions, but also under rather complicated reaction conditions, i.e., a non-stoichiometric ratio of the reactants, solvent evaporation and multi-step addition of the reactants. We introduced a novel network analysis, similar to size-exclusion chromatography based on graph theory, to study the growth of the network during the polymerization process. By combining the reaction simulations with these analysis methods, a set of correlations between the reaction conditions, reaction mechanisms and final network structure and properties is revealed. For instance, a two-step addition of materials for the reaction, i.e., first the dangling chain to the tri-isocyanate and then the diol, leads to the highest integrated network structure. We observed that different reaction conditions lead to different glass transition temperatures (T_g) of the network due to the distinct differences in the final network structures obtained. For example, by addition of dangling chains to the network, the T_g decreases as compared to the network without dangling chains, as also is commonly observed experimentally.

Interfacial Characteristics of Boron Nitride Nanosheet/Epoxy Resin Nanocomposites: A Molecular Dynamics Simulation

The interface between nanofillers and matrix plays a key role in determining the properties of nanocomposites, but the interfacial characteristics of nanocomposites such as molecular structure and interaction strength are not fully understood yet. In this work, the interfacial features of a typical nanocomposite, namely epoxy resin (EP) filled with boron nitride nanosheet (BNNS) are investigated by utilizing molecular dynamics simulation, and the effect of surface functionalization is analyzed. The radial distribution density (RDD) and interfacial binding energy (IBE) are used to explore the structure and bonding strength of nanocomposites interface. Besides, the interface compatibility and molecular chain mobility (MCM) of BNNS/EP nanocomposites are analyzed by cohesive energy density (CED), free volume fraction (FFV), and radial mean square displacement (RMSD). The results indicate that the interface region of BNNS/EP is composed of three regions including compact region, buffer region, and normal region. The structure at the interfacial region of nanocomposite is more compact, and the chain mobility is significantly lower than that of the EP away from the interface. Moreover, the interfacial interaction strength and compatibility increase with the functional density of BNNS functionalized by CH3–(CH2)4–O– radicals. These results adequately illustrate interfacial characteristics of nanocomposites from atomic level.

Model of the DGEBA-EDA Epoxy Polymer: Experiments and Simulation Using Classical Molecular Dynamics

Polyepoxy samples are synthesized from diglycidylether of bisphenol A (DGEBA) and ethylene diamine (EDA) monomers at a stoichiometric ratio of 2 DGEBA : 1 EDA in model conditions in order to promote a high degree of polymerization and a low density of defects and to try to approach the ideal models obtained by simulation. A slow polymerization (>24 h at ambient temperature) and a postcuring achieved in an inert atmosphere lead to a conversion degree of 92±2% and a midpoint glass transition temperature of 391±1 K. In parallel, a model is created with a multistep cross-linking procedure. In this work, all-atom molecular dynamics (MD) simulations are performed with LAMMPS and the GAFF 1.8 force field. In the initial liquid mixture of reactants (600 molecules), proper mixing is demonstrated by the calculation of the partial radial distribution functions (RDF), which show a minimum intermolecular distance of 2.8 Å and similar distributions for EDA-EDA, DGEBA-DGEBA, and DGEBA-EDA molecules in the simulation boxes. Then, in alternation with MD equilibrations, cross-linking is performed on frozen configurations by creating covalent bonds between reactive pairs within a reaction radius of 3 Å. The resulting boxes show conversion rates of 90-93% and densities close to the experimental value. Finally, a cooling ramp from 700 K to 25 K is applied in order to monitor the specific volume and the coefficient of volumetric thermal expansion (CVTE) of the polymer and to derive the glass transition temperature. Experimental thermomechanical analyses (TMA) compares well with simulations for both the specific volume and the CVTE evolutions with temperature. Whereas the uncertainty remains high with the fitting procedure used, we calculate a glass transition temperature of 390±8 K which compares very well with the experimental values (391±1 K from DSC and 380 K from TMA).

Thermoset Polymer Matrix Structure and Properties: Coarse-Grained Simulations

The formation of a thermoset polymer network is a complex process with great variability. In this study, we used dissipative particle dynamics and graph theory tools to investigate the curing process and network topology of a phthalonitrile thermoset to reveal the influence of initiator and plasticizer concentration on its properties. We also propose a novel way to characterize the network topology on the basis of two independent characteristics: simple cycle length (which is mainly affected by the initiator amount) and the number of simple cycles passing through a single covalent bond (which is determined primarily by plasticizer concentration). These values can be treated in the more familiar terms of network "mesh size" and "sponginess", correspondingly. The combination of these two topological parameters allows one to characterize any given network in an implicit but precise way and predict the resulting network properties, including the mechanical modulus. We believe that the same approach could be useful for other polymer networks as well, including rubbers and gels.

Parameter-free predictions of the viscoelastic response of glassy polymers from non-affine lattice dynamics

We study the viscoelastic response of amorphous polymers using theory and simulations. By accounting for internal stresses and considering instantaneous normal modes (INMs) within athermal non-affine theory, we make parameter-free predictions of the dynamic viscoelastic moduli obtained in coarse-grained simulations of polymer glasses at non-zero temperatures. The theoretical results show very good correspondence with rheology data collected from molecular dynamics simulations over five orders of magnitude in frequency, with some instabilities that accumulate in the low-frequency part on approach to the glass transition. These results provide evidence that the mechanical glass transition itself is continuous and thus represents a crossover rather than a true phase transition. The relatively sharp drop of the low-frequency storage modulus across the glass transition temperature can be explained mechanistically within the proposed theory: the proliferation of low-eigenfrequency vibrational excitations (boson peak and nearly-zero energy excitations) is directly responsible for the rapid growth of a negative non-affine contribution to the storage modulus.

Dynamic and Static Mechanical Properties of Crosslinked Polymer Matrices: Multiscale Simulations and Experiments

We studied the static and dynamic mechanical properties of crosslinked polymer matrices using multiscale simulations and experiments. We continued to develop the multiscale methodology for generating atomistic polymer networks, and applied it to the case of phthalonitrile resin. The mechanical properties of the resulting networks were analyzed using atomistic molecular dynamics (MD) and dissipative particle dynamics (DPD). The Young’s and storage moduli increased with conversion, due both to the appearance of a network of covalent bonds, and to freezing of degrees of freedom and lowering of the glass transition temperature during crosslinking. The simulations’ data showed good quantitative agreement with experimental dynamic mechanical analysis measurements at temperatures below the glass transition. The data obtained in MD and DPD simulations at elevated temperatures were conformable. This makes it possible to use the suggested approach for the prediction of mechanical properties of a broad range of polymer matrices, including ones with high structural heterogeneity.