计算化学公社

标题: gaussian计算Pt化学位移求助 [打印本页]

作者
Author: 最瘦二师兄 时间: 2021-9-7 19:32
标题: gaussian计算Pt化学位移求助
各位大神，根据sob老师的帖子http://sobereva.com/354，我用高斯计算Pt化合物的化学位移（DFT方法），选择K2PtCl4作为参考物。可是，我计算的K2PtCl4的Isotropic = -345.6452 ppm，而根据http://chem.ch.huji.ac.il/nmr/techniques/1d/row6/pt.html，K2PtCl4的化学位移在1618 ppm附近，我把我的输入和输出文件贴在这里了，请各位老师帮忙看一下哪里出错了，谢谢大家！

作者
Author: wzkchem5 时间: 2021-9-7 19:48
实验测的是相对化学位移，高斯给出来的是绝对化学位移。实验只有一个数，但高斯必须做两个计算，除了算你的目标物质外还要在相同理论级别下算参比物质，两者相减才能和实验比。
此外核磁不建议用赝势，因为内层电子对化学位移贡献很大，且相减的时候远不如能量抵消得那么好，除非你用的赝势有文献证明确实结果可信。有条件的话尽量做全电子相对论计算

作者
Author: 最瘦二师兄 时间: 2021-9-7 20:43

wzkchem5 发表于 2021-9-7 19:48
实验测的是相对化学位移，高斯给出来的是绝对化学位移。实验只有一个数，但高斯必须做两个计算，除了算你的 ...

谢谢老师！1）其实我是用K2PtCl4做参考物计算其他Pt化合物的化学位移，但根据两者的Isotropic计算出来化学位移只有负一百多，但实验值是负3000多，相差十分大。所以我才单独考虑K2PtCl4的化学位移，想与实验值对比，看看哪出错了。K2PtCl4已经是计算Pt化合物的化学位移常用的参考物了，我尚未找到计算K2PtCl4化学位移的参考物；2）我在关键词里加了int=dkh2考虑全电子相对论效应，可是输出结果显示“No NMR with relativistic Hamiltonian.”

作者
Author: wzkchem5 时间: 2021-9-7 21:11

最瘦二师兄发表于 2021-9-7 13:43
谢谢老师！1）其实我是用K2PtCl4做参考物计算其他Pt化合物的化学位移，但根据两者的Isotropic计算出来化 ...

那说明高斯没有实现考虑显式相对论效应的核磁计算，起码你用的版本没有。不过我不知道最新版有没有。
可以考虑用ORCA或BDF（http://182.92.69.169:7226/Introduction）算，其中BDF因为支持X2C，对相对论效应的处理比高斯、ORCA都准，但需要发邮件申请才能获取

作者
Author: sobereva 时间: 2021-9-8 03:50
Gaussian不支持相对论哈密顿下的NMR计算
对于重原子计算，相对论效应对磁屏蔽值的影响巨大，不考虑结果完全离谱，通常需要做二分量相对论DFT做NMR计算，开源免费的可以用Dirac（http://www.diracprogram.org/doku.php），我记得官方教程里专门有NMR计算部分

作者
Author: kimariyb 时间: 2025-4-3 17:04

sobereva 发表于 2021-9-8 03:50
Gaussian不支持相对论哈密顿下的NMR计算
对于重原子计算，相对论效应对磁屏蔽值的影响巨大，不考虑结果完 ...

卢老师，请问如果是用 Dirac 开二分量相对论 DFT 计算 NMR 以及自旋-自旋耦合常数，那么选择基组的时候还是使用平常推荐的 pcJ-2 这种还是使用专门为相对论准备的的全电子基组，例如 x2c-TZVPPall-2c 这种

作者
Author: sobereva 时间: 2025-4-4 16:37

kimariyb 发表于 2025-4-3 17:04
卢老师，请问如果是用 Dirac 开二分量相对论 DFT 计算 NMR 以及自旋-自旋耦合常数，那么选择基组的时候还 ...

NMR用后者。如果算J，对后者还得再做去部分收缩才能足够描述好内核电子

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