REACTER: A Heuristic Method for Reactive Molecular Dynamics

Atomistic Modeling of Cross-Linking in Epoxy-Amine Resins: An Open-Source Protocol

Atomistic modeling has become an extensively used method for studying thermosetting polymers, particularly in the analysis and development of high-performance composite materials. Despite extensive research on the topic, a widely accepted, standardized, flexible, and open-source approach for simulating the cross-linking process from precursor molecules has yet to be established. This study proposes, tests, and validates a Molecular Dynamics (MD) protocol to simulate the cross-linking process of epoxy resins. We developed an in-house code based on Python and LAMMPS, enabling the generation of epoxy resin structures with high degrees of cross-linking. In our work, the epoxy network is dynamically formed within the MD simulations, modeling the chemical bonding process with constraints based on the distance between the reactive sites. To validate our model against experimental data from the literature, we then computed the density, thermal conductivity, and elastic response. The results show that the produced structures align well with experimental evidence, validating our method and confirming its feasibility for further analyses and in silico experiments. Beyond the case study presented in this work, focusing on bisphenol A diglycidyl ether (DGEBA) epoxy resin and diethylenetriamine (DETA) as curing agents in a 5:2 ratio, our approach can be easily adapted to investigate different epoxy resins.

Optimizing Epoxy Nanocomposites with Oxidized Graphene Quantum Dots for Superior Mechanical Performance: A Molecular Dynamics Approach

Due to their excellent mechanical properties, epoxy composites are widely used in low-density applications. However, the brittle epoxy matrix often serves as the principal failure point. Matrix enhancements can be achieved by optimizing polymer combinations to maximize intermolecular interactions or by introducing fillers. While nanofillers such as clay, rubber, carbon nanotubes, and nanoplatelets enhance mechanical properties, they can lead to issues like agglomeration, voids, and poor load transfer. Quantum dots, being the smallest nanofillers, offer higher dispersion and the potential to promote intermolecular interactions, enhancing stiffness, strength, and toughness simultaneously. This study employed molecular dynamics simulations to design graphene quantum dot (GQD) reinforced epoxy nanocomposites. By functionalizing GQDs with oxygen-based groups-hydroxyl, epoxide, carboxyl, and mixed chemistries-their effects on the mechanical properties of nanocomposites were systematically evaluated. Results show that hydroxyl-functionalized GQDs provide optimal performance, increasing stiffness and yield strength by 18.4 and 56.1%, respectively. Structural analysis reveals that these GQDs promote a closely packed molecular configuration, resulting in reduced free volume.

Thermodynamics and kinetics of early stages of carbon dot formation: a case of citric acid and ethylenediamine reaction

Owing to their extraordinary photophysical properties, carbon dots (CDs) have found applications across various fields, including bioimaging, sensing, and environmental research. Despite huge application potential, the fabrication of CDs still lacks the desired control at the molecular level, and precise structural regulation towards property-tailored CDs remains elusive. The mechanistic details of nucleation, growth, and carbonization processes leading to CDs are still unknown, with key thermodynamic and kinetic parameters yet to be revealed. Herein, we performed quantum chemical calculations of explicitly micro-hydrated reaction systems to thoroughly explore the mechanism of a prototypical reaction of citric acid and ethylenediamine. The theoretical results showed activation barriers and thermodynamics along the reaction pathway, thus helping identify key heterocyclic intermediates and cyclization products. The cyclization and condensation reactions were further simulated via a reactive molecular dynamics protocol, suggesting potential growth scenarios and generating plausible structures for further exploration of the polymerization and carbonization processes. The theoretical calculations were cross-validated with NMR and MALDI-TOF measurements. The data obtained provide a comprehensive deterministic insight into the initial stages of CD formation, revealing new reaction intermediates and pathways, and rationally predicting the formation of specific structural arrangements of premature CDs. The presented deterministic approach represents an important step towards rational bottom-up design of these unique fluorescence systems.

Macroscopic Fracture Properties of Glassy Nanocomposites from Molecular Dynamics Simulations and Empirical Force Fields

In this study, we extend a previously developed methodology for calculating macroscopic fracture properties in glassy polymers using molecular dynamics (MD) and empirical force fields to glassy nanocomposites. We apply this approach to epoxy nanocomposites with randomly dispersed carbon nanotubes (CNTs), modeling four system types: two with pristine CNTs (1 and 2% by weight) and two with hydroxyl (OH)-functionalized CNTs at similar concentrations. Using a macroscopic analytical model, we calculate the fracture energy and the stress intensity factor for each system and examine how interfacial adhesion influences load transfer and failure mechanisms. Moreover, the cohesive energies and the mean squared displacements for the different systems are calculated, to further analyze the load transfer mechanisms between CNTs and the matrix. Our results demonstrate that functionalized CNTs significantly enhance fracture properties compared to pristine CNTs due to improved interfacial adhesion, enabling better load transfer and delaying crack propagation. This study offers a computationally efficient approach for exploring fracture characteristics in CNT-epoxy nanocomposites. Using empirical force fields, we get faster calculations enabling us to model bigger and more complex systems.

Molecularly Woven Polymer Aerogels

