24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

应助: 46 (小学生)
贵宾: 8.842
金币: 7440.4
散金: 9654
红花: 118
沙发: 23
帖子: 5306
在线: 1961.4小时
虫号: 537452
注册: 2008-04-01
性别: GG
专业: 理论和计算化学
管辖: 仿真模拟

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

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

Gaussian怎么算过渡态呢？是不是一样的原理？有什么差别？

回复此楼

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

第一性原理

» 猜你喜欢

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

1楼 2012-01-10 14:06:44

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

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫

QC强帖: 17
应助: 226 (大学生)
贵宾: 5.4
金币: 2180.5
红花: 217
帖子: 1241
在线: 366.8小时
虫号: 281271
注册: 2006-09-23
专业: 宇宙学

【答案】应助回帖

感谢参与，应助指数 +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里。