计算化学公社

标题: CP2K各项同性极化率计算问题请教 [打印本页]

作者
Author: xjtuglr 时间: 2024-12-5 14:06
标题: CP2K各项同性极化率计算问题请教
各位老师大家好，最近想通过CP2K对材料的极化率进行计算，参考的方法为2018年发表在PRB上一篇论文的方法DOI：10.1103/PhysRevB.98.064102，在文章的最后写了部分计算方法，但是有一些问题不是很清楚，想请教一下大家。计算过程的截图放在了帖子最后，在计算程序中，我通过PERIODIC_EFIELD字段添加了电场，分别计算了添加方向为100 010 001电场的结果，并且在PRINT字段里选择输出了MOMENTS文件，从而可以获得α，想请问大家如何通过三个不同的α值按照A2公式进行下一步计算？公式中的E(0)，E(|Fi|)是不同方向和不同电场上的Dipole moment输出值吗？还是计算单点后的能量？
(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 3)

作者
Author: sobereva 时间: 2024-12-5 16:18
CP2K直接可以计算解析的极化率，没必要用有限差分的方式
按下文用Multiwfn创建CP2K输入文件的时候把任务选成极化率就完了
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

作者
Author: xjtuglr 时间: 2024-12-5 17:12
本帖最后由 xjtuglr 于 2024-12-5 17:22 编辑

好的，谢谢老师的回答，但是我想继续问一下，通过CP2K获得了POLARIZABILITY TENSOR(stmoic units)和(Angstrom^3)以后，直接将xx，yy，zz三个相加取平均值就是各项同性极化了吗？两个不同大小晶胞的体系算完以后该怎么比较谁的极化大小呢？同时在计算极化的时候可以通过开OT对体系加周期性电场吗老师？

作者
Author: sobereva 时间: 2024-12-5 20:41

xjtuglr 发表于 2024-12-5 17:12
好的，谢谢老师的回答，但是我想继续问一下，通过CP2K获得了POLARIZABILITY TENSOR(stmoic units)和(Angstr ...

极化率（polarizability）不能简称为极化（polarization）
是各向同性极化率
可以根据晶胞体积折算成单位体积的极化率
解析地算极化率加电场干嘛

作者
Author: xjtuglr 时间: 2024-12-6 10:24

sobereva 发表于 2024-12-5 20:41
极化率（polarizability）不能简称为极化（polarization）
是各向同性极化率
可以根据晶胞体积折算成单 ...

好的谢谢卢老师，加外电场是我昨天理解错了，我昨天看那篇PRB里说计算极化率之前有两种方式，一种是固定原子位置，一种是加外电场允许原子弛豫，这个应该是结构优化时候加的外电场，然后用优化之后的结构可以来算极化率，另外请问一下不同材料的单位体积极化率比较有意义吗老师？

作者
Author: sobereva 时间: 2024-12-7 02:59

xjtuglr 发表于 2024-12-6 10:24
好的谢谢卢老师，加外电场是我昨天理解错了，我昨天看那篇PRB里说计算极化率之前有两种方式，一种是固定 ...

有

作者
Author: YwzY 时间: 2025-5-30 15:50

sobereva 发表于 2024-12-7 02:59
有

卢老师，我想问通过cp2k比较一批分子的极化率，来预测对电极极化程度较高的分子时。
我将xx或yy或zz大于700的分子作为预测结果，但是由于极化率张量（xx,yy,zz）的方向与体系的绝对坐标一致，这会不会导致同一个分子在不同的方向，计算出的极化率张量（xx,yy,zz）有较大差异，从而影响我的预测结果

作者
Author: sobereva 时间: 2025-5-30 20:14

YwzY 发表于 2025-5-30 15:50
卢老师，我想问通过cp2k比较一批分子的极化率，来预测对电极极化程度较高的分子时。
我将xx或yy或zz大于 ...

朝向影响极化率分量，不影响极化率总大小（看Multiwfn手册3.27.1节的公式）

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