Multiple Objective NSGA-II-Based Optimization Program and Its Application in Reactive Force Field for 2,4,6-Trinitrotoluene Diffusion in the Aqueous Phase

Guan Zhang,Jin Li,Zongkuan Liu
DOI: https://doi.org/10.1021/acs.jpcc.9b05036
2019-01-01
Abstract:The reactive force field (ReaxFF) molecular dynamics simulation is an effective method to study the behavior of energetic materials in the aqueous phase. In order to optimize the ReaxFF of energetic materials in the aqueous phase, we have developed a multiple objective nondominated sorting genetic algorithm II (NSGA-II)-based optimization program (MONOP). The MONOP was used by adopting the CHON-2017_weak as the initial force field to re-optimize the force field of a 2,4,6 trinitrotoluene (TNT)/water system. In addition, the ReaxFF-lg (low-gradient dispersion correction), CHON-2017_weak, and TNT/water were used to simulate TNT crystals correspondingly. Ultimately, the TNT diffusion coefficient in the aqueous phase was also calculated by the mean squared displacement (MSD) method. The optimization results showed that the MONOP optimization has a fast global convergence, in about 11 generations in this study. By comparing the lattice constants a, b, c, density, and the radial distribution functions of the TNT element pairs, the TNT/water system was found to describe TNT crystals well. Relatively, the TNT/water system had a slightly larger value of liquid water density than that of the experimental measurement. The MSD results showed that the migration rate of TNT molecules in water is much slower than that of OH anions in water; thus, it is reasonable to conclude that the TNT hydrolysis reaction is dominated by OH anions attacking TNT molecules, and the reaction mostly occurs at the liquid-solid interface. This study indicates that the MONOP optimizes the ReaxFF and has broad application prospects; the diffusion coefficient of TNT in water can provide theoretical support for the aqueous phase degradation of energetic materials.
What problem does this paper attempt to address?