Development of dissociative force field for all‐atomistic molecular dynamics calculation of fracture of polymers

Kazushi Fujimoto,Rajadeep Singh Payal,Tomonori Hattori,Wataru Shinoda,Masayuki Nakagaki,Shigeyoshi Sakaki,Susumu Okazaki

DOI: https://doi.org/10.1002/jcc.26034

2019-07-19

Journal of Computational Chemistry

Abstract:<p>A dissociative force field for all‐atomistic molecular dynamics calculations has been developed to investigate impact fracture of polymers accompanying dissociation of chemical bonds of polymer main chain. Energy of dimer molecules was evaluated as a function of both bond‐length <i>b</i> and bond‐angle <i>θ</i> by CASPT2 calculations, whose quality is enough to describe dissociation of chemical bonds. Because we found that the bond dissociation energy <i>D</i> decreases with increasing bond‐angle, we employed the Morse‐type function <span><i>V</i>_Bond(<i>b</i>, <i>θ</i>) = {<i>D</i> − <i>V</i>_Angle(<i>θ</i>)}[1 − exp{−<i>α</i>(<i>b</i> − <i>b</i>₀) − <i>β</i>(<i>b</i> − <i>b</i>₀)²}]</span> where a quartic function <span><i>V</i>_Angle(<i>θ</i>) = <i>k</i>₁(<i>θ</i> − <i>θ</i>₀) + <i>k</i>₂(<i>θ</i> − <i>θ</i>₀)² + <i>k</i>₃(<i>θ</i> − <i>θ</i>₀)³ + <i>k</i>₄(<i>θ</i> − <i>θ</i>₀)⁴</span>. This function reproduced well the CASPT2 potential energy surface in a wide range of <i>b</i> and <i>θ</i>. The parameters have been obtained for four popular glassy polymers, polyethylene, poly(methyl methacrylate), poly(styrene), and polycarbonate. © 2019 Wiley Periodicals, Inc.</p>

chemistry, multidisciplinary