高斯软件的介绍摘要.ppt

题 目：高斯软件的简介指导老师：徐畅老师姓 名：马乃宇学 号：专 业：高分子化学与物理专业 Gaussian是一个量子化学软件包，它是目前应用最广泛的计算化学软件之一，Gaussian软件的出现降低了量子化学计算的门槛，使得从头计算方法可以广泛使用，从而极大地推动了其在方法学上的进展。 (3) G98输出结果: 下面以H2O能量计算的输出结果为例: L1输出G98版权信息 %部分内容 关键词部分 计算所需调用的模块 Title内容 电荷, S以及分子构型 各原子的直角坐标 距离矩阵 需核查! 注意:计算结果是以该坐标系为准 需核查! MO初始猜测 自洽场迭代求解部分 收敛指标 迭代次数 (通常接近2.0) (布居分析结果) MO对称性及能级(a.u.) Mulliken键级和电荷 偶极矩及多极矩 对主要计算结果进行总结 (4) 能量计算的几个问题： 基组的选取： 理论上而言，计算中所选取的基组越大，计算结果越准确，但由 于受到硬件上的限制，需根据实际情况选择基组。此外，对于环状 分子或存在?作用的体系，通常需要考虑使用极化基组；对于带有 较多电荷的体系，或考察弱相互作用的体系，通常需考虑使用弥 散基组； 以O2为例(O-O为1.208A)： 6-31G -150.26755a.u. 6-31G* -150.31998a.u. dE=0.05243a.u. 6-31+G -150.27678a.u. dE=0.00923a.u. 由于O2中存在?作用，故极化函数的影响要比弥散函数来得大。 以O22-为例： 6-31G -149.78721a.u. 6-31G* -149.82396a.u. dE=0.03675a.u. 6-31+G -149.98117a.u. dE=0.19396a.u. 6-31+G* -150.02156a.u. dE=0.23435a.u. 此时，因带较多负电荷，故弥散函数的影响较极化函数来得大。 体系多重度的选取： 对于复杂体系，特别是含有过渡金属原子的体系，多重度的选取 凭经验不易得到，此时，可通过比较不同多重度时体系的能量来 确定基态能量。 以O2为例(基组为6-31G*)： S=1 -150.25734a.u. S=3 -150.31998a.u. S=5 -149.79200a.u. S=7 -148.80784a.u. 由上结果可见，O2的基态应为三重态。 键能的计算： 以O2为例(O-O键长为1.208A，计算O—O键能： 步骤如下： 首先确定O2基态的能量： 假设基组为6-31G*，则总能量为： -150.31998a.u.(基态为：3?g) 其次，计算O原子基态能量，结果为：-75.06061a.u. 则O2中O—O键的键能为： 2?-75.06061+150.31998 =0.19876a.u.=521.8KJ/mol 该结果略高于实验值（494KJ/mol）。 计算机在化学中的应用 课程作业 G98的功能和程序结构 1. 主要功能： 分子构型的优化 基态(Ground state) 激发态(Excited state) 反应过渡态(Transition state) 能量计算 基态和激发态能量 化学键的键能 电子亲合能和电离能 化学反应途径和势能面 光谱计算 IR光谱 Raman光谱 电子光谱 NMR 其它功能 电荷和电荷密度 偶极矩和超极矩 热力学参数 适用体系：气相和溶液 2.程序结构： 常用模块的功能： L0—初始化模块； L1—读入输入，根据所给关键词确定将要使用的模块； L101,102,…—与构型优化和反应过渡态相关的模块； L202—输出距离矩阵、判断化合物点群及确定新的坐标系； L301,302…309—与基组和赝势有关模块； L310,…319—计算单及双电子积分模块； L401,402—SCF初始猜测模块； L502,503,508—SCF模块； L601,608—Mulliken布居以及自然键轨道分析模块； L701,702…—计算能量一阶和二阶导数模块； L9999—进程结束模块； 说明： 1.根据不同的任务，某些模块需重复调用多次； 2.通常耗时较多的模块有：L5,L7,L8,L9,L10,L11等，此外， L8~L11这些模块的执行对内存和硬盘的需求较大； 3.若L9999未能正常执行完毕，则表明计算过程存在问题，需 检查之； 4.可根据各个模块的功能，对g98程序进行简化，例如如果用 户通常只用g98进行能量计算，则可只保留L1~6