线性响应测定U值，遇到报错“The HP code ...... lda_plus_u=.true.”

Freeman

大家好。

我试图用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
    
46. 
    
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赝势 mn.upf (3.87 MB, 下载次数 Times of downloads: 0)

2024-8-7 12:30 上传 Uploaded

点击下载
Click to download

dojo赝势的原子轨道部分貌似没有label（2P、3D之类的）。上面文件的label是我自己写上去的。

谢谢大佬解答！