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.

Evolution of Glassy Carbon Derived from Pyrolysis of Furan Resin.

Glassy carbon (GC) material derived from pyrolyzed furan resin was modeled by using reactive molecular dynamics (MD) simulations. The MD polymerization simulation protocols to cure the furan resin precursor material are validated via comparison of the predicted density and Young's modulus with experimental values. The MD pyrolysis simulations protocols to pyrolyze the furan resin precursor is validated by comparison of calculated density, Young's modulus, carbon content, sp2 carbon content, the in-plane crystallite size, out-of-plane crystallite stacking height, and interplanar crystallite spacing with experimental results from the literature for furan resin derived GC. The modeling methodology established in this work can provide a powerful tool for the modeling-driven design of next-generation carbon-carbon composite precursor chemistries for thermal protection systems and other high-temperature applications.

Internal circulation and mixing within tight-squeezing deformable droplets.

The internal flow and mixing properties inside deformable droplets, after reaching the steady state within two types of passive droplet traps, are visualized and analyzed as dynamical systems. The first droplet trap (constriction) is formed by three spheres arranged in an equilateral triangle, while the second consists of two parallel spherocylinders (capsules). The systems are assumed to be embedded in a uniform far-field flow at low Reynolds number, and the steady shapes and interfacial velocities on the drops are generated using the boundary-integral method. The internal velocity field is recovered by solving the internal Dirichlet problem, also via a desingularized boundary-integral method. Calculation of 2D streamlines within planes of symmetry reveals the internal equilibria of the flow. The type of each equilibrium is classified in 3D and their interactions probed using passive tracers and their Poincaré maps. For the two-capsule droplet, saddle points located on orthogonal symmetry planes influence the regular flow within the drop. For the three-sphere droplet, large regions of chaos are observed, embedded with simple periodic orbits. Flow is visualized via passive dyes, using material lines and surfaces. In 2D, solely the interface between two passive interior fluids is advected using an adaptive number of linked tracer particles. The reduction in dimension decreases the number of required tracer points, and also resolves arbitrarily thin filaments, in contrast to backward cell-mapping methods. In 3D, the advection of a material surface, bounded by the droplet interface, is enabled using an adaptive mesh scheme. Off-lattice 3D contour advection allows for highly resolved visualizations of the internal flow and quantification of the associated degree of mixing. Analysis of the time-dependent growth of material surfaces and 3D mixing numbers suggests the three-sphere droplet exhibits superior mixing properties compared to the two-capsule droplet.

Dynamics of carbohydrate strands in water and interactions with clay minerals: influence of pH, surface chemistry, and electrolytes.

Carbohydrate hydrogels are extensively used in pharmaceuticals and engineered biomaterials. Molecular conformations, assembly, and interactions of the carbohydrate strands with stabilizers such as clay minerals in aqueous solution are difficult to quantify in experiments and the hydrogel properties remain largely a result of trial-and-error studies. We analyzed the assembly of gellan gum in aqueous solution and interactions with dispersed clay minerals in all-atomic detail using molecular dynamics simulation, atomic force microscopy (AFM), and comparisons to earlier measurements. Gellan strands associate at low pH values of 2 and gradually disassemble to double strands with weak association of -0.4 kcal per mole carbohydrate ring as the pH values increases to 9. Ionization of the carbonic acid side groups in the backbone extends the chains and accelerates the conformational dynamics via rapidly changing intramolecular ion bridges. Gellan interactions with clay minerals depend on the strength of electric triple layers between clay, cations, and anionic polymer strands, as well as weaker hydrogen bonds along the edges, which are tunable as a function of the clay surface chemistry, local ionic strength, and pH values. Interaction energies range from -4 to +6 kcal per mol carbohydrate ring and were most favorable for electric triple layers with high charge mobility, which can be achieved for intermediate cation exchange capacity of the clay mineral and high pH values to increase ionization of the clay edges and of the polymer. The findings provide understanding and help control the dynamics and stabilization of carbohydrate hydrogels by clay minerals.

Nanoscale Structure-Property Relationships of Polyacrylonitrile/CNT Composites as a Function of Polymer Crystallinity and CNT Diameter.

