| 5 | 1/1 | 返回列表 |
| 查看: 2708 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
didi5158金蟲 (小有名氣)
|
[交流]
【求助】優(yōu)化的結(jié)果,scf計算該結(jié)構(gòu),force為什么仍然很大? 已有7人參與
|
||
|
vc-relax計算晶體的晶格常數(shù),用計算所得到的結(jié)構(gòu)再進行自洽計算,而所得結(jié)果中force,stress,total energy都差別較大,這個是為什么?請大家指教。 vc-relax計算結(jié)果: --------------------------------------------------------------------------------------------------------- total energy = -588.70418379 Ry Harris-Foulkes estimate = -588.70418379 Ry estimated scf accuracy < 3.9E-10 Ry Total force = 0.000201 Total SCF correction = 0.000059 SCF correction compared to forces is too large, reduce conv_thr entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= 0.06 0.00000076 0.00000000 0.00000000 0.11 0.00 0.00 0.00000000 0.00000076 0.00000000 0.00 0.11 0.00 0.00000000 0.00000000 -0.00000024 0.00 0.00 -0.03 bfgs converged in 14 scf cycles and 12 bfgs steps End of BFGS Geometry Optimization Final enthalpy = -588.7041837859 Ry new unit-cell volume = 657.99457 a.u.^3 ( 97.50473 Ang^3 ) CELL_PARAMETERS (alat) 0.781116722 0.000000000 0.000000000 -0.390558361 0.676466924 0.000000000 0.000000000 0.000000000 1.486696131 ATOMIC_POSITIONS (crystal) Li 0.333333333 0.666666667 0.250000000 Li -0.333333333 -0.666666667 -0.250000000 F 0.420499391 -0.010919477 0.122948838 F 0.010919477 0.431418866 0.122948838 F -0.431418866 -0.420499390 0.122948838 F 0.010919477 -0.420499390 0.377051161 F -0.431418866 -0.010919477 0.377051161 F 0.420499390 0.431418866 0.377051161 F -0.420499390 0.010919478 -0.122948838 F -0.010919477 -0.431418866 -0.122948838 F 0.431418866 0.420499390 -0.122948838 F -0.010919477 0.420499390 0.622948838 F 0.431418866 0.010919477 0.622948838 F -0.420499390 -0.431418866 0.622948838 Al 0.666666667 0.333333333 0.250000000 Al -0.666666667 -0.333333333 -0.250000000 Ca 0.000000000 0.000000000 0.000000000 Ca 0.000000000 0.000000000 0.500000000 ---------------------------------------------------------------------------------------------------------- 用上面所得到結(jié)構(gòu),scf計算結(jié)果: ---------------------------------------------------------------------------------------------------------- total energy = -571.49242692 Ry Harris-Foulkes estimate = -571.49242692 Ry estimated scf accuracy < 2.4E-11 Ry Total force = 0.262955 Total SCF correction = 0.000025 entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P=-4722.80 -0.03210576 0.00000000 0.00000000 -4722.92 0.00 0.00 0.00000000 -0.03210576 0.00000000 0.00 -4722.92 0.00 0.00000000 0.00000000 -0.03210340 0.00 0.00 -4722.57 ---------------------------------------------------------------------------------------------------------- |
至尊木蟲 (著名寫手)
| One possibility is that a fixed set of G-vector is selected using the initial cell based on the specified `ecutwfc`. This G-vector set is kept fixed during a vc-relax run (which can modify effectively the stored G-vectors through changes in the unit cell). As such, at the end of vc-relax, the planewave cutoff can differ from direct scf calculation if the unit cell undergoes sizable strain. One solution to this problem is to use the approach by M. Bernasconi et al, [J. Phys. Chem. Solids 56, 501 (1995), doi:10.1016/0022-3697(94)00228-2]. Relevant QE input parameters include: ecfixed, qcutz, and q2sigma see: https://www.quantum-espresso.org/Doc/INPUT_PW.html#idm406 |
金蟲 (著名寫手)
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗: +83 |
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 282求調(diào)劑 +18 | ycy1201 2026-04-01 | 20/1000 |
|
|---|---|---|---|---|
|
[考研] 290求調(diào)劑085701 +3 | 1314捧花 2026-04-02 | 3/150 |
|
|
[考研] 一志愿廈門大學(xué)材料工程專碩354找調(diào)劑。! +7 | 貝唄鋇鋇 2026-03-30 | 7/350 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +13 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 材料考研調(diào)劑 +9 | Gs大王 2026-04-02 | 9/450 |
|
|
[考研] 考研調(diào)劑0702 +3 | 梅子味晚霞 2026-03-28 | 3/150 |
|
|
[考研] 304求調(diào)劑 +12 | 素年祭語 2026-03-31 | 15/750 |
|
|
[考研] 286求調(diào)劑 +5 | lim0922 2026-03-26 | 5/250 |
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
[考研] 086000生物與醫(yī)藥298調(diào)劑求助 +4 | 元元青青 2026-03-31 | 6/300 |
|
|
[考研] 土木304求調(diào)劑 +5 | 頂級擦擦 2026-03-31 | 5/250 |
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
[考研] 江蘇蘇北高校誠邀調(diào)劑同學(xué) +3 | zzll406 2026-03-31 | 3/150 |
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
[考研] 一志愿雙一流機械285分求調(diào)劑 +4 | 幸運的三木 2026-03-29 | 5/250 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|