计算化学公社

标题: 高斯可以代替MOLCAS做文中这部分计算吗 [打印本页]

作者
Author: 天天121 时间: 2022-3-8 02:32
标题: 高斯可以代替MOLCAS做文中这部分计算吗
请问各位老师高斯可以代替Molcas做文中的QM/MM计算吗？这两个软件各自有什么样的优势呢？

作者
Author: sobereva 时间: 2022-3-8 02:39
不能
Gaussian不支持CASPT2。但如果不是非得用CASPT2的话，可以用Gaussian中的TDDFT或更好的EOM-CCSD来算。

另外，你提到的文献的做法不当，MM区域原子用Mulliken电荷表现十分不合适，对静电势重现性不理想。应当用拟合静电势电荷，看下文
基于背景电荷计算分子在晶体环境中的吸收光谱
http://sobereva.com/579（http://bbs.keinsci.com/thread-20960-1-1.html）

作者
Author: 天天121 时间: 2022-3-9 10:02

sobereva 发表于 2022-3-8 02:39
不能
Gaussian不支持CASPT2。但如果不是非得用CASPT2的话，可以用Gaussian中的TDDFT或更好的EOM-CCSD来算 ...

sob老师
1.我在网上下载好了cif文件然后扩胞到大概500个原子需要先用cp2k做个结构优化再截取15个分子左右的团簇
2.还是我直接用VMD截取大概15个分子的团簇然后再用ONIOM做结构优化单点能之类的任务呢

这两种您觉得哪种比较精确呢

作者
Author: 天天121 时间: 2022-3-9 17:23

sobereva 发表于 2022-3-8 02:39
不能
Gaussian不支持CASPT2。但如果不是非得用CASPT2的话，可以用Gaussian中的TDDFT或更好的EOM-CCSD来算 ...

Sob 老师我已经学会了您用背景电荷法绘制吸收光谱请问您有什么好方法绘制有机晶体的荧光吗
我现在激发态300多个原子已经用ONIOM进行了优化，关键词如下：
# opt oniom(td=(nstates=3,root=1)/b3lyp/6-311g(d,p):pm7) geom=connectivity

请问我可以基于这个结构直接跑荧光光谱吗

作者
Author: sobereva 时间: 2022-3-9 20:05

天天121 发表于 2022-3-9 10:02
sob老师
1.我在网上下载好了cif文件然后扩胞到大概500个原子需要先用cp2k做个结构优化再截取15个分 ...

看下文5.2节怎么弄分子团簇结构
使用Multiwfn做IGMH分析非常清晰直观地展现化学体系中的相互作用
http://sobereva.com/621（http://bbs.keinsci.com/thread-28147-1-1.html）

作者
Author: sobereva 时间: 2022-3-9 20:07

天天121 发表于 2022-3-9 17:23
Sob 老师我已经学会了您用背景电荷法绘制吸收光谱请问您有什么好方法绘制有机晶体的荧光吗
我现在激发 ...

这种方式计算没法体现low层对high层电子结构的极化作用

作者
Author: 喵星大佬 时间: 2022-3-10 02:04
原文里用6-31G*/LanL2DZ做CASPT2也是不咋靠谱的样子呢

作者
Author: 天天121 时间: 2022-3-10 08:19

喵星大佬发表于 2022-3-10 02:04
原文里用6-31G*/LanL2DZ做CASPT2也是不咋靠谱的样子呢

老师你好比方说我现在想扫描一条质子转移的势能曲线也是像上文这种有机晶体请问Low层的原子需要冻结吗？我看过很多做这方面文章的人都说优化过程Low层都是固定不动的，那像势能曲线扫描、过渡态结构优化这种任务也需要冻结Low层吗？

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