minmin_0082003

金蟲 (正式寫手)

[交流] 【求助完結(jié)】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 ]

回復此樓

» 猜你喜歡

工科材料085601 279求調(diào)劑已經(jīng)有3人回復
301求調(diào)劑已經(jīng)有3人回復
材料專碩306英一數(shù)二已經(jīng)有8人回復
材料與化工專碩調(diào)劑已經(jīng)有5人回復
085600材料與化工已經(jīng)有4人回復
085601材料工程專碩求調(diào)劑已經(jīng)有3人回復
293求調(diào)劑已經(jīng)有10人回復
有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～已經(jīng)有6人回復
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有3人回復
0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 已經(jīng)有3人回復

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求各種原子吸收標準曲線的配置方法啊或者有書推薦已經(jīng)有18人回復
C原子分布在α-Fe 鐵素體的什么間隙中？八面體還是四面體間隙？已經(jīng)有6人回復
各位大俠，我可以優(yōu)化結(jié)構(gòu)，但只要一計算頻率就出現(xiàn)下面的問題？已經(jīng)有4人回復
紫外分光光度計和原子吸收分光光度計單色器位置為什么不同？已經(jīng)有7人回復
討論：Au或Ag NPs均勻摻雜的微納米粒子的制備方法已經(jīng)有3人回復
計算一個簇模型時，被鈷原子難住了～～求教。已經(jīng)有7人回復
請問能用DFT方法計算一個體系中原子由于所處位置不同而勢能不同，能輸出原子勢能嗎已經(jīng)有4人回復
入門生，處方單因素考察急切求解答非常感謝已經(jīng)有8人回復
CASTEP在計算是出現(xiàn)問題啦，求大俠幫忙啊。。。已經(jīng)有13人回復
如何把cif文件中有統(tǒng)計的占有率不是1的原子在畫圖時變?yōu)橐粋€原子已經(jīng)有6人回復
請教有機分子中O原子和H原子電荷數(shù)怎么計算，O-H鍵鍵級怎么計算！謝謝！已經(jīng)有4人回復
【求助】Ms計算磁性要不要設(shè)置每個原子的磁矩？？已經(jīng)有10人回復
【求助】請問用cc-pvtz-pp基組算一個C兩個Au原子，怎么寫輸入文件，謝謝。已經(jīng)有12人回復

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

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

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

goldenfisher

金蟲 (著名寫手)

★ ★
zzy870720z(金幣+2):謝謝指點 2011-01-13 21:18:18

引用回帖:

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

非常感謝，我測試過的，確實比較大，我以為截斷能不會影響到spin的計算結(jié)果，至少siesta的計算中是這樣的，大截斷能太耗時，所以想用小的截斷，希望能做的快點.....
您意思是因為截斷能的原因才使得spin的結(jié) ...

很有可能的，階段能大小直接影響到所有的結(jié)果，包括各種構(gòu)型、電子結(jié)構(gòu)等等，所以至少要測試到能量收斂時得到的才是可靠的。
你可以逐步加大試試，不行就用diagnization = 'cg' ，這個雖然慢，但是內(nèi)存占得少。此外，你算一個單獨的Au原子，個人認為最多15A的supercell足夠了。

贊一下(1人)

回復此樓

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

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 10 個回答

goldenfisher

金蟲 (著名寫手)

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

ecutwfc太小了。你需要測試，norm-conserving需要的截斷往往很大

贊一下(2人)

回復此樓

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

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

minmin_0082003

金蟲 (正式寫手)

引用回帖:

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

非常感謝，我測試過的，確實比較大，我以為截斷能不會影響到spin的計算結(jié)果，至少siesta的計算中是這樣的，大截斷能太耗時，所以想用小的截斷，希望能做的快點.....
您意思是因為截斷能的原因才使得spin的結(jié)果不正確嗎？

贊一下

回復此樓

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

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

銳利的碎片

木蟲 (正式寫手)

1ST強帖: 1
應助: 136 (高中生)
金幣: 3637.1
帖子: 988
在線: 1224.9小時
蟲號: 961933

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

Au的磁矩應該很小，我用lapw算了一個10Bohr和15Bohr邊長盒子里的單獨Au原子，結(jié)果是 -0.2377196275E-02玻爾磁子左右。
如果這樣的話，感覺贗勢法算準應該比較困難。

贊一下(1人)

回復此樓

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

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考研] 085601材料工程專碩求調(diào)劑 +3	慕寒mio 2026-03-16	3/150	2026-03-17 13:46 by houyaoxu
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +3	share_joy 2026-03-16	3/150	2026-03-17 13:41 by houyaoxu
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 344求調(diào)劑 +3	knight344 2026-03-16	3/150	2026-03-16 09:42 by 無際的草原
[考研] 機械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 0805，333求調(diào)劑 +3	112253525 2026-03-10	3/150	2026-03-13 23:42 by JourneyLucky
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 【0856】化學工程（085602）313 分，本科學科評估A類院�；瘜W工程與工藝，誠求調(diào)劑 +7	小劉快快上岸 2026-03-11	7/350	2026-03-13 16:06 by ruiyingmiao
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 化工學碩306求調(diào)劑 +9	42838695 2026-03-12	9/450	2026-03-13 10:16 by houyaoxu
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321
[考研] 279求調(diào)劑 +3	莫xiao 2026-03-10	4/200	2026-03-11 08:06 by 斬魂滴兔子！