| 查看: 1676 | 回复: 16 | |||
| 当前主题已经存档。 | |||
cometring木虫 (著名写手)
爱你就等于爱自己
|
[交流]
【求助】请教一个关于LST/QST搜索过渡态的问题
|
||
|
请教这个LST/QST的原理是什么？我只会用，但是不知道怎么来的。 下边我附了一张图，不知道这张图上的路径说明什么？怎么来的？ 我是菜鸟，大家多帮帮我  [ Last edited by freshgirl on 2009-6-26 at 16:56 ] |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有10人回复
-
求推荐英文EI期刊
已经有5人回复
-
存款400万可以在学校里躺平吗
已经有27人回复
-
请教限项目规定
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
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