计算化学公社

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

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

作者
Author: sobereva 时间: 2023-12-7 03:21
CP2K直接就能算拉曼光谱，然后用Multiwfn绘制，特意输出这个没意义

作者
Author: wangknano 时间: 2023-12-7 12:34

sobereva 发表于 2023-12-7 03:21
CP2K直接就能算拉曼光谱，然后用Multiwfn绘制，特意输出这个没意义

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

作者
Author: sobereva 时间: 2023-12-7 12:39
没关键词
自己想办法做有限差分算

作者
Author: wangknano 时间: 2023-12-7 14:30

sobereva 发表于 2023-12-7 12:39
没关键词
自己想办法做有限差分算

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

作者
Author: wangknano 时间: 2023-12-8 15:54

sobereva 发表于 2023-12-7 12:39
没关键词
自己想办法做有限差分算

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

作者
Author: sobereva 时间: 2023-12-8 18:23

wangknano 发表于 2023-12-8 15:54
老师您好，我看了一下振动分析的输出文件，发现X+DX和X-DX的极化率差距很大，请问我该用哪个来计算这个偏 ...

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

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

作者
Author: wangknano 时间: 2023-12-8 22:43

sobereva 发表于 2023-12-8 18:23
显然用双向差分，没有这种矛盾

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

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

作者
Author: sobereva 时间: 2023-12-9 05:37

wangknano 发表于 2023-12-8 22:43
老师您好，再打扰您一下，振动分析输出文件中某一原子X+DX极化率某一项是3700多，但X-DX的极化率对应项只 ...

尝试和自己手动做有限差分得到的极化率进行对比（即两个结构下分别跑极化率计算任务），并且考察极化率随有限差分步长的影响

作者
Author: Superginger 时间: 2023-12-18 14:44

sobereva 发表于 2023-12-9 05:37
尝试和自己手动做有限差分得到的极化率进行对比（即两个结构下分别跑极化率计算任务），并且考察极化率随 ...

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

作者
Author: sobereva 时间: 2023-12-18 23:33

Superginger 发表于 2023-12-18 14:44
老师好，想咨询一下CP2K软件计算的极化率是在什么波长下的数值呢？这个可以自己另外设置吗？

静态的，相当于波长无限大

作者
Author: Superginger 时间: 2023-12-19 11:22

sobereva 发表于 2023-12-18 23:33
静态的，相当于波长无限大

谢谢老师的解答

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