计算化学-4.从头计算法原理.ppt

袁淑萍国家重点实验室培育基地408室E-mail:spyuanqdu电话：85950690，13589208235,计算化学,1.闭壳层,分子轨道法的基本方程,计算化学,2.开壳层,RHF,UHF,ROHF,根据量子力学基本假定4，电子总能量表达式：,3.基于分子轨道的电子总能量表达式,4.LCAO-MO近似下的电子总能量表达式,定义“单电子哈密顿矩阵”h和“密度矩阵”D。,双电子作用矩阵G,LCAO-MO近似下的电子总能量表达式,为简化Fock矩阵元的表达，双电子积分常采用密度函数标记,在这种标记方式下，Fock矩阵元简记为,RoothaanEquation:Fc=Sc(2.2-3),闭壳层体系Roothaan方程,FC=SCE，whereE=,定义MO系数矩阵,Spin-UnrestrictedRoothaanEquation,UHFEquations:,因，UHF方程须采用交叉叠代法求解。每次叠代的计算量比RHF方程约增加1倍，其自洽叠代的收敛速度一般也较前者慢。,开壳层体系自旋非限制Roothaan方程,4.从头计算法原理及应用PrincipleandApplicationsofSCFAbInitioMethod,计算化学,从头计算的定义在采用MO理论在物理模型上的三个基本近似，以及LCAO-MO的数学近似后，Schrdinger方程具体表达形式为Hatree-Fock-Roothaan方程(简称为Roothaan方程)。在此基础上严格求解分子的Roothaan方程，以获得MO波函数及其能级，并利用波函数进一步计算分子的其它性质。这样的处理称为从头算法，又称“量子力学第一原理计算”（firstprinciplecalculation），或“Hatree-Fock计算”。从头算法是分子轨道法的主流，较之更精确的“Post-Hatree-Fock”处理，另外各种较粗略的半经验量子化学计算方法，都是由从头算法派生出来的。,从头算程序的结构与输入StructureandInputofAbInitioProgram,以Gaussian程序为例，说明程序的结构特点：由于从头算程序计算详细，功能复杂，一个程序达几十万条之长，为使用户一目了然，了解程序的功能结构，也为了程序编写扩充的方便，Gaussian程序从早期就采用了一种特殊的结构。即程序按几个主要方面分成几大块，每一块称为一层（Overlay），每层之内又分为许多小块，一小块称为Link。整个程序就象一座大楼，一层一层叠上去，每层里的一个个Link象一个个房间。例如自恰场迭代放在第5层，里面又分为闭壳层（RHF）、开壳层（UHF）、限制开壳层（ROHF）等分别位于Link501，502，504，可在不同计算中调用，若要补充其他功能，只需再编写一个Link补上即可。,Gaussian,Overlay1,Overlay9,10,11,99,L101,L102,L122,Overlay0,L0,L001,link0:初始化程序，控制overlaylink1:读入并处理RouteSection，建立要执行的link列表link9999:终止计算,Overlay99,L9999,Gaussian程序的结构,Gaussian程序这种结构，早期还为了克服计算机内存有限的缺点，在计算中可实行覆盖技术，后面执行的Overlay可覆盖前面计算过的Overlay，达到节省内存的目的。现在Gaussian公司提供二种版本：微机版（适合PC机，其中就采用覆盖技术）并行版（用于计算机集群系统、并行机）,下面是Link功能一览表：,Link功能一览表：,Link功能一览表：,Gaussian程序界面和输入文件的构造,Gaussian03程序界面,Preferences:对Gaussian程序进行初始化设置,自定义外部文字编辑器，用来打开.out文件link.exe所在的文件夹临时文件存放文件夹,Gaussian03图形工具栏,开始作业暂停当前作业当前link后暂停终止当前作业和批处理恢复当前作业在当前作业完成后终止批处理终止当前作业编辑或建立批处理,批处理专用,打开外部编辑器编辑输出文件,将.chk文件转换为.fch文件，这种文件可以使用图形软件打开将.fch文件还原为.chk文件从指定.chk文件中显示作业的routesection和title将.chk文件转换成文本格式将以前版本的Gaussian产生的.chk文件转换为G03的.chk文件,Gaussian03程序工具,编辑批处理作业文件转换不同格式的分子结构文件读取.fch文件中的数据并生成三维空间网格图利用.chk文件中的分子轨道，生成电子密度和静电势的空间分布网格图从.chk文件中打印出频率和热化学数据,NewZmat工具界面,H2O分子单点计算后的.chk文件转换成.fch文件后，在Chem3D中所作分子轨道(HOMO)，电荷密度和静电势空间分布图,Gaussian程序的输入文件,%Section设定作业运行的环境变量RouteSection设定作业的控制项Title作业题目电荷与自旋多重度MoleculeSpecification分子说明,%chk=water.chk%Section：行首以%开始，段后无空行%rwf=water.rwf#phf/6-31gscfcyc=250scfcon=8RouteSection：行首以#开始，段后加空行Watertitle：作业的简要描述，段后加空行01MolecularSpecification：O分子说明部分，段后通常加空行H1R1H1R12a1R1=1.04a1=104.0,Gaussian作业的格式,%Section(link0),定义计算过程中的临时文件%chk=name.chk.chk文件在计算中记录分子几何构型，分子轨道，力常数矩阵等信息%rwf=name.rwf.rwf文件主要在作业重起时使用，当计算量比较大时，.rwf文件通常会非常大，此时需要将之分割保存（%rwf=a,245MW,b,245MW,c,245MW,d,245MW,e,245MW,f,245MW）%int=name.int,%d2e=name.d2e.int文件在计算过程中存储双电子积分,.d2e文件在计算过程中存储双电子积分的二阶导数,内存使用控制%mem=n控制运行过程中使用内存的大小,可以以W或者MB，GB为单位default：6MW48MB（Gaussian程序使用的内存单位W是双精度字，相当于8字节)综合考虑到计算的需要和硬件水平，内存并非给得越多越好，最有效率的方法是根据作业类型估算所需要内存的大小,建立检查点文件（CHK）有几个重要功能：运行大分子几何优化等大作业时如意外原因（如停电、死机等）导致作业不正常中断，其后可从最接近中断点前最后一次记录处重新启动作业，不致前功尽弃。对正常结束的作业，主要的计算结果如优化后的分子几何参数、电子光谱、IR和Raman光谱的频率和强度等都可随时从保存的检查点文件中精确提取，并利用主窗口“Utilities”中选项的功能转换为其它软件可读的格式，以便进行图形显示和分析。有GVIEW的用户，可以直接将检查点文件转换为图，但检查点文件的存在是必须的）。,对同一分子进行新的项目或更高精度的计算（CI或MP）时，分子的几何参数可以直接从检查点文件中精确读入。,因此，建议用户从一开始就养成建立和充分利用检查点文件的习惯。但鉴于大分子作业的有名检查点文件可能很大，为节省硬盘空间，用户应经常整理Scratch，将有用的检查点文件用CD-RW或外接硬盘备份备份，把无用的删除。,RouteSection,RouteSection以#开始，#控制作业的输出#N正常输出；默认(没有计算时间的信息)#P输出更多信息。包括每一执行模块在开始和结束时与计算机系统有关的各种信息(包括执行时间数据，以及SCF计算的收敛信息)#T精简输出：只打印重要的信息和结果。RouteSection主要由方法，基组，任务类型三部分组成方法与基组后续课程专门介绍Gaussian程序能完成的任务类型：,Gaussian程序能完成的任务类型,RouteSection采用自由格式，大小写不敏感同一行不同项之间可以使用空格，逗号和“/”连接；例：#phf/6-31gscfcyc=230scfcon=8#p,hf/6-31g,scfcyc=230,scfcon=8关键词可以通过keyword=option；keyword(option)；keyword=(option1,option2,.)，keyword(option1,option2,.)指定例：#pHF/STO-3Gopt#pHF/STO-3Gopt=(TS,readfc),Title部分必须输入，但是程序并不执行，起标识和说明作用,TitleSection,Charge&Multipl.,输入分子的电荷和自旋多重度例：电荷多重度2s+1H2O01H3O+11NO02电荷多重度部分通常也算作分子说明,MolecularSpecification,分子说明部分主要用来定义分子核相对位置分子核相对位置可以用笛卡尔坐标，内坐标(Z-matrix)，或者是二者混合表示笛卡尔坐标是内坐标的一种特殊形式分子坐标的格式为:元素符号,x,y,z元素符号(n)原子1键长原子2键角原子3二面角格式代码(0,180)以原子在分子中的序数表示通过右手规则确定,几何优化用内坐标比较容易到达终点。,长度单位一般默认为，个别软件使用nm角度单位一般为度(),原子和分子在体系中的位置,分子的构建,描述,能量,笛卡儿坐标内坐标,分子结构数据的获得,简单分子,大分子(如蛋白质分子),通过标准键长、键角、二面角等数据直接写出通过分子图形软件构建,复杂分子,有单晶结构数据的，直接转化为分子结构数据采用软件构建，必须通过适当的构象搜索方法方可得到合理的结构；,其结构的构建本身就是专门的计算化学研究的课题，只有通过特定的方法(如同源模建等)方可得到可能的合理结构。,ENT，PDBBrookhavenProteinDataBankMOL(MolecularDesign,Ltd.MDL)HINHyperchemGJFGaussianINPGamessZMTMOPACZ-MatrixCARAccelrys,InsightIICSDCambridgeStructuralDatabaseCrystallographicfilesMOL2TriposMSFAccelrys,QUANTAMSIAccelrys,Cerius2CPDAccelrys,CatalystCHMqueryAccelrys,CatalystqueryCIFCrystallographicInformationfileXYZgenericgenericSKC,ISISMDL,ISISSMISESfileSMILESfile,常见结构数据类型,笛卡儿坐标输入,C-0.25-0.470.0.00O-0.250.960.0.00H0.77-0.830.00H-0.77-0.830.89H-0.77-0.83-0.89H-1.161.280.00,用文本编辑器直接输入笛卡儿坐标原子符号x坐标y坐标z坐标,C0.00.00.00.00.00.0O0.00.01.221.220.00.0H0.90.940.540.0H-0.940.540.0,原子符号可以是元素符号，也可以是原子序数(核电荷数)，也可以采用在元素符号后加上原子的序列编号。,C10.00.00.0O20.01.220.0H30.940.540.0H4-0.940.540.0,60.00.00.080.01.220.010.940.540.01-0.940.540.0,内坐标输入(Z-Matrix),当没有构建分子的图像软件可用时，与直角坐标相比使用内坐标要方便得多。对具有较高对称性的分子，输入内坐标Z-矩阵还可减少优化变量、节省计算时间并保证优化结果与分子的对称性严格相符。,CC11.54H11.092109.47H11.092109.473-120.H11.092109.473120.H21.091109.473180.H21.091109.476-120.H21.091109.4761120.,作业,写出分子H2C=C=O分子输入坐标其中C=C:1.35C=O:1.20C-H:1.09写出苯分子内坐标C=C:1.30C-H:1.09,