Development of Chemical Kinetics Models from Atomistic Reactive Molecular Dynamics Simulations: Application to Iso-octane Combustion and Rubber Ablative Degradation.

Modeling the complex chemical phenomena resulting from multiple active species and long-chain polymers is limited by uncertainties in the reaction rate parameters, which increase rapidly with the number of active species and/or reaction pathways. Reactive molecular dynamics simulations have the potential to help obtain in-depth information on the chemical reactions that occur between the polymer (e.g.

Modeling high-temperature diffusion of gases in micro and mesoporous amorphous carbon.

In this work, we study diffusion of gases in porous amorphous carbon at high temperatures using equilibrium molecular dynamics simulations. Microporous and mesoporous carbon structures are computationally generated using liquid quench method and reactive force fields. Motivated by the need to understand high temperature diffusivity of light weight gases like H2, O2, H2O, and CO in amorphous carbon, we investigate the diffusion behavior as function of two important parameters: (a) the pore size and (b) the concentration of diffusing gases.

Passivation of aluminum nanoparticles by plasma-enhanced chemical vapor deposition for energetic nanomaterials.

We have produced passivating coatings on 80-nm aluminum particles by plasma-enhanced chemical vapor deposition (PECVD). Three organic precursors--isopropyl alcohol, toluene, and perfluorodecalin--were used to fabricate thin films with thicknesses ranging from 5 nm to 30 nm. The coated samples and one untreated sample were exposed to 85% humidity at 25 °C for two months, and the active Al content was determined by thermogravimetric analysis (TGA) in the presence of oxygen.

Comparison of ReaxFF, DFTB, and DFT for phenolic pyrolysis. 1. Molecular dynamics simulations.

A systematic comparison of atomistic modeling methods including density functional theory (DFT), the self-consistent charge density-functional tight-binding (SCC-DFTB), and ReaxFF is presented for simulating the initial stages of phenolic polymer pyrolysis. A phenolic polymer system is simulated for several hundred picoseconds within a temperature range of 2500 to 3500 K. The time evolution of major pyrolysis products including small-molecule species and char is examined.

Comparison of ReaxFF, DFTB, and DFT for phenolic pyrolysis. 2. Elementary reaction paths.

Reaction paths for the loss of CO, H2, and H2O from atomistic models of phenolic resin are determined using the hybrid B3LYP approach. B3LYP energetics are confirmed using CCSD(T). The energetics along the B3LYP paths are also evaluated using the PW91 generalized gradient approximation (GGA), the more approximate self-consistent charge density functional tight binding (SCC-DFTB), and the reactive force field (ReaxFF).

Polymer chain dynamics at interfaces: role of boundary conditions at solid interface.

Using classical molecular dynamics simulations, we study the dynamical properties of a single polymer chain dissolved in an explicit solvent and strongly adsorbed at solid-liquid interface. To circumvent a serious challenge posed by finite size effects due to long-range hydrodynamic effects, we developed a correction procedure that substantially limits the finite size effects. Concurrently, we provide an analysis of distinctly different size effects in the directions transverse and normal to the interface.

Modeling diffusion of adsorbed polymer with explicit solvent.

Computer simulations of a polymer chain of length N strongly adsorbed at the solid-liquid interface in the presence of explicit solvent are used to delineate the factors affecting the N dependence of the polymer lateral diffusion coefficient, D(||). We find that surface roughness has a large influence, and D(||) scales as D(||) approximately N(-x), with x approximately 3/4 and x approximately 1 for ideal smooth and corrugated surfaces, respectively. The first result is consistent with the hydrodynamics of a "particle" of radius of gyration R(G) approximately N(nu) (nu=0.

Molecular-dynamics simulations of the transport properties of a single polymer chain in two dimensions.

Molecular-dynamics simulations are conducted to elucidate the critical factors affecting the transport properties of isolated polymer chains in strictly two dimensions. The relevance of surface inhomogeneity is critically examined. We unequivocally find that surface inhomogeneity is critical in obtaining transport behavior consistent with the recent measurements of surface diffusion for polymers adsorbed at the solid-liquid interface.

Molecular dynamics simulations of polymer transport in nanocomposites.

Molecular dynamics simulations on the Kremer-Grest bead-spring model of polymer melts are used to study the effect of spherical nanoparticles on chain diffusion. We find that chain diffusivity is enhanced relative to its bulk value when polymer-particle interactions are repulsive and is reduced when polymer-particle interactions are strongly attractive. In both cases chain diffusivity assumes its bulk value when the chain center of mass is about one radius of gyration R(g) away from the particle surface.