计算化学公社

 找回密码 Forget password
 注册 Register

求助:有机晶体计算红外和Raman光谱,使用带弥散的基组遇到的SCF不收敛问题

查看数: 846 | 评论数: 4 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2024-1-1 21:05

正文摘要:

各位老师好,我的目标是模拟计算有机晶体的红外和拉曼振动光谱,单个晶胞约200个原子,使用的是CP2K 2023.1版本。(在结构优化阶段就出现的问题)输入文件使用的是Multiwfn生成的,由于要计算拉曼强度,打算使用带弥 ...

回复 Reply

Spadek 发表于 Post on 2024-1-2 12:00:54
sobereva 发表于 2024-1-2 04:14
量子化学计算带弥散函数的基组通常都没法用于固体计算使用。哪怕一些不带弥散的3-zeta基组,都由于最外层的 ...

好的好的,谢谢老师解答
sobereva 发表于 Post on 2024-1-2 04:14:28
量子化学计算带弥散函数的基组通常都没法用于固体计算使用。哪怕一些不带弥散的3-zeta基组,都由于最外层的基函数的指数偏小,可能带来固体计算时SCF不稳定性问题。

专门为固体计算优化的pob-TZVP-rev2基组远比量子化学常用的各种3-zeta基组适合周期性计算
CP2K常用的MOLOPT系列虽然有指数较小的高斯函数,但由于收缩时和指数较大的绑在一起,也没一般带弥散函数的基组算固体时的数值不稳定性问题

另外,本身固体计算时就没有分子体系计算的那种在分子外围有电子密度很低的尾巴区域,所以算分子体系对弥散函数的要求也根本不能简单套用在固体计算上,很多算分子需要带弥散函数的情况对于算固体并不需要。
Spadek 发表于 Post on 2024-1-1 21:33:56
嗷嗷,好的,谢谢老师。那就是说目前计算固体体系,加弥散的话大概率都是有一些问题的,我说之前用Crystal算咋一直不行呢
zjxitcc 发表于 Post on 2024-1-1 21:12:23
固体计算与分子计算毕竟不一样,我不建议你把加弥散的经验套用过去,否则基函数线性相关巨严重,周期性计算程序离稳健的SCF还有很长的路要走。建议你还是找个合适的没有弥散的triple-zeta基组算一算(例如pob-TZVP)。

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
卡开发发 + 3 我很赞同

查看全部评分 View all ratings

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

GMT+8, 2024-11-25 13:00 , Processed in 0.199819 second(s), 31 queries , Gzip On.

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