PSI4能量分解计算索烃结构621原子出错

dengyh

计算任务是索烃结构621原子，含W.Cu金属，基组def2-svp，130G内存，50核并行,分两个片段，分别带+1，+4电荷，结构是xtb优化的，还请各位老师的指点，谢谢
输入文件：
memory 130 GB
..........
set {
    scf_type DF
    freeze_core True
    basis def2-svp
}

energy('sapt0')

输出文件出错：
Fatal Error: READ failed. Error description from the OS: Bad address             !
! Error reading the first partial page, unit 97.                                   !
! Error occurred in file: /home/conda/feedstock_root/build_artifacts/psi4_16893433 !
!     50370/work/psi4/src/psi4/libpsio/error.cc on line: 135                       !
! The most recent 5 function calls were:                                           !
! psi::PsiException::PsiException(std::__cxx11::basic_string<char,                 !
!     std::char_traits<char>, std::allocator<char> >, char const*, int)            !
! psi::PSIO::rw(unsigned long, char*, psi::psio_address, unsigned long, int)       !
! psi::DiskDFJK::manage_JK_disk()                                                  !
! psi::DiskDFJK::compute_JK()                                                      !
!                                                                                  !
!----------------------------------------------------------------------------------!

    Psi4 stopped on: Friday, 28 July 2023 09:50PM
    Psi4 wall time for execution: 1:31:55.85

*** Psi4 encountered an error. Buy a developer more coffee!
*** Resources and help at github.com/psi4/psi4.

sobereva

621原子SAPT就别想了

用sobEDA/sobEDAw做能量分解倒是有戏。只要DFT算得动单点的级别就能做得动sobEDA/sobEDAw。看原文，以及里面提到的教程
Tian Lu, Qinxue Chen, A simple, efficient and universal energy decomposition analysis method based on dispersion-corrected density functional theory
ChemRxiv (2023) DOI: 10.26434/chemrxiv-2023-n79rz
文章已基本被JPCA接收了，最近会写博文专题介绍

dengyh

sobereva 发表于 2023-7-29 12:30
621原子SAPT就别想了

用sobEDA/sobEDAw做能量分解倒是有戏。只要DFT算得动单点的级别就能做得动sobEDA/s ...

好的，谢谢Sob老师一直来的指导，我认真拜读一下