分子的几何构型优化计算文档格式.docx

分子的几何构型优化计算文档格式.docx

《分子的几何构型优化计算文档格式.docx》由会员分享，可在线阅读，更多相关《分子的几何构型优化计算文档格式.docx（10页珍藏版）》请在冰豆网上搜索。

＞计算体系设置

2.1势能面

在不分解的前提下，分子可以有很多个可能的构型，每个构型都有一个能量值，所有这些可能的结构所对应的能量值的图形表示就是一个势能面，势能面描述的是分子结构和其能量之间的关系，以能量和坐标作图。

根据分子中的原子数和相互作用形式，有可能是二维的，也有可能是多维的。

势能面上的每一个点对应一个具有一个能量的结构。

能量最低的点叫全局最小点，局域最小点是在势能面上某一区域内能量最小的点，一般对应着可能存在的异构体。

鞍点是势能面上在一个方向有极大值而在其他方向上有极小值的点，通常对应的都是过渡态。

优化的目的就是找到势能面上的最小点，因为这个点所对应的构型能量最低，是最稳定的。

22确定能量最小值

构型优化就是找体系的最小点或鞍点。

能量的一阶导（也就是梯度，注意在数学中，一阶导表示着函数的变化趋势，一阶导为零就表明找到了极值点，这是确定最小值的数学基础）是零，这表明在这个点上的力也是零（因为梯度的负值是力）。

我们把势能面上这样的点称为静态点（也就是上面所说的极小点）。

所有成功的优化都会找到一个静态点，虽然有时找到的静态点并不是想要的静态点。

程序从输入的分子构型开始沿势能面进行优化计算，其目的是要找到一个梯度为零的点。

计算过程中，程序根据上一个点的能量和梯度来确定下一步计算的方向和步幅。

梯度其实就是我们所说的斜率，表示从当前点开始能量下降最快的方向。

以这种方式，程序始终沿能量下降最快的方向进行计算，只至找到梯度为零的点。

而梯度为零表明能量已是最小，所以这个点就是我们所要找的具有最小能量的结构。

很多程序还可以计算能量的二阶导，所以很多和能量的二阶导相关的性质（如频率计算）也可以得到。

2.3计算收敛的标准

优化计算不能无限制的进行下去，判定是否可以结束优化计算的判据就是我们现在所要了解的收敛标准。

需要特别强调的是这个标准规定的是两个SCF计算结果间的差距，当计算出的这两个能量值的差落在程序默认的标准之内时，程序就认为收敛达到，优化结束。

而单点能计算中也有一个收敛标准，这个收敛标准是用于判定CF计算是否完成。

SCF计算是一个迭代过程，假定一个解，带入到方程中，求出一个解，再将这个解带入到方程中，如此循环，直到两次解的差落在程序默认的范围之内，SCF计算完成。

在这里强调这二者区别的原因是：

优化是高斯计算中最易出错的地方，有时OPT=TIGH■可以帮助我们解决这个问题，所以我们要注意这个命令和SCF=TIGHE勺区别。

四个收敛标准：

（真正的0难以达到，程序给出了

个判断标准）

Item

Value

Thtreshold

Converged?

Maximujn.

Force

0.001037

0.000450

NO

RMS

0.000132

0.000300

YES

Masimujn

Displacejnerit

0.116214

0.001800

RHS

Displacement

Qt01S2S3

0.001200

MaximwnForce:

力的收斂标准是0,00045jRMSF皿口：

力的均方根的收敛标准0.0003

MaximumDisplacement位移的收數标准：

0.0018

RMSDisplacement&

根的收敛标淮0.0012*

在优化过程中，有时会出现只有前两项收敛（YES表示已收敛，NO表示不收敛），这种结

果是可以接受的。

这是因为高斯程序默认：

当计算所得的力已比收敛指标小两个数量极

时，即使nt值仍大于收敛指标，也认为整个计算已收敛。

这种情况对大分子

（具有较平缓的势能面）比较常见。

3.基本输入格式和输出解释

ExampleEL1EthyleneOptmiazation

EthyleneCeojnetiyOptijaization

01

C

C1CC

H1CH2HCC

H1CH2HO23180.

H2CH1HCC3180.

H2CH1HOC4180.

Variables：

CC=1.31

