计算化学公社

标题: PSI4能量分解计算索烃结构621原子出错 [打印本页]

作者
Author: dengyh 时间: 2023-7-28 23:31
标题: PSI4能量分解计算索烃结构621原子出错
计算任务是索烃结构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.

作者
Author: sobereva 时间: 2023-7-29 12:30
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接收了，最近会写博文专题介绍

作者
Author: dengyh 时间: 2023-7-29 14:10

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

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

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

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