 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (886)
- >導(dǎo)師招生 (89)
- >考研 (87)
- >蟲(chóng)友互識(shí) (36)
- >考博 (22)
- >論文道賀祈福 (20)
- >教師之家 (17)
- >碩博家園 (16)
- >無(wú)機(jī)/物化 (15)
- >文獻(xiàn)求助 (11)
- >論文投稿 (11)
- >休閑灌水 (8)
- >公派出國(guó) (7)
- >基金申請(qǐng) (6)
- >找工作 (5)
- >留學(xué)生活 (4)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3100 | 回復(fù): 23 | |||
dhc198611木蟲(chóng) (小有名氣)
|
[交流]
【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計(jì)算 已有8人參與
|
| 我用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計(jì)算，采用nelec=10,nelup =5,neldw =5, 設(shè)置反鐵磁有序的。但是輸出的結(jié)果到100步還是不收斂，而且奇怪的是輸出結(jié)果居然沒(méi)有顯示受力。不知是怎么回事？希望哪位大俠幫忙解決。小弟非常感謝！ |
回復(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)格式了，假專家開(kāi)始管格式了
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 平陽(yáng)霉素與Mn_2_絡(luò)合物結(jié)構(gòu)的_1H_NMR弛豫研究
已經(jīng)有3人回復(fù)
- 靜態(tài)計(jì)算能收斂，同樣的參數(shù)弛豫不能收斂，這是為什么？
已經(jīng)有4人回復(fù)
- 誰(shuí)做過(guò)224個(gè)LiFePO的弛豫啊。。。求救
已經(jīng)有6人回復(fù)
- 什么是多弛豫時(shí)間呀？介電方面的
已經(jīng)有6人回復(fù)
- 弛豫鐵電體的幾種模型
已經(jīng)有10人回復(fù)
- 為什么做表面弛豫時(shí)，真空層位置不同會(huì)對(duì)結(jié)果有顯著影響？
已經(jīng)有13人回復(fù)
- 【求助】非晶合金中的弛豫現(xiàn)象
已經(jīng)有9人回復(fù)
- 【求助】A-Fe晶體中加H無(wú)法弛豫
已經(jīng)有4人回復(fù)
157806400
金蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 822.4
- 帖子: 27
- 在線: 6.5小時(shí)
- 蟲(chóng)號(hào): 796630
- 注冊(cè): 2009-06-18
- 性別: GG
- 專業(yè): 功能材料計(jì)算
2樓2010-05-13 11:16:54
dhc198611
木蟲(chóng) (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2488.4
- 散金: 20
- 紅花: 2
- 帖子: 280
- 在線: 573.9小時(shí)
- 蟲(chóng)號(hào): 780288
- 注冊(cè): 2009-05-26
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-05-13 14:11:01
goldenfisher
金蟲(chóng) (著名寫(xiě)手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 55 (初中生)
- 貴賓: 0.153
- 金幣: 5779.1
- 散金: 2908
- 紅花: 21
- 帖子: 1552
- 在線: 351.5小時(shí)
- 蟲(chóng)號(hào): 117764
- 注冊(cè): 2005-11-26
- 性別: GG
- 專業(yè): 原子和分子物理
★ ★
dhc198611(金幣+2):感謝你的寶貴意見(jiàn)，但是nstep最大也只能設(shè)置成100。 2010-05-13 21:19:05
zzy870720z(金幣+2):感謝指導(dǎo) 2010-05-13 22:25:49
dhc198611(金幣+2):感謝你的寶貴意見(jiàn)，但是nstep最大也只能設(shè)置成100。 2010-05-13 21:19:05
zzy870720z(金幣+2):感謝指導(dǎo) 2010-05-13 22:25:49
|
我懷疑你的100步不收斂這里指的是scf計(jì)算不收斂，這是電子結(jié)構(gòu)計(jì)算的問(wèn)題，有些參數(shù)沒(méi)有設(shè)置準(zhǔn)確。你的scf計(jì)算都沒(méi)有收斂，當(dāng)然不會(huì)有力的輸出。 relax迭代步數(shù)的設(shè)置是nstep，你可以設(shè)置大一些，計(jì)算收斂后會(huì)自動(dòng)停止的。 |
4樓2010-05-13 16:44:41
goldenfisher
金蟲(chóng) (著名寫(xiě)手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 55 (初中生)
- 貴賓: 0.153
- 金幣: 5779.1
- 散金: 2908
- 紅花: 21
- 帖子: 1552
- 在線: 351.5小時(shí)
- 蟲(chóng)號(hào): 117764
- 注冊(cè): 2005-11-26
- 性別: GG
- 專業(yè): 原子和分子物理
5樓2010-05-13 16:46:41
dhc198611
木蟲(chóng) (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2488.4
- 散金: 20
- 紅花: 2
- 帖子: 280
- 在線: 573.9小時(shí)
- 蟲(chóng)號(hào): 780288
- 注冊(cè): 2009-05-26
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
starting_ns_eigenvalue這個(gè)是+U的設(shè)置，即使默認(rèn)頁(yè)還是不收斂。迭代100步的時(shí)候，不收斂就不會(huì)輸出力，還有PWSCF的默認(rèn)迭代次數(shù)為100，即使你設(shè)置nstep=200，程序也會(huì)在100步的時(shí)候自動(dòng)停掉。 我的主要問(wèn)題在于，當(dāng)加入 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è)置有問(wèn)題。 請(qǐng)問(wèn)有沒(méi)有人知道怎么設(shè)置 nelec，nelup ， neldw 這3個(gè)參數(shù)，這3個(gè)參數(shù)怎么樣設(shè)置成反鐵磁有序。對(duì)于這3個(gè)參數(shù)我不太了解，但是按計(jì)算的提示是要加的。 |
6樓2010-05-13 21:17:43
goldenfisher
金蟲(chóng) (著名寫(xiě)手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 55 (初中生)
- 貴賓: 0.153
- 金幣: 5779.1
- 散金: 2908
- 紅花: 21
- 帖子: 1552
- 在線: 351.5小時(shí)
- 蟲(chóng)號(hào): 117764
- 注冊(cè): 2005-11-26
- 性別: GG
- 專業(yè): 原子和分子物理
★ ★
dhc198611(金幣+2):感謝你的指點(diǎn) 2010-05-13 22:44:39
gavinliu7390(金幣+2):謝謝提示！ 2010-05-14 07:28:52
dhc198611(金幣+2):感謝你的指點(diǎn) 2010-05-13 22:44:39
gavinliu7390(金幣+2):謝謝提示！ 2010-05-14 07:28:52
|
需要說(shuō)明的是nstep不是電子自洽場(chǎng)迭代的步數(shù)，電子迭代步數(shù)應(yīng)該在&electron這個(gè)模塊設(shè)置，參數(shù)是electron_maxstep,你可以把這個(gè)參數(shù)設(shè)置為任意值。nstep指的是relax，也就是馳譽(yù)的步數(shù)。 另外，設(shè)置反鐵磁的例子在程序包里面有，你可以看一下。一般來(lái)說(shuō)不是設(shè)置電子的個(gè)數(shù)，而是設(shè)置原子的初始磁矩大小和方向。你可以把input貼出來(lái)大家看看，如果你對(duì)大家有防御之心，就沒(méi)必要在這里問(wèn)了。 |
7樓2010-05-13 22:29:34
dhc198611
木蟲(chóng) (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2488.4
- 散金: 20
- 紅花: 2
- 帖子: 280
- 在線: 573.9小時(shí)
- 蟲(chóng)號(hào): 780288
- 注冊(cè): 2009-05-26
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
&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)過(guò)VASP弛豫的，在PWSCF里應(yīng)該不會(huì)差很多的。 |
8樓2010-05-13 22:44:20
goldenfisher
金蟲(chóng) (著名寫(xiě)手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 55 (初中生)
- 貴賓: 0.153
- 金幣: 5779.1
- 散金: 2908
- 紅花: 21
- 帖子: 1552
- 在線: 351.5小時(shí)
- 蟲(chóng)號(hào): 117764
- 注冊(cè): 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弛豫時(shí)的對(duì)稱性） 2010-05-14 11:17:08
gavinliu7390(金幣+2):謝謝提示！ 2010-05-14 07:29:06
dhc198611(金幣+4):非常感謝 2010-05-14 07:44:26
dhc198611(金幣+1):期待你的回復(fù)（關(guān)于如何控制PWSCF弛豫時(shí)的對(duì)稱性） 2010-05-14 11:17:08
|
ecutwfc = 35.0可能小了點(diǎn)，另外應(yīng)該設(shè)置ecutrho=8*ecutwfc.另外，你根本沒(méi)有設(shè)置體系的自旋初始值，初始是默認(rèn)為0的。start_magnetization這個(gè)參數(shù)你應(yīng)該設(shè)置好。 我覺(jué)得你應(yīng)該設(shè)置Fe1的自旋為1，F(xiàn)e2的自旋為-1，把體系的總磁矩置零。 我覺(jué)得你應(yīng)該仔細(xì)看input_pwscf文件，然后再看例子。例子里有計(jì)算FeO的，和你這個(gè)很相似。 |
9樓2010-05-13 23:45:30
lono75
鐵桿木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6963.1
- 散金: 100
- 紅花: 2
- 帖子: 2433
- 在線: 536.4小時(shí)
- 蟲(chóng)號(hào): 932458
- 注冊(cè): 2009-12-23
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★
dhc198611(金幣+1):謝謝參與
dhc198611(金幣+1):謝謝參與
|
http://valenhou.blog.edu.cn/2007/133214.html [ Last edited by lono75 on 2010-5-14 at 10:49 ] |
10樓2010-05-14 10:26:26
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 dhc198611 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[論文投稿]
投稿問(wèn)題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 化工學(xué)碩306求調(diào)劑 +9 | 42838695 2026-03-12 | 9/450 |
|