REACTER: A Heuristic Method for Reactive Molecular Dynamics

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.