计算化学公社

 找回密码 Forget password
 注册 Register
Views: 11284|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[CASTEP/Dmol3/MS] 求助:金属团簇的磁矩计算

[复制链接 Copy URL]

141

帖子

0

威望

1606

eV
积分
1747

Level 5 (御坂)

最近测试了Aln(n=2-7)小团簇的磁矩,发现了比较奇怪的结果:n=偶数的团簇总磁矩为0 ub,而n=奇数的团簇总磁矩为1ub。而参考文献中说:n=偶数的团簇总磁矩为2ub,而n=奇数的团簇总磁矩为1ub。
这是为什么呢?计算在20A*20A*20A的盒子中,参数:CASTEP,GGA-PBE, cutoff=300eV, k-point=1。

2.png (4.32 KB, 下载次数 Times of downloads: 61)

偶数

偶数

1.png (4.31 KB, 下载次数 Times of downloads: 69)

1.png

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

2#
发表于 Post on 2019-3-14 11:16:16 | 只看该作者 Only view this author
这种孤立体系,而且还是小团簇,用量化程序诸如Gaussian算比用CASTEP强得多得多,而且还能充分利用对称性,结合Multiwfn考察电子结构特征
那两种磁矩值相当于计算时设的自旋多重度不同,只有其中能量最低的才是基态,才对应真实情况
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

3754

帖子

3

威望

1万

eV
积分
19678

Level 6 (一方通行)

围观吃瓜群众

3#
发表于 Post on 2019-3-14 14:38:23 | 只看该作者 Only view this author
在Modify-Electronic Configuration里面把每个原子的磁矩都猜+1看看,Al默认可能是0就会算成反铁磁的,还有动能截断最好测试下。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

141

帖子

0

威望

1606

eV
积分
1747

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2019-3-14 19:48:49 | 只看该作者 Only view this author
sobereva 发表于 2019-3-14 11:16
这种孤立体系,而且还是小团簇,用量化程序诸如Gaussian算比用CASTEP强得多得多,而且还能充分利用对称性, ...

社长,Gaussian计算磁矩是用 命令NBO=npa,然后通过alpha与beta电子做差吗?

141

帖子

0

威望

1606

eV
积分
1747

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2019-3-14 19:49:59 | 只看该作者 Only view this author
卡开发发 发表于 2019-3-14 14:38
在Modify-Electronic Configuration里面把每个原子的磁矩都猜+1看看,Al默认可能是0就会算成反铁磁的,还有 ...

谢谢老师,我给Al都赋予磁矩+1之后,计算出来的磁矩终于正常了。通过这次测试,我感觉以后再计算类似磁性团簇的时候都要先赋予磁矩才行。

3754

帖子

3

威望

1万

eV
积分
19678

Level 6 (一方通行)

围观吃瓜群众

6#
发表于 Post on 2019-3-14 21:00:37 | 只看该作者 Only view this author
本帖最后由 卡开发发 于 2019-3-15 04:41 编辑
yiranfengbai 发表于 2019-3-14 19:49
谢谢老师,我给Al都赋予磁矩+1之后,计算出来的磁矩终于正常了。通过这次测试,我感觉以后再计算类似磁性 ...

磁性系统吃不准默认给出的磁矩是啥样的情况建议都手工做初猜。这里CASTEP一部分原子默认是0,很容易算成反铁磁或无磁性,而VASP开自旋极化默认是1,很容易算成铁磁。很多人觉得开了自旋极化就能算出来想要的结果,这个还不一定。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

6万

帖子

99

威望

5万

eV
积分
120159

管理员

公社社长

7#
发表于 Post on 2019-3-15 06:44:02 | 只看该作者 Only view this author
yiranfengbai 发表于 2019-3-14 19:48
社长,Gaussian计算磁矩是用 命令NBO=npa,然后通过alpha与beta电子做差吗?

跟pop=NPA毫无直接关系

仔细看了此文就懂了,专门说了磁矩
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353http://bbs.keinsci.com/thread-4422-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

141

帖子

0

威望

1606

eV
积分
1747

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2019-3-15 09:32:26 | 只看该作者 Only view this author
卡开发发 发表于 2019-3-14 21:00
磁性系统吃不准默认给出的磁矩是啥样的情况建议都手工做初猜。这里CASTEP一部分原子默认是0,很容易算成 ...

明白了,谢谢老师

141

帖子

0

威望

1606

eV
积分
1747

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2019-3-15 09:32:42 | 只看该作者 Only view this author
sobereva 发表于 2019-3-15 06:44
跟pop=NPA毫无直接关系

仔细看了此文就懂了,专门说了磁矩

谢谢社长了

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 13:24 , Processed in 1.700927 second(s), 23 queries , Gzip On.

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