版塊導航: 正在加載中...

調(diào)劑小程序

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (1664)
>考研 (121)
>導師招生 (113)
>蟲友互識 (68)
>教師之家 (28)
>考博 (25)
>論文道賀祈福 (24)
>碩博家園 (19)
>無機/物化 (19)
>論文投稿 (17)
>休閑灌水 (17)
>文獻求助 (13)
>公派出國 (13)
>基金申請 (9)
>留學生活 (7)
>找工作 (5)

返回列表

dhc198611

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2488.4
散金: 20
紅花: 2
帖子: 280
在線: 573.9小時
蟲號: 780288
注冊: 2009-05-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計算已有8人參與

我用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計算，采用nelec=10,nelup =5,neldw =5, 設(shè)置反鐵磁有序的。但是輸出的結(jié)果到100步還是不收斂，而且奇怪的是輸出結(jié)果居然沒有顯示受力。不知是怎么回事？希望哪位大俠幫忙解決。小弟非常感謝！

回復(fù)此樓

» 猜你喜歡

293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(fù)
281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請已經(jīng)有3人回復(fù)

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

平陽霉素與Mn_2_絡(luò)合物結(jié)構(gòu)的_1H_NMR弛豫研究已經(jīng)有3人回復(fù)
靜態(tài)計算能收斂，同樣的參數(shù)弛豫不能收斂，這是為什么？已經(jīng)有4人回復(fù)
誰做過224個LiFePO的弛豫啊。。。求救已經(jīng)有6人回復(fù)
什么是多弛豫時間呀？介電方面的已經(jīng)有6人回復(fù)
弛豫鐵電體的幾種模型已經(jīng)有10人回復(fù)
為什么做表面弛豫時，真空層位置不同會對結(jié)果有顯著影響？已經(jīng)有13人回復(fù)
【求助】非晶合金中的弛豫現(xiàn)象已經(jīng)有9人回復(fù)
【求助】A-Fe晶體中加H無法弛豫已經(jīng)有4人回復(fù)

科研無止境

1樓 2010-05-13 08:19:19

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

157806400

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 822.4
帖子: 27
在線: 6.5小時
蟲號: 796630
注冊: 2009-06-18
性別: GG
專業(yè): 功能材料計算

至少應(yīng)該有輸入輸出文件吧？

dhc198611(金幣+2):謝謝參與 2010-05-13 14:11:20

你這樣描述很抽象。

贊一下

回復(fù)此樓

2樓2010-05-13 11:16:54

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

dhc198611

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2488.4
散金: 20
紅花: 2
帖子: 280
在線: 573.9小時
蟲號: 780288
注冊: 2009-05-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by 157806400 at 2010-05-13 11:16:54:
你這樣描述很抽象。

這是我的INPIUT
nelec= 10,
nelup = 5,
neldw = 5,
occupations='fixed', degauss=0.0
lda_plus_u=.true., Hubbard_U(2)=2.0, Hubbard_U(3)=2.0,
starting_ns_eigenvalue(3,2,2) = 1.d0
starting_ns_eigenvalue(3,1,3) = 1.d0

這是加磁性和強關(guān)聯(lián)勢，這樣設(shè)置反鐵磁對不對

贊一下

回復(fù)此樓

科研無止境

3樓2010-05-13 14:11:01

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

★ ★
dhc198611(金幣+2):感謝你的寶貴意見，但是nstep最大也只能設(shè)置成100。 2010-05-13 21:19:05
zzy870720z(金幣+2):感謝指導 2010-05-13 22:25:49

我懷疑你的100步不收斂這里指的是scf計算不收斂，這是電子結(jié)構(gòu)計算的問題，有些參數(shù)沒有設(shè)置準確。你的scf計算都沒有收斂，當然不會有力的輸出。
relax迭代步數(shù)的設(shè)置是nstep，你可以設(shè)置大一些，計算收斂后會自動停止的。

贊一下(1人)

回復(fù)此樓

4樓2010-05-13 16:44:41

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

dhc198611(金幣+2): 2010-05-13 21:17:54

另外，我步明白你為什么要設(shè)置starting_ns_eigenvalue這個參數(shù)，你確定你的設(shè)置是對的嗎？如果不確定，建議你用默認值

贊一下

回復(fù)此樓

5樓2010-05-13 16:46:41

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

dhc198611

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2488.4
散金: 20
紅花: 2
帖子: 280
在線: 573.9小時
蟲號: 780288
注冊: 2009-05-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by goldenfisher at 2010-05-13 16:46:41:
另外，我步明白你為什么要設(shè)置starting_ns_eigenvalue這個參數(shù)，你確定你的設(shè)置是對的嗎？如果不確定，建議你用默認值

starting_ns_eigenvalue這個是+U的設(shè)置，即使默認頁還是不收斂。迭代100步的時候，不收斂就不會輸出力，還有PWSCF的默認迭代次數(shù)為100，即使你設(shè)置nstep=200，程序也會在100步的時候自動停掉。
我的主要問題在于，當加入
  nelec= 10,
nelup = 5,
neldw = 5,
occupations='fixed',  degauss=0.0
lda_plus_u=.true., Hubbard_U(2)=2.0, Hubbard_U(3)=2.0,
starting_ns_eigenvalue(3,2,2) = 1.d0
starting_ns_eigenvalue(3,1,3) = 1.d0
這幾句的時候，其體系就是不收斂，去掉的話就能收斂了，所以可能是我的反鐵磁有序設(shè)置有問題。
請問有沒有人知道怎么設(shè)置  nelec，nelup ， neldw 這3個參數(shù)，這3個參數(shù)怎么樣設(shè)置成反鐵磁有序。對于這3個參數(shù)我不太了解，但是按計算的提示是要加的。

贊一下

回復(fù)此樓

科研無止境

6樓2010-05-13 21:17:43

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

★ ★
dhc198611(金幣+2):感謝你的指點 2010-05-13 22:44:39
gavinliu7390(金幣+2):謝謝提示！ 2010-05-14 07:28:52

需要說明的是nstep不是電子自洽場迭代的步數(shù)，電子迭代步數(shù)應(yīng)該在&electron這個模塊設(shè)置，參數(shù)是electron_maxstep,你可以把這個參數(shù)設(shè)置為任意值。nstep指的是relax，也就是馳譽的步數(shù)。
另外，設(shè)置反鐵磁的例子在程序包里面有，你可以看一下。一般來說不是設(shè)置電子的個數(shù)，而是設(shè)置原子的初始磁矩大小和方向。你可以把input貼出來大家看看，如果你對大家有防御之心，就沒必要在這里問了。

贊一下(1人)

回復(fù)此樓

7樓2010-05-13 22:29:34

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

dhc198611

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2488.4
散金: 20
紅花: 2
帖子: 280
在線: 573.9小時
蟲號: 780288
注冊: 2009-05-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by goldenfisher at 2010-05-13 22:29:34:
需要說明的是nstep不是電子自洽場迭代的步數(shù)，電子迭代步數(shù)應(yīng)該在&electron這個模塊設(shè)置，參數(shù)是electron_maxstep,你可以把這個參數(shù)設(shè)置為任意值。nstep指的是relax，也就是馳譽的步數(shù)。
另外，設(shè)置反鐵磁的 ...

&control
calculation='relax',
restart_mode='from_scratch',
pseudo_dir = '~/pseudo/',
outdir='.',
prefix='BFO',
tprnfor = .true.
tstress = .true.
etot_conv_thr = 1e-5
forc_conv_thr = 5e-5
nstep=200
/
&system
ibrav=  5,
celldm(1) =10.7241882,
celldm(4) =0.51099,
nat=10,
ntyp= 4,
ecutwfc    = 35.0
nspin = 2
nelec= 10,
nelup = 5,
neldw = 5,
occupations='fixed',  degauss=0.0
lda_plus_u=.true., Hubbard_U(2)=2.0, Hubbard_U(3)=2.0,
/
&electrons
mixing_beta = 0.3,
/
&IONS
/
ATOMIC_SPECIES
Bi  208.98  Bi.pbe-d-mt.UPF
Fe1  55.845 Fe.pbe-sp-van_ak.UPF
Fe2  55.845 Fe.pbe-sp-van_ak.UPF
O  15.9994  O.pbe-van_ak.UPF
ATOMIC_POSITIONS {crystal}
Bi  0.0000000000000000  0.0000000000000000  0.0000000000000000 0  0  0
Bi  0.5000000000000000  0.5000000000000000  0.5000000000000000 0  0  0
Fe1  0.2225464345999172  0.2225464345999172  0.2225464345999172
Fe2  0.7225464345999169  0.7225464345999169  0.7225464345999169
O  0.5367314264732763 -0.0624437787918008  0.3877084883892429
O -0.0624437787918008  0.3877084883892430  0.5367314264732763
O  0.3877084883892429  0.5367314264732763 -0.0624437787918008
O  0.0367314264732762  0.8877084883892430  0.4375562212081991
O  0.8877084883892431  0.4375562212081992  0.0367314264732762
O  0.4375562212081991  0.0367314264732761  0.8877084883892430
K_POINTS {automatic}
  4 4 4 0 0 0

