基于ReaxFF力场的对硝基苯酚臭氧氧化分子动力学模拟

Molecular Dynamics Simulation of Ozonation of p-Nitrophenol at Room Temperature with ReaxFF Force Field

本文采用基于ReaxFF反应力场的分子动力学方法(ReaxFF MD),利用自主研发的国际首个基于GPU加速的ReaxFF MD程序系统GMD-Reax和独特的化学反应分析工具VARx MD,探索臭氧氧化对硝基苯酚的反应机理。通过模拟考察了300 K恒温条件下臭氧氧化水中对硝基苯酚的行为,获得了酚结构开环、CO_2生成、主要自由基(·OH、·O_2、·O)及团簇型自由基的数量演变趋势,并可定性描述六元环开环和CO_2生成均遵循伪一级反应动力学规律。反应机理分析表明酚类分子在水溶液中被臭氧氧化的路径主要经过攫氢、六元环开环、碳链的氧化分解三个阶段,也揭示了自由基和团簇型自由基在臭氧降解对硝基苯酚时所发挥的重要作用。本工作是应用ReaxFF MD分子模拟方法对常温水环境下臭氧降解酚类污染物反应机理研究的一个尝试,可为深入认识该机理及相关的实验、理论研究提供一定的参考。

Understanding the reaction mechanism of phenol ozonation in coking wastewater is very important for industrial applications of the ozonation process. Ozonation of p-nitrophenol in water at 300 K was simulated by ReaxFF force field molecular dynamics （ReaxFF MD） employing the GPU-enabled high-performance code of GMD-Reax and a unique code of VARxMD developed in authors＇ group. Evolution trends of aromatic ring opening, CO2 generation, dominant radicals （.OH, .O2, O）, and H20 clusters were obtained. The simulated CO2 generation and reduction of aromatic ring could be described with pseudo-first-order kinetics. Moreover, the reaction pathways for the ozonation of p-nitrophenol can be divided into three stages： hydrogen abstraction, opening of the six-membered-ring structure, and the breaking of C-C bonds. The simulations revealed the important role of radicals and water clusters in the ozonation of p-nitrophenol. This work is an attempt to investigate the ozonation mechanism of phenols in aqueous solutions at room temperature using ReaxFF MD, which should be helpful for further experimental or theoretical investigation of the mechanism.