Aerogels with abundant nanopores and large specific surface areas have extensive potential in various applications but are constrained by fragility and difficulty in degradation. Currently, the exploration of adaptive and reprocessing aerogels has become increasingly urgent, as the demand for intelligent and sustainable materials intensifies. Here, we present a molecular weaving strategy to construct molecularly woven polymer aerogels (WPAs) via catalyst-free aldimine condensation between prewoven aldehyde-functionalized Cu(I) bisphenanthroline (Cu(PBD)2) and flexible 4,4'-diaminodibenzyl (DB). The key feature of this system consists entirely of dense woven nodes that can be readily activated by external stimuli, where Cu(I) ions can also be reversibly removed as needed, while preserving porous structures. Consequently, we achieve adjustable mechanical properties of WPAs, with a 10-fold enhancement in elasticity after removing Cu(I) ions. Moreover, the destroyed WPAs demonstrate a straightforward reprocessing capacity rather than tedious monomer recovery due to the dissociation of Cu(I)-coordination bonds, the activation of sequential polymer thread motions, and the accelerated imine bond exchange enabled by adjacent Cu(I) ions. This work offers a new perspective on designing customizable and sustainable aerogels and verifies the feasibility of the emergent molecularly woven technique in a more complex functional material system.

High-performance, multi-component epoxy resin simulation for predicting thermo-mechanical property evolution during curing

High-performance epoxy systems are extensively used in structural polymer‒matrix composites for aerospace vehicles. The evolution of the thermomechanical properties of these epoxies significantly impacts the evolution of process-induced residual stresses. The corresponding process parameters need to be optimized via multiscale process modeling to minimize the residual stresses and maximize the composite strength and durability. In this study, the thermomechanical properties of a multicomponent epoxy system are predicted via molecular dynamics (MD) simulation as a function of the degree of cure to provide critical property evolution data for process modeling. In addition, the experimentally validated results of this study provide critical insight into MD modeling protocols. Among these insights, harmonic- and Morse-bond-based force fields predict similar mechanical properties. However, simulations with the Morse-bond potential fail at intermediate strain values because of cross-term energy dominance. Additionally, crosslinking simulations should be conducted at the corresponding processing temperature, because the simulation temperature impacts shrinkage evolution significantly. Multiple analysis methods are utilized to process MD heating/cooling data for glass transition temperature prediction, and the results indicate that neither method has a significant advantage. These results are important for effective and comprehensive process modeling within the ICME (Integrated Computational Materials Engineering) and Materials Genome Initiative frameworks.

Polymer dynamics under tension: Mean first passage time for looping

This study deals with polymer looping, an important process in many chemical and biological systems. We investigate basic questions on the looping dynamics of a polymer under tension using the freely jointed chain (FJC) model. Previous theoretical approaches to polymer looping under tension have relied on barrier escape methods, which assume local equilibrium, an assumption that may not always hold. As a starting point, we use an analytical expression for the equilibrium looping probability as a function of the number of monomers and applied force, predicting an inverse relationship between looping time and looping probability. Using molecular dynamics simulations, the predictions of this theoretical approach are validated within the numerical precision achieved. We compare our predictions to those of the barrier escape approach by way of a calculation of the mean first passage time (MFPT) for the ends of a polymer to cross. For this purpose, we derive the exact free energy landscape, but the resulting temporal predictions do not agree with the observed inverse scaling. We conclude that the traditional barrier escape approach does not provide satisfactory predictions for polymer looping dynamics and that the inverse scaling with looping probability offers a more reliable alternative.

Investigating the structure-property correlations of pyrolyzed phenolic resin as a function of degree of carbonization

Carbon-carbon (C/C) composites are attractive materials for high-speed flights and terrestrial atmospheric reentry applications due to their insulating thermal properties, thermal resistance, and high strength-to-weight ratio. It is important to understand the evolving structure-property correlations in these materials during pyrolysis, but the extreme laboratory conditions required to produce C/C composites make it difficult to quantify the properties in situ. This work presents an atomistic modeling methodology to pyrolyze a crosslinked phenolic resin network and track the evolving thermomechanical properties of the skeletal matrix during simulated pyrolysis. First, the crosslinked resin is pyrolyzed and the resulting char yield and mass density are verified to match experimental values, establishing the model's powerful predictive capabilities. Young's modulus, yield stress, Poisson's ratio, and thermal conductivity are calculated for the polymerized structure, intermediate pyrolyzed structures, and fully pyrolyzed structure to reveal structure-property correlations, and the evolution of properties are linked to observed structural features. It is determined that reduction in fractional free volume and densification of the resin during pyrolysis contribute significantly to the increase in thermomechanical properties of the skeletal phenolic matrix. A complex interplay of the formation of six-membered carbon rings at the expense of five and seven-membered carbon rings is revealed to affect thermal conductivity. Increased anisotropy was observed in the latter stages of pyrolysis due to the development of aligned aromatic structures. Experimentally validated predictive atomistic models are a key first step to multiscale process modeling of C/C composites to optimize next-generation materials.

Probing the Nanostructure and Reactivity of Epoxy-Amine Interphases

