计算化学公社

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

[新手求助] 如何进行ESR顺磁谱图的模拟?

[复制链接 Copy URL]

14

帖子

0

威望

279

eV
积分
293

Level 3 能力者

16#
发表于 Post on 2023-5-19 15:35:19 | 只看该作者 Only view this author
wzkchem5 发表于 2023-5-18 18:31
手动删掉那些很小的、你不用计算就知道不会对谱图有显著影响的A值,再算。easyspin计算慢一般都是因为考 ...

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

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

17#
发表于 Post on 2023-5-19 16:17:14 | 只看该作者 Only view this author
自由风 发表于 2023-5-19 08:35
我算的A值大的有130,小的有15,最小的只有个位数。那手动删除的话15左右的算小吗?大小是指的绝对值吗?

结合你所用的lwpp,判断这个较小的A值有或者没有会不会显著影响峰形。最简单的方法是,另模拟一个只有一个A值的信号,其A值等于比如说15,lwpp和你当前模拟相同。然后看看去掉这个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
本团队长期招收研究生,有意者可私信联系

210

帖子

1

威望

4352

eV
积分
4582

Level 6 (一方通行)

18#
发表于 Post on 2023-5-19 16:17:19 | 只看该作者 Only view this author
之前用ORCA计算的时候用BP/EPR-II计算的结果和实验值最为吻合,不过应该是case by case,仅供参考。
Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

120

帖子

0

威望

2376

eV
积分
2496

Level 5 (御坂)

#未来可寄

19#
发表于 Post on 2023-6-13 10:24:20 | 只看该作者 Only view this author
考虑的原子数太多Matlab也会崩掉,我的建议是先把自旋密度图画出来,然后根据自旋密度分布大的位置,去选择计算那些磁性核。总的来说条件就是:磁性核、该磁性核周围自旋密度大

120

帖子

0

威望

2376

eV
积分
2496

Level 5 (御坂)

#未来可寄

20#
发表于 Post on 2023-7-6 16:23:40 | 只看该作者 Only view this author
wangzhe 发表于 2023-5-19 16:17
之前用ORCA计算的时候用BP/EPR-II计算的结果和实验值最为吻合,不过应该是case by case,仅供参考。

是BP86吗?

210

帖子

1

威望

4352

eV
积分
4582

Level 6 (一方通行)

21#
发表于 Post on 2023-7-6 16:49:21 | 只看该作者 Only view this author

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

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
QuantumicGuy + 2 thx

查看全部评分 View all ratings

Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

22#
发表于 Post on 2023-7-6 17:19:52 | 只看该作者 Only view this author
wangzhe 发表于 2023-7-6 09:49
BP86和BP计算的D和E张量结果是一样的。不知道这两个是否是同一个泛函。

是一个泛函
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
本团队长期招收研究生,有意者可私信联系

210

帖子

1

威望

4352

eV
积分
4582

Level 6 (一方通行)

23#
发表于 Post on 2023-7-6 17:27:22 | 只看该作者 Only view this author

感谢
Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

8

帖子

0

威望

92

eV
积分
100

Level 2 能力者

24#
发表于 Post on 2024-12-4 17:20:45 | 只看该作者 Only view this author
请问一下老师们,用orca+easyspin模拟EPR曲线,可以把数据导出用origin作图吗

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

25#
发表于 Post on 2024-12-4 18:04:57 | 只看该作者 Only view this author
青花 发表于 2024-12-4 10:20
请问一下老师们,用orca+easyspin模拟EPR曲线,可以把数据导出用origin作图吗

可以,matlab有一个函数叫做dlmwrite,可以把matlab里的任意向量或矩阵导出成csv格式,具体在matlab帮助窗口里搜
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
本团队长期招收研究生,有意者可私信联系

8

帖子

0

威望

92

eV
积分
100

Level 2 能力者

26#
发表于 Post on 2024-12-4 18:36:04 | 只看该作者 Only view this author
wzkchem5 发表于 2024-12-4 18:04
可以,matlab有一个函数叫做dlmwrite,可以把matlab里的任意向量或矩阵导出成csv格式,具体在matlab帮助 ...

谢谢老师,成功导出了。我看之前楼里之前的回复,可不可以理解为H数量过大时,允许增大lwpp(比如lwpp=1 → lwpp=2)用于抵消一部分未计算的较小A值的峰。

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

27#
发表于 Post on 2024-12-4 18:39:42 | 只看该作者 Only view this author
青花 发表于 2024-12-4 11:36
谢谢老师,成功导出了。我看之前楼里之前的回复,可不可以理解为H数量过大时,允许增大lwpp(比如lwpp=1  ...

对,调节lwpp直至和实验谱图大致相符,然后明确说明这里lwpp是调出来的,不是算出来的,因此峰宽和实验相符不说明算得准
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
本团队长期招收研究生,有意者可私信联系

8

帖子

0

威望

92

eV
积分
100

Level 2 能力者

28#
发表于 Post on 2024-12-4 19:44:16 | 只看该作者 Only view this author
wzkchem5 发表于 2024-12-4 18:39
对,调节lwpp直至和实验谱图大致相符,然后明确说明这里lwpp是调出来的,不是算出来的,因此峰宽和实验相 ...

感谢,再麻烦您一下。仪器上设置的Modulation Amplitude是我们用的lwpp对吗,然后实验谱图出来之后展宽与Modulation Amplitude设定值会有一些出入,我们再以实际为准再来计算,是这样?麻烦了

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

29#
发表于 Post on 2024-12-4 20:34:38 | 只看该作者 Only view this author
青花 发表于 2024-12-4 12:44
感谢,再麻烦您一下。仪器上设置的Modulation Amplitude是我们用的lwpp对吗,然后实验谱图出来之后展宽与 ...

不是。EPR的纵轴是EPR响应信号对磁场的导数(dX/dB,X其实是希腊字母chi),仪器测量的时候通过让B小幅度振荡,测X振荡的幅度,这样推出dX/dB,Modulation Amplitude就是B小幅度振荡的幅度。根据导数的定义知,Modulation Amplitude越小,结果的系统误差越大,但实验噪音的影响也越大(也就是随机误差越大)。Modulation Amplitude越大,lwpp越大,但Modulation Amplitude趋于零的时候lwpp不趋于0,因为自然展宽、振动展宽等等也会导致非零lwpp,Modulation Amplitude只是总展宽里面仪器展宽的那部分,甚至可能都不是仪器展宽的全部。所以lwpp没有任何理由应该接近或等于Modulation Amplitude
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
本团队长期招收研究生,有意者可私信联系

8

帖子

0

威望

92

eV
积分
100

Level 2 能力者

30#
发表于 Post on 2024-12-4 21:10:34 | 只看该作者 Only view this author
wzkchem5 发表于 2024-12-4 20:34
不是。EPR的纵轴是EPR响应信号对磁场的导数(dX/dB,X其实是希腊字母chi),仪器测量的时候通过让B小幅度 ...

太感谢了,明白了,谢谢您

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 13:39 , Processed in 0.251871 second(s), 23 queries , Gzip On.

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