24小时热门版块排行榜    

查看: 1404  |  回复: 5

daixing2011

铜虫 (初入文坛)

[求助] Dmol3优化周期结构,如何实现离域内坐标?

看MS说明书,以及网上的讨论,Dmol3用离域内坐标优化(delocalized internal coordinate),效率超级快。

但手册上的例子,只是在MS面板上设置任务等等,并没有体现出哪里设置了离域内坐标。

测试了几个手册中的关键词,如Opt_Del_Test, Opt_Bend_Lin,也都无效,程序提示输入错误。


请问,Dmol在周期性优化中怎样实现离域内坐标,还是默认就已经是了?不像吧,输出文件中都是笛卡尔坐标。

另外,做表面吸附的研究,物理吸附,适合用离域内坐标么。

Thanks in advance
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

默认是采用离域内坐标,但有些体系并不是这样,个人测试有些体系直角坐标反而比较快,甚至需要追加GDIIS。

物理吸附最好考虑一下DFT-D(建议Grimme06)
不一定挂在论坛,计算问题问题欢迎留言。
2楼2015-08-14 03:38:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daixing2011

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by 卡开发发 at 2015-08-14 03:38:08
默认是采用离域内坐标,但有些体系并不是这样,个人测试有些体系直角坐标反而比较快,甚至需要追加GDIIS。

物理吸附最好考虑一下DFT-D(建议Grimme06)

感谢回复!

仔细读了输出文件,确实默认就是离域内坐标。

同时发现一个问题,dmol好像不能固定晶格Z轴,只优化XY。如果您知道请指点。

已经考虑了dft-d   grimme色散,感谢提醒。

请您点应助,要不给不了金币,别嫌少,心意一下,哈哈。
3楼2015-08-15 01:41:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by daixing2011 at 2015-08-15 01:41:46
感谢回复!

仔细读了输出文件,确实默认就是离域内坐标。

同时发现一个问题,dmol好像不能固定晶格Z轴,只优化XY。如果您知道请指点。

已经考虑了dft-d   grimme色散,感谢提醒。

请您点应助,要不给不 ...

DMol3本身优化晶格也不准,建议不要那样做,可以的话可以对XY方向晶格常数“扫描”一下。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2015-08-15 01:46:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daixing2011

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by 卡开发发 at 2015-08-15 01:46:08
DMol3本身优化晶格也不准,建议不要那样做,可以的话可以对XY方向晶格常数“扫描”一下。...

感谢您再次回复。

您指的“扫描”,是否是逐个做单点能,比如X=Y=10,X=Y=10.5,.....,这样,然后要一个能量最低的?

dmol3优化一个周期表面,我这里一直都出现这样的情况:
opt==  Cycle    Total Energy   Energy change   Max Gradient   Max Displacement
opt==      tolerance:.......      0.0000100        0.001058       0.009449
opt==    1    -32105.2859592      0.0000000        0.000066       0.000089


                      Step    1                                       <-- DELOC
                                                                      <-- DELOC
+-----------+-----------------+-----------------+------------+-----+ <-- DELOC
| Parameter |      value      |    tolerance    |    units   | OK? | <-- DELOC
+-----------+-----------------+-----------------+------------+-----+ <-- DELOC
|  delta E  |   -0.204437E+02 |    0.100000E-04 |      Ha    | No  | <-- DELOC
|  |F|max   |    0.660822E-04 |    0.105835E-02 |      au    | Yes | <-- DELOC
|  |dR|max  |    0.889212E-04 |    0.944863E-02 |      au    | Yes | <-- DELOC
+-----------+-----------------+-----------------+------------+-----+ <-- DELOC

  Geometry optimization completed successfully in    1 steps.


然后就正常结束了。
很奇怪,delta E 并没有收敛,这是不是几何结构收敛了但晶格没有收敛?为什么还能正常结束呢。

请指点,多谢。
5楼2015-08-15 20:31:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
5楼: Originally posted by daixing2011 at 2015-08-15 20:31:43
感谢您再次回复。

您指的“扫描”,是否是逐个做单点能,比如X=Y=10,X=Y=10.5,.....,这样,然后要一个能量最低的?

dmol3优化一个周期表面,我这里一直都出现这样的情况:
opt==  Cycle    Total Energy ...

对,而且每一步都要进行非opt cell的geom opt,如果两个边长不相等或是夹角不确定,这个做起来就相当复杂了。

第一步的delta其实没啥意义(第0步的能量从何而来?),估计一般判定第一步如果有两个任意指标收敛就算是收敛了。
不一定挂在论坛,计算问题问题欢迎留言。
6楼2015-08-15 21:14:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 daixing2011 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见