zjxitcc 发表于 2025-8-19 20:25 非常感谢！ |
lycheeho 发表于 2025-8-5 09:37 “输出结果为S2=0.318475” 最好不要简写为S2，而采用更严谨的<S^2>或<S**2>符号，以免造成歧义。 “但S2偏离1较多，说明不是自旋纯态” 注意，<S^2>偏离0就说明不是自旋纯态了，不是偏离1.0的问题。若<S^2>接近1.0说明体系可能有较强的双自由基特征（但不一定，还需结合UNO占据数判断）。 “UM062X/DEF2TZVP guess=mix pop=NO输出UNO看轨道占据数” 这不是一个合理的计算。对于电子结构较为复杂的体系，不可以每次抛弃之前的计算结果，重新独立计算SCF。因为那样SCF需要从程序默认初猜开始迭代，到收敛的时候 不仅浪费了很多机时，也不方便确定收敛到了什么结果；万一收敛到了错误的SCF解，一切白搭。合理的做法有两种：（1）读取你之前um062x/def2tzvp stable=opt guess=mix em=gd3辛苦算出来的UM062X稳定波函数，关键词例如UM062X chkbasis nosymm guess=read geom=allcheck pop=NO；（2）不另做Gaussian计算，直接使用UM062X的fch文件，利用Multiwfn或MOKIT直接产生UNO轨道及占据数。 |
|
本帖最后由 lycheeho 于 2025-8-5 11:23 编辑 大家好，我首先对闭壳层用m062x/def2tzvp stable em=gd3结果输出为The wavefunction has an RHF -> UHF instability. 之后尝试计算开壳层单重态，um062x/def2tzvp stable=opt guess=mix em=gd3，输出结果为S2=0.318475，且能量比闭壳层低13 meV，但S2偏离1较多，说明不是自旋纯态。接着我想计算双自由基特征值，UM062X/DEF2TZVP guess=mix pop=NO输出UNO看轨道占据数，根据如图方法计算出y0=0.016，太小了，这些结果能说明体系是开壳层双自由基单重态基态吗？谢谢大家！ |
UNO计算y0.png (22.85 KB, 下载次数 Times of downloads: 55)
influenza 发表于 2025-3-13 20:41 有对称破缺单重态、自旋极化单重态、开壳层单重态，都是一个意思，但没有“自旋破缺态”。 “发现S2算符是0，自旋密度计算验证了确实是破缺态”说明您看错了数据。可以上传输出文件，介绍一下计算步骤和计算目的，大家可以帮你检查。 |
| 请问各位，计算了一个环辛三烯结构的对称破缺态，发现S2算符是0，自旋密度计算验证了确实是破缺态。这正常吗？同样是自旋破缺态，前面有个中间体是类似于环丁烷单键打开的那种双自由基，破缺态的S2算符就接近1。很奇怪。 |
zjxitcc 发表于 2024-5-12 12:55 好的，谢谢 |
yangjunfang 发表于 2024-5-12 09:50 表示自旋平方算符S^2在当前波函数下的期望值<Psi|S^2|Psi>。单重态的<S^2>理想值为0，三重态理想值为2，偏离理想值表示体系有自旋污染，波函数不是自旋纯态。如果不理解这句话，需要阅读Levin的Quantum Chemsitry一书，或者Szabo的Quantum Chmeistry一书前4章，了解基础知识。 |
zjxitcc 发表于 2023-2-28 12:49 请问一下，<S**2>的数值代表什么啊 |
wzkchem5 发表于 2023-3-1 18:26 老师您好，您说的方法试过了，还是不行，依然会报错 Read checkpoint file 3.chk Write formatted file 3.fchk GetFrg NYI. Error termination via Lnk1e at Wed Mar 1 18:42:09 2023. Segmentation fault 输入文件 %chk=/home-ys/users/nsys_DSY/zxy/3/3.chk %nproc=24 %mem=1GB #6-31g* em=gd3 m062x guess=read nosymm opt freq Title Card Required 0 1 C -1.28366186 4.60411310 3.10975128 H -0.88662223 5.33330949 3.78468250 H -1.33430593 3.57606799 3.40208401 C -1.78221248 5.03423181 1.71761437 H -2.82052506 5.28671752 1.77291320 C -0.98055090 6.25983159 1.24133305 H 0.05776167 6.00734588 1.18603422 C -1.47910153 6.68995030 -0.15080386 H -1.87614116 5.96075391 -0.82573508 H -1.42845746 7.71799541 -0.44313660 H -1.11546929 7.06502168 1.93297309 H -1.64729409 4.22904172 1.02597432 |
zhibisuyusheng 发表于 2023-3-1 11:04 我感觉可能是formchk不支持片段初猜计算。试试读取这个波函数，做一个不分片段的unrestricted计算（应该一步就收敛了），再做formchk |
wzkchem5 发表于 2023-3-1 18:02 这是在超算上的报错信息，只产生了0b的fchk Read checkpoint file 3.chk Write formatted file 3.fchk GetFrg NYI. Error termination via Lnk1e at Wed Mar 1 17:46:00 2023. Segmentation fault |
zhibisuyusheng 发表于 2023-3-1 10:53 但凡问报错问题，必须提供详细报错输出信息 |
wzkchem5 发表于 2023-2-28 22:44 老师您好，感谢您的解答，我明白您说的意思了，之后我通过sob老师片段组合波函数计算自旋极化单重态的方法练习了一下双自由基的计算，我分了三步 第一步：#6-31g* em=gd3 m062x guess(fragment=2) 第二步：#6-31g* em=gd3 m062x guess=read nosymm stable=opt 第三步: #6-31g* em=gd3 m062x guess=read nosymm opt freq 我确实得到了对应的自旋极化单重态， S**2 before annihilation 1.0085, after 0.0823 但在我想要绘制自旋密度等值面图时发现执行formchk xxx.chk时总是报错，这是我计算过程有问题才导致无法转换吗，期待您的回复 |
zhibisuyusheng 发表于 2023-2-28 15:37 因为你没有指定unrestricted计算，单重态体系默认做restricted计算，<S**2>必然是0，因此没有显示的必要。如果怀疑是单重态双自由基，必须显式指定做unrestricted计算，并用之前几楼说的几种方法准备自旋对称性破缺初猜 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-24 12:03 , Processed in 0.173859 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites
