版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

汕头大学海洋科学、生物学、生物与医药等3个专业接受调剂

返回列表

minmin_0082003

金虫 (正式写手)

[交流] 【求助完结】Au原子的磁性计算为0，为什么不是1呢？求大侠解答

想用spin-polarized 去计算金原子，用的模守恒赝式，configuration 5d106s1
但是为什么计算得到的“
   total magnetization    =    0.00 Bohr mag/cell
   absolute magnetization =    0.00 Bohr mag/cell
”
总觉得是1才对。求大侠解答
输入文件如下
“&CONTROL
  calculation = "relax",
  prefix       = "Au-30"
  nstep       = 500,
  pseudo_dir = "/home/Au",
  outdir    = "/scr",
/
&SYSTEM
  ibrav                = 4,
  celldm(1)             = 33.123112,
  celldm(3)             = 1.42628937, !a=25Ang
  nat                   = 1,
  ntyp                = 1,
  ecutwfc             = 30.0D0,
  nspin                = 2,
  starting_magnetization(1) = 0.2,
  occupations    = "smearing",
  degauss       = 0.05,
/
&ELECTRONS
  conv_thr = 1.0D-8,
  mixing_beta = 0.3D0,
  electron_maxstep = 400,
/
&IONS
  ds             = 1.D-8,
/
ATOMIC_SPECIES
Au 196.96655 Au-wc.UPF
ATOMIC_POSITIONS { angstrom }
Au       3.755923000    7.951333000    7.475759000
K_POINTS {gamma}

”

[ Last edited by minmin_0082003 on 2011-1-20 at 11:26 ]

回复此楼

» 猜你喜欢

材料专硕调剂已经有3人回复
08生物与医药调剂302 已经有7人回复
085601材料工程求调剂已经有3人回复
269求调剂已经有3人回复
298求调剂已经有5人回复
302求调剂已经有4人回复
材料科学与工程290求调剂已经有15人回复
调剂已经有4人回复
求调剂！一志愿南京理工大学材料工程，总分344，有1篇SCI，1个专利已经有10人回复
344求调剂材料工程，一志愿西北工业大学已经有7人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求各种原子吸收标准曲线的配置方法啊或者有书推荐已经有18人回复
C原子分布在α-Fe 铁素体的什么间隙中？八面体还是四面体间隙？已经有6人回复
各位大侠，我可以优化结构，但只要一计算频率就出现下面的问题？已经有4人回复
紫外分光光度计和原子吸收分光光度计单色器位置为什么不同？已经有7人回复
vasp 计算单原子能量失败！请教各位！已经有8人回复
讨论：Au或Ag NPs均匀掺杂的微纳米粒子的制备方法已经有3人回复
计算一个簇模型时，被钴原子难住了～～求教。已经有7人回复
请问能用DFT方法计算一个体系中原子由于所处位置不同而势能不同，能输出原子势能吗已经有4人回复
入门生，处方单因素考察急切求解答非常感谢已经有8人回复
CASTEP在计算是出现问题啦，求大侠帮忙啊。。。已经有13人回复
如何把cif文件中有统计的占有率不是1的原子在画图时变为一个原子已经有6人回复
请教有机分子中O原子和H原子电荷数怎么计算，O-H键键级怎么计算！谢谢！已经有4人回复
【求助】Ms计算磁性要不要设置每个原子的磁矩？？已经有10人回复
【求助】用DFT+U(U=7eV)计算单个铜原子的能量为正值，是不是错了？已经有14人回复
【求助】请问用cc-pvtz-pp基组算一个C两个Au原子，怎么写输入文件，谢谢。已经有12人回复
【求助】紧急求助，sieasta优化后原子重合问题，忘大侠们帮忙已经有4人回复

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2011-01-13 11:15:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
zzy870720z(金币+2):谢谢指教 2011-01-13 16:32:01

ecutwfc太小了。你需要测试，norm-conserving需要的截断往往很大

赞一下(2人)

回复此楼

2楼2011-01-13 15:53:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

引用回帖:

Originally posted by goldenfisher at 2011-01-13 15:53:01:
ecutwfc太小了。你需要测试，norm-conserving需要的截断往往很大

非常感谢，我测试过的，确实比较大，我以为截断能不会影响到spin的计算结果，至少siesta的计算中是这样的，大截断能太耗时，所以想用小的截断，希望能做的快点.....
您意思是因为截断能的原因才使得spin的结果不正确吗？

赞一下

回复此楼

3楼2011-01-13 16:27:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
zzy870720z(金币+2):谢谢指点 2011-01-13 21:18:18

引用回帖:

Originally posted by minmin_0082003 at 2011-01-13 16:27:04:

非常感谢，我测试过的，确实比较大，我以为截断能不会影响到spin的计算结果，至少siesta的计算中是这样的，大截断能太耗时，所以想用小的截断，希望能做的快点.....
您意思是因为截断能的原因才使得spin的结 ...

很有可能的，阶段能大小直接影响到所有的结果，包括各种构型、电子结构等等，所以至少要测试到能量收敛时得到的才是可靠的。
你可以逐步加大试试，不行就用diagnization = 'cg' ，这个虽然慢，但是内存占得少。此外，你算一个单独的Au原子，个人认为最多15A的supercell足够了。

赞一下(1人)

回复此楼

4楼2011-01-13 16:36:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

1ST强帖: 1
应助: 136 (高中生)
金币: 3717.1
帖子: 988
在线: 1224.9小时
虫号: 961933

★ ★
zzy870720z(金币+2):谢谢指教 2011-01-13 21:18:38

Au的磁矩应该很小，我用lapw算了一个10Bohr和15Bohr边长盒子里的单独Au原子，结果是 -0.2377196275E-02玻尔磁子左右。
如果这样的话，感觉赝势法算准应该比较困难。

赞一下(1人)

回复此楼

5楼2011-01-13 18:21:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
youzhizhe(金币+2):谢谢交流。 2011-01-13 23:59:44

还有一个问题，也许算Au的时候应该考虑强关联效应。可以查查文献，肯定是有描述的。

赞一下(1人)

回复此楼

6楼2011-01-13 22:47:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

★
zzy870720z(金币+1): 谢谢分享 2011-01-20 13:12:01

引用回帖:

Originally posted by goldenfisher at 2011-01-13 22:47:39:
还有一个问题，也许算Au的时候应该考虑强关联效应。可以查查文献，肯定是有描述的。

找到一个不知道算不算得上办法的办法，设置 tot_magnetization，因为我可以确定单电子的数量，所以，我直接给定住得了

赞一下(1人)

回复此楼

7楼2011-01-14 09:48:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

★
zzy870720z(金币+1): 谢谢交流观点，哈哈 2011-01-20 13:12:24

这是因为DFT计算方法的缺陷

赞一下(1人)

回复此楼

8楼2011-01-20 11:27:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

newton3915

铁杆木虫 (正式写手)

应助: 4 (幼儿园)
金币: 4742.3
帖子: 418
在线: 391.3小时
虫号: 609189

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

244989楼: Originally posted by minmin_0082003 at 2011-01-13 11:15:12
想用spin-polarized 去计算金原子，用的模守恒赝式，configuration 5d106s1
但是为什么计算得到的“
total magnetization = 0.00 Bohr mag/cell
absolute magnetization = 0.00 ...

lz的金原子模守恒赝势能发一个吗，谢谢 544761094@qq.com

赞一下

回复此楼

9楼2012-06-21 20:25:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

应助: 21 (小学生)
金币: 4993.8
帖子: 381
在线: 363.2小时
虫号: 629899

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

有单电子的体系不一定有磁性吧，如果体系只有一个原子，我认为应该用dmol这样的软件算，这样单原子的磁性就不会因为周期性平均掉了。

[ Last edited by huazhorg on 2012-6-22 at 20:16 ]

赞一下

回复此楼

10楼2012-06-22 20:11:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 minmin_0082003 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 302求调剂 +4	想要上学??? 2025-04-05	4/200	2025-04-05 08:40 by 川哥v5
[考研] 一志愿哈尔滨工程大学356求调剂化工相关专业 +4	若月孤倚 2025-04-04	5/250	2025-04-05 07:38 by lfhuang
[考研] 302求调剂 +8	金鑫02 2025-04-04	8/400	2025-04-04 23:54 by 褚迎松小小
[考研] 314分一志愿统考西农植物保护学硕，求调剂 +3	洋葱蛋堡 2025-04-04	7/350	2025-04-04 22:03 by cbzhang
[基金申请] 面上项目一般是几位通讯评审专家？ +10	443933501 2025-04-01	11/550	2025-04-04 18:46 by arthas_007
[考研] 290求调剂 +8	赐予（求调剂） 2025-04-04	9/450	2025-04-04 11:56 by 子訡千羽
[考研] 材料与化工334调剂求助 +10	DPQ123 2025-04-03	10/500	2025-04-04 10:46 by 小郑….
[考研] 08工学 300分，可调剂上岸 +3	tzsbaqjs 2025-04-03	4/200	2025-04-03 23:32 by tjx66
[考研] 济南大学机器学习及其应用课题组拟接收计算机硕士 +3	成剑火火 2025-04-02	3/150	2025-04-03 15:53 by 18cuu
[考研] 302分求调剂 +4	绵绵。。。。 2025-04-03	4/200	2025-04-03 10:28 by Mornach1988
[考研] 总分288分材料与化工英语没过a区线求B区院校调剂 +3	咕噜咕噜oh 2025-04-02	3/150	2025-04-02 20:33 by 雨晴月夜
[考研] 304求调剂 +5	zqj0909 2025-03-30	5/250	2025-04-02 20:02 by 绿色港湾
[考研] 0703化学调剂 +3	哈哈哈@. 2025-04-02	3/150	2025-04-02 19:27 by 噬血摄魂
[考研] 接受跨专业调剂 +4	没学上咋整 2025-04-02	5/250	2025-04-02 17:18 by 我爱我家娃
[有机交流] 终产物纯化 100+3	hanweiguo 2025-04-01	5/250	2025-04-02 14:16 by xiaobenpu
[考研] 本科河南农业大学，一志愿苏州大学071000生物学289分求调剂 +6	潮潮玫瑰少年 2025-04-01	6/300	2025-04-02 09:45 by Cloris-y
[考研] 307求调剂 +5	王世宁 2025-04-01	5/250	2025-04-01 19:10 by nwpu315
[考研] 307求调剂 +4	梦想碰过 2025-03-30	4/200	2025-04-01 15:23 by lj8886
[考研] 08工科一志愿杭州电子科技大学学硕280分 +5	xxyyyybb 2025-03-31	5/250	2025-03-31 19:43 by panzt
[考研] 材料318找调剂！！！！！！！ +4	一颗糖。 2025-03-29	4/200	2025-03-30 09:28 by 娄娄虫

24小时热门版块排行榜

[交流] 【求助完结】Au原子的磁性计算为0，为什么不是1呢？求大侠解答

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:查看全部散金贴

第一性原理

» 抢金币啦！回帖就可以得到:
查看全部散金贴