24小时热门版块排行榜    

查看: 1406  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[基金申请] fileCode有新解读? +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见