qnclqq

木蟲 (正式寫手)

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

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

對TiO2進(jìn)行+U計算，晶格弛豫（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ù)此樓

» 猜你喜歡

一志愿070300浙大化學(xué)358分，求調(diào)劑！已經(jīng)有4人回復(fù)
一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑已經(jīng)有4人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有5人回復(fù)
323求調(diào)劑已經(jīng)有6人回復(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ù)

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小時
蟲號: 1419918
注冊: 2011-09-27
性別: GG
專業(yè): 理論和計算化學(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小時
蟲號: 1835611
注冊: 2012-05-26
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

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

你確定你的計算中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小時
蟲號: 1464347
注冊: 2011-10-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

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

自旋極化只是加上看看效果，如果沒有極化并不會影響最后的效果。問題在于+U之后無法實現(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小時
蟲號: 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 感覺沒有這么多

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

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

Hubbard_U(1)=1, Hubbard_U(2)=2，這樣設(shè)置，我的本意是試一下給O加1eV，給Ti加2eV。請問錯在哪里，應(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小時
蟲號: 1419918
注冊: 2011-09-27
性別: GG
專業(yè): 理論和計算化學(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小時
蟲號: 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，是想看看效果。因為TiO2是電荷轉(zhuǎn)移絕緣體，價帶是氧的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小時
蟲號: 1419918
注冊: 2011-09-27
性別: GG
專業(yè): 理論和計算化學(xué)

引用回帖:

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

先不加U，relax一下，得到坐標(biāo)。然后把這個優(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小時
蟲號: 2396394
注冊: 2013-04-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

我用的別的一個軟件也是，加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小時
蟲號: 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)。然后把這個優(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！
[考研] 招08考數(shù)學(xué) +4	laoshidan 2026-03-20	8/400	2026-03-22 19:56 by 小皮蛋醬
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 化學(xué)求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-21 00:39 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11	adaie 2026-03-19	11/550	2026-03-20 21:30 by laoshidan
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong