minmin_0082003

金蟲 (正式寫手)

[交流] 【求助完結(jié)】Au原子的磁性計(jì)算為0，為什么不是1呢？求大俠解答

想用spin-polarized 去計(jì)算金原子，用的模守恒贗式，configuration 5d106s1
但是為什么計(jì)算得到的“
   total magnetization    =    0.00 Bohr mag/cell
   absolute magnetization =    0.00 Bohr mag/cell
”
總覺得是1才對(duì)。求大俠解答
輸入文件如下
“&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 ]

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請(qǐng) 已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)

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

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

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

查看全部散金貼

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

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

goldenfisher

金蟲 (著名寫手)

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

ecutwfc太小了。你需要測(cè)試，norm-conserving需要的截?cái)嗤艽?

贊一下(2人)

回復(fù)此樓

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

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

minmin_0082003

金蟲 (正式寫手)

引用回帖:

Originally posted by goldenfisher at 2011-01-13 15:53:01:
ecutwfc太小了。你需要測(cè)試，norm-conserving需要的截?cái)嗤艽?

非常感謝，我測(cè)試過的，確實(shí)比較大，我以為截?cái)嗄懿粫?huì)影響到spin的計(jì)算結(jié)果，至少siesta的計(jì)算中是這樣的，大截?cái)嗄芴臅r(shí)，所以想用小的截?cái)�，希望能做的快點(diǎn).....
您意思是因?yàn)榻財(cái)嗄艿脑虿攀沟胹pin的結(jié)果不正確嗎？

贊一下

回復(fù)此樓

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

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

goldenfisher

金蟲 (著名寫手)

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

引用回帖:

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

非常感謝，我測(cè)試過的，確實(shí)比較大，我以為截?cái)嗄懿粫?huì)影響到spin的計(jì)算結(jié)果，至少siesta的計(jì)算中是這樣的，大截?cái)嗄芴臅r(shí)，所以想用小的截?cái)啵Ｍ茏龅目禳c(diǎn).....
您意思是因?yàn)榻財(cái)嗄艿脑虿攀沟胹pin的結(jié) ...

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

贊一下(1人)

回復(fù)此樓

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

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

銳利的碎片

木蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 136 (高中生)
金幣: 3637.1
帖子: 988
在線: 1224.9小時(shí)
蟲號(hào): 961933

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

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

贊一下(1人)

回復(fù)此樓

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

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

goldenfisher

金蟲 (著名寫手)

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

還有一個(gè)問題，也許算Au的時(shí)候應(yīng)該考慮強(qiáng)關(guān)聯(lián)效應(yīng)�？梢圆椴槲墨I(xiàn)，肯定是有描述的。

贊一下(1人)

回復(fù)此樓

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

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

minmin_0082003

金蟲 (正式寫手)

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

引用回帖:

Originally posted by goldenfisher at 2011-01-13 22:47:39:
還有一個(gè)問題，也許算Au的時(shí)候應(yīng)該考慮強(qiáng)關(guān)聯(lián)效應(yīng)�？梢圆椴槲墨I(xiàn)，肯定是有描述的。

找到一個(gè)不知道算不算得上辦法的辦法，設(shè)置 tot_magnetization，因?yàn)槲铱梢源_定單電子的數(shù)量，所以，我直接給定住得了

贊一下(1人)

回復(fù)此樓

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

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

minmin_0082003

金蟲 (正式寫手)

★
zzy870720z(金幣+1): 謝謝交流觀點(diǎn)，哈哈 2011-01-20 13:12:24

這是因?yàn)镈FT計(jì)算方法的缺陷

贊一下(1人)

回復(fù)此樓

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

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

newton3915

鐵桿木蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 4742.3
帖子: 418
在線: 391.3小時(shí)
蟲號(hào): 609189

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

引用回帖:

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

lz的金原子模守恒贗勢(shì)能發(fā)一個(gè)嗎，謝謝 544761094@qq.com

贊一下

回復(fù)此樓

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

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

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
帖子: 381
在線: 363.2小時(shí)
蟲號(hào): 629899

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

有單電子的體系不一定有磁性吧，如果體系只有一個(gè)原子，我認(rèn)為應(yīng)該用dmol這樣的軟件算，這樣單原子的磁性就不會(huì)因?yàn)橹芷谛云骄袅恕?br />
[ Last edited by huazhorg on 2012-6-22 at 20:16 ]

贊一下

回復(fù)此樓

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

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

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

[考研] 301求調(diào)劑 +6	yy要上岸呀 2026-03-17	6/300	2026-03-17 23:58 by 星空星月
[考研] 274求調(diào)劑0856材料化工 +13	z2839474511 2026-03-11	14/700	2026-03-17 16:51 by share_joy
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 0703化學(xué)調(diào)劑 +6	妮妮ninicgb 2026-03-15	9/450	2026-03-16 16:40 by houyaoxu
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5	Becho359 2026-03-11	5/250	2026-03-14 11:35 by 哦哦123
[考研] 材料專碩288分求調(diào)劑一志愿211 +4	在家想你 2026-03-11	4/200	2026-03-13 22:49 by JourneyLucky
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 化工學(xué)碩306求調(diào)劑 +9	42838695 2026-03-12	9/450	2026-03-13 10:16 by houyaoxu