| 查看: 2943 | 回复: 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 |
回复此楼» 猜你喜欢
国自然评审意见
已经有12人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有12人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
-
2026年叶企孙基金
已经有6人回复
2楼2021-04-13 15:57:35
30