计算化学公社

标题: dmol3如何添加电场？添加电场后所对应的自由能的值是取哪一个？非常感谢 [打印本页]

作者
Author: lzc 时间: 2017-9-4 14:27
标题: dmol3如何添加电场？添加电场后所对应的自由能的值是取哪一个？非常感谢
各位老师，最近在做电催化方面的计算，请问dmol3做催化方面计算的时候如何添加电场？所添加的电池所对应的自由能值是哪一个？非常感谢

作者
Author: zyj19831206 时间: 2017-9-4 19:28
在input文件里面手动加入电场的关键词和方向。

作者
Author: lzc 时间: 2017-9-5 16:03

zyj19831206 发表于 2017-9-4 19:28
在input文件里面手动加入电场的关键词和方向。

非常感谢，如果加电场，最好是三个方向全都加是吗？

作者
Author: 卡开发发 时间: 2017-9-5 17:02

lzc 发表于 2017-9-5 16:03
非常感谢，如果加电场，最好是三个方向全都加是吗？

手册没有具体指明加电场的方向到底是按照怎样的规则，我测试下来感觉电场的三个方向是根据直角坐标来考虑的，建议你测试一下。一般催化外加的电场是垂直于表面的。DMol3加电场做有个问题就是没办法同时修正体系的偶极，另外DMol3没有合适的溶剂化模型来描述过渡金属，那套COSMO模型使用当中你可以在输出看到对很多元素没有定义，算出来自然也就没有合理性了。

加电位更流行的做法不是通过电场来做，而是给体系外加电子，再根据功函数响应来计算实际到底加了多少电势（体系功函数一般与氢标电势有特定关系），但这里面也有一些问题，比如多加了的电子需要对能量进行修正。刘智攀老师和Mattew Neurock都做过一些相关的工作，你可以仔细研究他们的文章。

周期性计算包含溶剂模型比较成熟的程序如VASP的VASPsol，QE的Enviroment提供SCCS0，CP2K内置SCCS0等，这些溶剂化模型都不依赖溶剂化半径这项参数，并且声称能够用于界面的溶剂问题。

作者
Author: lzc 时间: 2017-9-6 08:51

卡开发发发表于 2017-9-5 17:02
手册没有具体指明加电场的方向到底是按照怎样的规则，我测试下来感觉电场的三个方向是根据直角坐标来考虑 ...

非常感谢您的耐心解答，您指的给体系外加电子是指利用计算公式-neU来计算是吗？非常感谢您

作者
Author: 卡开发发 时间: 2017-9-6 12:10

lzc 发表于 2017-9-6 08:51
非常感谢您的耐心解答，您指的给体系外加电子是指利用计算公式-neU来计算是吗？非常感谢您

外加电荷之后会响应出一个work function，这个work function去与电势U建立关系，一般差一个常数4.4~4.8左右。

我的理解是，这些外加电荷可以看成是充放电引起的。但多加上的电荷所提供的能量需要经过修正，修正方法有，比如Angew. Chem. 2006, 118, 416 –420里面的公式(3)，如果程序自身能做Markov-Payne修正之类的也可以（比如Quantum Espresso）。

而-neU所指的电荷是电极氧化还原若涉及到ne的电荷转移，这些电子是从电极发射（吸收）电极上反应物的，自由能需要额外-neU。

作者
Author: lzc 时间: 2017-9-8 08:09

卡开发发发表于 2017-9-6 12:10
外加电荷之后会响应出一个work function，这个work function去与电势U建立关系，一般差一个常数4.4~4.8左 ...

非常感谢您的耐心解答，非常感谢，我仔细学习一下

作者
Author: liutianyang0531 时间: 2020-5-21 18:43

卡开发发发表于 2017-9-6 12:10
外加电荷之后会响应出一个work function，这个work function去与电势U建立关系，一般差一个常数4.4~4.8左 ...

您好，想再请教一下，我这里测试功函和外加电子数关系的时候，在DMOL3 setup标签改完charge之后计算功函，发现直接报错显示带电体系不能计算功函，那我该怎么办呢，是因为设置电荷的地方不对吗

作者
Author: 卡开发发 时间: 2020-5-21 18:51

liutianyang0531 发表于 2020-5-21 18:43
您好，想再请教一下，我这里测试功函和外加电子数关系的时候，在DMOL3 setup标签改完charge之后计算功函 ...

只是界面下面有这样的限制，手工改input可以做，但因为偶极问题其实最好还是不要这么做。

作者
Author: liutianyang0531 时间: 2020-5-21 18:58

卡开发发发表于 2020-5-21 18:51
只是界面下面有这样的限制，手工改input可以做，但因为偶极问题其实最好还是不要这么做。

好的谢谢，刚刚接触这一块，我准备试试。

作者
Author: liutianyang0531 时间: 2020-5-21 19:01

卡开发发发表于 2020-5-21 18:51
只是界面下面有这样的限制，手工改input可以做，但因为偶极问题其实最好还是不要这么做。

前辈还想追问一下，我之前用VASP尝试的时候，发现加电子之后，真空能级不是一条直线了，而考虑偶极校正的时候，也是直接报错体系不是立方超胞（实际体系是一个四方结构）。那如果加电子真空能级不平了，功函该怎么确定呢

作者
Author: 卡开发发 时间: 2020-5-21 19:21

liutianyang0531 发表于 2020-5-21 19:01
前辈还想追问一下，我之前用VASP尝试的时候，发现加电子之后，真空能级不是一条直线了，而考虑偶极校正的 ...

其实带过剩电荷的情况光修正偶极也是不行的，上面提到VASP、QE这些专门做溶剂化的插件有专门修正的方式。

作者
Author: liutianyang0531 时间: 2020-5-21 20:47

卡开发发发表于 2020-5-21 19:21
其实带过剩电荷的情况光修正偶极也是不行的，上面提到VASP、QE这些专门做溶剂化的插件有专门修正的方式。

好的，谢谢前辈我再查查

作者
Author: mimicry 时间: 2024-5-22 10:21

zyj19831206 发表于 2017-9-4 19:28
在input文件里面手动加入电场的关键词和方向。

请问具体该怎么操作啊，我在input文件里添加了电场的关键词和方向后，运行总是失败

作者
Author: zhang_hj0206 时间: 2025-7-24 20:53

zyj19831206 发表于 2017-9-4 19:28
在input文件里面手动加入电场的关键词和方向。

您好，我在input文件添加关键词后run file,计算报错，出现”INPUT_DMOL error： unrecognized keyword(s): Electric_filed “这个怎么解决呢？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)