Determination of interaction parameters in a bottom-up approach employed in reactive dissipative particle dynamics simulation for thermosetting polymers
Kaiwen Li,Gota Kikugawa,Yoshiaki Kawagoe,Yinbo Zhao,Tomonaga Okabe
DOI: https://doi.org/10.1039/d3sm01743e
IF: 4.046
2024-05-15
Soft Matter
Abstract:The limitations in previous dissipative particle dynamics (DPD) studies confined simulations to a narrow resin range. This study refines DPD parameter calculation methodology, extending its application to diverse polymer materials. Using a bottom-up approach with molecular dynamics (MD) simulations, we evaluated solubility parameters and bead number density governing nonbonded interactions via the Flory-Huggins parameter and covalent-bonded interactions. Two solubility parameter methods, Hildebrand and Krevelen–Hoftyzer, were compared for DPD simulations. The Hildebrand method, utilizing MD simulations, demonstrates higher consistency and broader applicability in determining solubility parameters for all DPD particles. The DPD/MD curing reaction process was examined in three epoxy systems: DGEBA/4,4′-DDS, DGEBA/MPDA and DGEBA/DETA. Calculations for curing profile, gelation point, radial distribution function and branch ratio were performed. Compared to MD data for DGEBA/4,4′-DDS, the maximum deviation in secondary reactions between epoxy and amine groups by DPD simulation with Krevelen–Hoftyzer was 14.8%, while with Hildebrand, it was 1.7%. The accuracy of DPD curing reaction in reproducing structural properties verifies its expanded application to general polymeric material simulation. The proposed curing DPD simulation, with a short run time and minimal computational resources, contributes to high-throughput screening for optimal resins and investigates mesoscopic inhomogeneous structures in large resin systems.
polymer science,chemistry, physical,materials science, multidisciplinary,physics