使用Gaussian做镧系金属配合物的量子化学计算

helpme · 发表于 Post on 2021-2-7 16:01:32

非常好的相对论效应体系理论计算的教材。

bottle · 发表于 Post on 2021-3-12 15:43:22

sob老师您好，我想问，在做单点计算和波函数稳定性测试时，小核赝势怎么都不能收敛，是否能用大核赝势来算呢？

sobereva · 发表于 Post on 2021-3-13 06:00:15

bottle 发表于 2021-3-12 15:43
sob老师您好，我想问，在做单点计算和波函数稳定性测试时，小核赝势怎么都不能收敛，是否能用大核赝势来算 ...

用恰当的大核赝势可以

luzujia · 发表于 Post on 2022-5-3 09:24:47

社长，有一个问题想请教。在这个体系中，在确定体系多重度的时候常规做法是在不同多重度下优化再比较能量，您在此用大核赝势好像可以规避这个过程，不确定的情况再对优化好的结构在不同多重度下算单点能量来比较就好了。我的问题是这是否可以作为通用的确定多重度的做法，因为单点比优化便宜得多，这种做法似乎很有优势。

kekexili_08 · 发表于 Post on 2022-5-3 10:02:33

但对于G16来说（至少截止到撰文时最新的G16 C.01而言），收敛不到默认的tight收敛限（相当于SCF=conver=8）的话就不给做涉及能量导数的任务、TDDFT，甚至连我们当前想做的波函数稳定性测试也没法做，这个要求是G16很恶心的一点。

这一点极其恶心。
IOP(5/13=1)也弄得不能用，也是很恶心的事。
有些关键词，或iop，可以提示使用风险，但是直接禁了，使得本来可以算好的都不能往下算了，是不应该的。

sobereva · 发表于 Post on 2022-5-3 11:38:01

kekexili_08 发表于 2022-5-3 10:02
但对于G16来说（至少截止到撰文时最新的G16 C.01而言），收敛不到默认的tight收敛限（相当于SCF=conver=8） ...

有特殊技巧可以用
令Gaussian 16中SCF未收敛到默认收敛限也能继续做后续计算的方法
http://sobereva.com/625（http://bbs.keinsci.com/thread-26427-1-1.html）

sobereva · 发表于 Post on 2022-5-3 11:39:14

luzujia 发表于 2022-5-3 09:24
社长，有一个问题想请教。在这个体系中，在确定体系多重度的时候常规做法是在不同多重度下优化再比较能量， ...

判断基态自旋多重度应当用小核赝势。图省事算单点对比，追求严格分别优化再对比能量

护卫天使 · 发表于 Post on 2022-8-7 23:13:36

老师好，我研究的体系这个博文中的很类似。我的想法是：1、用AIM和IRI定量定性分析稀土和配位原子的相互作用强弱，是否成键。2、用Hirshfeld方法分析成键轨道是哪些。3、用CDA和ETS-NOCV分析不同片段的化学键作用。4、SAPT能量分解。不知道这些分析思路是否可以，是否合适，拜托老师！

sobereva · 发表于 Post on 2022-8-8 04:00:02

护卫天使发表于 2022-8-7 23:13
老师好，我研究的体系这个博文中的很类似。我的想法是：1、用AIM和IRI定量定性分析稀土和配位原子的相互作 ...

SAPT不适合，SAPT主要用于非共价相互作用，配位键已经脱离了这个范畴
其它的可以

护卫天使 · 发表于 Post on 2022-8-18 12:44:06

sobereva 发表于 2022-8-8 04:00
SAPT不适合，SAPT主要用于非共价相互作用，配位键已经脱离了这个范畴
其它的可以

老师好，对于CDA和ETS-NOCV分析，每个配位的水分子都当成一个片段吗？还是只分析配体与稀土的电荷和能量转移？

sobereva · 发表于 Post on 2022-8-18 17:00:23

护卫天使发表于 2022-8-18 12:44
老师好，对于CDA和ETS-NOCV分析，每个配位的水分子都当成一个片段吗？还是只分析配体与稀土的电荷和能量 ...

如果你对不同配体间的相互作用不感兴趣，可以把所有配体一起当个片段用Multiwfn做CDA和ETS-NOCV分析。

护卫天使 · 发表于 Post on 2022-8-30 15:56:28

sobereva 发表于 2022-8-18 17:00
如果你对不同配体间的相互作用不感兴趣，可以把所有配体一起当个片段用Multiwfn做CDA和ETS-NOCV分析。

老师好，像您给的例子，想计算相互结合能。也是所有配体当做一个片段吗？还是每个水分子都当做一个片段？几何优化在pbe1pbe/SDD，计算单点能在pbe1pbe/def2TZVP,计算的泛函及基组是否可行？

sobereva · 发表于 Post on 2022-8-30 18:13:38

护卫天使发表于 2022-8-30 15:56
老师好，像您给的例子，想计算相互结合能。也是所有配体当做一个片段吗？还是每个水分子都当做一个片段？ ...

看你是考察结合一个配体的结合能，还是和所有配体同时的结合能

级别没问题

Quantum198907 · 发表于 Post on 2022-9-20 02:10:33

本帖最后由 Quantum198907 于 2022-9-20 02:43 编辑

sob老师，受益匪浅，学习了

护卫天使 · 发表于 Post on 2022-9-24 18:30:29

老师好，对于12-冠-4-五水合钬(III)和12-冠-4-五水合镧(III)两个配合物。想简单比较一下配合物稳定性，我的想法主要是前面两种:一、比较一下相互结合能（如图中的△E）。二、NBO分析，有很多文献都用NBO分析里对应二阶微扰校正能比较稳定性。三、能否用CDA分析，用转移电荷来比较，越稳定的是不是相应转移电荷越多用于电子共享？四、或者用原子间键径上ELF的最大点的数值，而且还能作图。

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

[Gaussian/gview] 使用Gaussian做镧系金属配合物的量子化学计算

评分 Rate