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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 其他 » 关于MS DMOL3模块算过渡态的问题，急！！！！
81/1返回列表
查看: 2393  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhuhy1116

新虫 (初入文坛)

[求助] 关于MS DMOL3模块算过渡态的问题，急！！！！

本人最近正在研究金属催化剂对煤焦还原NO的影响，但是采用的是dmol模块，我构建了金属切面，选择了焦炭模型和NO分子同时吸附在切面上，但是算出来的过渡态虚频很多，而且振动强度都为零，我想问一下是我的方法有问题么？用castep来算的话，怎么保证你算得过渡态的正确性，因为不像dmol模块可以通过频率分析来确认过渡态是否准确啊。
Frequency (1/cm)          Intensity (km/mol)
-699.22000000        0.00000000
-497.61000000        0.00000000
-44.83000000        0.00000000
-40.25000000        0.00000000
-38.84000000        0.00000000
-34.48000000        0.00000000
-32.98000000        0.00000000
-27.25000000        0.00000000
-24.72000000        0.00000000
-22.35000000        0.00000000
-21.64000000        0.00000000
-18.19000000        0.00000000
-14.31000000        0.00000000
-7.72000000        0.00000000
-1.51000000        0.00000000
11.69000000        0.00000000
16.69000000        0.00000000
18.17000000        0.00000000
22.62000000        0.00000000
27.23000000        0.00000000
28.68000000        0.00000000
29.27000000        0.00000000
32.47000000        0.00000000
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-12-15 15:35:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★ ★
gmy1990: 金币+2, 感谢应助 2014-12-19 18:54:38
DMol3你想看到频率的详细信息必须得进行性质中频率计算，这里给出的频率我估计应该是几何优化过程采用的粗糙的Hessian估计的。MS下CASTEP做反应确实稍弱一些。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
2楼2014-12-17 14:16:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuhy1116

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 卡开发发 at 2014-12-17 14:16:13
DMol3你想看到频率的详细信息必须得进行性质中频率计算，这里给出的频率我估计应该是几何优化过程采用的粗糙的Hessian估计的。MS下CASTEP做反应确实稍弱一些。

我现在是一直卡在这里进行不下去，进行过渡态优化也是收敛不了，怎么办呐···
赞一下 回复此楼
3楼2014-12-17 15:01:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
3楼: Originally posted by zhuhy1116 at 2014-12-17 15:01:53
我现在是一直卡在这里进行不下去，进行过渡态优化也是收敛不了，怎么办呐···...

我不知道你指的不收敛是？到底是频率算不出，还是过渡态搜索不出？或是过渡态优化结构的过程不能收敛？
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
4楼2014-12-19 09:50:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuhy1116

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by 卡开发发 at 2014-12-19 09:50:03
我不知道你指的不收敛是？到底是频率算不出，还是过渡态搜索不出？或是过渡态优化结构的过程不能收敛？...

频率不是算出来了呢，但是振动强度都是零，不收敛是在过渡态优化的过程中，而且这个算出来的虚频也太多了···
赞一下 回复此楼
5楼2014-12-22 15:34:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
5楼: Originally posted by zhuhy1116 at 2014-12-22 15:34:56
频率不是算出来了呢，但是振动强度都是零，不收敛是在过渡态优化的过程中，而且这个算出来的虚频也太多了···...

不是算的频率是看不到振动模的（振动强度），而且Hessian矩阵不够精确，没啥参考意义。补充频率计算吧。

如果过渡态不收敛，只能通过修改那些DIIS和charge mix等来增加收敛的可能。金属切面的话要保证k和smearing都非常充分、合适。一些带磁性的金属可能处理起来会比较棘手，如果是这样就没办法了。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
6楼2014-12-22 16:20:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuhy1116

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by 卡开发发 at 2014-12-22 16:20:00
不是算的频率是看不到振动模的（振动强度），而且Hessian矩阵不够精确，没啥参考意义。补充频率计算吧。

如果过渡态不收敛，只能通过修改那些DIIS和charge mix等来增加收敛的可能。金属切面的话要保证k和smeari ...

确实 做fe切面的话连结构优化都成功不了，更别提过渡态了，我现在用的是CAO切面。您的意思是我这么计算出来的虚频和振动强度没有意义是么？那应该怎么补充频率计算？过渡态搜索的时候outmol里面已经有反应能垒的信息了，根据我这种情况算出来的能量可信么？另外调节参数的时候，一直在用smearing来加速收敛，你说的这两个倒是很少使用，您能详细说说这两个参数一般怎么调才能使过渡态收敛呢？
赞一下 回复此楼
7楼2014-12-23 09:46:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
7楼: Originally posted by zhuhy1116 at 2014-12-23 09:46:55
确实 做fe切面的话连结构优化都成功不了，更别提过渡态了，我现在用的是CAO切面。您的意思是我这么计算出来的虚频和振动强度没有意义是么？那应该怎么补充频率计算？过渡态搜索的时候outmol里面已经有反应能垒的信 ...

不计算频率根本是得不到振动强度的，补充频率需要在Properties计算Freq。如果能够搜索到过渡态，outmol中的能垒一般稍高估些，可通过频率计算后再进行TS opt。

Smearing很大的结果是不可靠的，至少要保证输出文件中-KTeS项<1meV/atom。

如果很难收敛的情况一般建议先提升精度，如果Fine精度还很困难，减小charge mix和DIIS。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
8楼2014-12-23 15:10:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhuhy1116 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 272分材料子求调剂 +28 Loy0361 2026-04-10 36/1800 2026-04-11 00:45 by success5733
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+3 chenny174 2026-04-10 17/850 2026-04-11 00:26 by 骑牛渡寒江
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 314求调剂 +23 wakeluofu 2026-04-09 24/1200 2026-04-10 15:31 by MOF_Catal
[考研] 282，电气工程专业，求调剂，不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[基金申请] 有爆料，一个青年教师卖房得400万，然后换了一个四青帽子 +9 babu2015 2026-04-08 9/450 2026-04-10 11:43 by 苏东坡二世
[考研] 求调剂 +9 翩翩一书生 2026-04-09 9/450 2026-04-10 11:43 by wp06
[考研] 材料工程302分求调剂 +18 zyx上岸！ 2026-04-04 18/900 2026-04-10 10:07 by 314126402
[考研] 284求调剂 +7 archer.. 2026-04-09 8/400 2026-04-09 21:30 by bljnqdcc
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH（等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 280求调剂 +5 兮兮夜夜 2026-04-09 8/400 2026-04-09 11:15 by 兮兮夜夜
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +3 慕绝cc 2026-04-09 3/150 2026-04-09 09:57 by liuhuiying09
[考研] 296求调剂 +3 汪！？！ 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 生物学学硕，初试351分，求调剂 +4 …～、王…～ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 材料科学与工程320求调剂，080500 +12 黄瓜味薯片 2026-04-06 12/600 2026-04-08 16:26 by luoyongfeng
[考研] 一志愿南昌大学，085600，344分求调剂 +11 调剂上岸玘 2026-04-05 12/600 2026-04-08 16:17 by luoyongfeng
[考研] 363求调剂 +9 zh096 2026-04-04 9/450 2026-04-07 21:51 by 418490947
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
信息提示
关闭
请填处理意见
关闭