作者nunuaan (nunuaan)
看板Chemistry
标题[学习] gaussian计算问题
时间Tue May 29 20:20:25 2018
近期接触这一套软体
但接触时间不长所以有些地方还不太了解
希望有熟悉这套软体的人帮忙解惑
我在做二氧化矽的晶格结构
大略是排出来了 但计算到最後总是会出现L502
按过clean优化但会让整个结构会扭曲变形
不知道设定上是不是有问题
另外我被交代算出光谱和刚性系数
光谱我大概还有点概念 但刚性系数查不太到相关资料
不知道有没有熟悉的人能够告知一下 感谢
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 218.161.83.103
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Chemistry/M.1527596428.A.6C7.html
1F:推 powernba: 你是使用什麽Method跟Basis-set ? 05/30 16:31
2F:→ nunuaan: 我是用DFT 3-21G 因为不是很熟所以其实大多地方都没动 05/30 19:05
3F:→ nunuaan: 另外我没用过水球所以不太了解怎麽用 很抱歉... 05/30 19:05
4F:推 powernba: DFT是指B3LYP吗? 05/30 21:00
5F:→ nunuaan: 是的 05/30 21:00
6F:→ powernba: 你可以先试试看半经验方法的AM1、PM3、PM6 05/30 21:01
7F:→ powernba: 搭配你之前用的3-21G进行opt计算 05/30 21:02
8F:→ powernba: L502通常是指SCF自洽场没办法收敛 05/30 21:05
9F:→ nunuaan: 好的 我试试看 感谢 05/30 21:05
10F:→ powernba: 如果你会看关键字的话 或许可以找一下SCF Done那一栏 05/30 21:07
11F:→ powernba: 看一下Convg = D是否快接近-8 05/30 21:08
12F:→ powernba: 如果是-6 -7的话 可以试着输入SCF=(MaxCycle=N) 05/30 21:09
13F:→ powernba: 我记得预设的N=64 所以你可以尝试输入100以上 05/30 21:09
14F:→ powernba: 或许就有可能成功。 05/30 21:10
15F:→ nunuaan: 这个还真的没有接触到 等等会注意的 感谢! 05/30 21:15
16F:推 powernba: 话说原PO是要算什麽的光谱? 05/30 21:17
17F:→ nunuaan: 其实只交代说去跑出二氧化矽(石英)的光谱跟刚性系数来讨 05/30 21:31
18F:→ nunuaan: 论 我也不是很确定他明确指哪个 因为非本科系所以使用上 05/30 21:32
19F:→ nunuaan: 其实算是有点障碍... 05/30 21:32
20F:推 powernba: 阿....就我的直觉刚性系数好像是用查的? 05/30 21:34
21F:→ powernba: 至於光谱的话我觉得问清楚比较保险 05/30 21:36
22F:→ nunuaan: 我当初也是这麽说 但就被说不 要用算的所以有点混乱... 05/30 21:37
23F:→ nunuaan: 好的 总之我先试着让它能够收敛 感谢 05/30 21:37
24F:→ fauna: 贴个input吧 05/31 00:28
26F:推 sunev: 原来不是真的infinite system 05/31 09:59
27F:→ kkarthur: 算这种系统不该用Gaussian, 真要用请先去了解一下 05/31 16:12
28F:→ kkarthur: periodic boundary conditions (PBC) 如何在gaussian上 05/31 16:13
29F:→ kkarthur: 使用, 相关资讯网路上应该都不难查到 05/31 16:13
30F:→ kkarthur: 拿分子团簇去模拟infinite system的晶格,起始点就错了 05/31 16:19
31F:→ kkarthur: 即使真让你收敛,算出来的东西也是乱七八糟 05/31 16:19
32F:→ kkarthur: 这样的系统建议还是去找找专门计算solid state的软体 05/31 16:25
33F:→ kkarthur: 像是VASP,CASTEP...以plane wave为基础开发的计算软体 05/31 16:27
34F:推 powernba: 我算出来真的是乱七八糟的XDDDD 05/31 17:03
35F:→ nunuaan: 了解 我试着查察查BC的使用也看看授权软体中是否有上述 06/01 00:37
36F:→ nunuaan: 软体能提供使用 感谢解惑 06/01 00:37