GMXPolymer: a generated polymerization algorithm based on GROMACS

Jianchuan Liu,Haiyan Lin,Xun Li

2024-04-03

Abstract:This work introduces a method for generating generalized structures of amorphous polymers using simulated polymerization and molecular dynamics equilibration, with a particular focus on amorphous polymers. The techniques and algorithms used in this method are described in the main text, and example input scripts are provided for the GMXPolymer code, which is based on the GROMACS molecular dynamics package. To demonstrate the efficacy of our method, we apply it to different glassy polymers exhibiting varying degrees of functionality, polarity, and rigidity. The reliability of the method is validated by comparing simulation results to experimental data for various structural and thermal properties, all of which show excellent agreement.

Soft Condensed Matter,Computational Physics

What problem does this paper attempt to address?

This paper introduces a simulation aggregation algorithm called GMXPolymer based on the GROMACS molecular dynamics software, which aims to generate generalized structures of amorphous polymers. The research focuses on amorphous polymers, which have unique mechanical and thermoelectric properties due to their disordered molecular structure. Currently, understanding the structure-property relationship of amorphous polymers is a challenge because experimental studies are difficult to reveal their behavior at the molecular level.

In the paper, the authors propose a method for simulating polymerization and molecular dynamics equilibrium, which can generate representative structures of amorphous polymers. The algorithm includes initialization, polymerization steps, and iteration processes that involve the use of the GROMACS program to gradually form new chemical bonds and perform energy minimization and molecular dynamics simulations. The algorithm selects the nearest reactive atom pairs to form new bonds and adds relevant parameters in the ITP file.

To demonstrate the effectiveness of the method, the researchers apply it to different types of polymer systems and compare the simulation results with experimental data, showing good consistency. In addition, the paper also presents application examples in polymers such as polythiophene, trimethylsilyl chloride, and 5,5',6,6'-tetrahydroxy-3,3,3',3'-tetramethyl-1,1'-spirobiindane.

Finally, the paper generates different numbers of cross-linking units through simulation, and analyzes properties such as density and glass transition temperature (Tg) to verify the generality and accuracy of the method. The conclusion points out that the GMXPolymer code provides an effective tool for studying the structure generation of various amorphous polymers, contributing to a deeper understanding of their properties.