计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: Jerryluo
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 电子顺磁共振谱(EPR)的计算问题

[复制链接 Copy URL]

1245

帖子

1

威望

3860

eV
积分
5125

Level 6 (一方通行)

16#
发表于 Post on 2021-7-17 10:34:43 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-16 22:54
构象平均也很重要。如果体系只有一个优势构象,DFT就算不准,误差起码应该还比较系统。如果有很多个优势 ...

VPT2算振动平均的EPR参数这个应该也能提高不少精度吧

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

17#
发表于 Post on 2021-7-17 14:55:14 | 只看该作者 Only view this author
biogon 发表于 2021-7-17 03:34
VPT2算振动平均的EPR参数这个应该也能提高不少精度吧

嗯对,ORCA 5发布会讲过这个
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

82

帖子

0

威望

723

eV
积分
805

Level 4 (黑子)

18#
发表于 Post on 2022-1-17 17:18:31 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-8 15:09
g值,A值,D值,E值;各向同性,各向异性的,都可以算。算出来可以用easyspin画图

请问D值和g值不一样的地方只是比g值多了一个spin-flip项吗?另外,E值是什么呢?实验获取的数据只有g值和A值,ADF也只能计算出g值和A值,这个D值和E值是为了在理论上更准确描述吗?

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

19#
发表于 Post on 2022-1-17 18:01:25 | 只看该作者 Only view this author
Scarlett-ww 发表于 2022-1-17 10:18
请问D值和g值不一样的地方只是比g值多了一个spin-flip项吗?另外,E值是什么呢?实验获取的数据只有g值和 ...

D值和E值是零场分裂参数,只对S>1/2的体系存在,如果你的体系是S=1/2的,那么D和E都是0;如果你的体系是S>1/2的,那么即便实验上也应该把D值和E值拟合出来,不然属于数据没报全。参见https://en.wikipedia.org/wiki/Zero_field_splitting
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

82

帖子

0

威望

723

eV
积分
805

Level 4 (黑子)

20#
发表于 Post on 2022-1-17 22:30:29 | 只看该作者 Only view this author
wzkchem5 发表于 2022-1-17 18:01
D值和E值是零场分裂参数,只对S>1/2的体系存在,如果你的体系是S=1/2的,那么D和E都是0;如果你的体系是S ...

好的,谢谢,我学习一下

182

帖子

0

威望

2233

eV
积分
2415

Level 5 (御坂)

21#
发表于 Post on 2023-5-20 14:48:32 | 只看该作者 Only view this author
请问下我20几个原子的体系 matlab算了十几天还没有结束 这个正常吗

202305201447136104..png (39.51 KB, 下载次数 Times of downloads: 45)

202305201447136104..png

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

22#
发表于 Post on 2023-5-20 15:47:56 | 只看该作者 Only view this author
nianbin 发表于 2023-5-20 07:48
请问下我20几个原子的体系 matlab算了十几天还没有结束 这个正常吗

手动去掉所有的A值里面较小的、确定不会影响谱图的那些,再试
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

182

帖子

0

威望

2233

eV
积分
2415

Level 5 (御坂)

23#
发表于 Post on 2023-5-20 16:10:47 | 只看该作者 Only view this author
您的意思是把这些去掉是吧?

202305201610198960..png (123.36 KB, 下载次数 Times of downloads: 39)

202305201610198960..png

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

24#
发表于 Post on 2023-5-20 16:36:17 | 只看该作者 Only view this author
nianbin 发表于 2023-5-20 09:10
您的意思是把这些去掉是吧?

嗯,从这里去掉也可以,把数据读进matlab再去掉也可以。但是只去掉较小的A值,例如那个-1.77MHz的A值应该是可以去掉的
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

182

帖子

0

威望

2233

eV
积分
2415

Level 5 (御坂)

25#
发表于 Post on 2023-5-20 20:53:09 | 只看该作者 Only view this author
wzkchem5 发表于 2023-5-20 16:36
嗯,从这里去掉也可以,把数据读进matlab再去掉也可以。但是只去掉较小的A值,例如那个-1.77MHz的A值应该 ...

请问下 是就删掉我划线那部分还是说下面的A(iso)连同后面的数值也要删掉

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

26#
发表于 Post on 2023-5-21 05:33:29 | 只看该作者 Only view this author
nianbin 发表于 2023-5-20 13:53
请问下 是就删掉我划线那部分还是说下面的A(iso)连同后面的数值也要删掉

我没删过property文件,但是总的宗旨就是,把你自己代入easyspin的视角,假如删完以后剩余的部分的格式和原文件格式相同,easyspin看这个文件觉得是orca少算了一些A值,看不出是你手动改的,就行
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

84

帖子

0

威望

372

eV
积分
456

Level 3 能力者

27#
发表于 Post on 2023-8-3 11:27:05 | 只看该作者 Only view this author
上周看了几篇Neese组的文献,而且Pantazis曾分享有一个视频,讲orca算epr很透彻,建议看看。

提取码:p55h


总体感觉:

1、只要不是固态体系,首推orca,固态体系目前可以尝试cp2k;

2、easyspin计算跟实验可比拟的谱图之前,需要计算g、A、D张量。



本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 03:44 , Processed in 0.219190 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list