Understanding and controlling the structure of interphase regions in epoxy resins have been a long-standing goal in high-performance composite and coating development, since these are widely considered to be weak points in the microstructure of these materials, determining key properties such as fracture strength and barrier performance. These buried nanoscale regions are, however, inaccessible to conventional analytical techniques, and little is understood about their underlying formation mechanism. Here, we combine molecular dynamics (MD) simulation with nanoscale infrared chemical mapping to develop new understanding of the interphase using model epoxy-amine binders composed of diglycidyl ether of bisphenol A (DGEBA) cross-linked using m-xylylenediamine (MXDA). Iron oxide powders are used as exemplary surfaces, where we demonstrate that the electrostatic binding energies between the amine cross-linker and particles range from repulsive (magnetite, Fe3O4) to weakly attractive (hematite, Fe2O3) to strong immobilization (goethite, FeOOH). We find that interfacial binding occurs upon mixing and determines the overall level of residual amine content in the bulk matrix but does not correlate with a detectable amine depletion in the vicinity of particles. In all cases, an excess of both epoxy and amine functionality is detected close to particles, and the extent of matrix undercuring is found to be dependent on the entropic segregation of the unreacted material during the ambient cure. Detailed MD simulations demonstrate that spatial segregation of the unreacted precursors is expected in the interphase, leading to the experimental observation that, even after extensive postcure heating, individual particles remain embedded in a nanoscale underdeveloped environment.

Molecular structure data and modelling roadmap for optimized oxidized graphene quantum dot and epoxy interface and mechanical properties

Hybrid epoxy composites are highly considered for low-density applications due to the excellent specific strength and specific stiffness. Enhancements made to the epoxy matrix by addition of nanofillers like carbon nanotubes (CNTs) and graphene (GNPs) have been studied in detail over the course of few decades. These enhancements not only help elevate the material properties of the matrix but also activate different failure mitigating mechanisms in the composite. Although highly beneficial, there are few shortcomings due to the challenging fabrication process of integrating such. Common problems like filler agglomeration, formation of voids, wrinkling and more can result in poor load transfer within the composite. Graphene quantum dots (GQDs), on the other hand are the smallest carbon-based filler which are known to promote more intimate contact with the matrix. Their small size enables simultaneous enhancement of stiffness, strength and toughness. In addition, functionalization of these materials enables other supramolecular interactions like hydrogen bonding which improve the interfacial interaction with the epoxy. This study provides a molecular dynamics (MD) workflow to model a single functionalized GQD embedded in an epoxy matrix and the effective mechanical response of the nanocomposite. Ten chemistries were developed with different oxygen-based functional groups which capture the effect of GQD on the mechanical properties of the nanocomposite. Uniaxial strain simulations revealed that a maximum strength gain of 56 % and stiffness gain of 18 % was computed by the oxidized GQD-epoxy nanocomposite.

Interfacial Polymerization of Aromatic Polyamide Reverse Osmosis Membranes

Polyamide membranes are widely used in reverse osmosis (RO) water treatment, yet the mechanism of interfacial polymerization during membrane formation is not fully understood. In this work, we perform atomistic molecular dynamics simulations to explore the cross-linking of trimesoyl chloride (TMC) and m-phenylenediamine (MPD) monomers at the aqueous-organic interface. Our studies show that the solution interface provides a function of "concentration and dispersion" of monomers for cross-linking. The process starts with rapid cross-linking, followed by slower kinetics. Initially, amphiphilic MPD monomers diffuse in water and accumulate at the solution interface to interact with TMC monomers from the organic phase. As cross-linking progresses, a precross-linked thin film forms, reducing monomer diffusion and reaction rates. However, the structural flexibility of the amphiphilic film, influenced by interfacial fluctuations and mixed interactions with water and the organic solvent at the solution interface, promotes further cross-linking. The solubility of MPD and TMC monomers in different organic solvents (cyclohexane versus n-hexane) affects the cross-linking rate and surface homogeneity, leading to slight variations in the structure and size distribution of subnanopores. Our study of the interfacial polymerization process in explicit solvents is essential for understanding membrane formation in various solvents, which will be crucial for optimal polyamide membrane design.

Mechanically Interlocked [an]Daisy Chain Adhesives with Simultaneously Enhanced Interfacial Adhesion and Cohesion

Adhesives have been widely used to splice and repair materials to meet practical needs of humanity for thousands of years. However, developing robust adhesives with balanced adhesive and cohesive properties still remains a challenging task. Herein, we report the design and preparation of a robust mechanically interlocked [an]daisy chain network (DCMIN) adhesive by orthogonal integration of mechanical bonds and 2-ureido-4[1H]-pyrimidone (UPy) H-bonding in a single system. Specifically, the UPy moiety plays a dual role: it allows the formation of a cross-linked network and engages in multivalent interactions with the substrate for strong interfacial bonding. The mechanically interlocked [an]daisy chain, serving as the polymeric backbone of the adhesive, is able to effectively alleviate applied stress and uphold network integrity through synergistic intramolecular motions, and thus significantly improves the cohesive performance. Comparative analysis with the control made of the same quadruple H-bonding network but with non-interlocked [an]daisy chain backbones demonstrates that our DCMIN possesses superior adhesion properties over a wide temperature range. These findings not only contribute to a deep understanding of the structure-property relationship between microscopic mechanical bond motions and macroscopic adhesive properties but also provide a valuable guide for optimizing design principles of robust adhesives.

