计算化学公社

标题: 求助,NAdO分析时内存不够 [打印本页]

作者
Author:
计算消白    时间: 2025-11-25 22:29
标题: 求助,NAdO分析时内存不够
本帖最后由 计算消白 于 2025-11-25 23:26 编辑

用windows工作站先结构优化并计算了单点能,得到了.molden文件,将molden文件导入Multiwfn做NAdO分析(参考http://sobereva.com/719的2.3部分,社长的案例我试过了,在工作站上可以正常运行)。我用的命令窗口,发现问题好像是内存不够(如下图),如何解决这个问题?molden文件太大了没法上传。我上传了计算单点能的输入文件。
(, 下载次数 Times of downloads: 7)
体系:slab是Zn(256个原子),表面吸附一个COOH。
电脑配置:cpu 8375C 64核 内存128G。


作者
Author:
sobereva    时间: 2025-11-26 02:09
找个内存大的机子在Linux下算,或者尝试把晶胞弄小点
还一个可尝试的办法是先用Multiwfn主功能6里的38 Make orbital occupations integer and satisfy Aufbau principle,让轨道占据数都重排成整数(单行列式波函数),然后再做NAdO分析,这样计算AOM的时候就不需要考虑空轨道了,能大幅节约内存

作者
Author:
计算消白    时间: 2025-12-3 17:42
本帖最后由 计算消白 于 2025-12-3 17:45 编辑

谢谢站长!我把体系缩小了,现在可以得到FOM,但在做BOD和NAdO分析时,出现一下错误:
Error: Only restricted and unresticted single determinant wavefunction are supported!
我把 http://sobereva.com/719 的Pd(100)表面吸附苯案例跑完了,没出现这个问题。
我修改了波函数,按照轨道能量由低到高重排占据数,还是不行。
我把吸附构型(.cif)和单点能计算文件上传了。

作者
Author:
sobereva    时间: 2025-12-4 11:05
计算消白 发表于 2025-12-3 17:42
谢谢站长!我把体系缩小了,现在可以得到FOM,但在做BOD和NAdO分析时,出现一下错误:
Error: Only restri ...

Multiwfn里重排轨道占据数后没理由出现此问题。本身此体系的类型和Pd(100)表面吸附苯没有本质区别
至少不通过网盘给出压缩后的molden文件我没法测试
如果不点我帖子下方回复按钮进行回复,99%概率我看不到

作者
Author:
计算消白    时间: 2025-12-4 13:11
sobereva 发表于 2025-12-4 11:05
Multiwfn里重排轨道占据数后没理由出现此问题。本身此体系的类型和Pd(100)表面吸附苯没有本质区别
至少 ...

谢谢站长!我把 压缩后的 .molden 文件用网盘上传了:
通过网盘分享的文件:Zn_COOH_energy_smear-MOS-1_0.7z
链接: https://pan.baidu.com/s/1-Vvf0XTXgCeYcm5XIa__Iw 提取码: p89s
--来自百度网盘超级会员v6的分享

作者
Author:
sobereva    时间: 2025-12-4 15:32
计算消白 发表于 2025-12-4 13:11
谢谢站长!我把 压缩后的 .molden 文件用网盘上传了:
通过网盘分享的文件:Zn_COOH_energy_smear-MOS-1 ...

&SMEAR里加上FIXED_MAGNETIC_MOMENT 1让alpha和beta电子数差值保持为1可解决
作者
Author:
计算消白    时间: 2025-12-4 17:19
本帖最后由 计算消白 于 2025-12-4 17:22 编辑
sobereva 发表于 2025-12-4 15:32
&SMEAR里加上FIXED_MAGNETIC_MOMENT 1让alpha和beta电子数差值保持为1可解决

1. 谢谢社长!在&SMEAR里加上FIXED_MAGNETIC_MOMENT 1后,上述问题解决了。
2. 得到 .molden文件后,我在用FOM做NAdO分析时,按照http://sobereva.com/719 Pd(100)表面吸附苯的步骤:设置片段1为COOH,片段2为和C相连的Zn(编号84),发现Zn表面吸附COOH后,1号NAdO轨道成分主要由Zn和C的s和p轨道贡献,Zn的d轨道贡献很小。感觉结果不太对?Zn的d轨道怎么贡献这么小?

我把吸附构型、1号NAdO轨道、轨道成分的截图放上了。
作者
Author:
sobereva    时间: 2025-12-4 21:55
计算消白 发表于 2025-12-4 17:19
1. 谢谢社长!在&SMEAR里加上FIXED_MAGNETIC_MOMENT 1后,上述问题解决了。
2. 得到 .molden文件后,我 ...

本来Zn就是d壳层满占据,和Pd的成键明显不同




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3