计算化学公社

标题: 想问下cp2k计算分子晶体的极性是否可以？ [打印本页]

作者
Author: xuanzi 时间: 2021-1-4 20:35
标题: 想问下cp2k计算分子晶体的极性是否可以？
要计算极化率，第一超级化率，但用没有找到用cp2k计算的相关的例子....谢谢各位老师

作者
Author: doublezhang 时间: 2021-1-4 22:28
https://manual.cp2k.org/trunk/CP ... VAL/DFT/EFIELD.html
数值求解吧

作者
Author: biogon 时间: 2021-1-5 16:28

doublezhang 发表于 2021-1-4 22:28
https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/EFIELD.html
数值求解吧

数值计算的话算极化率还有超极化率太昂贵了

作者
Author: doublezhang 时间: 2021-1-5 17:01

biogon 发表于 2021-1-5 16:28
数值计算的话算极化率还有超极化率太昂贵了

那不如买个crystal?
或者免费的abinit也可以算。
浙江师大的蓝老师自己写了一个从CASTEP输出算NLO的小程序，可以发邮件索要。

作者
Author: biogon 时间: 2021-1-5 17:13

doublezhang 发表于 2021-1-5 17:01
那不如买个crystal?
或者免费的abinit也可以算。
浙江师大的蓝老师自己写了一个从CASTEP输出算NLO的小 ...

那些程序都有高阶解析导数嘛？我刚发现没有用簇模型算固体极化率的，在想能不能用簇模型算极化率这些东西

作者
Author: xuanzi 时间: 2021-1-5 19:17
你们的对话表示看不太懂

大佬们

作者
Author: doublezhang 时间: 2021-1-5 19:25

biogon 发表于 2021-1-5 17:13
那些程序都有高阶解析导数嘛？我刚发现没有用簇模型算固体极化率的，在想能不能用簇模型算极化率这些东西

crystal绝对是解析求导
abinit好像是根据dfpt？我不太懂。
蓝老师的程序应该是基于SOS态求和方法的

作者
Author: biogon 时间: 2021-1-5 19:46

doublezhang 发表于 2021-1-5 19:25
crystal绝对是解析求导
abinit好像是根据dfpt？我不太懂。
蓝老师的程序应该是基于SOS态求和方法的

SOS的话也是个办法，不用依赖于高阶导数了

作者
Author: lhd 时间: 2022-1-10 21:47
请问您找到计算的方法啦吗，想咨询一下

作者
Author: sobereva 时间: 2022-1-22 17:20

lhd 发表于 2022-1-10 21:47
请问您找到计算的方法啦吗，想咨询一下

目前最新版Multiwfn创建CP2K输入文件功能里直接就能选择计算极化率任务
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

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

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