QuantumEspresso拉曼光谱计算 ph.in运行报错The Hessian file: ./tmp/CBP-Alph...

Jerry_Wang

各位老师好，我在用QuantumEspresso计算一个分子晶体的拉曼光谱，在运行ph.in文件时，遇到如下报错，请问该如何解决呢？     



Error in routine d2ionq_dispd3 (1):
The Hessian file: ./tmp/CBP-Alpha.hess is missing.


以下为ph.in的输入文件
&INPUTPH
  tr2_ph = 1.0d-14
  prefix = 'CBP-Alpha'
  outdir = './tmp'
  fildyn = 'CBP-Alpha'
  ldisp = .false.
  epsil = .true.
  trans = .true.
  qplot = .false.
  nq1 = 1, nq2 = 1, nq3 = 1
/
0.0 0.0 0.0

以下为上一步scf计算的pw.in输入文件
&CONTROL
  title='CBP-Alpha-scf',
  calculation='scf',
  prefix = 'CBP-Alpha'
  pseudo_dir='./',
  outdir='./tmp',
  verbosity='high',
/
&SYSTEM
  ibrav= 0, nat= 124, ntyp= 3,
  occupations = 'fixed',
  ecutwfc = 80,
  ecutrho = 800,
  vdw_corr = 'grimme-d3',
/
&ELECTRONS
  conv_thr = 1.0d-10
  mixing_beta = 0.5d0
/
&IONS
/
&CELL
/
ATOMIC_SPECIES
  C 12.0107 C.pbe-n-rrkjus_psl.1.0.0.UPF
  H 1.00794 H.pbe-rrkjus_psl.1.0.0.UPF
  N 14.00674 N.pbe-n-rrkjus_psl.1.0.0.UPF
CELL_PARAMETERS (angstrom)
   8.459549526   0.000000000  -0.013055016
   0.000000000  15.662496878   0.000000000
  -3.735405692   0.000000000   9.404073167
