计算化学公社

标题: 线性响应测定U值，遇到报错“The HP code ...... lda_plus_u=.true.” [打印本页]
作者
Author: Freeman    时间: 2024-8-7 12:25
标题: 线性响应测定U值，遇到报错“The HP code ...... lda_plus_u=.true.”
本帖最后由 Freeman 于 2024-8-7 12:36 编辑

大家好。

我试图用hp.x计算U值，pw和hp的输入文件如下（请不要在意我给Ni和Zn都用了Pb的赝势，那是故意的。）

  1. &CONTROL
  2.   calculation = 'scf'
  3.   pseudo_dir = '../paw'
  4.   prefix='mn'
  5.   max_seconds = 86000
  6.   nstep = 75
  7.   verbosity = 'high'
  8. /
  9. &SYSTEM
  10.   ecutwfc =   40
  11.   ibrav = 1
  12.   A =  17.376089024
  13.   nat = 135
  14.   ntyp = 7
  15.   nspin = 2
  16.   !tot_charge = 5
  17.   tot_magnetization = 5
  18.   !occupations = 'smearing'
  19.   degauss = 0.01
  20. /
  21. &ELECTRONS
  22.   !startingpot = 'file'
  23.   startingwfc = 'file'
  24.   !electron_maxstep = 200
  25.   conv_thr = 1.0d-10
  26. /
  27. &IONS
  28. /
  29. &CELL
  30. cell_dofree = 'ibrav'
  31. /
  32. ATOMIC_SPECIES
  33. Mn 1 mn.upf
  34. Ni 1 pb.upf
  35. Zn 1 pb.upf
  36. Pb 1 pb.upf
  37. Cl 1 cl.upf
  38. Cs 1 cs.upf
  39. Br 1 br.upf
  40. ATOMIC_POSITIONS angstrom
  41. ......
  42. K_POINTS automatic
  43. 2 2 2 0 0 0
  44. HUBBARD (ortho-atomic)
  45. U Mn-3d 1.d-8

  47. &inputhp
  48.   prefix = 'mn_u'
  49.   nq1 = 2
  50.   nq2 = 2
  51.   nq3 = 2
  52.   perturb_only_atom(1) = .true.
  53.   iverbosity = 2
  54. /
复制代码
用到hp.x时，遇到了这个报错“The HP code can be used only when lda_plus_u=.true.”。


猜测的原因：

我现在用的是dojo的paw，而dojo的paw文件的特点是，他不光给了价电子的原子轨道，还给了空轨道，比如mn.upf不仅有半占据的3d轨道，还有空的4d轨道。我怀疑这就是问题所在。这个贴子（https://lists.quantum-espresso.o ... 22-July/049197.html）提到说“However, the hp.x currently supports the calculations of U and V for one Hubbard channel per atomic type.”，也就是说hp.x恐怕不能处理一个原子两个d轨道的情况吧？

猜测的解决方案：
我可否偷偷删掉mn.upf的4d轨道信息，或者将4d的角动量从2写成1，让hp.x识别不到？


QE版本7.1
Mn赝势 (, 下载次数 Times of downloads: 0)
dojo赝势的原子轨道部分貌似没有label（2P、3D之类的）。上面文件的label是我自己写上去的。

谢谢大佬解答！






欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3