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

Local structure of Amorphous carbon investigated by X-ray total scattering and RMC modeling

Amorphous carbon is a promising candidate as an energy storage material. In this paper, we performed an X-ray total scattering measurement, RMC modeling, and persistent homology analysis for amorphous carbon samples fabricated at two different heat treatment temperatures. According to the analysis of the nearest-neighbor carbon atoms and their angular histogram, the sample treated at higher temperature shows higher connectivity between carbon atoms than that treated at lower temperature. Furthermore, topological data analysis (persistent homology, PH) reveals quantitative results that relate ring structure and the connectivity between carbon atoms.

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., ablative material) and the multiple active species in an aggressive environment. In this work, we demonstrate that molecular dynamics (MD) simulations using the ReaxFF interatomic potential can be used to obtain the reaction kinetics of complex reaction pathways at high temperatures. We report two recently developed tools, namely, MolfrACT and KinACT, designed to extract chemical kinetic pathways by postprocessing reactive MD simulation data. The pathway extraction is based on a new algorithm, Consistent Reaction Stoichiometry via an Iterative Scheme (CReSIS), for the automated extraction of reactions and kinetics from MD trajectories. As a validation of the methodology, we first report the kinetic analysis and mechanisms for the high-temperature combustion of iso-octane. The observed reaction pathways are consistent with experimental models. In addition, we compare the activation energies of select iso-octane combustion pathways with experimental data and show that nanosecond timescale molecular dynamics simulations are sufficient to obtain realistic estimates of activation energies for different fuel consumption reaction pathways at high temperatures. The framework developed here can potentially be combined with time-series forecasting and machine learning methods to further reduce the computational complexity of transient molecular dynamics simulations, yet yielding realistic chemical kinetics information. Subsequently, the CReSIS scheme applied to ethylene-propylene-diene-monomer (EPDM) rubber ablative reveals that H₂O, C₂H₄, and HCHO are the major products during the initial stages of the polymer degradation in high-temperature oxidative environments. While prior work involving ReaxFF simulations is restricted to overall rates of formation of any species, we extract kinetic information for individual reaction pathways. In this paper, we present several reaction pathways observed during the EPDM rubber degradation into the dominant products and report the pathway-specific reaction rates. Arrhenius analysis reveals that the dominant reaction pathway activation energies for the formation of water, ethylene, and formaldehyde are 34.42, 27.26, and 6.37 kcal/mol, respectively. In contrast, the activation energies for the overall formation (across all reaction pathways) of these products are in the 40-50 kcal/mol range.

The changing state of porous materials

Porous materials contain regions of empty space into which guest molecules can be selectively adsorbed and sometimes chemically transformed. This has made them useful in both industrial and domestic applications, ranging from gas separation, energy storage and ion exchange to heterogeneous catalysis and green chemistry. Porous materials are often ordered (crystalline) solids. Order-or uniformity-is frequently held to be advantageous, or even pivotal, to our ability to engineer useful properties in a rational way. Here we highlight the growing evidence that topological disorder can be useful in creating alternative properties in porous materials. In particular, we highlight here several concepts for the creation of novel porous liquids, rationalize routes to porous glasses and provide perspectives on applications for porous liquids and glasses.

Characterizing Structural Complexity in Disordered Carbons: From the Slit Pore to Atomistic Models

The reliable characterization of nanoporous carbons is critical to the design and optimization of their numerous applications; however, the vast majority of carbons in industrial use are highly disordered, with complex structures whose understanding has long challenged researchers. The idealized slit pore model represents the most commonly used approximation to a carbon nanopore; nevertheless, it has been only partially successful in predicting adsorption isotherms and fails significantly in predicting transport properties because of its inability to capture structural disorder and its effect on fluid accessibility. Atomistic modeling of the structure has much potential for overcoming this limitation, and among such approaches, hybrid reverse Monte Carlo simulation has emerged as the most attractive. This method reconstructs the structure of a carbon based on the fitting of its experimentally measured pair distribution function and appropriate properties such as porosity while minimizing the energy. The method is shown to be best implemented using a multistage strategy, with the first stage used to attain a deep minimum of the energy and subsequent stages to refine the structure based on the fitting of specific properties. Methods to determine the accessibility of gases based on the atomistic structure are outlined, and it is shown that energy barriers are very sensitive to small differences in the sizes of constrictions and pore entries. The ability to accurately predict macroscopic transport coefficients of adsorbates in nanoporous carbons appears to be the greatest limitation of such models. Overcoming this will require the fitting of properties more sensitive to long-range disorder than the currently used pair distribution and the use of a suitable multiscaling strategy, which is suggested as a future direction for advancing atomistic models. The inclusion of heteroatoms in the structure is also an important area requiring further attention, particularly in the development of computationally efficient force fields incorporating their interactions.