Sato Ryuhei, Kawakami Shinya, Ejima Hirotaka, Ujii Takahiro, Sato Koichi, Ichiki Takanori, Shibuta Yasushi
Department of Materials Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan.
Technology Division Research and Development Department, NIPPON PAINT CO. LTD., 4-1-15, Minamishinagawa, Shinagawa-ku, Tokyo 140-8675, Japan.
J Chem Theory Comput. 2024 Dec 24;20(24):10798-10806. doi: 10.1021/acs.jctc.4c01418. Epub 2024 Dec 2.
Coarse-grained molecular dynamics (CG-MD) simulations and subsequent persistent homology (PH) analysis were performed to correlate the structure and stress-strain behavior of polymer films. During uniaxial tensile MD simulations, the first principal component of the persistence diagram obtained by principal component analysis (PCA) was in good agreement with the stress-strain curve. This indicates that PH + PCA can identify critical ring structures relevant to the dynamic changes in MD simulations without requiring any prior knowledge. Inverse analysis of the persistence diagram revealed that smaller rings with ten or fewer CG beads mainly contribute to changes in the first principal component of the persistence diagram. This is due to the properties of the poly(ethylene oxide) chain, which favors the formation of a seven-membered helical structure during the self-entanglement process. The PH + PCA approach successfully reproduced the stress-strain curves for polymers with different nonbonding interactions and bond lengths. Moreover, the changes in the yield stress of each polymer film were qualitatively explained by the ring distribution in the persistence diagram. These results suggest that persistent homology analysis followed by PCA provides a versatile and powerful framework for correlating structural features with physical properties, such as ring distribution and stress-strain behavior in polymer films.
进行了粗粒度分子动力学(CG-MD)模拟以及后续的持久同调(PH)分析,以关联聚合物薄膜的结构与应力-应变行为。在单轴拉伸MD模拟过程中,通过主成分分析(PCA)获得的持久图的第一主成分与应力-应变曲线吻合良好。这表明PH + PCA能够识别与MD模拟中的动态变化相关的关键环结构,而无需任何先验知识。对持久图的逆分析表明,具有十个或更少CG珠子的较小环主要对持久图的第一主成分的变化有贡献。这是由于聚环氧乙烷链的性质,其在自缠结过程中有利于形成七元螺旋结构。PH + PCA方法成功再现了具有不同非键相互作用和键长的聚合物的应力-应变曲线。此外,通过持久图中的环分布定性地解释了每个聚合物薄膜屈服应力的变化。这些结果表明,随后进行PCA的持久同调分析为将结构特征与物理性质(如聚合物薄膜中的环分布和应力-应变行为)相关联提供了一个通用且强大的框架。