想问下cp2k计算分子晶体的极性是否可以？

xuanzi · 发表于 Post on 2021-1-4 20:35:52

要计算极化率，第一超级化率，但用没有找到用cp2k计算的相关的例子....谢谢各位老师

doublezhang · 发表于 Post on 2021-1-4 22:28:52

biogon · 发表于 Post on 2021-1-5 16:28:33

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

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

doublezhang · 发表于 Post on 2021-1-5 17:01:45

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

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

biogon · 发表于 Post on 2021-1-5 17:13:55

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

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

xuanzi · 发表于 Post on 2021-1-5 19:17:55

你们的对话表示看不太懂

大佬们

doublezhang · 发表于 Post on 2021-1-5 19:25:35

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

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

biogon · 发表于 Post on 2021-1-5 19:46:18

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

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

lhd · 发表于 Post on 2022-1-10 21:47:05

请问您找到计算的方法啦吗，想咨询一下

sobereva · 发表于 Post on 2022-1-22 17:20:51

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 想问下cp2k计算分子晶体的极性是否可以？