24小时热门版块排行榜    

查看: 1909  |  回复: 30

lingtao

铜虫 (正式写手)

[求助] Dmol3 中经行周期性结构TS过程, Vibrational Analysis中出现可虚频但是没有强度值 已有1人参与

如题:小弟利用Dmol3 模块进行周期性分子筛结构下的催化反应TS search 过程,经过40几个小时的计算,算出来了过渡态,之后再 Vibrational Analysis中出现了2个虚频(-264和-25)其他的都不是虚频,但是没有对应的强度值(Intensity),自己琢磨的很久,不知道哪里出现了问题,想请教坛里的大神这方面的经验。下面是我的两张图:

该数据应楼主要求删除。

[ Last edited by franch on 2014-8-12 at 15:46 ]
回复此楼

» 猜你喜欢

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

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★ ★
lingtao(franch代发): 金币+5, 谢谢回帖交流。。。 2014-08-11 17:48:54
引用回帖:
5楼: Originally posted by lingtao at 2014-08-10 17:27:14
我昨天就开始用TS optimization那个选项,对搜到的过渡态经行了优化,但是那三个指标(Energy change, Max Displacement,Force )一直没有收敛,而且看了那个在计算过程中一直变动的.xsd文件,感觉优化的这个过渡态 ...

要不然考虑加上OPT_GDIIS的关键字,并且使用直角坐标。
不一定挂在论坛,计算问题问题欢迎留言。
6楼2014-08-10 18:12:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huxiao.gl

木虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★
franch: 金币+8, 鼓励交流。。。鼓励DIY。。。呵呵。。 2014-08-11 17:44:59
引用回帖:
15楼: Originally posted by 卡开发发 at 2014-08-10 20:43:23
不,所有平面波方法必须是有周期边界的,这样的话结果不会有什么不同,这是相同理论方法之间的比较,统计的方式也是根据Debye的理论。

Dmol3这个程序用的原子轨道方法,不必要求有边界条件。这样就可以做比较, ...

我刚才用Dmol3算了H2O分子的频率,的确没有区别,你要不信自己再试试。
实验值:1594.7, 3651.1,3755.9
Dmol3非周期模型: 1614.7, 3730.1 , 3845.4
DMol3周期模型,边长1nm: 1615.4,  3729.9 ,3844.4
对应的ZPE分别是0.56eV,0.57eV,0.57eV
good luck
17楼2014-08-10 21:01:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huxiao.gl

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
franch: 金币+5, 鼓励交流。。。 2014-08-11 17:47:14
引用回帖:
22楼: Originally posted by 卡开发发 at 2014-08-10 22:33:39
我没遇到到过表面处理热化学的,最多通过这样修正ZPE。

可以想象个分子,一段吸附在金属表面,整体虽然不能像原来那样转动,但是上侧的支链上的烷基转动其实丝毫不影响。

问题在于,要这样识别转动的片段我觉 ...

表面的热力学可以参考以下文献:
David Loffreda,Surface Science 600 (2006) 2103–2112
其中,引言里有这样一段话:
The vibrational energy and entropy contributions of the adsorbed phase should also be significant since the reactant loses translational and rotational energy during the chemisorption process.

烷基的确可以转动,但是存在能垒。我认为大部分情况下,这样的能垒是很大的。因此我们不需要考虑吸附态的转动能。
如果你坚持认为吸附态有转动能,请给出参考文献或其它证据。

此外,计算转动能用一个很简单的公式就能搞定。
good luck
24楼2014-08-10 23:02:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★
franch: 金币+4, 谢谢回帖交流。。。 2014-08-11 17:48:22
引用回帖:
29楼: Originally posted by huxiao.gl at 2014-08-10 23:25:59
直接忽略,因为多出的频率就不存在。...

对啊,所以表面体系直接计算频率也就很难确定多出的频率了,因为没有开放体系做对照,这样一来,如果有2~3个虚频,其实很难判断。TS Opt做出来似乎就不可靠了。
不一定挂在论坛,计算问题问题欢迎留言。
30楼2014-08-10 23:28:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-08-11 17:38:54
要对TS的构型计算了频率才有频率的强度。你可以在得到的构型上计算能量及性质中的频率。

另外表面体系dmol3计算频率尤其小心,仅可作为参考,零点能及热化学性质不具备明确,一般也不会按照这样的方式来处理。做个过渡态确定还是很有必要的。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2014-08-10 11:48:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingtao

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by 卡开发发 at 2014-08-10 11:48:11
要对TS的构型计算了频率才有频率的强度。你可以在得到的构型上计算能量及性质中的频率。

另外表面体系dmol3计算频率尤其小心,仅可作为参考,零点能及热化学性质不具备明确,一般也不会按照这样的方式来处理。做 ...

