 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3096 | 回復(fù): 23 | |||
dhc198611木蟲 (小有名氣)
|
[交流]
【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計(jì)算 已有8人參與
|
| 我用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計(jì)算，采用nelec=10,nelup =5,neldw =5, 設(shè)置反鐵磁有序的。但是輸出的結(jié)果到100步還是不收斂，而且奇怪的是輸出結(jié)果居然沒有顯示受力。不知是怎么回事？希望哪位大俠幫忙解決。小弟非常感謝！ |
回復(fù)此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國重點(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博士申請
已經(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)價值貼推薦，對您同樣有幫助:
- 平陽霉素與Mn_2_絡(luò)合物結(jié)構(gòu)的_1H_NMR弛豫研究
已經(jīng)有3人回復(fù)
- 靜態(tài)計(jì)算能收斂，同樣的參數(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ù)
157806400
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 822.4
- 帖子: 27
- 在線: 6.5小時
- 蟲號: 796630
- 注冊: 2009-06-18
- 性別: GG
- 專業(yè): 功能材料計(jì)算
2樓2010-05-13 11:16:54
dhc198611
木蟲 (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2488.4
- 散金: 20
- 紅花: 2
- 帖子: 280
- 在線: 573.9小時
- 蟲號: 780288
- 注冊: 2009-05-26
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-05-13 14:11:01
goldenfisher
金蟲 (著名寫手)
- 1ST強(qiáng)帖: 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):感謝指導(dǎo) 2010-05-13 22:25:49
dhc198611(金幣+2):感謝你的寶貴意見，但是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ì)算的問題，有些參數(shù)沒有設(shè)置準(zhǔn)確。你的scf計(jì)算都沒有收斂，當(dāng)然不會有力的輸出。 relax迭代步數(shù)的設(shè)置是nstep，你可以設(shè)置大一些，計(jì)算收斂后會自動停止的。 |
4樓2010-05-13 16:44:41
goldenfisher
金蟲 (著名寫手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 55 (初中生)
- 貴賓: 0.153
- 金幣: 5779.1
- 散金: 2908
- 紅花: 21
- 帖子: 1552
- 在線: 351.5小時
- 蟲號: 117764
- 注冊: 2005-11-26
- 性別: GG
- 專業(yè): 原子和分子物理
5樓2010-05-13 16:46:41
dhc198611
木蟲 (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2488.4
- 散金: 20
- 紅花: 2
- 帖子: 280
- 在線: 573.9小時
- 蟲號: 780288
- 注冊: 2009-05-26
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
starting_ns_eigenvalue這個是+U的設(shè)置，即使默認(rèn)頁還是不收斂。迭代100步的時候，不收斂就不會輸出力，還有PWSCF的默認(rèn)迭代次數(shù)為100，即使你設(shè)置nstep=200，程序也會在100步的時候自動停掉。 我的主要問題在于，當(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è)置 nelec，nelup ， neldw 這3個參數(shù)，這3個參數(shù)怎么樣設(shè)置成反鐵磁有序。對于這3個參數(shù)我不太了解，但是按計(jì)算的提示是要加的。 |
6樓2010-05-13 21:17:43
goldenfisher
金蟲 (著名寫手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 55 (初中生)
- 貴賓: 0.153
- 金幣: 5779.1
- 散金: 2908
- 紅花: 21
- 帖子: 1552
- 在線: 351.5小時
- 蟲號: 117764
- 注冊: 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
|
需要說明的是nstep不是電子自洽場迭代的步數(shù)，電子迭代步數(shù)應(yīng)該在&electron這個模塊設(shè)置，參數(shù)是electron_maxstep,你可以把這個參數(shù)設(shè)置為任意值。nstep指的是relax，也就是馳譽(yù)的步數(shù)。 另外，設(shè)置反鐵磁的例子在程序包里面有，你可以看一下。一般來說不是設(shè)置電子的個數(shù)，而是設(shè)置原子的初始磁矩大小和方向。你可以把input貼出來大家看看，如果你對大家有防御之心，就沒必要在這里問了。 |
7樓2010-05-13 22:29:34
dhc198611
木蟲 (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2488.4
- 散金: 20
- 紅花: 2
- 帖子: 280
- 在線: 573.9小時
- 蟲號: 780288
- 注冊: 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)過VASP弛豫的，在PWSCF里應(yīng)該不會差很多的。 |
8樓2010-05-13 22:44:20
goldenfisher
金蟲 (著名寫手)
- 1ST強(qiáng)帖: 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
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可能小了點(diǎn)，另外應(yīng)該設(shè)置ecutrho=8*ecutwfc.另外，你根本沒有設(shè)置體系的自旋初始值，初始是默認(rèn)為0的。start_magnetization這個參數(shù)你應(yīng)該設(shè)置好。 我覺得你應(yīng)該設(shè)置Fe1的自旋為1，F(xiàn)e2的自旋為-1，把體系的總磁矩置零。 我覺得你應(yīng)該仔細(xì)看input_pwscf文件，然后再看例子。例子里有計(jì)算FeO的，和你這個很相似。 |
9樓2010-05-13 23:45:30
lono75
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6963.1
- 散金: 100
- 紅花: 2
- 帖子: 2433
- 在線: 536.4小時
- 蟲號: 932458
- 注冊: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +7 | yangfz 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +8 | 曦熙兮 2026-03-15 | 8/400 |
|
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +4 | zbzihdhd 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|