| 2 | 1/1 | 返回列表 |
| 查看: 2725 | 回復: 1 | ||
a_lin916新蟲 (初入文壇)
|
[求助]
求助! QE中如何加范德瓦爾斯力修正
|
|
如題。 QE 中如何加范德瓦爾斯力修正? 看文獻中,vasp是用vdW-DF functional. 但是QE官網(wǎng)上并沒有相關(guān)贗勢。 嘗試過計算mos2層狀結(jié)構(gòu), 設(shè)置input_dft='vdW-DF',但是這種情況下結(jié)構(gòu)優(yōu)化的結(jié)果: 面內(nèi)晶格常數(shù)a b沒有變化, 只有晶格常數(shù)c發(fā)生變化。 具體輸入輸出參數(shù)如下: ###########scf.in################ &control calculation = 'vc-relax' restart_mode='from_scratch', prefix='mos2', tstress = .true. tprnfor = .true. pseudo_dir = '$PSEUDO_DIR', outdir='$TMP_DIR' etot_conv_thr = 1.0D-6 forc_conv_thr = 1.0D-4 / &system ibrav = 4 celldm(1) = 5.9722904441d0 celldm(3) =3.8903303379d0 nat = 6 ntyp = 2 occupations = 'smearing' smearing = 'methfessel-paxton' degauss = 0.005 ecutwfc = 50.0 ecutrho = 200.0 input_dft = 'vdW-DF' / &electrons electron_maxstep = 250 conv_thr = 1.0d-8 mixing_mode = 'plain' mixing_beta = 0.7 mixing_ndim = 8 diagonalization = 'david' diago_david_ndim = 4 / &ions / &cell press_conv_thr = 0.5D0 press = 0.D0 cell_dynamics = 'bfgs' cell_dofree = 'z' / ATOMIC_SPECIES Mo 95.940000d0 Mo.pbe-spn-kjpaw_psl.0.2.UPF S 32.065000d0 S.pbe-n-kjpaw_psl.0.1.UPF ATOMIC_POSITIONS {crystal} Mo 0.3333333333d0 0.6666666667d0 0.2500000000d0 Mo 0.6666666667d0 0.3333333333d0 0.7500000000d0 S 0.3333333333d0 0.6666666667d0 0.6210000000d0 S 0.6666666667d0 0.3333333333d0 0.1210000000d0 S 0.3333333333d0 0.6666666667d0 0.8790000000d0 S 0.6666666667d0 0.3333333333d0 0.3790000000d0 K_POINTS automatic 8 8 3 0 0 0 #########scf.out########## entering subroutine stress ... total stress (Ry/bohr**3) (kbar) P= 44.97 0.00045861 -0.00000000 0.00000000 67.46 -0.00 0.00 -0.00000000 0.00045861 0.00000000 -0.00 67.46 0.00 0.00000000 0.00000000 -0.00000012 0.00 0.00 -0.02 bfgs converged in 9 scf cycles and 8 bfgs steps (criteria: energy < 1.0E-04, force < 1.0E-04, cell < 5.0E-01) End of BFGS Geometry Optimization Final enthalpy = -1005.7174148516 Ry Begin final coordinates new unit-cell volume = 772.34851 a.u.^3 ( 114.45024 Ang^3 ) CELL_PARAMETERS (alat= 5.97229044) 1.000000000 0.000000000 0.000000000 -0.500000000 0.866025404 0.000000000 0.000000000 0.000000000 4.186585056 ATOMIC_POSITIONS (crystal) Mo 0.333333333 0.666666667 0.250000000 Mo 0.666666667 0.333333333 0.750000000 S 0.333333333 0.666666667 0.628788852 S 0.666666667 0.333333333 0.128788852 S 0.333333333 0.666666667 0.871211148 S 0.666666667 0.333333333 0.371211148 End final coordinates |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 290求調(diào)劑 +10 | ADT 2026-03-11 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 材料,紡織,生物(0856、0710),化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 297求調(diào)劑 +4 | 戲精丹丹丹 2026-03-17 | 4/200 |
|
|
[考研] 材料專碩306英一數(shù)二 +9 | z1z2z3879 2026-03-16 | 12/600 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +9 | 負心者當誅 2026-03-11 | 9/450 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 工科,求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 土木第一志愿276求調(diào)劑,科研和技能十分豐富,求新興方向的導師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|