qnclqq

木蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 3862.6
紅花: 4
帖子: 467
在線: 132.8小時(shí)
蟲號: 1464347
注冊: 2011-10-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] +U后，晶格弛豫無法實(shí)現(xiàn)力收斂已有1人參與

對TiO2進(jìn)行+U計(jì)算，晶格弛豫（vc-relax）過程力（stress）不收斂，求指導(dǎo)。
軟件版本：QE 5.4
操作環(huán)境：Cent OS
代碼如下：
&CONTROL
  calculation  = "vc-relax",
!  restart_mode = 'restart',
  prefix    = "rutile",
  pseudo_dir = "/root/Desktop/qe/pseudo",
  outdir    = "./",
/
&SYSTEM
  ibrav = 6,  A = 4.594,  C = 2.959,  nat = 6,  ntyp  = 2,
  ecutwfc = 60,  ecutrho = 500, nbnd = 100,
  nspin = 2,  starting_magnetization(1)=0.1,  starting_magnetization(2)=0.1,
  occupations= 'smearing',  smearing = 'gaussian',  degauss = 0.02
  lda_plus_u=.true. Hubbard_U(1)=1, Hubbard_U(2)=2,
/
&ELECTRONS
  conv_thr = 1.D-6,
  mixing_beta = 0.7D0,
/
&IONS
/
&CELL
/
ATOMIC_SPECIES
O 15.9989995956  O.pz-van_ak.UPF
Ti  47.9000015259  Ti.pz-sp-van.UPF

ATOMIC_POSITIONS {crystal}
O 0.30479974  0.30479974  0.0000000000000000
O 0.69520026  0.69520026 0.0000000000000000
O 0.1952002612102742 0.8047997387897258 0.5000000000000000
O 0.8047997387897258 0.1952002612102742 0.5000000000000000
Ti 0.0000000000000000 0.0000000000000000 0.0000000000000000
Ti 0.5000000000000000 0.5000000000000000 0.5000000000000000

K_POINTS {automatic}
  3 3 5 0 0 0

回復(fù)此樓

» 猜你喜歡

323求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有4人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

1樓 2016-07-18 11:46:51

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

回帖支持 ( 顯示支持度最高的前 50 名 )

張玲海

木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 3027.1
散金: 766
紅花: 3
帖子: 399
在線: 1577.4小時(shí)
蟲號: 1419918
注冊: 2011-09-27
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

★ ★
漫天飄雪: 金幣+2, 謝謝交流 2016-07-21 22:54:38

Hubbard_U(1)=1, Hubbard_U(2)=2 設(shè)置不對。mixing_beta 太大了。 ecutwfc = 80。nbnd = 100 感覺沒有這么多

贊一下(1人)

回復(fù)此樓

簡單的生活

3樓2016-07-20 16:10:11

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

普通回帖

hexin1982

新蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 546
帖子: 148
在線: 259.6小時(shí)
蟲號: 1835611
注冊: 2012-05-26
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

你確定你的計(jì)算中TiO2要加自旋極化嗎？

贊一下

回復(fù)此樓

2樓2016-07-19 20:51:20

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

qnclqq

木蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 3862.6
紅花: 4
帖子: 467
在線: 132.8小時(shí)
蟲號: 1464347
注冊: 2011-10-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by hexin1982 at 2016-07-19 20:51:20
你確定你的計(jì)算中TiO2要加自旋極化嗎？

自旋極化只是加上看看效果，如果沒有極化并不會影響最后的效果。問題在于+U之后無法實(shí)現(xiàn)vc-relax，結(jié)果顯示力是不收斂的。

贊一下

回復(fù)此樓

4樓2016-07-21 17:00:20

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

qnclqq

木蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 3862.6
紅花: 4
帖子: 467
在線: 132.8小時(shí)
蟲號: 1464347
注冊: 2011-10-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

3樓: Originally posted by 張玲海 at 2016-07-20 16:10:11
Hubbard_U(1)=1, Hubbard_U(2)=2 設(shè)置不對。mixing_beta 太大了。 ecutwfc = 80。nbnd = 100 感覺沒有這么多

對于金紅石單胞，這個(gè)能量確實(shí)設(shè)置比較高，這樣能夠達(dá)到我的收斂標(biāo)準(zhǔn)。能帶數(shù)目，本來就有48個(gè)價(jià)電子，增大了一倍隨手設(shè)置了100。理論上，能量和能帶數(shù)也只是影響運(yùn)算速度，應(yīng)該不會導(dǎo)致不收斂。

mixing_beta應(yīng)該是對收斂速度有影響，這里在不加U的時(shí)候是可以的，+U之后我沒有改。修改成為0.2后，結(jié)果類似，力不收斂。

