文档介绍：计算化学理论和应用 -第八讲
putational Chemistry laboratory
Beijing Normal university
分子几何结构优化(2)
Peter Pulay, "Analytical Derivative Methods in Quantum Chemistry", Adv. Chem. Phys. 69, 241 (1987) (Ab Initio Methods in Quantum Chemistry II, ed. . Lawley (Wiley, 1987))
. Schlegel, "Optimization of Equilibrium Geometries and Transition Structures" Adv. Chem. Phys., 67, 249 (1987). (Ab Initio Methods in Quantum Chemistry I, ed. . Lawley (Wiley 1987))
. Schlegel, "Geometry optimization on Potential Energy Surfaces" in Modern Electronic Structure Theory, ed. . Yarkony (World Scientific Press, 1995)
. Schlegel, “Geometry Optimization” in Encyclopedia putational Chemistry, ed. PvR Schleyer, NL Allinger, T Clark, J Gasteiger, P Kollman, HF Schaefer PR Schreiner, (Wiley, Chichester, 1998)
关于分子几何结构优化方面的重要文献
初猜分子几何结构和Hessian矩阵
计算分子能量和梯度
根据Hessian矩阵和梯度进行优化
更新Hessian矩阵
根据Hessian矩阵继续优化
根据力和位移判断收敛性
更新几何结构
是
完成
否
分子结构优化过程
Gaussian中的分子结构优化过程 与输出文件示例 乙烯分子优化
#P RHF/6-31G(d) Opt Test
Ethylene Geometry Optimization
0 1
C
C
H 1 CH 2 HCC
H 1 CH 2 HCC 3 180.
H 2 CH 1 HCC 3 180.
H 2 CH 1 HCC 4 180.
Variables:
CC=
CH=
HCC=
输入文件
优化作业执行过程
1/18=20,38=1/1,3;
2/9=110,17=6,18=5,40=1/2;
3/5=1,6=6,7=1,11=1,16=1,25=1,30=1/1,2,3;
4/7=1/1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7//1,2,3,16;
1/18=20/3(1);
99//99;
2/9=110/2;
3/5=1,6=6,7=1,11=1,16=1,25=1,30=1/1,2,3;
4/5=5,7=1,16=3/1;
5/5=2,38=5/2;
7//1,2,3,16;
1/18=20/3;
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
L101
L103
L202
L301
L302
L303
L401
L502
L601
L701
L702
L703
L716
L9999
优化开始的标志
第一次做 L103
初猜Hessian矩阵
优化结束的标志