计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8312|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[其它量化程序] 过渡态搜索中频率部分出错问题

[复制链接 Copy URL]

9

帖子

0

威望

89

eV
积分
98

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 CJM1111 于 2017-9-20 15:58 编辑

我在用Dmol3搜索过渡态时,老是出错。outmol 文件部分已显示过渡态结果已找到,但是在频率计算部分出错,显示没有收敛。如果用过渡态搜索的xsd文件重新提交频率计算任务,还是会出错。精度已经是medium了,那这该怎么解决啊。求各位大神指导!!!附件是outmol,  input 文件。

TS7-5.outmol

1.08 MB, 下载次数 Times of downloads: 7

TS7-5.err

55 Bytes, 下载次数 Times of downloads: 0

TS7-5.input

2.46 KB, 下载次数 Times of downloads: 9

3809

帖子

3

威望

1万

eV
积分
20344

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2017-9-21 09:19:24 | 只看该作者 Only view this author
磁矩初猜的时候只需要猜Fe就行,但几个Fe的磁序到底如何你最好确定一下,或者都试试,非周期性的情况下并不推荐使用smearing。结构本身比较奇怪,我对你的研究体系不很了解,所以结构的可靠性还需要你自己把握。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

9

帖子

0

威望

89

eV
积分
98

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2017-9-21 10:16:29 | 只看该作者 Only view this author
卡开发发 发表于 2017-9-21 09:19
磁矩初猜的时候只需要猜Fe就行,但几个Fe的磁序到底如何你最好确定一下,或者都试试,非周期性的情况下并不 ...

老师,请问您第一句话什么意思啊?我设置的时候没有涉及到磁矩啊。还有我想问一下,如果过渡态任务正常结束,但是结果却没有频率计算结果,该怎么办呐?谢谢您,求指导啊

TS-H.outmol

88.36 KB, 下载次数 Times of downloads: 1

3809

帖子

3

威望

1万

eV
积分
20344

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2017-9-21 10:37:41 | 只看该作者 Only view this author
CJM1111 发表于 2017-9-21 10:16
老师,请问您第一句话什么意思啊?我设置的时候没有涉及到磁矩啊。还有我想问一下,如果过渡态任务正常结 ...
  1. ……
  2. Start_Spin_Populations on <--
  3. 1 5.0000 <-s
  4. 2 5.0000 <-s
  5. 3 5.0000 <-s
  6. 4 5.0000 <-s
  7. 5 5.0000 <-s
  8. ……
复制代码

这里你可以调整一下。另外我觉得你给的结构挺古怪的,我不知道你的H是本来就有还是后来手工添加用于饱和边缘的。过渡态如果真的能正确找到的话,原则上可以在找到的过渡态上重新做频率计算,没有太多影响。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

9

帖子

0

威望

89

eV
积分
98

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2017-9-21 10:57:25 | 只看该作者 Only view this author
卡开发发 发表于 2017-9-21 10:37
这里你可以调整一下。另外我觉得你给的结构挺古怪的,我不知道你的H是本来就有还是后来手工添加用于饱 ...

结构上的H是后来加上去用于饱和边缘的,看文献上做团簇加了。对于Spin,结构里只有一个Fe, 我只是在modify--Electronic Configuration 里面设定了Fe的电荷。谢谢老师

9

帖子

0

威望

89

eV
积分
98

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2017-9-21 11:08:41 | 只看该作者 Only view this author
卡开发发 发表于 2017-9-21 10:37
这里你可以调整一下。另外我觉得你给的结构挺古怪的,我不知道你的H是本来就有还是后来手工添加用于饱 ...

搜索出来的过渡态能量怎么比产区还低,但是看结构和搜索路径应该是合理的呀,求老师指导!

1505962998(1).jpg (7.69 KB, 下载次数 Times of downloads: 53)

反应物,产物,过渡态能量

反应物,产物,过渡态能量

1505963206(1).jpg (19.49 KB, 下载次数 Times of downloads: 54)

搜索路径

搜索路径

3809

帖子

3

威望

1万

eV
积分
20344

Level 6 (一方通行)

围观吃瓜群众

7#
发表于 Post on 2017-9-21 11:26:38 | 只看该作者 Only view this author
CJM1111 发表于 2017-9-21 10:57
结构上的H是后来加上去用于饱和边缘的,看文献上做团簇加了。对于Spin,结构里只有一个Fe, 我只是在modif ...

你可能没设置好初始磁矩,把其它原子全加上+5的磁矩了。我不明白为啥会有个H挤在固体结构当中((14.90625 26.83078 15.81995)这个H),非边缘的、暴露在表面的O原子不用饱和。

你的过渡态比产物能量低有可能是由于固定不合理导致的。当然如上所说模型本身也有诸多不合理之处,这样的体系确实不好处理,宁可稍微降低点精度按照PBC做。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

9

帖子

0

威望

89

eV
积分
98

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2017-9-21 14:30:01 | 只看该作者 Only view this author
卡开发发 发表于 2017-9-21 11:26
你可能没设置好初始磁矩,把其它原子全加上+5的磁矩了。我不明白为啥会有个H挤在固体结构当中((14.90625 ...

老师,我是从叶腊石抠出来一个团簇,然后把本来连接其他原子的位置,添加了H。然后替换了一些原子,模拟了Fe-蒙脱石。您能帮我看看那里有不合理的地方吗,谢谢您!

pyrophyllite.xsd

210.85 KB, 下载次数 Times of downloads: 6

叶腊石

pyrophyllite-5.xsd

340.87 KB, 下载次数 Times of downloads: 3

团簇

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 23:19 , Processed in 0.290539 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list