ORCA计算DMPO捕获自由基EPR计算报错

阿社社社

各位老师好，我用ORCA模拟了一个含DMPO捕获自由基结构，做了EPR模拟（输入文件如下），计算正常结束。我用Easyspin 6.0版本读取property.txt文件，并加上easyspin官网给的一些参数（我想先试试能不能模拟出来），在Matlab上命令如下，结果报错内存不足。我用笔记本算的，笔记本是24G。请问计算这个是真的需要很大的内存吗，还是说我的参数哪些设置有误？

ORCA输入文件：
! B3LYP D3 def2-TZVP def2/J AUTOAUX
%maxcore  8000
%pal nprocs   16 end
* xyz   1   2
原子坐标
*

%EPRNMR
        GTENSOR   TRUE
        NUCLEI    = ALL H,25 {AISO, ADIP, AORB}
END


easyspin读取结果；

xyz: [59×3 double]
          S: 0.5000
          g: [2.0027 2.0097 2.0184]
     gFrame: [2.9261 0.3086 5.9390]
       Nucs: 'H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H,H'
          A: [30×3 double]
     AFrame: [30×3 double]
       data: [1×1 struct]
    NucsIdx: [1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 … ]


Matlab命令：
clear,clc,clf
Sys = orca2easyspin('ssp-ome-dmpo_epr_property.txt');
Exp.mwFreq = 9.5;   
Exp.Range = [315 350];
[x,y] = pepper(Sys,Exp);
plot (x,y);

TANGYIHAN

只是计算量太大算不出来。你要是想拟合EPR，选一些靠近单电子的磁性核做模拟就行了。距离太远的HFC本身就很小，全纳入没什么意义。