应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求QST2方法输入文件
61/1返回列表
查看: 1052  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者默默地别动将赠送您 10 个金币

匿名

本帖仅楼主可见

» 猜你喜欢
1楼2020-05-07 13:33:16
已阅   同方向广播   申请QC强帖   回复此楼   编辑   查看我的主页

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhou2009: 金币+5 2020-05-08 15:23:23
首先GAUSSIAN的QST方法其实算法很差的。稍微复杂一点儿的过渡态根本就根据反应物和产物猜不出来过渡态初始结构。导致各种会报错。
建议还是使用 OPT=TS方法直接优化过渡态。而过渡态初始结构 根据反应物结构 用柔性扫描等方法去确定成功率要高很多。
不过既然你要例子我就给你发一个算例。但是我个人很不推荐用QST 方法。或者说至少不是首选方法。

算例 La系金属有机物。 对La用SDD 而对其它原子使用6-311+g(d,p), 输入文件中的反应物和产物结构是已经优化过的. 如果你用的g09 就把关键词里的 recalc=3 去掉就可以。

%nprocshared=20
%mem=16GB
%chk=QST2_example.chk
# opt(qst2,recalc=3) ub3lyp/genecp

reactant

0 2
La                 0.28580300   -0.00811400   -0.00021600
N                 -2.26630400    0.06638300    0.00171300
H                 -2.63322500   -0.40365100   -0.83140800
H                 -2.63760600    1.02030300    0.01176100
H                 -2.63194400   -0.42064000    0.82558500

product

0 2
La                 0.23245900    0.00480200   -0.00001500
N                 -1.98036400   -0.17390000   -0.00087500
H                 -2.58317900   -0.14448700   -0.82076700
H                  0.64272200    2.06348700    0.01531700
H                 -2.58298400   -0.16460400    0.81966700

La 0
SDD
****
N H 0
6-311+G(d,p)
****

La 0
SDD
赞一下 回复此楼
2楼2020-05-07 16:04:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

本帖仅楼主可见
赞一下
3楼2020-05-08 12:23:14
 已阅   申请QC强帖   回复此楼   编辑   查看我的主页

paramecium86

版主 (著名写手)
引用回帖:
3楼: Originally posted by 默默地别动 at 2020-05-07 23:23:14
老师，我就按照这个输入，但是出现了错误。关键词里需要加freq吗？
Old curvilinear step not converged, using linear step:
SCX= 3.94D+01 DXMaxT= 6.18-317 SCLim= 0.00D+00 Fact= 0.00D+00
RedCar/ORedCr ...

这就是我上一个帖子提到的。QST2的算法很傻。不是你给什么样的初始结构 和 末态结构都能猜得出来的。这个提示就是高斯没猜出来TS结构所以出错。所以一般你自己都猜不出来大概过渡态结构的体系，QST算法也大概率猜不出来。
赞一下 回复此楼
4楼2020-05-08 13:40:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
3楼: Originally posted by 默默地别动 at 2020-05-07 23:23:14
老师，我就按照这个输入，但是出现了错误。关键词里需要加freq吗？
Old curvilinear step not converged, using linear step:
SCX= 3.94D+01 DXMaxT= 6.18-317 SCLim= 0.00D+00 Fact= 0.00D+00
RedCar/ORedCr ...

计算方法就是这样的，也不需要改什么。根据你给的结构 QST能猜出来，就能开始算。猜不出来非常快就会报错。所以我一直都建议用OPT=TS方法的
赞一下 回复此楼
5楼2020-05-08 13:43:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

本帖仅楼主可见
赞一下
6楼2020-05-08 15:57:19
 已阅   申请QC强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 默默地别动 的主题更新
61/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 求助大神，Fe和Al离子对MOF都有淬灭，当两种离子共存时，怎么区分两种离子？ 10+5 maoxiao 2024-06-02 6/300 2024-06-03 19:35 by maolC
[基金申请] 年年想打听，年年打听不到。。 +12 hdzw9071 2024-06-03 16/800 2024-06-03 19:00 by hdzw9071
[硕博家园] 讨厌酒桌文化 +10 hahamyid 2024-05-31 10/500 2024-06-03 17:50 by beefly
[基金申请] 前些天开会有个人见到人就搞关系，一查此人全是MDPI/Hindawi论文，鄙视! +19 zju2000 2024-06-02 20/1000 2024-06-03 16:25 by lkhjiayou
[硕博家园] 每到中夜，情难自抑 +37 sioc-sunj 2024-05-28 61/3050 2024-06-03 15:16 by EllieWU
[硕博家园] 实验室太吵闹，无法安静学习，怎么办？ +11 utahh 2024-05-31 22/1100 2024-06-03 11:18 by utahh
[基金申请] 化学B02口青基 代表作都是什么水平的？向大佬求助 +14 arthas_007 2024-06-01 19/950 2024-06-03 11:14 by prisonrabbit
[基金申请] 离职后国自然项目不能变更单位，在新单位还能申请新的国自然项目吗 +3 宜兰 2024-05-31 11/550 2024-06-03 10:26 by lilyxjw
[硕博家园] 材料博士毕业 +3 112233ssg 2024-06-02 3/150 2024-06-03 09:59 by Lwyy1106
[教师之家] 职能部门工作人员态度不好是普遍的吗？怎么让他们态度好一些？ +7 河西夜郎 2024-06-01 7/350 2024-06-03 08:37 by YY392619
[基金申请] 2024杰青和万人领军什么时候会评 +6 墨香琴韵 2024-06-02 6/300 2024-06-03 07:24 by llhljsy
[硕博家园] 论大家对6070后普通教授导师的看法 +8 SNaiL1995 2024-05-28 14/700 2024-06-03 05:07 by SNaiL1995
[微米和纳米] 银浆叠压烧结粘片的问题有没有好的解决方法 +3 微纳材料 2024-05-30 5/250 2024-06-02 21:17 by manmansboy
[教师之家] 选择 +8 459582015 2024-05-31 11/550 2024-06-02 18:30 by auvauv
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[考博] 申博求助?本硕双非一篇三区和四区去985工科非天坑专业有没有希望？ +4 LYS1200 2024-05-29 6/300 2024-06-01 08:50 by 学术渣渣神
[论文投稿] 纠结选哪一个期刊，电化学领域 50+8 Freya163 2024-05-28 10/500 2024-05-31 15:09 by wzykobe
[高分子] MMA预聚体光固化发雾问题求助 +3 惠亚金总 2024-05-29 10/500 2024-05-31 14:59 by 惠亚金总
[硕博家园] 求问endnote上Wiley引用格式在哪下载 +3 Elon. 2024-05-30 3/150 2024-05-31 10:13 by crazy peng
[论文投稿] 高手朋友推荐比较容易投稿和录用的SCI期刊，不投稿收费SC，对分区没有要求 5+3 xintangren 2024-05-28 4/200 2024-05-29 10:46 by xintangren
信息提示
关闭
请填处理意见
关闭