你好,先谢谢你了,我在搜索过渡态的时候,已经勾选上了Frequence这一项了,我的构型是附件这样的。想请问你一下,搜索过渡态的时候,频率的强度值是不是不重要?有虚频的存在就行了吧?因为我找到虚频了
Dmol3 中经行周期性结构TS过程, Vibrational Analysis中出现可虚频但是没有强度值
计算参数3.png


Dmol3 中经行周期性结构TS过程, Vibrational Analysis中出现可虚频但是没有强度值-1
过渡态文件.png

3楼2014-08-10 12:03:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★ ★ ★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-08-11 17:40:52
lingtao(franch代发): 金币+5, 谢谢回帖交流。。。 2014-08-11 17:50:30
引用回帖:
3楼: Originally posted by lingtao at 2014-08-10 12:03:51
你好,先谢谢你了,我在搜索过渡态的时候,已经勾选上了Frequence这一项了,我的构型是附件这样的。想请问你一下,搜索过渡态的时候,频率的强度值是不是不重要?有虚频的存在就行了吧?因为我找到虚频了

计算参数 ...

强度没那么重要,但正确的过渡态应该只有一个虚频。你这有两个。这种情况只能试试做个过渡态的优化。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2014-08-10 14:53:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingtao

铜虫 (正式写手)

引用回帖:
4楼: Originally posted by 卡开发发 at 2014-08-10 14:53:05
强度没那么重要,但正确的过渡态应该只有一个虚频。你这有两个。这种情况只能试试做个过渡态的优化。...

我昨天就开始用TS optimization那个选项,对搜到的过渡态经行了优化,但是那三个指标(Energy change, Max Displacement,Force )一直没有收敛,而且看了那个在计算过程中一直变动的.xsd文件,感觉优化的这个过渡态文件,不是我想要的那个结构了。是不是一定要进行TS optimization这个选项?因为第二个虚频比较小(-25),能不能忽略呢?
5楼2014-08-10 17:27:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huxiao.gl

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 卡开发发 at 2014-08-10 11:48:11
要对TS的构型计算了频率才有频率的强度。你可以在得到的构型上计算能量及性质中的频率。

另外表面体系dmol3计算频率尤其小心,仅可作为参考,零点能及热化学性质不具备明确,一般也不会按照这样的方式来处理。做 ...

用DMol3算表面的热力学性质,算法有问题吗?具体什么问题呢?
一般是用什么方法计算?
good luck
7楼2014-08-10 19:08:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-08-11 17:41:24
引用回帖:
7楼: Originally posted by huxiao.gl at 2014-08-10 19:08:17
用DMol3算表面的热力学性质,算法有问题吗?具体什么问题呢?
一般是用什么方法计算?...

你可以比较一下分子放在很大的格子(PBC & Gamma Points)和自由态(非PBC)的区别。

这有可能和统计方式有关,看了一下计算的振动模式似乎有转动混入。

否则根本不必做复杂的热化学循环http://muchong.com/bbs/viewthread.php?tid=3695520
不一定挂在论坛,计算问题问题欢迎留言。
8楼2014-08-10 19:14:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huxiao.gl

木虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
franch: 金币+3, 谢谢回帖交流。。。 2014-08-11 17:41:45
引用回帖:
8楼: Originally posted by 卡开发发 at 2014-08-10 19:14:08
你可以比较一下分子放在很大的格子(PBC & Gamma Points)和自由态(非PBC)的区别。

这有可能和统计方式有关,看了一下计算的振动模式似乎有转动混入。

否则根本不必做复杂的热化学循环http://muchong.com/b ...

Dmol3是把非周期性体系当做气态分子处理,吉布斯能包含了Translational, rotational, vibrational 及ZPE的贡献。
而对于周期性体系(即使你把气体分子放入周期性盒子),Dmol3是当成固态来处理,所以吉布斯能只包括了vibrational 及ZPE的贡献。
对于同一分子,两种模型算出来的吉布斯能当然不一样。

那个帖子的论文,没有这种方法算。而是设置DFT在0K下的总能为参考态,H和S的变化取实验值。
good luck
9楼2014-08-10 19:24:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
franch: 金币+2, 鼓励交流。。。。。。 2014-08-11 17:42:23
引用回帖:
9楼: Originally posted by huxiao.gl at 2014-08-10 19:24:57
Dmol3是把非周期性体系当做气态分子处理,吉布斯能包含了Translational, rotational, vibrational 及ZPE的贡献。
而对于周期性体系(即使你把气体分子放入周期性盒子),Dmol3是当成固态来处理,所以吉布斯能只包 ...

对于气态分子平均自由程相对大,振转平均存在,如果按照估计的统计方式显然物理意义缺失,并且振动能也不一样。。。同样表面体系分子个别基团仍然不丧失转动,这样物理意义也缺失了吧?

既然能够算的H和S,为何还要取实验值?
不一定挂在论坛,计算问题问题欢迎留言。
10楼2014-08-10 19:33:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Zeochem 的主题更新
信息提示
请填处理意见