Effect of Thickness on the Thermo-Mechanical Response of Free-Standing Thermoset Nanofilms from Molecular Dynamics

Relief of Residual Stress in Bulk Thermosets in the Glassy State by Local Bond Exchange

The covalently cross-linked network gives thermosets superior thermal, mechanical, and electrical properties, which, however, squarely makes the large residual stress that is inevitably induced during preparation hardly relieved in the glassy state. In this work, an incredible reduction in residual stress is successfully achieved in bulk thermosets in the glassy state through introducing highly dynamic thiocarbamate bonds by "click" reactions of thiols and isocyanates. Due to the excellent dynamic behaviors of thiocarbamate bonds, local network rearrangement is achieved through thermal stimulation, while the strong 3D cross-linked network is well maintained. Ultimately, a decrease by 44% in residual stress is detected by simply annealing samples at 30 °C below glass transition temperature (Tg), during which they could well maintain more than 98.4% of the storage modulus. After the annealing, more uniform residual stress distribution is also observed, showing a 32% decline in sample standard deviation. However, the residual stress of epoxy resin, a typical thermoset as a reference, changes little even after annealing at Tg. The results prove it a feasible strategy to reduce residual stress in bulk thermosets in the glassy state by introducing proper dynamic covalent bonds.

Uncovering the Mechanism of Size Effect on the Thermomechanical Properties of Highly Cross-Linked Epoxy Resins

Epoxy resins are widely used as matrix resins, especially for carbon-fiber-reinforced plastic, due to their outstanding physical and mechanical properties. To date, most research into cross-linking processes using simulation has considered only a distance-based criterion to judge the probability of reaction. In this work, a new algorithm was developed for use with the large-scale atomic/molecular massively parallel simulator (LAMMPS) simulation package to study the cross-linking process; this new approach combines both a distance-based criterion and several kinetic criteria to identify whether the reaction has occurred. Using this simulation framework, we investigated the effect of model size on predicted thermomechanical properties of three different structural systems: diglycidyl ether of bisphenol A (DGEBA)/4,4'-diaminodiphenyl sulfone (4,4'-DDS), DGEBA/diethylenetriamine (DETA), and tetraglycidyl diaminodiphenylmethane (TGDDM)/4,4'-DDS. Derived values of gel point, volume shrinkage, and cross-linked resin density were found to be insensitive to model size in these three systems. Other thermomechanical properties, i.e., glass-transition temperature, Young's modulus, and yield stress, were found to reach stable values for systems larger than ∼40 000 atoms for both DGEBA/4,4'-DDS and DGEBA/DETA. However, these same properties modeled for TGDDM/4,4'-DDS did not stabilize until the system size reached 50 000 atoms. Our results provide general guidelines for simulation system size and procedures to more accurately predict the thermomechanical properties of epoxy resins.

Identifying Structure-Property Relationships of Micro-Architectured Porous Scaffolds through 3D Printing and Finite Element Analysis

This study integrates 3D printing and finite element analysis (FEA) to investigate the effect of micro-architectural characteristics on the mechanical properties of porous scaffolds. The studied characteristics include the thickness of the scaffold walls and the number of domains at the cross-section. We use 3D printing to fabricate scaffolds of deliberately designed microstructures to enable strict control of the structures. The longitudinal compressive properties of different scaffolds are first analyzed through experimental testing. Then, FEA is conducted to investigate the mechanical properties and the deformation mechanisms of the scaffolds. We find that decreasing wall thickness leads to failure mechanism transition from wall compression failure to buckling instability. For scaffolds with different wall thicknesses, the failure mechanisms and the critical loads are evaluated using the theory of thin plate buckling. For the characteristic of the number of domains, both experimental and FEA results indicate increasing effective stiffness with increasing domains. Interestingly, we find that with the material properties extracted from a single wall scaffold, the computational models tend to overestimate the effective compression modulus of scaffolds with larger numbers of walls or domains than the experimental data. This observation indicates possible size-dependent material properties in 3D printed structs. Our study demonstrates that integrating experiments and computational modeling can provide fundamental insights into the mechanical properties and deformation mechanisms of micro-architectured scaffolds and unveil unique small-scale material behaviors.

Dynamic Mechanical Behaviors of Nacre-Inspired Graphene-Polymer Nanocomposites Depending on Internal Nanostructures