ATOMIC_POSITIONS (crystal)
C                0.9391414775        0.6511901948        0.0001336347
C                0.3784976529        0.9006437189        0.0155178630
H                0.7757116114        0.2124087035        0.0187583199
H                1.0019841321        0.5897111019        0.0306536179
C                0.9206504794        0.8052563194        0.0383402022
C                0.5506940455        0.0355539285        0.0435877664
C                0.7189543912        0.1652396522        0.0665723707
H                0.6325579656        0.4887538700        0.0548401428
H                0.8993127218        0.9906024777        0.0795285774
H                0.6023302268        0.3313104779        0.0765492092
C                0.0024835333        0.7250014214        0.0795842995
C                0.2261448565        0.1888536358        0.0848644218
H                0.3951836956        0.9024270494        0.1273945662
C                0.2063954886        0.3409479726        0.1226546038
H                0.2897961331        0.1277071056        0.1175317552
H                0.3012245132        0.6118530530        0.1146420810
C                0.0201530199        0.8672945206        0.1411095687
C                0.7292089644        0.4535095030        0.1395458206
H                0.2897961331        0.3722928944        0.6175317552
H                0.3012245132        0.8881469470        0.6146420810
C                0.0201530199        0.6327054794        0.6411095687
C                0.7292089644        0.0464904970        0.6395458206
C                0.0060365193        0.5444443856        0.6545881225
C                0.7123141847        0.1346406675        0.6496738235
C                0.2885512487        0.2379342436        0.6651907325
H                0.2546885844        0.1027339140        0.6882598666
C                0.5811894619        0.4575321578        0.6897656561
N                0.1472689194        0.7637812090        0.7040610279
C                0.7493496897        0.3299254200        0.7119607479
H                0.5258479272        0.5034909736        0.7409130304
C                0.3216422077        0.8918055436        0.7271974022
H                0.8773951383        0.9330155564        0.7239177372
C                0.1596653066        0.6772655983        0.7411508941
C                0.8685822085        0.0017679097        0.7353366102
H                0.4011250486        0.2417865996        0.7628008744
H                0.5988749514        0.7417865996        0.7371991256
C                0.1314177915        0.5017679097        0.7646633898
C                0.8403346934        0.1772655983        0.7588491059
H                0.1226048617        0.4330155564        0.7760822628
C                0.6783577923        0.3918055436        0.7728025978
H                0.4741520728        0.0034909736        0.7590869696
C                0.2506503103        0.8299254200        0.7880392521
N                0.8527310806        0.2637812090        0.7959389721
C                0.4188105381        0.9575321578        0.8102343439
H                0.7453114156        0.6027339140        0.8117401334
C                0.7114487513        0.7379342436        0.8348092675
C                0.2876858153        0.6346406675        0.8503261765
C                0.9939634807        0.0444443856        0.8454118775
C                0.2707910356        0.5464904970        0.8604541794
C                0.9798469801        0.1327054794        0.8588904313
H                0.6987754868        0.3881469470        0.8853579190
H                0.7102038669        0.8722928944        0.8824682448
C                0.7936045114        0.6590520274        0.8773453962
H                0.6048163044        0.0975729506        0.8726054338
C                0.7738551435        0.8111463642        0.9151355782
C                0.9975164667        0.2749985786        0.9204157005
H                0.3976697732        0.6686895221        0.9234507908
H                0.1006872782        0.0093975223        0.9204714226
H                0.3674420344        0.5112461300        0.9451598572
C                0.2810456088        0.8347603478        0.9334276293
C                0.4493059545        0.9644460715        0.9564122336
C                0.0793495206        0.1947436806        0.9616597978
H               -0.0019841321        0.4102888981        0.9693463821
H                0.2242883886        0.7875912965        0.9812416801
C                0.6215023471        0.0993562811        0.9844821370
C                0.0608585225        0.3488098052        0.9998663653
K_POINTS automatic
1 1 1 0 0 0

丁越

QE计算拉曼光谱时必须使用NC赝势并结合LDA泛函，你可以考虑使用ONCV的赝势。

另外ph.x的输入文件中需要添加lraman=.true.才能进行拉曼计算。

Jerry_Wang

丁越 发表于 2026-1-6 10:13
QE计算拉曼光谱时必须使用NC赝势并结合LDA泛函，你可以考虑使用ONCV的赝势。

另外ph.x的输入文件中需要 ...

多谢老师的解答！
我想再问一下， 我目前出现的报错，是不是正是因为您说的这几个问题引起的呢？
Error in routine d2ionq_dispd3 (1):
The Hessian file: ./tmp/CBP-Alpha.hess is missing.

丁越

Jerry_Wang 发表于 2026-1-6 22:26
多谢老师的解答！
我想再问一下， 我目前出现的报错，是不是正是因为您说的这几个问题引起的呢？
Error ...

应该是这个dftd3_hess关键词需要读取hessian矩阵，而你的相关文件正好缺失了。
File where the D3 dispersion hessian matrix is read. Set to
'automatic.hess' to enable automatic mode (experimental). In
this mode, D3 Hessian is computed if 'automatic.hess' file is
missing.

Jerry_Wang

丁越 发表于 2026-1-7 09:28
应该是这个dftd3_hess关键词需要读取hessian矩阵，而你的相关文件正好缺失了。
File where the D3 dispe ...

感谢老师的回复！！！
我想接着请教一下，这个文件缺失可能是由于什么原因引起的呢？

Jerry_Wang

丁越 发表于 2026-1-7 09:28
应该是这个dftd3_hess关键词需要读取hessian矩阵，而你的相关文件正好缺失了。
File where the D3 dispe ...

老师，您好！
我在vc-relax和scf 这两步计算中都写了vdw_corr = 'grimme-d3',
到了ph这步，输入文件是没有写vdw_corr = 'grimme-d3'的，那D3 dispersion hessian matrix应该怎么产生呢？需要我额外再写什么关键词吗？