Hubbard_U(1)=1, Hubbard_U(2)=2，這樣設(shè)置，我的本意是試一下給O加1eV，給Ti加2eV。請問錯(cuò)在哪里，應(yīng)如何改正呢？

贊一下

回復(fù)此樓

5樓2016-07-21 17:26:02

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

張玲海

木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 3027.1
散金: 766
紅花: 3
帖子: 399
在線: 1577.4小時(shí)
蟲號: 1419918
注冊: 2011-09-27
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

可以先不加U優(yōu)化一下。完了再加上U。氧沒有必要加U吧

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

簡單的生活

6樓2016-07-21 17:47:01

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

qnclqq

木蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 3862.6
紅花: 4
帖子: 467
在線: 132.8小時(shí)
蟲號: 1464347
注冊: 2011-10-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

6樓: Originally posted by 張玲海 at 2016-07-21 17:47:01
可以先不加U優(yōu)化一下。完了再加上U。氧沒有必要加U吧

不加U是可以的，很快就可以完成。單獨(dú)的+U+relax也可以完成。但是vc-relax+U就不行了。之所以給氧+U，是想看看效果。因?yàn)門iO2是電荷轉(zhuǎn)移絕緣體，價(jià)帶是氧的p電子，所以感覺單獨(dú)給Ti加U是不夠的。

贊一下

回復(fù)此樓

7樓2016-07-21 19:10:37

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

張玲海

木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 3027.1
散金: 766
紅花: 3
帖子: 399
在線: 1577.4小時(shí)
蟲號: 1419918
注冊: 2011-09-27
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

7樓: Originally posted by qnclqq at 2016-07-21 19:10:37
不加U是可以的，很快就可以完成。單獨(dú)的+U+relax也可以完成。但是vc-relax+U就不行了。之所以給氧+U，是想看看效果。因?yàn)門iO2是電荷轉(zhuǎn)移絕緣體，價(jià)帶是氧的p電子，所以感覺單獨(dú)給Ti加U是不夠的。...

先不加U，relax一下，得到坐標(biāo)。然后把這個(gè)優(yōu)化的坐標(biāo) 重新輸入input文件。再用 vc-relax+U。把conv_thr = 1.D-6 提高一點(diǎn) conv_thr = 1.D-8看看。

贊一下

回復(fù)此樓

簡單的生活

8樓2016-07-21 19:30:44

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

applezju

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 2723.5
散金: 90
紅花: 11
沙發(fā): 1
帖子: 365
在線: 343.5小時(shí)
蟲號: 2396394
注冊: 2013-04-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

我用的別的一個(gè)軟件也是，加U的話收斂比較難

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

混，混混

9樓2016-07-21 20:43:20

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

qnclqq

木蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 3862.6
紅花: 4
帖子: 467
在線: 132.8小時(shí)
蟲號: 1464347
注冊: 2011-10-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

8樓: Originally posted by 張玲海 at 2016-07-21 19:30:44
先不加U，relax一下，得到坐標(biāo)。然后把這個(gè)優(yōu)化的坐標(biāo) 重新輸入input文件。再用 vc-relax+U。把conv_thr = 1.D-6 提高一點(diǎn) conv_thr = 1.D-8看看。...

剛又試了幾組參數(shù)，發(fā)現(xiàn)如果只給Ti加U，是可以vc-relax的。但是一旦給O也加上，只能運(yùn)行3步bfgs，后面能量就可是開始不收斂了（力也不收斂）

贊一下

回復(fù)此樓

10樓2016-07-21 21:01:09

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

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

[考研] 323求調(diào)劑 +6	洼小桶 2026-03-18	6/300	2026-03-23 00:29 by king123！
[考研] 317求調(diào)劑 +12	申子申申 2026-03-19	18/900	2026-03-22 22:23 by luoyongfeng
[考研] 0854電子信息求調(diào)劑 +3	α____ 2026-03-22	3/150	2026-03-22 21:28 by zhq0425
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 287求調(diào)劑 +8	晨昏線與星海 2026-03-19	9/450	2026-03-22 17:01 by i_cooler
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 070300化學(xué)求調(diào)劑 +3	苑豆豆 2026-03-20	3/150	2026-03-22 14:48 by ColorlessPI
[考研] 291 求調(diào)劑 +3	化工2026屆畢業(yè)?/a> 2026-03-21	3/150	2026-03-22 14:26 by ColorlessPI
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +4	zxts12 2026-03-21	7/350	2026-03-22 09:57 by zxts12
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 化學(xué)求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-18	6/300	2026-03-21 00:32 by JourneyLucky
[考研] 材料專業(yè)求調(diào)劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan