本发明专利技术公开了一种利用遗传算法优化反应过渡态的方法，属于计算化学领域。本发明专利技术方法通过输入反应物的分子式自动转化为结构式，先设置目标原子的坐标为变量，随机生成多组初值，先通过常规算法计算得到满足条件的初值，在针对得到的初值利用利用遗传算法自动变异和交叉组合得到新的初值，重复计算得到符合条件的初值，直至符合条件的初值计算完毕，再根据满足条件的初值对应的反应路径，最终选出所需的过渡态及反应路径。本发明专利技术无需人为去调分子结构，计算过程快且结果更加准确。计算过程快且结果更加准确。计算过程快且结果更加准确。
【技术实现步骤摘要】
一种利用遗传算法优化反应过渡态的方法
[0001]本专利技术涉及一种利用遗传算法优化反应过渡态的方法，属于计算化学领域。
技术介绍
[0002]在分子构象转变或者化学反应中往往都会存在过渡态，过渡态结构对应着势能面上反应路径中能量最高点，通过最小能量路径连接反应物和产物。确定过渡态有助于了解反应机理，计算能垒推算反应速率。
[0003]但在化学反应过程中，过渡态不稳定，存在时间极短，无法通过实验的方法获得，计算化学方法是目前寻找反应过渡态最好的办法，但目前在量子化学计算上仍有一些困难。目前寻找化学反应过渡态主要有3种方法：TS、QST2以及QST3，目前TS法是常用方法。在高斯中选用opt＝ts方法来寻找过渡态，需要提供一个初猜结构，初猜结构对优化结果有很大影响，往往会存在如下几个问题：1、初猜结构错误；2、初猜结构不是本反应的过渡态；3、初猜结构的能量不是本反应路径中过渡态能量最小值。初猜结构的不准确容易导致寻找过渡态的失败，且初猜不正确往往需要重新改变初猜结构，耗时长且不一定能够获得最优过渡态。
[0004]遗传算法是模拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型，是一种通过模拟自然进化过程搜索最优解的方法。即“适者生存，不适者淘汰”。主要过程：从代表问题可能潜在的解集的一个种群(population)开始，种群又是由个体(individual)组成的，每个个体都是经由基因(gene)编码的，所以个体实际上是染色体chromosome)带有特征的实体。染色体是物质遗传的载体，多种基因的集合才决定了外在的表现。因此，一开始需要将所需的表现型映射到内容控制的基因型即编码工作。由于仿基因编码比较繁琐，所以进行简化采用二进制编码，初代种群产生之后，按照适者生存和优胜劣汰的原理，逐代(generation)演化产生出越来越好的近似解，在每一代，根据问题域中个体的适应度(fitness)大小选择(selection)个体，并借助于自然遗传学的遗传算子(genetic operators)进行组合交叉(crossover)和变异(mutation)，产生出代表新的解集的种群。这个过程将导致种群像自然进化一样的后生代种群比前代更加适应于环境，末代种群中的最优个体经过解码(decoding)，可以作为问题近似最优解。因此，寻找一种更简便且准确度更高的计算化学反应过渡态的方法对于化学反应过渡态的确定以及对后续化学工业的指导具有重要意义。
技术实现思路
[0005]针对现有的化学计算反应过渡态的过程中存在的初猜结构不易确定、随机性较强，且寻找反应过渡态的过程耗时长，计算结果往往错误或者不是理想结果，以及计算周期长的问题，本专利技术提供了一种利用遗传算法优化反应过渡态的方法。本专利技术方法通过输入反应物的分子式自动转化为结构式再利用遗传算法寻找出反应物可能的过渡态，计算出对应的反应路径，最终选出所需的过渡态及反应路径。本专利技术无需人为去调分子结构，计算过
程快且结果更加准确。
[0006]具体的，本专利技术的技术方案为，一种利用遗传算法优化反应过渡态的方法，所述方法包括：
[0007](1)输入反应物的分子式，通过RDKit将分子式转化为结构式并输出对应的.mol文件；
[0008](2)利用Multiwfn将步骤(1)得到的.mol文件转化为.gjf文件；
[0009](3)将.gjf文件中的目标原子的坐标设为变量，同时，设定坐标变量的范围，并随机生成多组目标原子坐标的初值；
[0010](4)输入目标原子坐标变量的初值，利用基组算法进行计算，计算得到.out文件，进行判断：当结果非正常结束，则删除结果并代入其他初值进行迭代计算；当结果为正常结束，但是频率计算结果中无虚频或者有超过一个虚频，则删除结果并代入其他初值进行迭代计算；当结果为正常结束，且频率计算结果中有且只有一个虚频，则将.out文件通过文件的读写方式写入并计算该.out文件中对应的过渡态所对应的反应路径，之后代入其他初值进行迭代计算；
[0011](5)初值迭代完成后，整理步骤(4)计算得到的符合条件的初值，并利用遗传算法自动变异和交叉组合得到新的初值，重复步骤(4)和(5)，直至将符合条件的初值计算完毕；
[0012](6)迭代完成后，整理所有计算得到的过渡态及其对应的反应路径，根据反应产物来判断和确定符合条件的过渡态。
[0013]在本专利技术的一种实施方式中，步骤(1)中，所述输入反应物的分子式应根据RDKit的规则进行输入，其中，所述RDKit为一款开源化学信息学软件。
[0014]在本专利技术的一种实施方式中，所述Multiwfn是一款量子化学波函数分析程序。
[0015]在本专利技术的一种实施方式中，步骤(2)中，需要根据实际情况将选用的基组以及注意事项(例如溶剂化效应)写入.gjf文件，与传统计算方法的输入一致。
[0016]在本专利技术的一种实施方式中，步骤(4)中基组算法为：只要能在高斯软件上计算的算法在此步均能写入文件并调用高斯软件；包括半经验法、从头算法、密度泛函方法、MPn方法中的任一种。
[0017]在本专利技术的一种实施方式中，所述坐标变量的范围优选为其余三个原子坐标范围的最大值或更大的范围，即其余三个原子的坐标的x、y、z的最大值和最小值分别作为坐标变量的x、y、z坐标的最大值和最小值，即获得目标原子的坐标范围，或大于前述坐标范围的坐标范围。
[0018]在本专利技术的一种实施方式中，所述坐标变量的初值是在坐标变量的范围内随机生成的目标原子的坐标值。
[0019]在本专利技术的一种实施方式中，步骤(3)中所述的多组初值可以根据实际情况选定，例如100组、200组、300组等。
[0020]在本专利技术的一种实施方式中，步骤(5)中所述遗传算法自动变异、交叉组合得到新的初值，即通过设定交叉率和变异率的参数值，然后利用遗传算法对符合条件的初值进行变异交叉组合，得到新的初值。
[0021]在本专利技术的一种实施方式中，所述交叉率为0.6～1，优选0.8。
[0022]在本专利技术的一种实施方式中，所述变异率小于0.1，优选0.003。
[0023]在本专利技术的一种实施方式中，步骤(4)中，当结果为正常结束，且有且只有一个虚频时，将此时得到的.out文件中的能量值与之前得到的所有的满足条件的.out文件中的能量值进行比对，当其与之前的能量值差值在
±
0.00002Hatree范围之内时(其中，1Hartree＝2625.5kJ mol-1
)，则将此次获得的能量值进行删除，且无需计算此时对应的反应路径；删除的原因是由于二者能量值差别过小，对应的反应路径一致。
[0024]在本专利技术的一种实施方式中，所述反应路径的计算方式为常规的反应路径的计算方式。
[0025]在本专利技术的一种实施方式中，迭代过程中的迭代次数可以根据实际需要进行自行设定。
[0026]本专利技术还提供了上述利用遗传算法优化反应过渡态本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种利用遗传算法优化反应过渡态的方法，其特征在于，所述方法包括：(1)输入反应物的分子式，通过RDKit将分子式转化为结构式并输出对应的.mol文件；(2)利用Multiwfn将步骤(1)得到的.mol文件转化为.gjf文件；(3)将.gjf文件中的目标原子的坐标设为变量，同时，设定坐标变量的范围，并随机生成多组目标原子坐标的初值；(4)输入目标原子坐标变量的初值，利用基组算法进行计算，计算得到.out文件，进行判断：当结果非正常结束，则删除结果并代入其他初值进行迭代计算；当结果为正常结束，但是频率计算结果中无虚频或者有超过一个虚频，则删除结果并代入其他初值进行迭代计算；当结果为正常结束，且频率计算结果中有且只有一个虚频，则将.out文件通过文件的读写方式写入并计算该.out文件中对应的过渡态所对应的反应路径，之后代入其他初值进行迭代计算；(5)初值迭代完成后，整理步骤(4)计算得到的符合条件的初值，并利用遗传算法自动变异和交叉组合得到新的初值，重复步骤(4)和(5)，直至将符合条件的初值计算完毕；(6)迭代完成后，整理所有计算得到的过渡态及其对应的反应路径，根据反应产物来判断和确定符合条件的过渡态。2.根据权利要求1所述的一种利用遗传算法优化反应过渡态的方法，其特征在于，步骤(4)中所述基组算法包括半经验法、从头算法、密度泛函方法、...
【专利技术属性】
技术研发人员：王海军，胡爱云，
申请(专利权)人：江南大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。