计算化学公社

标题: 如何进行ESR顺磁谱图的模拟？ [打印本页]

作者
Author: nnzhang 时间: 2017-10-17 14:49
标题: 如何进行ESR顺磁谱图的模拟？
找了一些这方面的文献，但是基本没有详细说明用的什么软件或者是什么基组什么的，有没有哪位大神做过这方面的模拟，提供一下资源或方法，不胜感激？

作者
Author: kyuu 时间: 2017-10-17 15:05
ADF 有钱你就上

作者
Author: liyuanhe211 时间: 2017-10-17 19:36
用高斯：基组使用EPR-N系列能够较好描述核处电子密度的基组，一般的泛函。用高斯计算各核的何处自旋密度（什么都不用写就会输出），计算出后乘以对应原子核的特定常数即可

作者
Author: nnzhang 时间: 2017-10-17 21:34

liyuanhe211 发表于 2017-10-17 19:36
用高斯：基组使用EPR-N系列能够较好描述核处电子密度的基组，一般的泛函。用高斯计算各核的何处自旋密度（ ...

“对应原子核的特定常数”是指的哪个常数啊？

作者
Author: sobereva 时间: 2017-10-17 22:14
计算超精细耦合在exploring第三版里就有例子
如果是要计算ESR谱图，用easyspin结合ORCA

作者
Author: sobereva 时间: 2017-10-17 22:32

nnzhang 发表于 2017-10-17 21:34
“对应原子核的特定常数”是指的哪个常数啊？

参看例如
J. Phys. Chem. A 1997, 101, 3174-3181

作者
Author: nnzhang 时间: 2017-10-17 22:48

sobereva 发表于 2017-10-17 22:14
计算超精细耦合在exploring第三版里就有例子
如果是要计算ESR谱图，用easyspin结合ORCA

嗯嗯谢谢老师！

作者
Author: nnzhang 时间: 2017-10-17 22:49

sobereva 发表于 2017-10-17 22:32
参看例如
J. Phys. Chem. A 1997, 101, 3174-3181

好的老师

作者
Author: dana19851128 时间: 2018-4-1 17:42
我是做电子顺磁计算的有什么问题？

作者
Author: nnzhang 时间: 2018-5-10 11:21

dana19851128 发表于 2018-4-1 17:42
我是做电子顺磁计算的有什么问题？

在用matlab结合easyspin进行谱图模拟的时候，原子数大于20个后就计算很慢很慢，这个是什么原因呢？有没有什么解决办法呢？

作者
Author: mengtian1234 时间: 2020-7-1 17:57

dana19851128 发表于 2018-4-1 17:42
我是做电子顺磁计算的有什么问题？

请问我要怎么做电子顺磁计算，今天刚装了一个ORCA，我要使用什么基组和关键词

作者
Author: fond 时间: 2021-1-4 11:26

dana19851128 发表于 2018-4-1 17:42
我是做电子顺磁计算的有什么问题？

通过实验电子自旋共振可以定性测出实验中产生自由基的类别，比如有DMPO-OH，DMPO-SO4.-。想问一下可以模拟计算DMPO-OH,DMPO-SO4.-的EPR谱图吗？这样的话实验和计算可以吻合，同时证明一下，确实是羟基自由基和硫酸根自由基。可以计算吗？如果可以的话，请问是怎么做呢？

作者
Author: wzkchem5 时间: 2021-1-4 17:25

fond 发表于 2021-1-4 11:26
通过实验电子自旋共振可以定性测出实验中产生自由基的类别，比如有DMPO-OH，DMPO-SO4.-。想问一下可以模 ...

可以，但没必要。这种自由基，尤其是DMPO-OH肯定有文献的实验谱图，跟实验比是更加靠谱的，没必要做计算。如果捕捉了一个比较奇怪的自由基，比如什么金属自由基、簇合物自由基之类的，那还值得算一下。

作者
Author: 自由风 时间: 2023-5-18 15:47

nnzhang 发表于 2018-5-10 11:21
在用matlab结合easyspin进行谱图模拟的时候，原子数大于20个后就计算很慢很慢，这个是什么原因呢？有没有 ...

是很慢，慢到我以为是软件卡死了。你有解决方法吗？

作者
Author: wzkchem5 时间: 2023-5-18 18:31

自由风发表于 2023-5-18 08:47
是很慢，慢到我以为是软件卡死了。你有解决方法吗？

手动删掉那些很小的、你不用计算就知道不会对谱图有显著影响的A值，再算。easyspin计算慢一般都是因为考虑了太多核的A值，实际上很多核的A值很小或者可以用展宽（lwpp）等效代替，没必要在计算中显式考虑。

作者
Author: 自由风 时间: 2023-5-19 15:35

wzkchem5 发表于 2023-5-18 18:31
手动删掉那些很小的、你不用计算就知道不会对谱图有显著影响的A值，再算。easyspin计算慢一般都是因为考 ...

我算的A值大的有130，小的有15，最小的只有个位数。那手动删除的话15左右的算小吗？大小是指的绝对值吗？

作者
Author: wzkchem5 时间: 2023-5-19 16:17

自由风发表于 2023-5-19 08:35
我算的A值大的有130，小的有15，最小的只有个位数。那手动删除的话15左右的算小吗？大小是指的绝对值吗？

结合你所用的lwpp，判断这个较小的A值有或者没有会不会显著影响峰形。最简单的方法是，另模拟一个只有一个A值的信号，其A值等于比如说15，lwpp和你当前模拟相同。然后看看去掉这个A值会不会显著影响峰形

作者
Author: wangzhe 时间: 2023-5-19 16:17
之前用ORCA计算的时候用BP/EPR-II计算的结果和实验值最为吻合，不过应该是case by case，仅供参考。

作者
Author: QuantumicGuy 时间: 2023-6-13 10:24
考虑的原子数太多Matlab也会崩掉，我的建议是先把自旋密度图画出来，然后根据自旋密度分布大的位置，去选择计算那些磁性核。总的来说条件就是：磁性核、该磁性核周围自旋密度大

作者
Author: QuantumicGuy 时间: 2023-7-6 16:23

wangzhe 发表于 2023-5-19 16:17
之前用ORCA计算的时候用BP/EPR-II计算的结果和实验值最为吻合，不过应该是case by case，仅供参考。

是BP86吗？

作者
Author: wangzhe 时间: 2023-7-6 16:49

QuantumicGuy 发表于 2023-7-6 17:23
是BP86吗？

BP86和BP计算的D和E张量结果是一样的。不知道这两个是否是同一个泛函。

作者
Author: wzkchem5 时间: 2023-7-6 17:19

wangzhe 发表于 2023-7-6 09:49
BP86和BP计算的D和E张量结果是一样的。不知道这两个是否是同一个泛函。

是一个泛函

作者
Author: wangzhe 时间: 2023-7-6 17:27

wzkchem5 发表于 2023-7-6 18:19
是一个泛函

感谢

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