简单量子化学问题答疑专帖

糖糖糖豆9988 · 发表于 Post on 2023-5-9 16:39:10

请教各位老师，G4计算结果如图，从数量级看 G4 Enthalpy= -10892.840833和E(CCSD(T))= -230.546023 的差异就大得离谱，是计算有什么bug吗？

KSeGaSn · 发表于 Post on 2023-5-9 22:28:15

各位老师好，请问我们有办法单独衡量某几个原子/片段对于整个分子旋轨耦合效应/系间窜跃速率的贡献呢？比如对于常见的包含过渡金属的配合物来说，能否单独比较两个配合物中金属原子对SOC/ISC的影响/贡献呢？

@CHQ啦 · 发表于 Post on 2023-5-10 14:52:06

本帖最后由 @CHQ啦于 2023-5-10 14:54 编辑

老师，请问在做IRC时，出现L123错误，如图，怎么办

sobereva · 发表于 Post on 2023-5-10 15:32:53

@CHQ啦发表于 2023-5-10 14:52
老师，请问在做IRC时，出现L123错误，如图，怎么办

老生常谈的问题。遇到问题别急着发帖，先搜在公社论坛首页google框搜了，若10分钟找不到答案再问
立马就能搜到此文
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

zgdx · 发表于 Post on 2023-5-10 16:17:04

请问老师，用DFT泛函优化结构时已经收敛了，将最后一帧的收敛结构保存为新的gjf，改用QCISD做单点计算，基组不变，结果不收敛。怎么解决呢？
还有一个问题是计算单点能的log文件中的能量收敛标准是下图标记的吗？

sobereva · 发表于 Post on 2023-5-10 16:39:37

zgdx 发表于 2023-5-10 16:17
请问老师，用DFT泛函优化结构时已经收敛了，将最后一帧的收敛结构保存为新的gjf，改用QCISD做单点计算，基 ...

解决SCF不收敛问题的方法
http://sobereva.com/61

是

sobereva · 发表于 Post on 2023-5-10 22:11:22

糖糖糖豆9988 发表于 2023-5-9 16:39
请教各位老师，G4计算结果如图，从数量级看 G4 Enthalpy= -10892.840833和E(CCSD(T))= ...

有可能中途有的计算收敛到不稳定波函数之类。不正常。没细节没法说

sobereva · 发表于 Post on 2023-5-10 22:11:49

qcn1211 发表于 2023-5-9 13:11
请问在windows系统中如何同时把多个log文件分别转换成gjf文件？

一键把所有gjf文件转成xyz文件、把所有Gaussian输出文件转成gjf文件的脚本
http://sobereva.com/530（http://bbs.keinsci.com/thread-16161-1-1.html）

sobereva · 发表于 Post on 2023-5-10 22:16:17

护卫天使发表于 2023-5-9 10:12
老师您好，第一次在CP2K程序跑AIMD，有几个问题请教一下您。1、如图您给的例子，假如实用Packmol构建体系 ...

xtb不适合跑周期性体系动力学
CP2K支持GFN1-xTB，适合周期性体系的动力学。无论哪个程序，GFN-xTB的并行效率都不怎么样，核数给得太多甚至反倒速度慢。自己测试便知。

CP2K里建议用PBE-D3(BJ)，GFN1-xTB跑这个未必能得到合理结果，虽然可以试试

护卫天使 · 发表于 Post on 2023-5-10 22:38:51

本帖最后由护卫天使于 2023-5-10 22:40 编辑

sobereva 发表于 2023-5-10 22:16
xtb不适合跑周期性体系动力学
CP2K支持GFN1-xTB，适合周期性体系的动力学。无论哪个程序，GFN-xTB的并行 ...

感谢老师的回答。另外，我在Multiwfn导出的CP2K输入文件，尝试了选择TZVP-GTH和TZVP-MOLOPT-GTH基组都表示不支持稀土元素，请问这21个基组有适用稀土的吗？

sobereva · 发表于 Post on 2023-5-10 23:10:23

护卫天使发表于 2023-5-10 22:38
感谢老师的回答。另外，我在Multiwfn导出的CP2K输入文件，尝试了选择TZVP-GTH和TZVP-MOLOPT-GTH基组都表 ...

北京科音CP2K第一性原理培训班ppt：

需要自己编辑输入文件

勇往直前 · 发表于 Post on 2023-5-11 15:32:54

如果分子中带一正电荷以及一负电荷，gjf文件中的电荷该如何设置？

sobereva · 发表于 Post on 2023-5-11 15:47:01

勇往直前发表于 2023-5-11 15:32
如果分子中带一正电荷以及一负电荷，gjf文件中的电荷该如何设置？

设0

吴双 · 发表于 Post on 2023-5-11 16:21:36

想跟大家请教一下，计算自旋极化单重态的具体方法是什么，我看了Sob老师的帖子http://sobereva.com/82 谈谈片段组合波函数与自旋极化单重态，文中老师的建议先尝试比如# UB3LYP/6-31G* guess=mix nosymm stable关键词，因为又快又省事，如果提示波函数是稳定的，而且自旋密度确实不是处处为0（怎么绘制自旋密度看http://sobereva.com/353），那就OK了。这块有点不理解，跟几何优化的前后顺序应该怎么处理呢。我现在的操作是UB3LYP opt freq guess=mix对几何结构进行优化，自旋多重度是1，看输出文件中自选密度有0.000002，然后在用UB3LYP/6-31G* guess=mix nosymm stable检验一下波函数的稳定性么？

sxh3231796 · 发表于 Post on 2023-5-11 20:32:13

wzkchem5 发表于 2023-5-4 20:14
一般来说，研究过渡金属配合物，应该把过渡金属的价层d轨道全部选进来。因为即便有些d轨道不参与成键，这 ...

想请问下老师，这种情况下即使轨道离homo很远也要囊括进去吗？之前查看轨道时，只看了homo，lumo前后10个轨道，是要筛选所有轨道，选取活性空间吗？

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 简单量子化学问题答疑专帖

IRC出现L123

Gaussian计算单点能不收敛怎么办？

电荷设置