Growth of two-dimensional covalent organic frameworks on substrates: insight from microsecond atomistic simulations

While growing two-dimensional covalent organic frameworks (2D COFs) on substrates holds promise for producing functional monolayers, the presence of many defects in the resulting crystals often hinders their practical applications. Achieving structural order while suppressing defect formation necessitates a detailed atomic-level understanding. The key lies in understanding the polymerization process with high nano-scale accuracy, which presents significant challenges. Here, we perform microsecond atomistic molecular dynamics simulations to describe the deposition and polymerization of cyclohexa-m-phenylene on metal substrates, closely mimicking experimental conditions. Our improved approach highlights that 2D polymerization occurs through monomer addition and island coalescence, with a pre-bonding stage allowing monomers/oligomers to dynamically adjust their configurations to the expanding island structures. Our results elucidate the mechanisms underlying the formation of vacancy and dislocation defects during 2D polymerization as well as their healing processes. Overall, our findings underscore the significant roles that high surface mobility, effective monomer-substrate anchoring, high framework rigidity, moderate monomer coordination, and low bonding rate play in forming large, extended 2D crystals while suppressing vacancy and dislocation defects. We demonstrate how these factors can be tuned through substrate selection, deposition rate modulation, and temperature control, thereby offering valuable insight for strategically optimizing on-surface 2D polymerizations.

Evolution of Physical, Thermal, and Mechanical Properties of Poly(methyl Methacrylate)-Based Elium Thermoplastic Polymer During Polymerization

Elium-based thermoplastic composites are a key material for future use in the marine, wind energy, and automotive industries because of their recyclability and ease of manufacture. To optimize the processing of the Elium composites to yield optimal structural properties, computational process modeling can be used to relate processing parameters to residual stresses and material durability. The key ingredient for reliable and accurate process modeling is the evolution of physical, thermal, and mechanical properties during polymerization. The objective of this study is to use molecular dynamics to predict the mass density, bulk modulus, shear modulus, Young's modulus, Poisson's ratio, glass transition temperature, and coefficient of thermal expansion as a function of the extent of reaction of the polymer. The predicted properties compare favorably to the experimentally measured values in the fully polymerized state. This data set of properties provides needed input data for process modeling of Elium-based composites for process parameter optimization and improved durability and performance.

Implementing reactivity in molecular dynamics simulations with harmonic force fields

The simulation of chemical reactions and mechanical properties including failure from atoms to the micrometer scale remains a longstanding challenge in chemistry and materials science. Bottlenecks include computational feasibility, reliability, and cost. We introduce a method for reactive molecular dynamics simulations using a clean replacement of non-reactive classical harmonic bond potentials with reactive, energy-conserving Morse potentials, called the Reactive INTERFACE Force Field (IFF-R). IFF-R is compatible with force fields for organic and inorganic compounds such as IFF, CHARMM, PCFF, OPLS-AA, and AMBER. Bond dissociation is enabled by three interpretable Morse parameters per bond type and zero energy upon disconnect. Use cases for bond breaking in molecules, failure of polymers, carbon nanostructures, proteins, composite materials, and metals are shown. The simulation of bond forming reactions is included via template-based methods. IFF-R maintains the accuracy of the corresponding non-reactive force fields and is about 30 times faster than prior reactive simulation methods.

Topological Programmability of Isomerizable Polymers

Topology isomerizable networks (TINs) can be programmed into numerous polymers exhibiting unique and spatially defined (thermo-) mechanical properties. However, capturing the dynamics in topological transformations and revealing the intrinsic mechanisms of mechanical property modulation at the microscopic level is a significant challenge. Here, we use a combination of coarse-grained molecular dynamics simulations and reaction kinetic theory to reveal the impact of dynamic bond exchange reactions on the topology of branched chains. We find that, the grafted units follow a geometric distribution with a converged uniformity, which depends solely on the average grafted units of branched chains. Furthermore, we demonstrate that the topological structure can lead to spontaneous modulation of mechanical properties. The theoretical framework provides a research paradigm for studying the topology and mechanical properties of TINs.

Fracture of Epoxy Networks Using Atomistic Simulations

Predicting fracture properties through all-atomistic simulations poses challenges due to classical force field limitations in breaking covalent bonds and the computational demands of reactive force fields like ReaxFF. In addressing this, we propose a scale-bridging method for forecasting the fracture behavior of highly cross-linked epoxy combining classical force fields, the LAMMPS package REACTER, and for bond breaking a parameter based on experimental distance criterion. In our analysis, we anticipate the macroscopic fracture energy GC of the epoxy network through the application of a continuum fracture mechanics model developed for fibrils. In addition, we extract the value of the stress intensity factor KI. This modeling approach is specifically implemented for a frequently used epoxy system that consists of bisphenol F and DETDA hardener. Notably, our results demonstrate a robust correlation with existing literature and experimental studies. Moreover, our approach boasts a substantial computational time advantage, facilitating calculations that are significantly faster compared to those performed using reactive force fields.

