| 查看: 2674 | 回复: 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 |
回复此楼» 猜你喜欢
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有9人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有5人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
自荐读博
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
2楼2021-04-13 15:57:35