计算化学公社

标题: 关于体系电场施加方式及电偶极矩计算的疑惑 [打印本页]

作者
Author: Iridium-191 时间: 2022-1-9 09:59
标题: 关于体系电场施加方式及电偶极矩计算的疑惑
本帖最后由 Iridium-191 于 2022-1-9 10:07 编辑

各位老师好，存在两点疑惑，首先是关于电场施加方式：在CP2K软件中，提供三种施加电场的方式，包括：
EFIELD、PERIODIC_EFIELD和EXTERNAL_POTENTIAL。个人不太清楚这三个函数的使用场合。以下是几点主要的疑惑：
1. 周期性体系AIMD计算施加电场应该使用哪个函数？个人查阅相关文献觉得应该使用PERIODIC_EFIELD，但是CP2K目前不支持在计算STRESS_TENSOR的同时施加PERIODIC_EFIELD。
2. 几何优化时如果想做电场下的几何优化，是需要选用哪个函数？个人查阅相关资料认为应该使用EXTERNAL_POTENTIAL，但是大电场下体系出现真空层，与实际严重不符。
3. 个人始终不太理解CP2K Manual中EFIELD板块下

parameters for finite, time dependent, nonperiodic electric fields. For static fields use EXTERNAL_POTENTIAL

的相关表述，他的意思是，施加与时间无关的匀强电场就一定要使用EXTERNAL_POTENTIAL来进行么？

第二点疑惑是关于Berry-Phase方法计算体系的电偶极矩（DFT-PRINT-MOMENT)
使用Berry-Phase方法计算体系每一步AIMD运算的电偶极矩，发生非常严重的波动的原因是什么（如图所示）？如果不采纳Berry-Phase方法计算得到的电矩，采用经典的计算轨迹文件正负电荷中心的方法得到的结果足够可信么？

作者
Author: sobereva 时间: 2022-1-9 10:07
一般DFT计算用EFIELD（限孤立体系）

2023-May-19补充：

(, 下载次数 Times of downloads: 12)

(, 下载次数 Times of downloads: 7)

作者
Author: Iridium-191 时间: 2022-1-9 10:08

sobereva 发表于 2022-1-9 10:07
一般DFT计算用EFIELD

好的，非常感谢老师

作者
Author: goldNAN 时间: 2023-5-18 00:08

Iridium-191 发表于 2022-1-9 10:08
好的，非常感谢老师

对于周期性体系要慎重，直接用efiled 软件就当非周期性体系出来了算偶极矩可能不对

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