计算化学公社

标题: gaussian优化871个原子的MOF体系总是出现问题,请老师们指点一下 [打印本页]

作者
Author: jiahui 时间: 2023-10-23 12:39
标题: gaussian优化871个原子的MOF体系总是出现问题,请老师们指点一下
在优化MOF材料时总是出现如下的问题，我的体系共871个原子，结果也没有报错，直接就停止计算了，请问老师们我应该怎么解决这个问题？输入文件如下：
%chk=four.chk
%mem=15600mb
%nprocs=36
# PBE1PBE/genecp opt SCF=NoVarAcc

Title Card Required

0 2
坐标

C H O  0
6-31G*
****
Zr  0
LANL2DZ
****

Zr  0
LANL2DZ

输出文件最后如下：（有附件）
There are 10356 symmetry adapted cartesian basis functions of A symmetry.
There are  9686 symmetry adapted basis functions of A symmetry.
  9686 basis functions, 19492 primitive gaussians, 10356 cartesian basis functions
  2401 alpha electrons    2400 beta electrons
   nuclear repulsion energy 467281.9936586973 Hartrees.
NAtoms=  871 NActive=  871 NUniq=  871 SFac= 1.00D+00 NAtFMM= 60 NAOKFM=T Big=T
Integral buffers will be 131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
One-electron integrals computed using PRISM.
PrsmSu:  requested number of processors reduced to:  14 ShMem 1 Linda.
PrsmSu:  requested number of processors reduced to:  28 ShMem 1 Linda.
1 Symmetry operations used in ECPInt.
ECPInt:  NShTT= 4900015 NPrTT=33227508 LenC2=  769695 LenP2D= 2069548.
LDataN:  DoStor=T MaxTD1= 5 Len=  102
Number of processors reduced to 5 by ecpmxn.

作者
Author: snljty2 时间: 2023-10-23 13:07
算不动的，不用试了。

作者
Author: jiahui 时间: 2023-10-23 13:41

snljty2 发表于 2023-10-23 13:07
算不动的，不用试了。

那我应该用哪个软件优化呢？或者把分子拆开分开计算行吗？

作者
Author: zjxitcc 时间: 2023-10-23 13:43
跟这个帖子类似http://bbs.keinsci.com/thread-40548-1-1.html
周期性体系强烈建议使用擅长周期性计算的程序（例如QE, VASP, CP2K等），人生短暂，时间宝贵。

作者
Author: muyueyan 时间: 2023-10-23 17:01
确实不好搞，Gaussian算了两个月也没算出来我的CuO

作者
Author: muyueyan 时间: 2023-10-23 17:02
就是不收敛

作者
Author: xsc6 时间: 2023-10-24 01:36
871个辛苦你的cpu了

作者
Author: sobereva 时间: 2023-10-24 04:39
Gaussian根本没戏，以目前的主流计算条件，400原子以上就根本别考虑了。

免费的CP2K是优化这个首选，没有之一。算这个是小儿科
用像样的双路服务器，纯泛函优化这个体系一个小时内准保得到想要的结果。创建输入文件还极其容易：
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）
北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里还有极为全面的使用讲解和各个方面的海量例子

一般计算条件下，CP2K用纯泛函优化两、三千原子的晶胞都是可行的。

作者
Author: cccxq 时间: 2024-10-17 21:08

sobereva 发表于 2023-10-24 04:39
Gaussian根本没戏，以目前的主流计算条件，400原子以上就根本别考虑了。

免费的CP2K是优化这个首选，没 ...

老师，请问假如我想用Orca算Mof的SOC，那我在高斯中应该选择不对称单元还是以金属为中心，生长配位的配体去进行进一步计算呢

作者
Author: sobereva 时间: 2024-10-20 10:11

cccxq 发表于 2024-10-17 21:08
老师，请问假如我想用Orca算Mof的SOC，那我在高斯中应该选择不对称单元还是以金属为中心，生长配位的配体 ...

如果你是指旋轨耦合的TDDFT，还不如直接用CP2K做，CP2K直接就支持周期性体系的SOC-TDDFT

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