计算化学公社

标题: 对单层GaN做振动分析out文件出现很多且中途自动停止算 [打印本页]

作者
Author: 风里禾 时间: 2026-3-12 17:41
标题: 对单层GaN做振动分析out文件出现很多且中途自动停止算
用gaussview做的GaN.gjf
优化：设置约束：Z方向保持不变：使用较严的阈值；选择基组：pob-DZVP-rev2；开启OT；产生输入文件GaN.inp
振动：选择基组：pob-DZVP-rev2；开启OT；产生输入文件GaN-1.inp
如果把RESTART.wfn里的波函数作为初猜的话
GaN-1.inp里面的
&DFT
BASIS_SET_FILE_NAME  BASIS_pob
POTENTIAL_FILE_NAME  POTENTIAL
# WFN_RESTART_FILE_NAME GaN-RESTART.wfn
CHARGE 0 #Net charge
MULTIPLICITY 1 #Spin multiplicity
还有
&SCF
   MAX_SCF 25 #Maximum number of steps of inner SCF
   EPS_SCF 1.0E-07 #Convergence threshold of density matrix of inner SCF
# SCF_GUESS RESTART #Use wavefunction from WFN_RESTART_FILE_NAME file as initial guess
# IGNORE_CONVERGENCE_FAILURE #Continue calculation even if SCF not converged, works for version >= 2024.1

这三个#都删了或者都不删或者保留一个，都会输出一样大小的64个48KB的out文件，并且会同时停止计算，想问一下该怎么解决呢？
也尝试过从头算，不要初猜，依旧输出64个48KB的out文件，为什么会这样呢
  谢谢！

作者
Author: sobereva 时间: 2026-3-14 08:15
一律自己看out文件内容或者提供out文件，光看out文件大小而不看内容，完全是迷惑行为
如果是按下文用的Multiwfn创建的CP2K输入文件（貌似是），必须帖子里明确说清楚，否则都不知道你在操作什么、怎么产生的。并且输入文件直接上传，别就光给出某个信息片段，否则有硬伤别人也发现不了、没法测试
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

作者
Author: northern_mar 时间: 2026-3-14 12:09
你把输入文件上传一下呢不然这一个小片段大家也没法帮你看的。我之前每个副本设置的核数太少也出过问题，仅供参考。

作者
Author: 杨小狗 时间: 2026-3-16 15:52

sobereva 发表于 2026-3-14 08:15
一律自己看out文件内容或者提供out文件，光看out文件大小而不看内容，完全是迷惑行为
如果是按下文用的Mul ...

(, 下载次数 Times of downloads: 1)
感谢回复，发帖人刚注册的，可发文件内存不太够，所以帮忙发一下
后面还试着把OT关键词由IRAC改成STRICT，也是很多out文件同时停止计算了

作者
Author: 杨小狗 时间: 2026-3-16 15:56

northern_mar 发表于 2026-3-14 12:09
你把输入文件上传一下呢不然这一个小片段大家也没法帮你看的。我之前每个副本设置的核数太少也出过问题， ...

(, 下载次数 Times of downloads: 2)
(, 下载次数 Times of downloads: 2)
感谢回复
GaN.inp是第一部优化的输入文件，GaN-1.inp是第二部振动分析的输入文件
我用的64核计算的

作者
Author: northern_mar 时间: 2026-3-16 16:06

杨小狗发表于 2026-3-16 15:56
感谢回复
GaN.inp是第一部优化的输入文件，GaN-1.inp是第二部振动分析的输入文件
我用的64 ...

NPROC_REP太小了，可以设为16或32试试

作者
Author: LittlePupil 时间: 2026-3-16 23:44
可能是内存问题，此时换个内存大的节点或者减少同时跑的副本数（增加“NPROC_REP”）应该可以解决。

作者
Author: 风里禾 时间: 2026-3-18 19:53

northern_mar 发表于 2026-3-16 16:06
NPROC_REP太小了，可以设为16或32试试

改成16跑成了，谢谢大佬们

欢迎光临计算化学公社 (http://bbs.keinsci.com/)