LUNAR: Automated Input Generation and Analysis for Reactive LAMMPS Simulations

Generating simulation-ready molecular models for the LAMMPS molecular dynamics (MD) simulation software package is a difficult task and impedes the more widespread and efficient use of MD in materials design and development. Fixed-bond force fields generally require manual assignment of atom types, bonded interactions, charges, and simulation domain sizes. A new LAMMPS pre- and postprocessing toolkit (LUNAR) is presented that efficiently builds molecular systems for LAMMPS. LUNAR automatically assigns atom types, generates bonded interactions, assigns charges, and provides initial configuration methods to generate large molecular systems. LUNAR can also incorporate chemical reactivity into simulations by facilitating the use of the REACTER protocol. Additionally, LUNAR provides postprocessing for free volume calculations, cure characterization calculations, and property predictions from LAMMPS thermodynamic outputs. LUNAR has been validated via building and simulation of pure epoxy and cyanate ester polymer systems with a comparison of the corresponding predicted structures and properties to benchmark values, including experimental results from the literature. LUNAR provides the tools for the computationally driven development of next-generation composite materials in the Integrated Computational Materials Engineering (ICME) and Materials Genome Initiative (MGI) frameworks. LUNAR is written in Python with the usage of NumPy and can be used via a graphical user interface, a command line interface, or an integrated design environment. LUNAR is freely available via GitHub.

Phase patterning of liquid crystal elastomers by laser-induced dynamic crosslinking

Liquid crystal elastomers hold promise in various fields due to their reversible transition of mechanical and optical properties across distinct phases. However, the lack of local phase patterning techniques and irreversible phase programming has hindered their broad implementation. Here we introduce laser-induced dynamic crosslinking, which leverages the precision and control offered by laser technology to achieve high-resolution multilevel patterning and transmittance modulation. Incorporation of allyl sulfide groups enables adaptive liquid crystal elastomers that can be reconfigured into desired phases or complex patterns. Laser-induced dynamic crosslinking is compatible with existing processing methods and allows the generation of thermo- and strain-responsive patterns that include isotropic, polydomain and monodomain phases within a single liquid crystal elastomer film. We show temporary information encryption at body temperature, expanding the functionality of liquid crystal elastomer devices in wearable applications.

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.

Dynamic Insights into the Growth Mechanisms of 2D Covalent Organic Frameworks on Graphene Surfaces

Producing high-quality two-dimensional (2D) covalent organic frameworks (COFs) is crucial for industrial applications. However, this remains significantly challenging with current synthetic techniques. A deep understanding of the intermolecular interactions, reaction temperature, and oligomers is essential to facilitate the growth of highly crystalline COF films. Herein, molecular dynamics simulations were employed to explore the growth of 2D COFs from monomer assemblies on graphene. Our results showed that chain growth reactions dominated the COF surface growth and that van der Waals (vdW) interactions were important in enhancing the crystallinity through monomer preorganization. Moreover, appropriately tuning the reaction temperature improved the COF crystallinity and minimized the effects of amorphous oligomers. Additionally, the strength of the interface between the COF and the graphene substrate indicated that the adhesion force was proportional to the crystallinity of the COF. This work reveals the mechanisms for nucleation and growth of COFs on surfaces and provides theoretical guidance for fabricating high-quality 2D polymer-based crystalline nanomaterials.

Type Label Framework for Bonded Force Fields in LAMMPS

New functionality is added to the LAMMPS molecular simulation package, which increases the versatility with which LAMMPS can interface with supporting software and manipulate information associated with bonded force fields. We introduce the "type label" framework that allows atom types and their higher-order interactions (bonds, angles, dihedrals, and impropers) to be represented in terms of the standard atom type strings of a bonded force field. Type labels increase the human readability of input files, enable bonded force fields to be supported by the OpenKIM repository, simplify the creation of reaction templates for the REACTER protocol, and increase compatibility with external visualization tools, such as VMD and OVITO. An introductory primer on the forms and use of bonded force fields is provided to motivate this new functionality and serve as an entry point for LAMMPS and OpenKIM users unfamiliar with bonded force fields. The type label framework has the potential to streamline modeling workflows that use LAMMPS by increasing the portability of software, files, and scripts for preprocessing, running, and postprocessing a molecular simulation.

Micromechanical Dilution of PLA/PETG-Glass/Iron Nanocomposites: A More Efficient Molecular Dynamics Approach