Polyacrylonitrile (PAN)/carbon nanotube (CNT) composites are used as precursors for ultrastrong and lightweight carbon fibers. However, insights into the structure at the nanoscale and the relationships to mechanical and thermal properties have remained difficult to obtain. In this study, molecular dynamics simulation with accurate potentials and available experimental data were used to describe the influence of different degrees of PAN preorientation and CNT diameter on the atomic-scale structure and properties of the composites. The inclusion of CNTs in the polymer matrix is favored for an intermediate degree of PAN orientation and small CNT diameter whereas high PAN crystallinity and larger CNT diameter disfavor CNT inclusion. The glass transition at the CNT/PAN interface involves the release of rotational degrees of freedom of the polymer backbone and increased mobility of the protruding nitrile side groups in contact with the carbon nanotubes. The glass-transition temperature of the composite increases in correlation with the amount of CNT/polymer interfacial area per unit volume, i.e., in the presence of CNTs, for higher CNT volume fraction,  and inversely with CNT diameter. The increase in glass-transition temperature upon CNT addition is larger for PAN of lower crystallinity than for PAN of higher crystallinity. Interfacial shear strengths of the composites are higher for CNTs of smaller diameter and for PAN with preorientation, in correlation with more favorable CNT inclusion energies. The lowest interfacial shear strength was observed in amorphous PAN for the same CNT diameter. PAN with ∼75% crystallinity exhibited hexagonal patterns of nitrile groups near and far from the CNT interface which could influence carbonization into regular graphitic structures. The results illustrate the feasibility of near-quantitative insights into macroscale properties of polymer/CNT composites from simulations of nanometer-scale composite domains. Guidance is most effective when key assumptions in experiment and simulation are closely aligned, such as exfoliation versus bundling of CNTs, size, type, potential defects of CNTs, and precise measures for polymer crystallinity.

Carbon Nanotube Dispersion in Solvents and Polymer Solutions: Mechanisms, Assembly, and Preferences.

Debundling and dispersion of carbon nanotubes (CNTs) in polymer solutions play a major role in the preparation of carbon nanofibers due to early effects on interfacial ordering and mechanical properties. A roadblock toward ultrastrong fibers is the difficulty to achieve homogeneous dispersions of CNTs in polyacrylonitrile (PAN) and poly(methyl methacrylate) (PMMA) precursor solutions in solvents such as dimethyl sulfoxide (DMSO), N,N-dimethylacetamide (DMAc), and N,N-dimethylformamide (DMF). In this contribution, molecular dynamics simulations with accurate interatomic potentials for graphitic materials that include virtual π electrons are reported to analyze the interaction of pristine single wall CNTs with the solvents and polymer solutions at 25 °C. The results explain the barriers toward dispersion of SWCNTs and quantify CNT-solvent, polymer-solvent, as well as CNT-polymer interactions in atomic detail. Debundling of CNTs is overall endothermic and unfavorable with dispersion energies of +20 to +30 mJ/m2 in the pure solvents, + 20 to +40 mJ/m2 in PAN solutions, and +20 to +60 mJ/m2 in PMMA solutions. Differences arise due to molecular geometry, polar, van der Waals, and CH-π interactions. Among the pure solvents, DMF restricts CNT dispersion less due to the planar geometry and stronger van der Waals interactions. PAN and PMMA interact favorably with the pure solvents with dissolution energies of -0.7 to -1.1 kcal per mole monomer and -1.5 to -2.2 kcal per mole monomer, respectively. Adsorption of PMMA onto CNTs is stronger than that of PAN in all solvents as the molecular geometry enables more van der Waals contacts between alkyl groups and the CNT surface. Polar side groups in both polymers prefer interactions with the polar solvents. Higher polymer concentrations in solution lead to polymer aggregation via alkyl groups and reduce adsorption onto CNTs. PAN and PMMA solutions in DMSO and dilute solutions in DMF support CNT dispersion more than other combinations whereby the polymers significantly adsorb onto CNTs in DMSO solution. The observations by molecular simulations are consistent with available experimental data and solubility parameters and aid in the design of carbon nanofibers. The methods can be applied to other multiphase graphitic materials.

Chemical Reactions in Classical Molecular Dynamics.

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.