|
发布时间: 2025-11-6 18:16
正文摘要:老师们好，我将同一个体系，用Gaussian软件进行了GKS结合int=DKHSO考虑旋轨耦合算了单点能。以及用ORCA软件用sf-X2C/SARC-DKH-TZVP/X2C-TZVPALL进行了标量相对论计算。输入文件已经上传，有一些问题想请教老师！ 1. ... |
zjxitcc 发表于 2025-11-10 13:06 感谢老师的指点！ |
|
本帖最后由 zjxitcc 于 2025-11-10 13:08 编辑 1. GKS计算目前难以进行键级、静电势、自由能等计算，如果非要这些信息，建议退一步采用普通KS-DFT结合标量相对论哈密顿做相应的计算和分析。 Gaussian目前不支持GKS波函数稳定性检验，不过免费、开源的PySCF程序支持，因此我们可以使用PySCF检测波函数稳定性。注意，没有必要在PySCF里重做一遍GKS计算，那样不仅费时费力，是否能收敛到目标SCF解 是一个不简单的问题。我们可以直接基于Gaussian计算结果，利用MOKIT把GKS波函数传给PySCF，PySCF会立即收敛，马上进入GKS波函数稳定性检验步骤，例子见 利用MOKIT从Gaussian向其他量化程序传轨道 内容比较长，与你问题相关的文字内容从“该功能还支持GHF...”开始。使用该功能建议在gjf文件里写上nosymm int=nobasistransform再做GKS计算。 2. 是这样写。但ORCA目前没有二分量或GKS方法。检验一下常规KS-DFT计算的波函数稳定性是可以的，标量相对论应该也可以检验。是的。 |
本帖最后由 wangzh 于 2025-11-10 09:45 编辑 sobereva 发表于 2025-11-9 02:53 感谢老师您的解答，您说提示波函数稳定也未必就真的稳定，这只是必要非充分条件，那我如何判定波函数是否真的稳定？如果不稳定我应该怎么让ORCA程序找到稳定的波函数？感谢老师指点！ |
|
1 用标量相对论计算产生的波函数做波函数分析就够了。本来Multiwfn也不支持二分量波函数 2 是 通常是。但提示波函数稳定也未必就真的稳定，这只是必要非充分条件 4 光写int=DKHSO是DKH4标量哈密顿，同时做GKS时对应DKH4二分量哈密顿 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-25 22:08 , Processed in 0.207289 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites