sobereva 发表于 2025-3-30 19:36
我忘了，如果只有这一处看起来是电子能量，那就只可能是这个

好的，谢谢sob老师

Xx_DerEnTh 发表于 2025-3-30 18:49
sob老师，电子能量是在这里(Total energy)读对吗？

我忘了，如果只有这一处看起来是电子能量，那就只可能是这个

sobereva 发表于 2025-3-29 23:09
内核电子能量很低，数量级自然很大

sob老师，电子能量是在这里(Total energy)读对吗？

Xx_DerEnTh 发表于 2025-3-29 21:18
sob老师您好！想再请教一下Dirac输出的电子能量具体是在哪里读？是下面这个吗？
Sum of all contributio ...

内核电子能量很低，数量级自然很大

sobereva 发表于 2025-3-28 16:30
是

不需要特意考虑不代表不能考虑，考虑了当然从精度来说有益，代价就是多花时间而已

sob老师您好！想再请教一下Dirac输出的电子能量具体是在哪里读？是下面这个吗？
Sum of all contributions to the energy
   Total energy                             :    -31332.519464472796
感觉数值好大

sobereva 发表于 2025-3-28 16:30
是

不需要特意考虑不代表不能考虑，考虑了当然从精度来说有益，代价就是多花时间而已

好的，谢谢sob老师！

Xx_DerEnTh 发表于 2025-3-28 15:49
sob老师您好！对于这篇博文算TlCl键解离能的例子，我如果用Dirac的X2C二分量哈密顿考虑耦合，单独算Cl的时 ...

是

不需要特意考虑不代表不能考虑，考虑了当然从精度来说有益，代价就是多花时间而已

sobereva 发表于 2025-3-14 03:51
不能用。对于基态DFT计算，此时需要做严格的二分量相对论计算体现旋轨耦合效应

sob老师您好！对于这篇博文算TlCl键解离能的例子，我如果用Dirac的X2C二分量哈密顿考虑耦合，单独算Cl的时候为了保证能量可相加减，是不是同样也得用X2C？如果是这样的话，对像Cl原子这样通常不需要考虑相对论的原子，用X2C这样处理相对论的方法会不会造成什么负面影响？

sobereva 发表于 2025-3-14 03:51
不能用。对于基态DFT计算，此时需要做严格的二分量相对论计算体现旋轨耦合效应

好的，谢谢sob老师！

Xx_DerEnTh 发表于 2025-3-13 09:54
感谢sob老师回复，请问这种没给旋轨势的情况是不是就不能用NWChem分析了？

不能用。对于基态DFT计算，此时需要做严格的二分量相对论计算体现旋轨耦合效应

sobereva 发表于 2025-3-13 04:33
这是赝势而非旋轨势

感谢sob老师回复，请问这种没给旋轨势的情况是不是就不能用NWChem分析了？

Xx_DerEnTh 发表于 2025-3-12 22:11
sob老师您好！我想用NWChem考虑锕系元素旋轨耦合，请教一下Stuttgart赝势的官网上蓝色的这一部分（图1）是 ...

这是赝势而非旋轨势

sob老师您好！我想用NWChem考虑锕系元素旋轨耦合，请教一下Stuttgart赝势的官网上蓝色的这一部分（图1）是旋轨势吗？手动加旋轨势是不是加这一部分？

冰释之川 发表于 2017-4-18 13:01
问题解决了，是Linear dependencies的锅：

Unfortunately, the spin-orbit code does not correctly h ...

nwchem判断基组线性依赖的默认阈值是1.0e-5。由于sodft不支持消除冗余积函数的算法，这么高的阈值显然不适用。可以用手册里没有讲的一个命令降低阈值。根据出错信息里的“Smallest S eigenvalue”，假如此值在1.0e-8以上，可以用

set lindep:tol 1.0e-8

但是如果出现更低的特征值，降低阈值的做法就不保险了，可能会带来其他问题，比如scf不收敛。