计算化学公社

 找回密码 Forget password
 注册 Register
Views: 591|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] CP2K各项同性极化率计算问题请教

[复制链接 Copy URL]

33

帖子

0

威望

323

eV
积分
356

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师大家好,最近想通过CP2K对材料的极化率进行计算,参考的方法为2018年发表在PRB上一篇论文的方法DOI:10.1103/PhysRevB.98.064102,在文章的最后写了部分计算方法,但是有一些问题不是很清楚,想请教一下大家。计算过程的截图放在了帖子最后,在计算程序中,我通过PERIODIC_EFIELD字段添加了电场,分别计算了添加方向为100 010 001电场的结果,并且在PRINT字段里选择输出了MOMENTS文件,从而可以获得α,想请问大家如何通过三个不同的α值按照A2公式进行下一步计算?公式中的E(0),E(|Fi|)是不同方向和不同电场上的Dipole moment输出值吗?还是计算单点后的能量?


西安交通大学等离子体与微波电子学研究所
THz无损检测与光谱模拟
gonglr@stu.xjtu.edu.cn

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2024-12-5 16:18:21 | 只看该作者 Only view this author
CP2K直接可以计算解析的极化率,没必要用有限差分的方式
按下文用Multiwfn创建CP2K输入文件的时候把任务选成极化率就完了
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587http://bbs.keinsci.com/thread-21668-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

33

帖子

0

威望

323

eV
积分
356

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-12-5 17:12:35 | 只看该作者 Only view this author
本帖最后由 xjtuglr 于 2024-12-5 17:22 编辑

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

202412051709211801..png (32.95 KB, 下载次数 Times of downloads: 2)

202412051709211801..png
西安交通大学等离子体与微波电子学研究所
THz无损检测与光谱模拟
gonglr@stu.xjtu.edu.cn

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

4#
发表于 Post on 2024-12-5 20:41:54 | 只看该作者 Only view this author
xjtuglr 发表于 2024-12-5 17:12
好的,谢谢老师的回答,但是我想继续问一下,通过CP2K获得了POLARIZABILITY TENSOR(stmoic units)和(Angstr ...

极化率(polarizability)不能简称为极化(polarization)
是各向同性极化率
可以根据晶胞体积折算成单位体积的极化率
解析地算极化率加电场干嘛
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

33

帖子

0

威望

323

eV
积分
356

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-12-6 10:24:44 | 只看该作者 Only view this author
sobereva 发表于 2024-12-5 20:41
极化率(polarizability)不能简称为极化(polarization)
是各向同性极化率
可以根据晶胞体积折算成单 ...

好的谢谢卢老师,加外电场是我昨天理解错了,我昨天看那篇PRB里说计算极化率之前有两种方式,一种是固定原子位置,一种是加外电场允许原子弛豫,这个应该是结构优化时候加的外电场,然后用优化之后的结构可以来算极化率,另外请问一下不同材料的单位体积极化率比较有意义吗老师?
西安交通大学等离子体与微波电子学研究所
THz无损检测与光谱模拟
gonglr@stu.xjtu.edu.cn

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

6#
发表于 Post on 2024-12-7 02:59:45 | 只看该作者 Only view this author
xjtuglr 发表于 2024-12-6 10:24
好的谢谢卢老师,加外电场是我昨天理解错了,我昨天看那篇PRB里说计算极化率之前有两种方式,一种是固定 ...

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

9

帖子

0

威望

135

eV
积分
144

Level 2 能力者

7#
发表于 Post on 2025-5-30 15:50:10 | 只看该作者 Only view this author

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

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

8#
发表于 Post on 2025-5-30 20:14:18 | 只看该作者 Only view this author
YwzY 发表于 2025-5-30 15:50
卢老师,我想问通过cp2k比较一批分子的极化率,来预测对电极极化程度较高的分子时。
我将xx或yy或zz大于 ...

朝向影响极化率分量,不影响极化率总大小(看Multiwfn手册3.27.1节的公式)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 06:59 , Processed in 0.150308 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list