Nacre, a natural nanocomposite with a brick-and-mortar structure existing in the inner layer of mollusk shells, has been shown to optimize strength and toughness along the laminae (in-plane) direction. However, such natural materials more often experience impact load in the direction perpendicular to the layers (i.e., out-of-plane direction) from predators. The dynamic responses and deformation mechanisms of layered structures under impact load in the out-of-plane direction have been much less analyzed. This study investigates the dynamic mechanical behaviors of nacre-inspired layered nanocomposite films using a model system that comprises alternating multi-layer graphene (MLG) and polymethyl methacrylate (PMMA) phases. With a validated coarse-grained molecular dynamics simulation approach, we systematically study the mechanical properties and impact resistance of the MLG-PMMA nanocomposite films with different internal nanostructures, which are characterized by the layer thickness and number of repetitions while keeping the total volume constant. We find that as the layer thickness decreases, the effective modulus of the polymer phase confined by the adjacent MLG phases increases. Using ballistic impact simulations to explore the dynamic responses of nanocomposite films in the out-of-plane direction, we find that the impact resistance and dynamic failure mechanisms of the films depend on the internal nanostructures. Specifically, when each layer is relatively thick, the nanocomposite is more prone to spalling-like failure induced by compressive stress waves from the projectile impact. Whereas, when there are more repetitions, and each layer becomes relatively thin, a high-velocity projectile sequentially penetrates the nanocomposite film. In the low projectile velocity regime, the film develops crazing-like deformation zones in PMMA phases. We also show that the position of the soft PMMA phase relative to the stiff graphene sheets plays a significant role in the ballistic impact performance of the investigated films. Our study provides insights into the effect of nanostructures on the dynamic mechanical behaviors of layered nanocomposites, which can lead to effective design strategies for impact-resistant films.

Novel Mode of Noncrystallographic Branching in the Initial Stages of Polymer Fibril Growth

Spherulites are the most ubiquitous of polycrystalline microstructure of polymers; they develop under a wide range of conditions by the subsequent branching of crystalline lamella that results in an overall spherical shape. Despite significant efforts over decades, the mechanisms behind branching remain unclear. Molecular dynamics simulations in polyethylene reveal the molecular-level origin of noncrystallographic branching and the initial formation of fibrils. We find that the growth of crystalline lamella by reeling in and folding of polymer chains causes surprisingly large local deformation which, in turn, aligns the chains in the neighboring undercooled liquid. Thus, subsidiary grains nucleate with preferred orientations resulting in fibril growth with branching at small angles, consistent with those observed experimentally.

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.

Theory of the spatial transfer of interface-nucleated changes of dynamical constraints and its consequences in glass-forming films

We formulate a new theory for how caging constraints in glass-forming liquids at a surface or interface are modified and then spatially transferred, in a layer-by-layer bootstrapped manner, into the film interior in the context of the dynamic free energy concept of the Nonlinear Langevin Equation (NLE) theory approach. The dynamic free energy at any mean location (cage center) involves contributions from two adjacent layers where confining forces are not the same. At the most fundamental level of the theory, the caging component of the dynamic free energy varies essentially exponentially with distance from the interface, saturating deep enough into the film with a correlation length of modest size and weak sensitivity to the thermodynamic state. This imparts a roughly exponential spatial variation of all the key features of the dynamic free energy required to compute gradients of dynamical quantities including the localization length, jump distance, cage barrier, collective elastic barrier, and alpha relaxation time. The spatial gradients are entirely of dynamical, not structural or thermodynamic, origin. The theory is implemented for the hard sphere fluid and diverse interfaces which can be a vapor, a rough pinned particle solid, a vibrating (softened) pinned particle solid, or a smooth hard wall. Their basic description at the level of the spatially heterogeneous dynamic free energy is identical, with the crucial difference arising from the first layer where dynamical constraints can be weakened, softened, or hardly changed depending on the specific interface. Numerical calculations establish the spatial dependence and fluid volume fraction sensitivity of the key dynamical property gradients for five different model interfaces. A comparison of the theoretical predictions for the dynamic localization length and glassy modulus with simulations and experiments for systems with a vapor interface reveals good agreement. The present advance sets the stage for using the Elastically Collective NLE theory to make quantitative predictions for the alpha relaxation time gradient, decoupling phenomena, T_g gradient, and many film-averaged properties of both model and experimental (colloids, molecules, and polymers) systems with diverse interfaces and chemical makeup.

Engineering curvature in graphene ribbons using ultrathin polymer films

We propose a method to induce curvature in graphene nanoribbons in a controlled manner using an ultrathin thermoset polymer in a bimaterial strip setup and test it via molecular dynamics (MD) simulations. Continuum mechanics shows that curvature develops to release the residual stress caused by the chemical and thermal shrinkage of the polymer during processing and that this curvature increases with decreasing film thickness; however, significant deformation is only achieved for ultrathin polymer films. Quite surprisingly, explicit MD simulations of the curing and annealing processes show that the predicted trend not just continues down to film thicknesses of 1-2 nm but that the curvature development is enhanced significantly in such ultrathin films due to surface tension effects. This combination of effects leads to very large curvatures of over 0.14 nm(-1) that can be tuned via film thickness. This provides a new avenue to engineer curvature and, thus, electromagnetic properties of graphene.

Coupled effects of substrate adhesion and intermolecular forces on polymer thin film glass-transition behavior

Intermolecular noncovalent forces between polymer chains influence the mobility and glass-transition temperature (Tg), where weaker interchain interactions, all else being the same, typically results in lower bulk polymer Tg. Using molecular dynamics simulations, here we show that this relation can become invalid for supported ultrathin films when the substrate-polymer interaction is extremely strong and the polymer-polymer interactions are much weaker. This contrasting trend is found to be due to a more pronounced substrate-induced appreciation of the film Tg for polymers with weaker intermolecular interactions and low bulk Tg. We show that optimizing this coupling between substrate adhesion and bulk Tg maximizes thin film Tg, paving the way for tuning film properties through interface nanoengineering.