计算化学公社

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

[CP2K] 用CP2K做金刚石单轴(z)应变的仿真求助

[复制链接 Copy URL]

5

帖子

0

威望

39

eV
积分
44

Level 2 能力者

大家好,我想咨询下,用CP2K做金刚石单轴(z)应变的仿真,应该是调整晶格参数,如做z轴拉伸金刚石supercell3 3 3的10*10*10的晶胞,是将晶胞参数调整为10*10*10.1,然后再跑GEO_ OPT吗?还是做CELL_OPT,然后逐渐调整EXTERNAL_PRESSURE来实现啊?

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

2#
发表于 Post on 2021-11-13 00:38:47 | 只看该作者 Only view this author
看你是想考察晶胞z轴变化特定距离时原子坐标怎么变化,还是想考察特定外压下结构怎么变化,取决于你
北京科音自然科学研究中心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

5

帖子

0

威望

39

eV
积分
44

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2021-11-13 14:22:02 | 只看该作者 Only view this author
sobereva 发表于 2021-11-13 00:38
看你是想考察晶胞z轴变化特定距离时原子坐标怎么变化,还是想考察特定外压下结构怎么变化,取决于你

我是想计算晶体受到外压力或外拉力作用下,晶体的能量变化,现在用的external_pressure,发现能量没有什么变化,所以在考虑是不是需要做应变

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

4#
发表于 Post on 2021-11-14 03:27:17 | 只看该作者 Only view this author
yb6231290 发表于 2021-11-13 14:22
我是想计算晶体受到外压力或外拉力作用下,晶体的能量变化,现在用的external_pressure,发现能量没有什 ...

数值设更大再试
北京科音自然科学研究中心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

148

帖子

0

威望

888

eV
积分
1036

Level 4 (黑子)

5#
发表于 Post on 2021-11-14 03:35:20 | 只看该作者 Only view this author
做之前最好做一下cell opt得到relaxed cell

5

帖子

0

威望

39

eV
积分
44

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2021-11-15 15:13:01 | 只看该作者 Only view this author
Aridea 发表于 2021-11-14 03:35
做之前最好做一下cell opt得到relaxed cell

谢谢回复,请问relaxed cell后再加geo,是可以直接上0.01 0.02 ... 0.05这样做吗?

5

帖子

0

威望

39

eV
积分
44

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2021-11-15 15:14:36 | 只看该作者 Only view this author
本帖最后由 yb6231290 于 2021-11-15 15:23 编辑

我做pressure已经加拉力到20GPa了,可以看到晶格常数的变化,但是里面不同位置,不同方向的缺陷(空位和掺杂)完全没有变化,不知道是什么原因啊,我用的PBE泛函,我看网上说这里应该把potential GTH改成ALL,您看是否需要呢?

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

8#
发表于 Post on 2021-11-16 01:11:46 | 只看该作者 Only view this author
yb6231290 发表于 2021-11-15 15:14
我做pressure已经加拉力到20GPa了,可以看到晶格常数的变化,但是里面不同位置,不同方向的缺陷(空位和掺 ...

跟potential完全没直接关系。
我不知道你说的没变化是怎么回事,晶格尺寸都变了,不可能原子坐标都一点不变,顶多是肉眼不好分辨的程度
北京科音自然科学研究中心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

5

帖子

0

威望

39

eV
积分
44

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2021-11-16 12:54:06 | 只看该作者 Only view this author
本帖最后由 yb6231290 于 2021-11-16 14:30 编辑
sobereva 发表于 2021-11-16 01:11
跟potential完全没直接关系。
我不知道你说的没变化是怎么回事,晶格尺寸都变了,不可能原子坐标都一点 ...

老师您好,是这样的,我想做111面的金刚石拉伸,同时观察金刚石内部不同方向NV的能量变化,这个在其他文献中的总结中可知:沿着111的NV方向和非111的NV方向随着z轴的拉力或压力,逐渐出现能量的分叉,文献结果如图所示。([img][/img]见附件图片),文献中使用vasp计算,计算中采用了672个原子,0应变时,二者的能量差在0.0018eV,应变到了0.01能量差大概在0.1eV,未使用k点。
但是我自己做的沿111面拉伸,原子数也为672个,未使用k点,cutoff为400,泛函和基组分别为:GTH-PBE和DZVP-MOLOPT-SR-GTH,同时添加了dftd3色散,OT收敛精度1e-6,计算cellopt,external_pressure从0-20GPa,做了10组,其中10GPa时金刚石z轴应变约为0.01,但发现这两个nv方向的能量差始终保持在0.0019-0.0023eV,所以我想请问下到底是什么原因导致了cp2k和vasp差异如此的大?

文献结果-应变与能量差.png (243.04 KB, 下载次数 Times of downloads: 10)

文献结果-应变与能量差.png

7

帖子

0

威望

29

eV
积分
36

Level 2 能力者

10#
发表于 Post on 2024-11-25 09:07:50 | 只看该作者 Only view this author
我想用cp2k做一个第一性原理分子动力学的拉伸,改变z轴方向的长度来进行拉伸,但是我不知道怎么写输入文件,能指教一下吗?

5

帖子

0

威望

31

eV
积分
36

Level 2 能力者

11#
发表于 Post on 2024-11-26 20:32:57 | 只看该作者 Only view this author
xjc 发表于 2024-11-25 09:07
我想用cp2k做一个第一性原理分子动力学的拉伸,改变z轴方向的长度来进行拉伸,但是我不知道怎么写输入文件, ...

您好,请问这个问题解决了嘛,我也想学习一下。

7

帖子

0

威望

29

eV
积分
36

Level 2 能力者

12#
发表于 Post on 2024-12-12 14:10:01 | 只看该作者 Only view this author
huangyh 发表于 2024-11-26 20:32
您好,请问这个问题解决了嘛,我也想学习一下。

还没有呢,你有什么想法吗

4

帖子

0

威望

61

eV
积分
65

Level 2 能力者

13#
发表于 Post on 2025-4-11 14:19:59 | 只看该作者 Only view this author
xjc 发表于 2024-12-12 14:10
还没有呢,你有什么想法吗

请问实现了吗

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 02:00 , Processed in 0.180126 second(s), 29 queries , Gzip On.

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