Polylactic acid (PLA) and poly(ethylene terephthalate glycol) (PETG) are popular thermoplastics used in additive manufacturing applications. The mechanical properties of PLA and PETG can be significantly improved by introducing fillers, such as glass and iron nanoparticles (NPs), into the polymer matrix. Molecular dynamics (MD) simulations with the reactive INTERFACE force field were used to predict the mechanical responses of neat PLA/PETG and PLA-glass/iron and PETG-glass/iron nanocomposites with relatively high loadings of glass/iron NPs. We found that the iron and glass NPs significantly increased the elastic moduli of the PLA matrix, while the PETG matrix exhibited modest increases in elastic moduli. This difference in reinforcement ability may be due to the slightly greater attraction between the glass/iron NP and PLA matrix. The NASA Multiscale Analysis Tool was used to predict the mechanical response across a range of volume percent glass/iron filler by using only the neat and highly loaded MD predictions as input. This provides a faster and more efficient approach than creating multiple MD models per volume percent per polymer/filler combination. To validate the micromechanics predictions, experimental samples incorporating hollow glass microspheres (MS) and carbonyl iron particles (CIP) into PLA/PETG were developed and tested for elastic modulus. The CIP produced a larger reinforcement in elastic modulus than the MS, with similar increases in elastic modulus between PLA/CIP and PETG/CIP at 7.77 vol % CIP. The micromechanics-based mechanical predictions compare excellently with the experimental values, validating the integrated micromechanical/MD simulation-based approach.

Simulation data for engineering graphene quantum dot epoxy nanocomposites using molecular dynamics

Graphene quantum dots (GQDs) were reported to fill the role of nanofillers that enhance composite properties. Detailed investigation of this nanofiller in composites is largely unexplored. To understand the fundamental mechanisms in play, this study uses molecular dynamics simulations to reveal the effects of GQDs on epoxy properties. Mechanical simulations were performed on three varying GQD chemistries which included a pristine GQD and 2 edge aminated GQDs with different degrees of functionalization (5.2 % and 7.6 %). These GQDs were separately inserted in a polymer matrix across five individual replicates. The nanocomposite mechanical properties were computed using uniaxial strain simulations to display the effect of embedded GQDs.

Noncovalent Interactions in Mechanical Response of Thermoset Epoxy Resin

Free volume in polymers is known to influence the mechanical response of the polymers. Noncovalent interactions such as hydrogen bonds, Coulombic electrostatic interactions, and van der Waals interactions are present within these free volume regions. The manuscript presents a comprehensive identification, characterization, and evolution of noncovalent interactions as a thermoset epoxy resin (typically used as an interfacial adhesive material) is subjected to uniaxial tension, shear, and shock loading. Even though noncovalent interactions dominate uniaxial tension and shear response (up to strain levels of 50% wherein covalent bond dissociation is not observed), both covalent and noncovalent interactions define response for shock loading. Van der Waals interactions dominate the response as the samples are subjected to strain levels of 50% in tension and shear. In contrast, hydrogen bonds influence shock response.

An elastomer with ultrahigh strain-induced crystallization

Strain-induced crystallization (SIC) prevalently strengthens, toughens, and enables an elastocaloric effect in elastomers. However, the crystallinity induced by mechanical stretching in common elastomers (e.g., natural rubber) is typically below 20%, and the stretchability plateaus due to trapped entanglements. We report a class of elastomers formed by end-linking and then deswelling star polymers with low defects and no trapped entanglements, which achieve strain-induced crystallinity of up to 50%. The deswollen end-linked star elastomer (DELSE) reaches an ultrahigh stretchability of 12.4 to 33.3, scaling beyond the saturated limit of common elastomers. The DELSE also exhibits a high fracture energy of 4.2 to 4.5 kJ m-2 while maintaining low hysteresis. The heightened SIC and stretchability synergistically promote a high elastocaloric effect with an adiabatic temperature change of 9.3°C.

Mesoscale Simulation of Polymer Pyrolysis by Coarse-Grained Molecular Dynamics: A Parametric Study

Full comprehension of the pyrolysis of polymer materials is crucial for the design and application of thermal protection systems; however, it involves complex phenomena at different spatial and temporal scales. To bridge the gap between the abundant atomistic simulations and continuum modeling in the literature, we perform a novel mesoscale study of the pyrolysis process using coarse-grained molecular dynamics (CG MD) simulations. Polyethylene (PE) consisting of united atoms including implicit hydrogen is considered a model polymer, and the configurational change of PE in thermal degradation is modeled by applying the bond-breaking phenomenon based on bond energy or bond length criteria. A cook-off simulation is implemented to optimize the heuristic protocol of bond dissociation by comparing the reaction products with a ReaxFF simulation. The aerobic hyperthermal pyrolysis under oxygen bombardment is simulated at a large scale of hundreds of nanometers to observe the intricate phenomena occurring from the surface to the depth inside the material. The intrinsic thermal durability of the model polymer at extreme conditions with and without oxygen environment can be effectively simulated from the proposed mesoscale simulation to predict important thermal degradation properties required for continuum-scale pyrolysis and ablation simulations. This work serves as an initial investigation of polymer pyrolysis at the mesoscale and helps understand the concept at a larger scale.

Modeling of a two-stage polymerization considering glass fibre sizing using molecular dynamics

