As the demand for PET plastic products continues to grow, developing effective processes to reduce their pollution is of critical importance. Pyrolysis, a promising technology to produce lighter and recyclable components from wasted plastic products, has therefore received considerable attention. In this work, the rapid pyrolysis of PET was studied by using reactive molecular dynamics (MD) simulations. Mechanisms for yielding gas species were unraveled, which involve the generation of ethylene and TPA radicals from ester oxygen−alkyl carbon bond dissociation and condensation reactions to consume TPA radicals with the products of long chains containing a phenyl benzoate structure and CO2. As atomistic simulations are typically conducted at the time scale of a few nanoseconds, a high temperature (i.e. >1000 K) is adopted for accelerated reaction events. To apply the results from MD simulations to practical pyrolysis processes, a kinetic model based on a set of ordinary differential equations was established, which is capable of describing the key products of PET pyrolysis as a function of time and temperature. It was further exploited to determine the optimal reaction conditions for low environmental impact. Overall, this study conducted a detailed mechanism study of PET pyrolysis and established an effective kinetic model for the main species. The approach presented herein to extract kinetic information such as detailed kinetic constants and activation energies from atomistic MD simulations can also be applied to related systems such as the pyrolysis of other polymers.
more »
« less
Multiply Accelerated ReaxFF Molecular Dynamics: Coupling Parallel Replica Dynamics with Collective Variable Hyper Dynamics
To tackle the time scales required to study complex chemical reactions, methods performing accelerated molecular dynamics are necessary even with the recent advancement in high-performance computing. A number of different acceleration techniques are available. Here we explore potential synergies between two popular acceleration methods – Parallel Replica Dynamics (PRD) and Collective Variable Hyperdynamics (CVHD), by analysing the speedup obtained for the pyrolysis of n-dodecane. We observe that PRD + CVHD provides additional speedup to CVHD by reaching the required time scales for the reaction at an earlier wall-clock time. Although some speedup is obtained with the additional replicas, we found that the effectiveness of the inclusion of PRD is depreciated for systems where there is a dramatic increase in reaction rates induced by CVHD. Similar observations were made in the simulation of ethylene-carbonate/Li system, which is inherently more reactive than pyrolysis, indicate that the speedup obtained via the combination of the two acceleration methods can be generalised to most practical chemical systems.
more »
« less
- Award ID(s):
- 1807740
- PAR ID:
- 10110480
- Date Published:
- Journal Name:
- Molecular simulation
- ISSN:
- 0892-7022
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Simulating the dynamics of ions near polarizable nanoparticles (NPs) using coarse-grained models is extremely challenging due to the need to solve the Poisson equation at every simulation timestep. Recently, a molecular dynamics (MD) method based on a dynamical optimization framework bypassed this obstacle by representing the polarization charge density as virtual dynamic variables and evolving them in parallel with the physical dynamics of ions. We highlight the computational gains accessible with the integration of machine learning (ML) methods for parameter prediction in MD simulations by demonstrating how they were realized in MD simulations of ions near polarizable NPs. An artificial neural network–based regression model was integrated with MD simulation and predicted the optimal simulation timestep and optimization parameters characterizing the virtual system with 94.3% success. The ML-enabled auto-tuning of parameters generated accurate dynamics of ions for ≈ 10 million steps while improving the stability of the simulation by over an order of magnitude. The integration of ML-enhanced framework with hybrid Open Multi-Processing / Message Passing Interface (OpenMP/MPI) parallelization techniques reduced the computational time of simulating systems with thousands of ions and induced charges from thousands of hours to tens of hours, yielding a maximum speedup of ≈ 3 from ML-only acceleration and a maximum speedup of ≈ 600 from the combination of ML and parallel computing methods. Extraction of ionic structure in concentrated electrolytes near oil–water emulsions demonstrates the success of the method. The approach can be generalized to select optimal parameters in other MD applications and energy minimization problems.more » « less
-
Abstract Dong et al. (2020,https://doi.org/10.1029/2019JD030691) employed a new compressible model to examine gravity wave (GW) self‐acceleration dynamics, instabilities, secondary gravity wave (SGW) generation, and mean forcing for GW packets localized in two dimensions (2D). This paper extends the exploration of self‐acceleration dynamics to a GW packet localized in three dimensions (3D) propagating into tidal winds in the mesosphere and thermosphere. As in the 2D packet responses, 3D GW self‐acceleration dynamics are found to be significant and include 3D GW phase distortions, stalled GW vertical propagation, local instabilities, and SGW and acoustic wave generation. Additional 3D responses described here include refraction by tidal winds, localized 3D instabilities, asymmetric SGW propagation, reduced SGW and acoustic wave responses at higher altitudes relative to 2D responses, and forcing of transient, large‐scale, 3D mean responses that may have implications for chemical and microphysical processes operating on longer time scales.more » « less
-
Molecular dynamics (MD) simulations with full-dimensional potential energy surfaces (PESs) obtained from high-level ab initio calculations are frequently used to model reaction dynamics of small molecules (i.e., molecules with up to 10 atoms). Construction of full-dimensional PESs for larger molecules is, however, not feasible since the number of ab initio calculations required grows rapidly with the increase of dimension. Only a small number of coordinates are often essential for describing the reactivity of even very large systems, and reduced-dimensional PESs with these coordinates can be built for reaction dynamics studies. While analytical methods based on transition-state theory framework are well established for analyzing the reduced-dimensionalPESs, MD simulation algorithms capable of generating trajectories on such surfaces are more rare. In this work, we present a new MD implementation that utilizes the relaxed reduced-dimensional PES for standard micro canonical (NVE) and canonical (NVT) MD simulations.The method is applied to the pyramidal inversion of a NH3molecule. The results from the MD simulations on a reduced, three-dimensional PES are validated against the ab initio MD simulations, as well as MD simulations on full-dimensional PES and experimental data.more » « less
-
null (Ed.)The rates of liquid-phase, acid-catalyzed reactions relevant to the upgrading of biomass into high-value chemicals are highly sensitive to solvent composition and identifying suitable solvent mixtures is theoretically and experimentally challenging. We show that the complex atomistic configurations of reactant–solvent environments generated by classical molecular dynamics simulations can be exploited by 3D convolutional neural networks to enable accurate predictions of Brønsted acid-catalyzed reaction rates for model biomass compounds. We develop a 3D convolutional neural network, which we call SolventNet, and train it to predict acid-catalyzed reaction rates using experimental reaction data and corresponding molecular dynamics simulation data for seven biomass-derived oxygenates in water–cosolvent mixtures. We show that SolventNet can predict reaction rates for additional reactants and solvent systems an order of magnitude faster than prior simulation methods. This combination of machine learning with molecular dynamics enables the rapid, high-throughput screening of solvent systems and identification of improved biomass conversion conditions.more » « less