| 查看: 2067 | 回复: 16 | |||
| 当前主题已经存档。 | |||
cometring木虫 (著名写手)
爱你就等于爱自己
|
[交流]
【求助】请教一个关于LST/QST搜索过渡态的问题
|
||
|
请教这个LST/QST的原理是什么？我只会用，但是不知道怎么来的。 下边我附了一张图，不知道这张图上的路径说明什么？怎么来的？ 我是菜鸟，大家多帮帮我  [ Last edited by freshgirl on 2009-6-26 at 16:56 ] |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有109人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
明天放榜？
已经有4人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
daniel2001
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 293.3
- 红花: 1
- 帖子: 81
- 在线: 17.7小时
- 虫号: 125677
- 注册: 2005-12-07
- 性别: GG
- 专业: 能源化工
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
cometring(金币+10,VIP+0):谢谢！！
cometring(金币+10,VIP+0):谢谢！！
|
有关过渡态的计算比较复杂，根据每个人的经验以及计算硬件的条件的不同，完成计算的时间相差很大。如果你是新手，使用两个至强CPU的服务器，40个原子左右的体系，估计得花一年左右的时间才能把所有过渡态计算完。这里我先不作过多的理论解释，只把具体的操作步骤简述如下： 1.建立并优化反应物和最终产物的分子结构，分析反应可能经过的路径，假设几个中间产物，然后建立中间产物的结构并优化。有过渡金属参与的话，选择GGA方法，RPBE泛函和DNP基组。先进行Quality为Coarse的计算，然后选Fine再优化一遍。为了加快计算速度，这一步可以在DMol的Calculation－Properties对话框里不选择任何的性质选项； 2.以你的贴图为例，假定这是从反应物到第一个中间产物的过渡态计算。我们姑且把反应物称为R0，中间产物称为P0，计算出来的过渡态称为T0，此时需要对T0进行频率分析，如果只有一个虚频，恭喜你，这就是你要找的过渡态。通常不会那么顺利，总会有几个虚频存在，这时需要具体问题具体分析，有经验者可以通过Tools－Vibrational Analysis就可以判断是否需要进行进一步的过渡态优化。初学者为了保险起见，只要发现T0有多个虚频，就应该进行TS Optimization，这一步可以保证得到只有一个虚频的过渡态。如果这一步无法完成，就表示你设想的产物结构不正确，R0与P0之间没有过渡态； 3.经过优化的过渡态还需要进行TS Confirmation，也就是从过渡态出发，反算与之相应的反应物和产物结构。计算出来的反应物和产物结构如果与你设想的反应物与产物的结构一致，表示R0与P0之间只有一个过渡态，反之则表示R0与P0之间存在多个过渡态； 4.假定通过TS Confirmation计算出的反应物结构为R01，产物结构为P01，那么以R0为反应物，R01为产物，重复步骤2、3就可以计算得到过渡态T01。以P01为反应物，P0为产物，就可以计算得到从P01到P0的过渡态T02，以此类推，可以计算得到R0到P0的所有过渡态。 5.重复以上个步骤，就可以算出从最初反应物到最终产物的所有过渡态。提取反应物和最终产物以及所有过渡态的能量，做二维的势能面图，这样才算完成了一个完整的工作。 |
9楼2008-06-06 10:13:51
10楼2008-06-06 10:26:54
cometring
木虫 (著名写手)
爱你就等于爱自己
- 应助: 0 (幼儿园)
- 贵宾: 1.03
- 金币: 4315.7
- 帖子: 2911
- 在线: 36.1小时
- 虫号: 195508
- 注册: 2006-02-23
- 性别: GG
- 专业: 催化化学
2楼2008-06-04 18:51:35
spur
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.21
- 金币: 4994.7
- 散金: 559
- 红花: 5
- 帖子: 866
- 在线: 374.9小时
- 虫号: 215683
- 注册: 2006-03-11
- 性别: GG
- 专业: 第一性原理计算
| LST/QST的方法与NEB的方法比较类似。NEB的相关信息相对要多些，可以去henkelman的主页看看。http://theory.cm.utexas.edu/henkelman/ |
3楼2008-06-04 21:31:24
cometring
木虫 (著名写手)
爱你就等于爱自己
- 应助: 0 (幼儿园)
- 贵宾: 1.03
- 金币: 4315.7
- 帖子: 2911
- 在线: 36.1小时
- 虫号: 195508
- 注册: 2006-02-23
- 性别: GG
- 专业: 催化化学
4楼2008-06-05 01:09:52
★ ★ ★ ★ ★ ★ ★ ★
spur(金币+3,VIP+0):3Q!感谢帮助！
cometring(金币+5,VIP+0):谢谢！
spur(金币+3,VIP+0):3Q!感谢帮助！
cometring(金币+5,VIP+0):谢谢！
| 我最近也在做TS SEARCH，工作中遇到了很多这样那样的问题，不过以我的经验来说楼主的过渡态搜索没有完成，就算完成也注定失败，因为你的路径太多，LST-QST是搜TS的两种方法，最终TS星点应该落在能量轨迹的最高出，单用LST搜的是最方便的，有时候会QST也会出面帮忙但是常常帮倒忙，如果QST和LST同时搜到一个吻合的TS点，那么这种搜索才是成功的，而且楼主的反应物怎么比产物能量还高啊？仅是个人理解未必正确啊！还是希望有高人回帖解答 |
5楼2008-06-05 09:14:32
cometring
木虫 (著名写手)
爱你就等于爱自己
- 应助: 0 (幼儿园)
- 贵宾: 1.03
- 金币: 4315.7
- 帖子: 2911
- 在线: 36.1小时
- 虫号: 195508
- 注册: 2006-02-23
- 性别: GG
- 专业: 催化化学
6楼2008-06-05 11:44:23
7楼2008-06-05 14:42:15
cometring
木虫 (著名写手)
爱你就等于爱自己
- 应助: 0 (幼儿园)
- 贵宾: 1.03
- 金币: 4315.7
- 帖子: 2911
- 在线: 36.1小时
- 虫号: 195508
- 注册: 2006-02-23
- 性别: GG
- 专业: 催化化学
8楼2008-06-05 16:35:56