CH=1.07

HCC=12L.5

输出解释

GradG-radGxadGradGradGradGradGradGzadGxadGradGradGradGradGzadGradGradGrad

优化计算的分隔符

Bernyoptijnization^

InternalFaices：

Mgx0.173679522RMS0,052518600

Searchforalocalminijiim・（优化目的是寻找极小値』对于过渡态是寻找鞍点）

Stepnumber1outofamaximumof197（197表示程序将对构型进行197次优化彳、、

注意这是程序默认要计算的次数，优化计算有可能提前完成也有可能在默认的次数内不能完成计算。

如是后者，通常是用iew打开输出文件，这时所得的结构对应着输出文件中第7

次计算的结果，在这个结构的基础上再接着进行优化计算；

1表示这是第一次优化计算）

Variable

OldX

-DE/DX

DeltaX

DeltaK

NewX

（Linear）

（Quad）

（Total）

R1

3.67961

-0.00491

0.00000

-0.04257

3.53704

R2

3.15269

-0,01965

-0.04830

-0.04830

3.10439

R3

2.35834

-0.01475

-0.02333

-0.02338

2.33496

Threshold

Maximujii

0,173680

RJIS

0.052519

0.000300

恥

MaKimum

CL302996

0.001800

EMS

0.089161

OldX结构旧的变量值,

四个收敛标准，NO表示还未收敛。

YES表示已收敛。

在这一项输出之后给出计算所得的分子结构参数。

然后是一个单点能计算，该计算会给该结构的能量（以ee为单位）。

收

敛后会有以下的输出。

Thxeshald

Maxinum

0.000022

0.000003

0.001476

0.000307

PredictedchangeinEnerg^-3^918833D-09

Optimizationcompleted,表乎1尤化结束

Stationarypointfound・

以下输出是以内坐标形式给出优化好的结构。

从中可得到需要的参数（键长，键角和二面角）

JNane

Definition

Derivativf

3Info*

■

!

Fl

R（l,23）

1.8726

-de/dx=

0.0

V

R2

R（2,26）

1.6244

-de/dx二

R3

R（3,15）

1.2297

-DE/DX=

JR4

E（4,5）

1.4038

-DE/DX=

以直角坐标表示的最终优化结枸

Standardorien丈盘1:

:

10]1：

CenterAtcunicAtomic

Coordinates（Angstroms）

Number

Humbex

Type

K

Y

Z

35

-5.62446S

0.619069

0.101S77

2

16

S.513621

-4.S11162

-0.440263

3

8

0,810636

1.135321

-1.663476

4

7

2.529167

1,684744

-0.074011

5

2.856067

1.017185

1.167681

0.41]0512-2.

646709P.

105241

剩下的输出的部分是优化结构的布局分析（用POP=FULl命令会有详细的输出），分子轨道,原子电荷和偶级距。

4高斯中自带的练习（通常都在ample文件夹和rcise文件夹中）

se3.1异构体优化练习

se3.2异构体优化练习考察结构和能量的关系

se3.3优化练习，考察取代基对分子键长，电荷等的影响

se3.4优化练习在这个练习中用了CC这个命令助其收敛。

其意思是：

SCF

计算一开始就使用正常级别的积分精度。

通常情况下是先用中等积分精度粗算，稍后再转换为正常积分精度计算。

se3.5优化核磁计算练习。

先优化，算核磁；

优化TMS算核磁。

两个值相减后就可以得到和实验值相比较的结果。

计算时通常都是先优化，然后在更高的计算方法上算性质（如核磁）。

但要注意：

频率计算的优化和频率计算必须在同一方法上进行，否则计算结果无意义。

本练习是LYP/6-31G（D）的基础上，用HF/6-311+G（2D,P）方法算核磁。

Exese3.6C60优化。

记得在做练习时，报错，如是，在命令行中加NOSYM即卩可。

因为体系默认在优化时保持对称性。

由于所算体系对称性非常高，优化过程中结构稍有变化，对称性就被破坏。

程序就认为不满足限制条件，计算终止。

加后，去掉对称性限制,允许

结构的变化，所以可计算。

se3.7过渡态优化计算。

这个计算会报错。

这个命令要求在每一个输入

的结构后都有一个指定特定输出的内容[也就是输入内容中的4.5]

