计算化学公社

标题: ORCA计算DMPO捕获自由基EPR计算报错 [打印本页]

作者
Author: 阿社社社 时间: 2023-5-10 15:06
标题: ORCA计算DMPO捕获自由基EPR计算报错
本帖最后由阿社社社于 2023-5-10 15:11 编辑

各位老师好，我用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);

作者
Author: TANGYIHAN 时间: 2023-5-31 13:44
只是计算量太大算不出来。你要是想拟合EPR，选一些靠近单电子的磁性核做模拟就行了。距离太远的HFC本身就很小，全纳入没什么意义。

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