24小时热门版块排行榜    

查看: 3984  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

[求助] Guassian和VASP计算过渡态的原理和算法上差别在哪里?

VASP是在周期性模型中,用string, elastic band等模型来描述反应路径的。实际上是分析势能面上的能量走势。

Gaussian怎么算过渡态呢?是不是一样的原理?有什么差别?
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

感谢参与,应助指数 +1
cenwanglai(金币+5): ★★★★★最佳答案 非常感谢sobereva MM 2012-01-11 09:24:05
VASP和Gaussian所用到的搜索过渡态、确定IRC的方法在这里详细讨论过
http://hi.baidu.com/sobereva/blo ... Stamp=1326195337907
http://hi.baidu.com/sobereva/blo ... ee62da572c841f.html

string、NEB这类方法,获得的是势能面上两个极小点之间的最小能量路径,以化学的角度来说,就是一次性将IRC和过渡态同时确定了。原理可以简单描述为将一串珠子的两端分别固定在势能面上反应物和产物的位置上,然后让这串珠子自行弛豫,最后这串珠子就描述出了最小能量路径(或称IRC),而这个路径也经过了过渡态位置。
实际上在Gaussian中path=N的方法也可以达到同样的目的,但是原理还是很不同的。

在Gaussian中用的方法都是人为指定(opt=ts)或者以自动方式(QST2,QST3)给出初猜位置,然后通过不同的迭代方法(berny、GEDIIS等)从这个初猜位置逐渐收敛到过渡态位置。若再想获得IRC则从过渡态位置开始顺着虚频矢量方向以不同方法做迭代。

无论是VASP还是Gaussian用的方法,研究的都是势能面结构,只不过算法的思路有本质不同。原理上讲,VASP的那些方法可以用在Gaussian上,Gaussian的方法也能用在VASP里。
6楼2012-01-10 19:50:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
zhou2009(金币+3): 2012-01-10 16:48:35
Gaussian是在零维模型中,用TS(实际上是berny算法,RFO/P-RFO法或牛顿法),QSTn(n=2,3,QST2是二次同步转变方法;而QST3是由用户自己为QST2提供TS的初猜)等模型来描述反应路径的。而不管是过渡态的计算还是IRC的计算,实际上仍旧是分析势能面上的能量走势。
2楼2012-01-10 16:33:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★
zhou2009(金币+2): 2012-01-10 16:49:04
其实从根子来说,VASP和Gaussian没区别,都是薛定谔的直系后裔,咱们才是一家人。
对咱们来说,做力场的那些人才是异端,统统绑上火刑柱,点火烧了
3楼2012-01-10 16:37:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

【答案】应助回帖

cenwanglai(金币+5): ★★★很有帮助 谢谢 2012-01-10 16:50:59
干脆把量化版和第一版合并吧,称为“正教版”
把分子模拟板改名为“邪教版”或者“异端版”
4楼2012-01-10 16:46:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:34 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:23 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 17:01 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 4/200 2026-11-07 16:28 by Xc2FxHsHkBNc
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:17 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-07 15:44 by Xc2FxHsHkBNc
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见