#TRHF/S-31G（d）Opt=（QST2,AddRedundant）Test

SiH2+H2—>

SiH4Reactants

0,1

Si

X11,0

H11.48255.0

H11.48255.03180.0

H1R2Al390.0

H1R5A22180.0

A1=SO,0

A2=22,O

这个垃蚤少了45（见下所以计算会报错.加上后就可以JG

S1H2+H2—>

S1H4Products

X11.0

H11.4S256.0

H11.45255.03130.0

H1R5A22130.0

E*48

Al=125.2

A2=109.5

进行优化。

如想用矩阵坐标或直角坐标优化，可用OP瑞令实现，参见OPT关键词说明。

Exse3.8比较矩阵坐标，直角坐标和内坐标的优化优势。

现在普遍认为内坐标在优化方面比矩阵坐标，直角坐标有优势。

所以在默认情况下程序会将输入的坐标自动转化为内坐标形式进行优化。

如想用矩阵坐标或直角坐标优化，可用OPT命令实现，参见OPT关键词说明。

5当优化遇阻时常用的解决办法

1）看所给的初始构型是否合理，这是初学者最易犯的错误。

解决办法：

检查初始构型的空

间构型，然后先用半经验方法或小基组（如STO-3G优化，然后再用大基组优化。

2）一般的问题用1的方法就可以解决，如还不行可以用OPT命令增大循环次数，减少步长和提高收敛精度来解决。

详见G98或G03手册PT关键词。

注意报错其实不是错误，而

只是在程序默认的次数内未完成优化任务，用LE命令增大循环次数即可解决这

类报错。

3）对于过渡态优化，由于分子构型需要手动调整，所以更难给出合理的初始构型，通常在命令行里加，该命令意思是在优化前先进行一个频率计算以获得用于指明优化

方向的力常数。

这个方法也可用于基态难以完成的优化。

对于基态的优化，还可用

命令来获得力常数，前提是在低一级别的计算水平上作了频率计算，且保留了检

查点文件。

L,最无奈的办法，但及耗时，对于我们所要处理的四五十个原子的体系，

用单机进行这样的计算太费时间了，已不具有任何实际的意义。

这个命令的意思是在每一步

优化前都要做一个频率计算获得指明下一步计算的力常数。

5•当计算因外因（如停电等）意外中止时，如果保留有检查点文件可用RESTAT命令继续这个计算，用ck命令后，分子的电荷，多重度和结构说明部分都不需要，如果

是用om=check命令，则需要有电荷和多重度的说明。

详细格式如下：

%cht=E：

\GO3W\Scratch\optqst3t3ng»

chk

^ftnproc=l

Krwf=匕:

VL,e：

\2,eA3,2gb,-1

#opt=（qst3,xestaxtjnoeig&

n）rb31yp^0_31g（d）nosyunguess=readaeom=allcheck

这是一个关于过渡态优化因停电终止而的例子。

t和

ck命令必须要连用。

rwf文件是用来指定临时交换文件空间大小的，通常情况下程序默认在RATCH文件中产生一个约2G的临时文件，对一般计算来说，这个存放临时文

件的空间够了。

但对于大的频率计算，大的从头算法（如MP2等）和CIS计算等一些会产生

大的临时交换数据的算法来说，可能会因F文件空间不够而报错，这时可通过这个命令格

式加大临时数据的存放空间，据说对于单机版格式的分区，程序最大可允许设定8个

临时文件，总空间是G24G（具体视操作系统而定，据说XP系统可达24G,但W9系统只

能达到G；

对于NTS格式，可产生的单个文件约为，总共可产生2G的临时文件空间。

［注：

关于RWF的大多介绍来自大话西游，量子化学研究小组和厦大论坛，自己未试过］Geom=（Check,Step=n）这个命令要求从已停止的优化计算的第步开始重新优化。

一般在失败的优化中有一步非常接近收敛值，但在这一步后又偏移了的情况下用这一命令，有时还要加上促使其收敛。

6实验内容

练习高斯中自带的几何构型优化例题，熟悉几何构型优化的基本方法

对乙烯醇异构体（或自己感兴趣的简单分子）进行几何构型优化，将理论计算结果与

文献实验结果分析比较，初步了优化方法的特点。