Yingzhe Du, Jun Li, Ning Wen, Zheng Zhang, Dan Song
{"title":"The pyrolysis of polyimide and epoxy resin by the ReaxFF molecular dynamics simulation","authors":"Yingzhe Du, Jun Li, Ning Wen, Zheng Zhang, Dan Song","doi":"10.1007/s00894-025-06492-8","DOIUrl":null,"url":null,"abstract":"<div><h3>Context</h3><p>Polyimide (PI) and epoxy resin will age by hot corrosion and long-term high temperature, losing the heat insulting property and forming tremendous potential dangers. In order to evaluate the thermal properties accurately and detect the potential damage of them, the pyrolysis processes of them were studied. The results show that the main products of PI are CO<sub>2</sub> and CN· at high temperature, and their formation are both associated with the break of C-N bond in the imide rings. With the increase of the temperature, the number of CN· increases, but the number of CO<sub>2</sub> decrease. Among several reaction path of the PI productions, the <i>p1</i> is the lowest activation energy and can form CO<sub>2</sub>. The bond dissociation energies of C-N bond in <i>p3</i> and <i>p4</i> are higher than <i>p1</i>. We also investigated the pyrolysis process of the epoxy resin. The results show that the main products of epoxy resin are H<sub>2</sub>, CH<sub>2</sub>O, H<sub>2</sub>O, and CH<sub>4</sub> at 1300 K, and the H<sub>2</sub> is generated by the collision of the hydrogen atoms (<i>p1</i> path), the CH<sub>2</sub>O is generated by the partial decomposition of the C<sub>2</sub> or C<sub>3</sub>, which can form the epoxy groups on the ends of the epoxy resin.</p><h3>Methods</h3><p>Reactive force field (ReaxFF) molecular dynamics simulations were used to study the pyrolysis of PI and epoxy resin. The initial structures of PI and epoxy resin were constructed using Material Studio software, followed by geometry optimization to achieve the most stable configuration. Pyrolysis simulations were performed using the large-scale atomic/molecular massively parallel simulator (LAMMPS). The simulation employed NPT ensemble (0.1 MPa, 298 K) to adjust the system density to 1.0 g/cm<sup>3</sup>, and NVT ensemble for pyrolysis Calculations with a time step of 1 fs and total simulation time of 1 ns. Temperature was controlled using the Bersenden method, with key simulation temperatures including 1300 K (epoxy resin) and 2800–3800 K (PI).</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 10","pages":""},"PeriodicalIF":2.5000,"publicationDate":"2025-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Molecular Modeling","FirstCategoryId":"92","ListUrlMain":"https://link.springer.com/article/10.1007/s00894-025-06492-8","RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q4","JCRName":"BIOCHEMISTRY & MOLECULAR BIOLOGY","Score":null,"Total":0}
引用次数: 0
Abstract
Context
Polyimide (PI) and epoxy resin will age by hot corrosion and long-term high temperature, losing the heat insulting property and forming tremendous potential dangers. In order to evaluate the thermal properties accurately and detect the potential damage of them, the pyrolysis processes of them were studied. The results show that the main products of PI are CO2 and CN· at high temperature, and their formation are both associated with the break of C-N bond in the imide rings. With the increase of the temperature, the number of CN· increases, but the number of CO2 decrease. Among several reaction path of the PI productions, the p1 is the lowest activation energy and can form CO2. The bond dissociation energies of C-N bond in p3 and p4 are higher than p1. We also investigated the pyrolysis process of the epoxy resin. The results show that the main products of epoxy resin are H2, CH2O, H2O, and CH4 at 1300 K, and the H2 is generated by the collision of the hydrogen atoms (p1 path), the CH2O is generated by the partial decomposition of the C2 or C3, which can form the epoxy groups on the ends of the epoxy resin.
Methods
Reactive force field (ReaxFF) molecular dynamics simulations were used to study the pyrolysis of PI and epoxy resin. The initial structures of PI and epoxy resin were constructed using Material Studio software, followed by geometry optimization to achieve the most stable configuration. Pyrolysis simulations were performed using the large-scale atomic/molecular massively parallel simulator (LAMMPS). The simulation employed NPT ensemble (0.1 MPa, 298 K) to adjust the system density to 1.0 g/cm3, and NVT ensemble for pyrolysis Calculations with a time step of 1 fs and total simulation time of 1 ns. Temperature was controlled using the Bersenden method, with key simulation temperatures including 1300 K (epoxy resin) and 2800–3800 K (PI).
期刊介绍:
The Journal of Molecular Modeling focuses on "hardcore" modeling, publishing high-quality research and reports. Founded in 1995 as a purely electronic journal, it has adapted its format to include a full-color print edition, and adjusted its aims and scope fit the fast-changing field of molecular modeling, with a particular focus on three-dimensional modeling.
Today, the journal covers all aspects of molecular modeling including life science modeling; materials modeling; new methods; and computational chemistry.
Topics include computer-aided molecular design; rational drug design, de novo ligand design, receptor modeling and docking; cheminformatics, data analysis, visualization and mining; computational medicinal chemistry; homology modeling; simulation of peptides, DNA and other biopolymers; quantitative structure-activity relationships (QSAR) and ADME-modeling; modeling of biological reaction mechanisms; and combined experimental and computational studies in which calculations play a major role.