| 2 | 1/1 | 返回列表 |
| 查看: 573 | 回復: 1 | |||
[交流]
【求助】晶格優(yōu)化崩潰?
|
|
各位朋友。本人馳豫原子位置的時候,坐標一直變大?偏離實驗數(shù)據(jù)跟文獻上的計算數(shù)據(jù)很大,麻煩大家給我看看,我的輸入文件有什么不妥的地方,萬分感激。。。 以下是我的輸入文件內(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/89
+1/85
+1/74
+1/63
+2/58
+1/49
+1/39
+1/17
+1/16
+1/5
+1/5
+1/5
+1/5
+1/4
+1/2
+1/2
+1/1
+1/1
+1/1
+1/1
|
以下是輸出: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工專碩調(diào)劑 +6 | heming3743 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +4 | 晨昏線與星海 2026-03-19 | 5/250 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +7 | Doleres 2026-03-19 | 7/350 |
|
|
[考研] 一志愿 南京航空航天大學大學 ,080500材料科學與工程學碩 +4 | @taotao 2026-03-20 | 4/200 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 0703化學調(diào)劑 ,六級已過,有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 一志愿福大288有機化學,求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 070300化學319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 復試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|