[发明专利]预测水相中有机化合物与羟基自由基反应速率常数的方法

权利要求书

1.一种预测水相中有机化合物与羟基自由基反应速率常数的方法，其特征在于，

首先收集526种有机化合物的水相k_OH数值，将其按照4:1随机拆分为训练集和验证集；训练集中的421种化合物用于构建模型，验证集中的105种化合物用于评估模型的外部预测能力；对上述526种有机化合物的分子结构进行优化，得到相应有机化合物的稳定构型并选取9个量子化学描述符；通过上述稳定构型得到2418个Dragon描述符；采用MLR回归分析方法筛选分子描述符和构建QSAR模型；

筛选出的最优QSAR模型如下：

logk_OH＝11.566+6.233×E_HOMO-0.074×HATS2s-0.183×Mor23u+0.238×GATS1e-0.099×N-075+0.107×nR＝Cp-0.230×nRCONH2-0.070×C-001+0.080×MLOGP+0.113×nS-0.265×nBR+0.651×q_H⁺+0.119×Eig03_EA(dm)；

其中，E_HOMO表示最高占据分子轨道能量，HATS2s是与内蕴状态相关的GETAWAY描述符，Mor23u表示未加权的3D分子结构描述符，GATS1e表示Sanderson电负性加权的lag 1的Geary自相关指数，N-075表示苯环上的N原子或与O、N、S、卤素等电负性原子形成离域键的N原子碎片数，nR＝Cp表示末端sp2杂化的主碳数目，nRCONH2表示分子中含RCONH₂结构的数目，C-001表示分子中-CH₃/CH₄结构信息，MLOGP表示Moriguchi辛醇-水分配系数，nS表示分子中含S原子的数目，nBR表示分子中含Br原子的数目，q_H⁺表示H原子最正净电荷，Eig03_EA(dm)表示偶极矩加权的本征值的边界邻接指数。

2.根据权利要求1所述的方法，其特征在于，所述有机化合物为烷烃类化合物、烯烃类化合物、炔烃类化合物、芳香烃类化合物、醇类化合物、醛类化合物、酮类化合物、醚类化合物、酸类化合物、脂类化合物、卤代类化合物、含氮化合物、含硫化合物或含磷化合物。