Fibre reinforced polymers are an important class of materials due to their light weight, high strength, and stiffness. However, there is a lack of knowledge about the interaction of fibre surface, sizing (fibre coating), and resin. Often only idealised academic systems are studied, and only rarely realistic systems that are used in an industrial context. Therefore, methods for studying the behaviour of complex sizing are highly desirable, especially as they play a crucial role in the performance of fibre reinforced polymers. Here, a simplified, yet industrially used resin system is extended using molecular dynamics simulations by adding a fibre surface and sizing layers. Furthermore, a common coupling agent was selected, and several additional assumptions were made about the structure of the sizing. Based on this, a systematic procedure for the development of a final cured system is introduced: a condensation reaction to form oligomers from coupling agent monomers is conducted. Subsequently, a two stage reaction, a polyurethane reaction and a radical polymerisation, is modelled based on an established approach. Using the final cured system, evaluations of averaged quantities during the reactions are carried out. Moreover, the system is evaluated along the normal direction of the fibre surface, which proves a spatial analysis of the fibre-sizing-resin interface.

Simulating Polymerization by Boltzmann Inversion Force Field Approach and Dynamical Nonequilibrium Reactive Molecular Dynamics

The radical polymerization process of acrylate compounds is, nowadays, numerically investigated using classical force fields and reactive molecular dynamics, with the aim to probe the gel-point transition as a function of the initial radical concentration. In the present paper, the gel-point transition of the 1,6-hexanediol dimethacrylate (HDDMA) is investigated by a coarser force field which grants a reduction in the computational costs, thereby allowing the simulation of larger system sizes and smaller radical concentrations. Hence, the polymerization is investigated using reactive classical molecular dynamics combined with a dynamical approach of the nonequilibrium molecular dynamics (D-NEMD). The network structures in the polymerization process are probed by cluster analysis tools, and the results are critically compared with the similar all-atom system, showing a good agreement.

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.

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.

Capturing Free-Radical Polymerization by Synergetic Ab Initio Calculations and Topological Reactive Molecular Dynamics

Photocurable polymers are used ubiquitously in 3D printing, coatings, adhesives, and composite fillers. In the present work, the free radical polymerization of photocurable compounds is studied using reactive classical molecular dynamics combined with a dynamical approach of the nonequilibrium molecular dynamics (D-NEMD). Different concentrations of radicals and reaction velocities are considered. The mechanical properties of the polymer resulting from 1,6-hexanediol dimethacrylate systems are characterized in terms of viscosity, diffusion constant, and activation energy, whereas the topological ones through the number of cycles (polymer loops) and cyclomatic complexity. Effects like volume shrinkage and delaying of the gel point for increasing monomer concentration are also predicted, as well as the stress-strain curve and Young's modulus. Combining ab initio, reactive molecular dynamics, and the D-NEMD method might lead to a novel and powerful tool to describe photopolymerization processes and to original routes to optimize additive manufacturing methods relying on photosensitive macromolecular systems.

Molecular Methods for Assessing the Morphology, Topology, and Performance of Polyamide Membranes

The molecular-scale morphology and topology of polyamide composite membranes determine the performance characteristics of these materials. However, molecular-scale simulations are computationally expensive and morphological and topological characterization of molecular structures are not well developed. Molecular dynamics simulation and analysis methods for the polymerization, hydration, and quantification of polyamide membrane structures were developed and compared to elucidate efficient approaches for producing and analyzing the polyamide structure. Polymerization simulations that omitted the reaction-phase solvent did not change the observed hydration, pore-size distribution, or water permeability, while improving the simulation efficiency. Pre-insertion of water into the aggregate pores (radius ≈ 4 Å) of dry domains enabled shorter hydration simulations and improved simulation scaling, without altering pore structure, properties, or performance. Medial axis and Minkowski functional methods were implemented to identify permeation pathways and quantify the polyamide morphology and topology, respectively. Better agreement between simulations and experimentally observed systems was accomplished by increasing the domain size rather than increasing the number of ensemble realizations of smaller systems. The largest domain hydrated was an order of magnitude larger by volume than the largest domain previously reported. This work identifies methods that can enable more efficient and meaningful fundamental modeling of membrane materials.

Solid with infused reactive liquid (SWIRL): A novel liquid-based separation approach for effective CO2 capture

Economical CO₂ capture demands low-energy separation strategies. We use a liquid-infused surface (LIS) approach to immobilize reactive liquids, such as amines, on a textured and thermally conductive solid substrate with high surface-area to volume ratio (A/V) continuum geometry. The infused, micrometer-thick liquid retains that high A/V and directly contacts the gas phase, alleviating mass transport resistance typically encountered in mesoporous solid adsorbents. We name this LIS class "solid with infused reactive liquid" (SWIRL). SWIRL-amine requires no water dilution or costly mixing unlike the current liquid-based commercial approach. SWIRL-tetraethylenepentamine (TEPA) shows stable, high capture capacities at power plant CO₂ concentrations near flue gas temperatures, preventing energy-intensive temperature swings needed for other approaches. Water vapor increases CO₂ capacity of SWIRL-TEPA without compromising stability.

Kinetically Corrected Monte Carlo-Molecular Dynamics Simulations of Solid Electrolyte Interphase Growth

