| 2 | 1/1 | 返回列表 |
| 查看: 572 | 回復(fù): 1 | |||
[交流]
【求助】晶格優(yōu)化崩潰?
|
|
各位朋友。本人馳豫原子位置的時(shí)候,坐標(biāo)一直變大?偏離實(shí)驗(yàn)數(shù)據(jù)跟文獻(xiàn)上的計(jì)算數(shù)據(jù)很大,麻煩大家給我看看,我的輸入文件有什么不妥的地方,萬(wàn)分感激。。。 以下是我的輸入文件內(nèi)容: &control calculation='vc-relax', restart_mode='from_scratch', prefix='Ba' tstress = .true. tprnfor = .true. pseudo_dir = './', outdir='./' nstep = 105 , etot_conv_thr = 1.0E-6 , forc_conv_thr = 1.0D-5 , dt = 50 , / &system ibrav = 0, nat = 5, ntyp = 3, ecutwfc = 40.0 , ecutrho = 400.0 , nosym = .false. , nbnd = 80, occupations = 'smearing' , degauss = 0.05, smearing = 'marzari-vanderbilt' , nspin = 1 , / &electrons mixing_beta = 0.7 conv_thr = 1.0d-8 / &IONS / &CELL cell_dynamics = 'bfgs' press = 1000 / 。。。 。。。 。。。 K_POINTS automatic 8 8 8 0 0 0 |
» 搶金幣啦!回帖就可以得到:
+1/241
+1/90
+1/74
+1/69
+1/63
+2/56
+1/52
+1/39
+1/36
+2/32
+1/32
+1/17
+1/16
+1/7
+1/7
+1/6
+1/5
+1/5
+1/5
+1/4
|
以下是輸出: Forces acting on atoms (Ry/au): atom 1 type 1 force = 0.00000000 0.00000000 0.00000000 atom 2 type 2 force = -0.00076425 0.00076425 0.00000000 atom 3 type 2 force = 0.00076425 -0.00076425 0.00000000 atom 4 type 3 force = -0.00079442 -0.00079442 0.43102097 atom 5 type 3 force = 0.00079442 0.00079442 -0.43102097 Total force = 0.000000 Total SCF correction = 0.000000 entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= 989.12 0.00670743 0.00000026 -0.00001451 986.70 0.04 -2.13 0.00000026 0.00670743 -0.00001451 0.04 986.70 -2.13 -0.00001451 -0.00001451 0.00675680 -2.13 -2.13 993.96 number of scf cycles = 11 number of bfgs steps = 10 enthalpy old = -666.0630713336 Ry enthalpy new = -666.0638799665 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.0033684758 bohr new conv_thr = 0.0000000010 Ry new unit-cell volume = 379.17935 a.u.^3 ( 56.18858 Ang^3 ) CELL_PARAMETERS (alat= 3.96250000) 1.627610836 -0.000000404 0.000086123 -0.000000404 1.627610836 0.000086123 0.813717998 0.813717998 2.300657444 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +4 | Charlieyq 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +4 | 西南交通專材355 2026-03-19 | 4/200 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 329求調(diào)劑 +6 | 想上學(xué)吖吖 2026-03-19 | 6/300 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|