| 查看: 2915 | 回复: 1 | ||
[求助]
求助！ QE中如何加范德瓦尔斯力修正
|
|
如题。 QE 中如何加范德瓦尔斯力修正？ 看文献中，vasp是用vdW-DF functional. 但是QE官网上并没有相关赝势。 尝试过计算mos2层状结构， 设置input_dft='vdW-DF',但是这种情况下结构优化的结果： 面内晶格常数a b没有变化， 只有晶格常数c发生变化。 具体输入输出参数如下： ###########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 |
回复此楼» 猜你喜欢
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有13人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有13人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
2楼2021-04-13 15:57:35
50