计算化学公社
标题: 电子顺磁共振谱（EPR）的计算问题 [打印本页]
作者Author: Jerryluo 时间: 2018-11-1 22:15
标题: 电子顺磁共振谱（EPR）的计算问题
大家好，请问金属配合物团簇的电子顺磁共振谱（EPR）应该用什么程序计算，用Gaussian能算吗，有相关的教程吗。谢谢！
作者Author: sobereva 时间: 2018-11-2 00:24
easyspin可以
原文：Journal of Magnetic Resonance 178 (2006) 42–55
一个应用例子Phys. Status Solidi B 2018, 1800026
可以结合ORCA使用
Gaussian没有现成的关键词
作者Author: dana19851128 时间: 2019-8-3 16:23
easyspin 只是画图软件，不建议用gaussian做顺磁谱的计算，ORCA软件可以，我是专业做顺磁谱图计算的，多看FrankNeese的工作。
作者Author: biogon 时间: 2020-7-11 12:48
你用过vms-tool画epr谱图吗
作者Author: 胡说 时间: 2020-8-17 17:27
本帖最后由 胡说 于 2020-8-17 17:35 编辑
老师好，请教个EPR计算的问题，我看基础班讲义中有涉及到高斯算g tensor和超精细耦合常数的内容，但g16 manul中却没有找到它们具体的表达式。ORCA手册中有给出了这两部分详细的定义，超精细耦合常数公式中还包含有SOC的贡献（除了fermi contact和spin-dipole）。高斯输出的g tensor也有四部分贡献，这对应ORCA定义中的四项，两者应当一致。但高斯算的超精细耦合常数不包含SOC贡献，所以更建议用ORCA来算，因为它可以包含了所有的贡献。不知道这样理解是否对？
作者Author: sobereva 时间: 2020-8-18 09:13
可能是这样，我没有仔细研究过ORCA这部分的信息
作者Author: 胡说 时间: 2020-8-18 12:33
好的，多谢老师
作者Author: abdoman 时间: 2021-6-8 08:44
您好，想问一下，orca计算EPR可以提供什么信息呢？
1.实验可以测得g值，那与理论计算预测的g值做比较吗？
2. 还能够提供哪些信息呢？
作者Author: wzkchem5 时间: 2021-6-8 15:09
g值，A值，D值，E值；各向同性，各向异性的，都可以算。算出来可以用easyspin画图
作者Author: LJF 时间: 2021-6-26 09:28
ORCA可以，团簇没有尝试过，分子配合物算得很快的。
作者Author: dana19851128 时间: 2021-7-15 15:20
理论计算G值的误差较大，目前主要是比较趋势
作者Author: biogon 时间: 2021-7-15 15:59
DFT算这种不过图个乐
作者Author: dana19851128 时间: 2021-7-15 18:25
DFT可以算的 见Frank Neese的工作
作者Author: biogon 时间: 2021-7-16 22:36
这我当然知道，但是准确度实在是一般，orca5新出的NEVPT2算EPR参数就准确了不少
作者Author: wzkchem5 时间: 2021-7-16 22:54
构象平均也很重要。如果体系只有一个优势构象，DFT就算不准，误差起码应该还比较系统。如果有很多个优势构象，但是只挑一个算，哪怕NEVPT2也无力回天
作者Author: biogon 时间: 2021-7-17 10:34
VPT2算振动平均的EPR参数这个应该也能提高不少精度吧
作者Author: wzkchem5 时间: 2021-7-17 14:55
嗯对，ORCA 5发布会讲过这个
作者Author: Scarlett-ww 时间: 2022-1-17 17:18
请问D值和g值不一样的地方只是比g值多了一个spin-flip项吗？另外，E值是什么呢？实验获取的数据只有g值和A值，ADF也只能计算出g值和A值，这个D值和E值是为了在理论上更准确描述吗？
作者Author: wzkchem5 时间: 2022-1-17 18:01
D值和E值是零场分裂参数，只对S>1/2的体系存在，如果你的体系是S=1/2的，那么D和E都是0；如果你的体系是S>1/2的，那么即便实验上也应该把D值和E值拟合出来，不然属于数据没报全。参见https://en.wikipedia.org/wiki/Zero_field_splitting
作者Author: Scarlett-ww 时间: 2022-1-17 22:30
好的，谢谢，我学习一下
作者Author: nianbin 时间: 2023-5-20 14:48
请问下我20几个原子的体系 matlab算了十几天还没有结束 这个正常吗
作者Author: wzkchem5 时间: 2023-5-20 15:47
手动去掉所有的A值里面较小的、确定不会影响谱图的那些，再试
作者Author: nianbin 时间: 2023-5-20 16:10
您的意思是把这些去掉是吧？
作者Author: wzkchem5 时间: 2023-5-20 16:36
嗯，从这里去掉也可以，把数据读进matlab再去掉也可以。但是只去掉较小的A值，例如那个-1.77MHz的A值应该是可以去掉的
作者Author: nianbin 时间: 2023-5-20 20:53
请问下 是就删掉我划线那部分还是说下面的A(iso)连同后面的数值也要删掉
作者Author: wzkchem5 时间: 2023-5-21 05:33
我没删过property文件，但是总的宗旨就是，把你自己代入easyspin的视角，假如删完以后剩余的部分的格式和原文件格式相同，easyspin看这个文件觉得是orca少算了一些A值，看不出是你手动改的，就行
作者Author: sigma-jlu 时间: 2023-8-3 11:27
上周看了几篇Neese组的文献，而且Pantazis曾分享有一个视频，讲orca算epr很透彻，建议看看。
提取码：p55h
总体感觉：
1、只要不是固态体系，首推orca，固态体系目前可以尝试cp2k；
2、easyspin计算跟实验可比拟的谱图之前，需要计算g、A、D张量。
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |