24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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 ]

回复此楼

» 猜你喜欢

广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有11人回复
闲聊已经有3人回复
江苏省自然基金什么时候出结果已经有4人回复
如何查询2026年获批项目名称？已经有3人回复
两块石头已经有10人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有19人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有15人回复
面上没中，邀请各位路过的虫友分析一下分数已经有13人回复
初秋的晨风已经有6人回复
Ei源刊怎么投已经有4人回复

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

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴

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

已阅回复此楼关注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足够了。