计算化学公社

 找回密码 Forget password
 注册 Register

在cp2k计算跃迁偶极时遇到报错

查看数: 106 | 评论数: 3 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-7-15 13:08

正文摘要:

在使用cp2k 计算周期性分子晶体时遇到如图中的报错 我希望计算并输出S1态的跃迁偶极矩,但报错显示cp2k好像无法输出周期性体系的跃迁偶极矩 报错信息:Periodic moments not implemented with TDDFT

回复 Reply

sobereva 发表于 Post on 2025-7-18 06:23:51
wangxc 发表于 2025-7-17 15:39
我在Gaussian的计算振动分辨电子光谱的fchk文件中找到ETran state values项,我想这应该是电子跃迁偶极矩 ...

跃迁偶极矩在Gaussian输出文件里直接就有,下文明确说了,根本没必要到fchk里去读,若不需要考虑振动耦合也完全不需要算振动电子光谱
使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱
http://sobereva.com/462http://bbs.keinsci.com/thread-12114-1-1.html

而且CP2K也不支持振动电子光谱的计算

如北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)详细讲CP2K电子激发计算的这页幻灯片所示,CP2K做TDDFT直接就输出了跃迁偶极矩,太好找了,完全没有藏着掖着


wangxc 发表于 Post on 2025-7-17 15:39:36
sobereva 发表于 2025-7-16 11:45
分清楚激发态的偶极矩,和基态到激发态的跃迁偶极矩,完全是两码事
做TDDFT自动就会计算后者,你写下面内 ...

我在Gaussian的计算振动分辨电子光谱的fchk文件中找到ETran state values项,我想这应该是电子跃迁偶极矩的数据。
现在换成周期性体系,我想在cp2k中计算中获取相应的电子跃迁偶极矩的数据,我应该添加什么关键词,或者在哪个文件可以找到这部分的输出
sobereva 发表于 Post on 2025-7-16 11:45:18
分清楚激发态的偶极矩,和基态到激发态的跃迁偶极矩,完全是两码事
做TDDFT自动就会计算后者,你写下面内容是多余的,自找这个麻烦
      &MOMENTS
        PERIODIC T #Use Berry phase formula (T) or simple operator (F), the latter normally applies to isolated systems
      &END MOMENTS

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

GMT+8, 2025-8-14 04:49 , Processed in 0.386013 second(s), 25 queries , Gzip On.

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