请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 95|回复: 8

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

[复制链接]

54

帖子

0

威望

339

eV
积分
393

Level 3 能力者

发表于 2019-3-14 10:37:21 | 显示全部楼层 |阅读模式
最近测试了Aln(n=2-7)小团簇的磁矩,发现了比较奇怪的结果:n=偶数的团簇总磁矩为0 ub,而n=奇数的团簇总磁矩为1ub。而参考文献中说:n=偶数的团簇总磁矩为2ub,而n=奇数的团簇总磁矩为1ub。
这是为什么呢?计算在20A*20A*20A的盒子中,参数:CASTEP,GGA-PBE, cutoff=300eV, k-point=1。

偶数

偶数
1.png

1万

帖子

25

威望

2万

eV
积分
39415

管理员

公社社长

发表于 2019-3-14 11:16:16 | 显示全部楼层
这种孤立体系,而且还是小团簇,用量化程序诸如Gaussian算比用CASTEP强得多得多,而且还能充分利用对称性,结合Multiwfn考察电子结构特征
那两种磁矩值相当于计算时设的自旋多重度不同,只有其中能量最低的才是基态,才对应真实情况
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

1460

帖子

3

威望

5236

eV
积分
6756

Level 6 (一方通行)

Ab Initio Amateur

发表于 2019-3-14 14:38:23 | 显示全部楼层
在Modify-Electronic Configuration里面把每个原子的磁矩都猜+1看看,Al默认可能是0就会算成反铁磁的,还有动能截断最好测试下。
满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

54

帖子

0

威望

339

eV
积分
393

Level 3 能力者

 楼主| 发表于 2019-3-14 19:48:49 | 显示全部楼层
sobereva 发表于 2019-3-14 11:16
这种孤立体系,而且还是小团簇,用量化程序诸如Gaussian算比用CASTEP强得多得多,而且还能充分利用对称性, ...

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

54

帖子

0

威望

339

eV
积分
393

Level 3 能力者

 楼主| 发表于 2019-3-14 19:49:59 | 显示全部楼层
卡开发发 发表于 2019-3-14 14:38
在Modify-Electronic Configuration里面把每个原子的磁矩都猜+1看看,Al默认可能是0就会算成反铁磁的,还有 ...

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

1460

帖子

3

威望

5236

eV
积分
6756

Level 6 (一方通行)

Ab Initio Amateur

发表于 2019-3-14 21:00:37 | 显示全部楼层
本帖最后由 卡开发发 于 2019-3-15 04:41 编辑
yiranfengbai 发表于 2019-3-14 19:49
谢谢老师,我给Al都赋予磁矩+1之后,计算出来的磁矩终于正常了。通过这次测试,我感觉以后再计算类似磁性 ...

磁性系统吃不准默认给出的磁矩是啥样的情况建议都手工做初猜。这里CASTEP一部分原子默认是0,很容易算成反铁磁或无磁性,而VASP开自旋极化默认是1,很容易算成铁磁。很多人觉得开了自旋极化就能算出来想要的结果,这个还不一定。
满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

1万

帖子

25

威望

2万

eV
积分
39415

管理员

公社社长

发表于 7 天前 | 显示全部楼层
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程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

54

帖子

0

威望

339

eV
积分
393

Level 3 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
卡开发发 发表于 2019-3-14 21:00
磁性系统吃不准默认给出的磁矩是啥样的情况建议都手工做初猜。这里CASTEP一部分原子默认是0,很容易算成 ...

明白了,谢谢老师

54

帖子

0

威望

339

eV
积分
393

Level 3 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
sobereva 发表于 2019-3-15 06:44
跟pop=NPA毫无直接关系

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

谢谢社长了
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-3-22 17:00 , Processed in 0.155098 second(s), 26 queries .

快速回复 返回顶部 返回列表