计算化学公社

标题: 关于TDDFT的各个态跃迁振子强度全为0的问题 [打印本页]

作者
Author: din5g 时间: 2022-8-11 19:49
标题: 关于TDDFT的各个态跃迁振子强度全为0的问题
本帖最后由 din5g 于 2022-8-16 14:12 编辑

各位老师好，我想用cp2k（9.1版本）计算WS2单层的tddft的光吸收谱，为此我构建了两个不同大小的结构，扩胞为3*3和5*5。3*3共27个原子，a=b=9.56 埃，c=18 埃，而5*5共75个原子，a=b=15.94埃。通过VASP（PBE）计算得到的带隙值均为1.81 eV，均为在K点的直接带隙。利用CP2K计算其TDDFT的光吸收谱，采用PBE/DZVP-MOLOPT-SR-GTH赝势基组。问题在于通过CP2K计算的TDDFT的振子强度不同，3*3体系看起来正常，激发能位置在1.90 eV 振子强度为3.81左右（如下图），而5*5体系的振子强度全为0（如下图），检查out文档也是如此，不知5*5体系的问题在哪，向各位老师请教，为什么会有这么大的区别？
按照sob老师的建议，已计算从2*2到7*7的tddft，见下图，但其规律非常怪异，不知问题在哪？

作者
Author: sobereva 时间: 2022-8-12 08:17
可以做个激发能、振子强度随扩胞倍数的收敛性测试，从2倍做到7倍
说不定3*3的时候定性错误

作者
Author: din5g 时间: 2022-8-14 12:43

sobereva 发表于 2022-8-12 08:17
可以做个激发能、振子强度随扩胞倍数的收敛性测试，从2倍做到7倍
说不定3*3的时候定性错误

多谢sob老师的回复，按照您的建议，已计算从2*2到7*7的tddft，但其规律非常怪异，只有6*6看起来才像正常的。

作者
Author: din5g 时间: 2022-8-16 14:14

sobereva 发表于 2022-8-12 08:17
可以做个激发能、振子强度随扩胞倍数的收敛性测试，从2倍做到7倍
说不定3*3的时候定性错误

sob老师，为啥收敛测试完全不收敛，没有呈现规律，操作步骤都是一样的，难道是模型的原因？

作者
Author: sobereva 时间: 2022-8-17 17:43

din5g 发表于 2022-8-16 14:14
sob老师，为啥收敛测试完全不收敛，没有呈现规律，操作步骤都是一样的，难道是模型的原因？

建议去CP2K官方论坛问问

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