| 查看: 2940 | 回复: 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 |
回复此楼» 猜你喜欢
今天务委会开完了，明天出结果吗
已经有21人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有14人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
2楼2021-04-13 15:57:35
28