Rescaled Coarse-Grained Potentials for Maximum Timestep in Simulation of Polyethylene

In order to better understand the physical properties of polyethylene, an extremely common plastic used mostly in packaging, many scientists and engineers use olecular dynamics. To reduce the computational expense associated with traditional atomistic molecular dynamics, coarse-grained molecular dynamics is often used. Coarse-grained molecular dynamics groups multiple atoms into single beads, reducing the number of degrees of freedom in a system and eliminating smaller atoms with faster kinematics. However, even coarse-grained methods have their limitations, one of which is timestep duration, which is limited by the maximum vibrational frequency in the coarse-grained system. To study this limitation, a coarse-grained model of polyethylene was created such that every C 2 H 4 unit was replaced with a bead. Coarse-grained potentials for bond-stretching, bond-bending, and non-bonded interaction were generated using the iterative Boltzmann inversion method, which matches coarse-grained distribution functions to atomistic distribution functions. After the creation of the model, the coarse-grained potentials were rescaled by a constant so that they were less stiff, decreasing the maximum vibrational frequency of the system. It is found that by diminishing the bond-stretching potential to 6.25% of its original value, the maximum stable timestep can be increased 85% over that of the unmodified potential functions. The results of this work suggest that it may be possible to simulate lengthy processes, such as the crystallization of polyethylene, in less time with adjusted coarse-grained potentials. Additionally, the large discrepancies in the speed of bond-stretching, bond-bending, and non- bonded interaction dynamics suggest that a multi-timestep method may be worth investigating in future work.