Guassian基本概念和用法

更新时间:2023-05-10 03:02:11 阅读: 评论:0

1. Guassian基本概念和用法 

Gaussian计算流程
210.34.15.126/i?forum=3&topic=1729

Gaussian菜鸟常见问题分析

210.34.15.126/i?forum=3&topic=2061

一些容易忽略的地方

210.34.15.126/i?forum=3&topic=1099

Gaussian计算出错和对应的解决方法 

210.34.15.126/i?forum=3&topic=1055

GAUSSIAN输出文件中出错信息大汇集

210.34.15.126/i?forum=3&topic=460

Gaussian98高级注释

210.34.15.126/i?forum=3&topic=299

A Brief Introduction to Gaussian

210.34.15.126/i?forum=3&topic=302

有关Gaussian 98的补充


210.34.15.126/i?forum=3&topic=33

Gaussian98 tips

210.34.15.126/i?forum=3&topic=31

How to start and run a simple calculation with Gaussian

210.34.15.126/i?forum=3&topic=30


Gaussian软件应用第一到第十章

210.34.15.126/i?action=list&forum=3&show=75


Gaussian98使用经验谈14

210.34.15.126/i?action=list&forum=3&show=100

Gaussian 98 – some hints for doing calculations  

210.34.15.126/i?forum=3&topic=300
Gaussian计算流程
 
此主题相关图片如下:
基本优化过程
Gaussian菜鸟常见问题分析
1 检查是否有初始文件错误
在命令行中加入 %kJob L301 or %kJob L302
如果通过则一般初始文件ok
常见初级错误:
a. 自旋多重度错误
b. 变量赋值为整数
c. 变量没有赋值或多重赋值
d. 键角小于等于0度,大于等于180
e. 分子描述后面没有空行
f. 二面角判断错误,造成两个原子距离过近
g. 分子描述一行内两次参考同一原子,或参考原子共线
2 SCF(自洽场)不收敛
则一般是L502错误
省却情况做64cycle迭代(G03缺省128 cycles)
a. 修改坐标,使之合理
b. 改变初始猜  Guess=Huckel 或其他的,看Guess关键词。
c. 增加叠代次数 SCFCYC=N (对小分子作计算时最好不要增加,很可能结构不合理)
d. iop5/13=1)这样忽略不收敛,继续往下做。
3 分子对称性改变
a. 修改坐标,强制高对称性或放松对称性
b. 给出精确的、对称性确定的角度和二面角。 CH4的角度给到109.47122
c. 放松对称性判据  Symm=loo
  d. 不做对称性检查 iop(2/16=1) (最好加这个选项)
     iop(2/16=2) 则保持新的对称性来计算
4 Opt时收敛的问题
a. 修改坐标,使之合理
b. 增加叠代次数optcyc=N
5 优化过渡态,若势能面太平缓,则不好找到。
iop(1/8=10) 默认30(下一个结构和该结构的差别0.3Å),可改成10。如果每一步都要用到小的步长,应该加opt(notrustupdate)
6 CI(组态)方法中
QCISD(T),CCSD(T),CID方法中,省却最大循环50,若出错(L913错误)
解决方法:#P QCISD(maxcyc=N)    注:N512
7 优化过渡态
opt=TS (给出过渡态
opt=qst2 (给出反应物和产物
opt=qst3 (给出反应物和产物和过渡态)
a.  G03时的出错  opt=ts 必须加FC (force constant)
写法:opt=(TS, calcFc)
or    opt=(TS,calchffc)计算HF力常数,对QCISDCCSD等方法用;
or    opt=(TS,modRedundant)  (最好写这个)
b. 如果计算采用QCISD计算(不好计算FC
则写为QCISD opt=(TS, calcHFFC) (HF计算FC)
8. 无法写大的Scratch文件RWF
a. 劈裂RWF文件
%rwf=loc1,size1,loc2,size2,……..,locN,-1
   b. 改变计算方法 MP2=Direct可以少占硬盘空间
   c. 限制最大硬盘 maxdisk=N GB****MB,有些系统写2GB会出错,可以写2000MB
9. FOPT出错
原因是变量数与分子自由度数不相等。 可用POPT 或直接用OPT
10. 优化过渡态只能做一个STEP
原因是负本征数目不对 添加 iop1/11=1或者noeigentest
一些容易忽略的地方
不挖坑。直接说了:
1
低频率振动模式时,力常数很小,必须使用opttight以确保适当的收敛和随后任务步骤中频率计算的可靠度。
比如我曾经算过一个体系,第一个频率振动的力常数只有1020,很不放心,加入tight后,优出来的结构却与前面大相径庭,而且还出现了一个虚频。
2
IRC计算的时候,如果步长很小,必须要用verytight。就曾有人步长设为1,而不加verytight,算出来的反应路径是v字形倒海鸥状的。
3
使用int指定用于数值积分的积分网格来用于消虚频的时候,必须注意,在比较能量的时候,对所有的计算要使用相同的积分网格。
Gaussian计算出错和对应的解决方法
有时Gaussian出错,没有经验的时候,不知道该从哪儿入手。
本贴的目的,在于指明一个方向,具体的解决方法,如果过于复杂,不在此讨论,以便于
查找。比如 L502出错,这里将仅仅指出是收敛问题。具体不收敛怎么办,我想就可以查看相应的其它帖子了。
大家可以帮忙总结一下。
supi指出:
G98的出错讯息都在用户手册中的‘Gaussian 98的链接’一节中就给出了,常见错误有:
1。输入有误,这没话说了。一般都是L1L101出错
2。基组有误。特别是刚开始用自定义基组的,或者不知道基组适合范围的,这是L301出错
3。无法收敛、优化和频率计算出错,这就很多了,比如L502L716L9999。这就只好具体问题具体分析了。

本文发布于:2023-05-10 03:02:11,感谢您对本站的认可!

本文链接:https://www.wtabcd.cn/fanwen/fan/89/877008.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:计算   出错   对称性   收敛   优化   频率   时候
相关文章
留言与评论(共有 0 条评论)
   
验证码:
推荐文章
排行榜
Copyright ©2019-2022 Comsenz Inc.Powered by © 专利检索| 网站地图