计算化学公社

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

[CP2K] CP2K计算二维过渡金属硫族化合物频率分析,怎么设置才能输出极化率随正则坐标的偏导

[复制链接 Copy URL]

9

帖子

0

威望

205

eV
积分
214

Level 3 能力者

各位老师好,最近希望用CP2K计算二维过渡金属硫族化合物结构的拉曼光谱,想请问一下怎么设置相关关键词才可以把极化率随正则坐标的偏导给输出出来?
还有就是想构建异质结加一个Au衬底,请问Au衬底的晶格常数如果和实际情况有所不同,对计算结果是否影响很大?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

2#
发表于 Post on 2023-12-7 03:21:44 | 只看该作者 Only view this author
CP2K直接就能算拉曼光谱,然后用Multiwfn绘制,特意输出这个没意义
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

9

帖子

0

威望

205

eV
积分
214

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-12-7 12:34:33 | 只看该作者 Only view this author
sobereva 发表于 2023-12-7 03:21
CP2K直接就能算拉曼光谱,然后用Multiwfn绘制,特意输出这个没意义

老师您好,是这样的,我想用极化率随正则坐标的偏导再另外计算加上选率的拉曼光谱,还是想咨询一下是否有相关输出关键词可以输出极化率随正则坐标的偏导,可以在inp里直接设置一下?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

4#
发表于 Post on 2023-12-7 12:39:32 | 只看该作者 Only view this author
没关键词
自己想办法做有限差分算
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

9

帖子

0

威望

205

eV
积分
214

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-12-7 14:30:34 | 只看该作者 Only view this author
sobereva 发表于 2023-12-7 12:39
没关键词
自己想办法做有限差分算

好的好的,谢谢老师您的及时回复

9

帖子

0

威望

205

eV
积分
214

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2023-12-8 15:54:01 | 只看该作者 Only view this author
sobereva 发表于 2023-12-7 12:39
没关键词
自己想办法做有限差分算

老师您好,我看了一下振动分析的输出文件,发现X+DX和X-DX的极化率差距很大,请问我该用哪个来计算这个偏导?或者CP2K是否可以另外写任务来做这个有限差分呢?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

7#
发表于 Post on 2023-12-8 18:23:27 | 只看该作者 Only view this author
wangknano 发表于 2023-12-8 15:54
老师您好,我看了一下振动分析的输出文件,发现X+DX和X-DX的极化率差距很大,请问我该用哪个来计算这个偏 ...

显然用双向差分,没有这种矛盾

没有对应的专门的任务。这种事自己写程序/脚本实现
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

9

帖子

0

威望

205

eV
积分
214

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2023-12-8 22:43:59 | 只看该作者 Only view this author
sobereva 发表于 2023-12-8 18:23
显然用双向差分,没有这种矛盾

没有对应的专门的任务。这种事自己写程序/脚本实现

老师您好,再打扰您一下,振动分析输出文件中某一原子X+DX极化率某一项是3700多,但X-DX的极化率对应项只有0.16,而每个副本的输出文件中显示的极化率对应项是3700多,X+DX感觉数值比较合理,但此处的X-DX极化率张量为什么会这么小?(相应图片和附件也已经上传)

X DX.png (219.73 KB, 下载次数 Times of downloads: 7)

X DX.png

X-DX.png (213.67 KB, 下载次数 Times of downloads: 7)

X-DX.png

raman-r-0.out.png (370.82 KB, 下载次数 Times of downloads: 7)

raman-r-0.out.png

raman.out.zip

221.65 KB, 下载次数 Times of downloads: 1

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

9#
发表于 Post on 2023-12-9 05:37:47 | 只看该作者 Only view this author
wangknano 发表于 2023-12-8 22:43
老师您好,再打扰您一下,振动分析输出文件中某一原子X+DX极化率某一项是3700多,但X-DX的极化率对应项只 ...

尝试和自己手动做有限差分得到的极化率进行对比(即两个结构下分别跑极化率计算任务),并且考察极化率随有限差分步长的影响
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

12

帖子

0

威望

197

eV
积分
209

Level 3 能力者

10#
发表于 Post on 2023-12-18 14:44:40 | 只看该作者 Only view this author
sobereva 发表于 2023-12-9 05:37
尝试和自己手动做有限差分得到的极化率进行对比(即两个结构下分别跑极化率计算任务),并且考察极化率随 ...

老师好,想咨询一下CP2K软件计算的极化率是在什么波长下的数值呢?这个可以自己另外设置吗?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

11#
发表于 Post on 2023-12-18 23:33:50 | 只看该作者 Only view this author
Superginger 发表于 2023-12-18 14:44
老师好,想咨询一下CP2K软件计算的极化率是在什么波长下的数值呢?这个可以自己另外设置吗?

静态的,相当于波长无限大
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

12

帖子

0

威望

197

eV
积分
209

Level 3 能力者

12#
发表于 Post on 2023-12-19 11:22:07 | 只看该作者 Only view this author
sobereva 发表于 2023-12-18 23:33
静态的,相当于波长无限大

谢谢老师的解答

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 11:43 , Processed in 0.209796 second(s), 23 queries , Gzip On.

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