這是我的INPUT，我想收斂的話，應(yīng)該能在100內(nèi)能搞定的，但是電子步就是不收斂。我的初始位置是經(jīng)過VASP弛豫的，在PWSCF里應(yīng)該不會差很多的。

贊一下

回復(fù)此樓

科研無止境

8樓2010-05-13 22:44:20

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

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

★ ★
gavinliu7390(金幣+2):謝謝提示！ 2010-05-14 07:29:06
dhc198611(金幣+4):非常感謝 2010-05-14 07:44:26
dhc198611(金幣+1):期待你的回復(fù)（關(guān)于如何控制PWSCF弛豫時的對稱性） 2010-05-14 11:17:08

ecutwfc = 35.0可能小了點，另外應(yīng)該設(shè)置ecutrho=8*ecutwfc.另外，你根本沒有設(shè)置體系的自旋初始值，初始是默認為0的。start_magnetization這個參數(shù)你應(yīng)該設(shè)置好。
我覺得你應(yīng)該設(shè)置Fe1的自旋為1，F(xiàn)e2的自旋為-1，把體系的總磁矩置零。
我覺得你應(yīng)該仔細看input_pwscf文件，然后再看例子。例子里有計算FeO的，和你這個很相似。

贊一下(1人)

回復(fù)此樓

9樓2010-05-13 23:45:30

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

lono75

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6963.1
散金: 100
紅花: 2
帖子: 2433
在線: 536.4小時
蟲號: 932458
注冊: 2009-12-23
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
dhc198611(金幣+1):謝謝參與

http://valenhou.blog.edu.cn/2007/133214.html

[ Last edited by lono75 on 2010-5-14 at 10:49 ]

贊一下(1人)

回復(fù)此樓

由天命,成人事！

10樓2010-05-14 10:26:26

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 dhc198611 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

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

[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3	QGZDSYS 2026-03-13	3/150	2026-03-18 00:13 by PY-Hzb
[考研] 302求調(diào)劑 +6	呼呼呼。。。。 2026-03-17	6/300	2026-03-17 23:56 by 星空星月
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3	Eember. 2026-03-17	7/350	2026-03-17 20:20 by 花125533
[考研] 268求調(diào)劑 +8	一定有學上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 考研化學學碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 293求調(diào)劑 +6	世界首富 2026-03-11	6/300	2026-03-17 17:04 by ruiyingmiao
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 工科材料085601 279求調(diào)劑 +3	困于星晨 2026-03-17	3/150	2026-03-17 14:08 by ms629
[考研] 材料與化工304求B區(qū)調(diào)劑 +7	邱gl 2026-03-11	8/400	2026-03-17 09:36 by 努力學習賺彩禮
[考研] 267一志愿南京工業(yè)大學0817化工求調(diào)劑 +6	SUICHILD 2026-03-12	6/300	2026-03-17 09:24 by 霧散后相遇lc
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 070300化學學碩求調(diào)劑 +6	太想進步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[基金申請] 有必要更換申報口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky

24小時熱門版塊排行榜

[交流] 【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計算 已有8人參與

» 猜你喜歡

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

至少應(yīng)該有輸入輸出文件吧？

[交流] 【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計算已有8人參與

至少應(yīng)該有輸入輸出文件吧？