主权项 |
1.一种基于约束正则模式的复杂函数极小值搜索方法,其特征在于,利用输入的原子坐标q<sub>i</sub>, i=1,2,…n,已知的势能面能量函数E=<i>f</i>(q<sub>i</sub>), i=1,2,…n,能量函数对坐标的一阶导数F=<i>g</i>(q<sub>i</sub>), i=1,2,…n,求解能量所有极小值对应的原子坐标,以及连接这些极小值的反应通道,具体步骤如下:步骤一、输入一个能量极小值坐标A,并随机产生一个模式<b>N0</b>;步骤二、在随机产生的模式<b>N0</b>基础上,通过约束双子方法优化得到一个约束的正则模式<b>N</b>;步骤三、在原有的势函数基础上,新添加该N方向的一个偏置高斯势函数,则总势能面<i>V</i><sub><i>tot</i></sub>可以表示为真实势能面<i>V</i><sub><i>real</i></sub>和偏置势函数<i>V</i><sub><i>G</i></sub>的和,即<img file="964079DEST_PATH_IMAGE001.GIF" wi="96" he="25" />,沿着该<b>N</b>方向外推当前坐标,到由偏置势函数和原有势函数叠加产生的能量极小值附近;步骤四、能量极小值优化当前坐标;步骤五、判定是否达到最多的偏置高斯势函数个数(NG)的上限,或达到预设的能量阈值E<sub>r</sub> ;达到则转到步骤六,否则重复步骤二、步骤三、步骤四;步骤六、去掉所有添加的偏置高斯势函数,能量极小值优化得到新的坐标体系B;步骤七、采用Metropolis 蒙特卡洛方法,判定B是否接受;如果接受,A被B替代;步骤八、重复步骤一到步骤七,寻找下一个能量极小值,或程序结束。 |