We present a kinetic approach to the Monte Carlo-molecular dynamics (MC-MD) method for simulating reactive liquids using nonreactive force fields. A graphical reaction representation allows definition of reactions of arbitrary complexity, including their local solvation environment. Reaction probabilities and molecular dynamics (MD) simulation times are derived from ab initio calculations. Detailed validation is followed by studying the development of the solid electrolyte interphase (SEI) in lithium-ion batteries. We reproduce the experimentally observed two-layered structure on graphite, with an inorganic layer close to the anode and an outer organic layer. This structure develops via a near-shore aggregation mechanism.

The role of interface on the toughening and failure mechanisms of thermoplastic nanocomposites reinforced with nanofibrillated rubber

The interface plays a crucial role in the physical and functional properties of polymer nanocomposites, yet its effects have not been fully recognized in the setting of classical continuum-based modeling. In the present study, we investigate the roles of interface and interfiber interactions on the toughening effects of rubber nanofibers embodied in thermoplastic-based materials. Emphasis is placed on establishing comprehensive theoretical and atomistic descriptions of the nanocomposite systems subjected to pull-out and uniaxial extension in the longitudinal and transverse directions. Using the framework of molecular dynamics, the annealed melt-drawn nanofibers were spontaneously formed via the proposed four-step methodology. The generated nanofibers were then crosslinked using the proposed robust topology-matching algorithm, through which the chemical reactions arising in the crosslinking were closely assimilated. The interfiber interactions were also examined with respect to separation distances and nanofiber radius via a nanofiber-pair atomistic scheme, and the obtained results were subsequently incorporated into the pull-out and uniaxial test simulations. The results indicate that the compatibilizer grafting results in enhanced interfacial shear strength by introducing extra chemical interactions at the interface. In particular, it was found that the compatibilizer restricts the formation and coalescence of nanovoids, resulting in enhanced toughening effects. Together, we have shown that the presence of a small amount of well-dispersed rubber nanofibrillar network whose surfaces are grafted with maleic anhydride compatibilizer can dramatically increase the toughness and alter the failure mechanisms of the nanocomposites without any deterioration in the stiffness, which is also consistent with the recent experimental observations in our lab. The interfacial failure mechanism was also investigated by monitoring the changes in the atomic concentration profiles, mean square displacement and fractional free volume. The results obtained may serve as a promising alternative for the continuum-based modeling and analysis of interfaces.

Characterization of Cure Behavior in Epoxy Using Molecular Dynamics Simulation Compared with Dielectric Analysis and DSC

The curing behavior of a thermosetting material that influences the properties of the material is a key issue for predicting the changes in material properties during processing. An empirical equation can describe the reaction kinetics of the curing behavior of an investigated material, which is usually estimated using experimental methods. In this study, the curing process of an epoxy resin, the polymer matrix in an epoxy molding compound, is computed concerning thermal influence using molecular dynamics. Furthermore, the accelerated reaction kinetics, which are influenced by an increased reaction cutoff distance, are investigated. As a result, the simulated crosslink density with various cutoff distances increases to plateau at a crosslink density of approx. 90% for the investigated temperatures during curing time. The reaction kinetics are derived according to the numerical results and compared with the results using experimental methods (dielectric analysis and differential scanning calorimetry), whereby the comparison shows a good agreement between experiment and simulation.

Thermoset resin curing simulation using quantum-chemical reaction path calculation and dissipative particle dynamics

Thermoset resin, which is commonly used as a matrix in carbon-fiber-reinforced plastic, requires curing procedures. We propose a curing simulation technique involving a dissipative particle dynamics (DPD) simulation, which can simulate a larger system and longer time period than those of conventional all-atom molecular dynamics (AA-MD) simulations. The proposed curing DPD simulation can represent the thermoset resin exothermic reaction process precisely by considering each reactivity according to the reaction types calculated via quantum-chemical reaction path calculations. The cure reaction process given by the curing DPD simulation agrees well with that given by a conventional curing AA-MD simulation, but with run-time and computational-resource reductions of 1/480 and 1/10 times, respectively. We also conduct reverse mapping, through which the AA-MD system can be reconstructed from the DPD system, to evaluate the structural and thermomechanical properties. The X-ray diffraction pattern and thermomechanical properties of the reconstructed system agree well with those of the systems derived from the curing AA-MD simulation and experimental setup. Therefore, a cured-resin AA-MD system can be obtained from a curing DPD simulation at an extremely low computational cost, and the thermomechanical properties can be evaluated precisely using this system. The proposed curing simulation technique can be applied in high-throughput screening for better materials properties and in large system calculations.

Dynamic behaviour of water molecules in heterogeneous free space formed in an epoxy resin

Although an epoxy resin is a stable material, it absorbs moisture over a long period of time, causing deterioration of its material properties. We here applied a full-atomistic molecular dynamics (MD) simulation to study where water molecules exist in an epoxy resin and how they dynamically behave. First, the curing reaction was simulated to obtain a network structure so that the time course of the density, and thereby the free space, in the resin were obtained. The results made it possible to discuss the formation and size distribution of the free spaces which were not connected to each other. Then, a few percent of water were inserted into the free space of the cured epoxy resin to examine the location and dynamics of their molecules. We found that several water molecules were clustered at a preferred site, where hydrogen bonds can be formed with hydroxy, ether and amino groups of the network, in the free space, and they heterogeneously moved from there to other sites.