丁越

Jerry_Wang 发表于 2026-1-8 05:33
老师，您好！
我在vc-relax和scf 这两步计算中都写了vdw_corr = 'grimme-d3',
到了ph这步，输入文件是 ...

这不是解决你问题的重点。

QE计算拉曼光谱坑挺多的，首先是之前说的你得用NC赝势，并且要使用LDA泛函。然后关于色散校正这块LDA泛函不支持DFT-D3,所以你也得换成DFT-D2。

另外ONCV赝势推荐使用修改版的（https://github.com/pipidog/ONCV*SP/tree/master），修改后的赝势与QE的各种计算更加兼容了。

Exodus

缺少的hess文件可以使用d3hess.x计算获得，DFT-D2则需不要

Jerry_Wang

丁越 发表于 2026-1-8 12:24
这不是解决你问题的重点。

QE计算拉曼光谱坑挺多的，首先是之前说的你得用NC赝势，并且要使用LDA泛函 ...

原来如此，谢谢老师的解答！感谢！
不过您发给我的网址我打开不了，显示404了，您能重新给我发一下吗？
我自己也在网上找了一下，http://quantum-simulation.org/potentials/sg15_oncv/
但是从这个网址里的下载的全是PBE的，请问哪里可以直接下载到LDA泛函的呢？

Jerry_Wang

Exodus 发表于 2026-1-8 12:38
缺少的hess文件可以使用d3hess.x计算获得，DFT-D2则需不要

我试了一下，确实可以！多谢啦！

丁越

Jerry_Wang 发表于 2026-1-9 06:21
原来如此，谢谢老师的解答！感谢！
不过您发给我的网址我打开不了，显示404了，您能重新给我发一下吗？ ...

应该是论坛的词屏蔽规则把*位置的P给屏蔽了,你粘贴URL时修改一下

1. https://github.com/pipidog/ONCV*SP/tree/master

复制代码

你给出的链接是适用于QE的原始ONCV赝势。用input_dft关键词指定理论方法，比如input_dft='pz'使用LDA泛函。

Jerry_Wang

丁越 发表于 2026-1-9 11:43
应该是论坛的词屏蔽规则把*位置的P给屏蔽了,你粘贴URL时修改一下
你给出的链接是适用于QE的原始ONCV赝势 ...

老师，您好！感谢您的回复！
您给的链接，我把*号换成P，就可以打开了。
但是打开后，abinit 和 sg15 这两个文件夹里的赝势，文件名全是带有pbe的，还是没找到LDA的，
是我找的地方不对吗？想请教一下该怎么找呢？

丁越

Jerry_Wang 发表于 2026-1-10 04:42
老师，您好！感谢您的回复！
您给的链接，我把*号换成P，就可以打开了。
但是打开后，abinit 和 sg15  ...

就是只有PBE的赝势，没有专门为LDA泛函搞一套赝势。直接借用PBE的赝势计算就行。

Jerry_Wang

丁越 发表于 2026-1-10 11:32
就是只有PBE的赝势，没有专门为LDA泛函搞一套赝势。直接借用PBE的赝势计算就行。

好嘞，谢谢老师！您的解答对我太有帮助了！感谢感谢！
我还想再请教一下，您上一条回复说的 “用input_dft关键词指定理论方法，比如input_dft='pz'使用LDA泛函”。是指我在做 pw.x 计算的时候，加上这个关键词吗？相当于就能用上LDA泛函了？

丁越

Jerry_Wang 发表于 2026-1-10 11:46
好嘞，谢谢老师！您的解答对我太有帮助了！感谢感谢！
我还想再请教一下，您上一条回复说的 “用input_d ...

对，这个关键词就是指定计算使用的泛函的，默认是直接从赝势文件读取。