求助QE模拟红外光谱出错

qczgzly

    我有一个氢氧化铜吸附有机小分子体系，首先用pw.x优化结构，得到正常收敛的结果。然后用ph.x计算声子，也能正常结束，但得不到dielec相关信息，通过fildyn关键词输出的文件中也没有任何信息。用于ph.x计算的输入文件如下：

1. Normal modes for Cu001-Val2
   
2. &inputph
   
3.   tr2_ph= 1.0d-12
   
4.   prefix= 'Cu001-Val2'
   
5.   epsil= .false.
   
6.   trans= .true.
   
7.   asr= .true.
   
8.   outdir= './'
   
9.   fildyn= 'Cu001-Val2_IR.dynG'
   
10. /
    
11. 0.0 0.0 0.0

复制代码

经排查，如果设置epsil= .true. ，则产生如下信息：

1.      Reading xml data from directory:
   
2. 
   
3.      ./Cu001-Val2.save/
   
4.      Message from routine read_pp_mesh:
   
5.      mesh size missing, using the one in header
   
6.      Message from routine read_pp_mesh:
   
7.      mesh size missing, using the one in header
   
8.      Message from routine read_pp_mesh:
   
9.      mesh size missing, using the one in header
   
10.      file O.pbe-n-kjpaw_psl.0.1.UPF: wavefunction(s)  2P renormalized
    
11.      file C.pbe-n-kjpaw_psl.1.0.0.UPF: wavefunction(s)  2S 2P renormalized
    
12. 
    
13.      IMPORTANT: XC functional enforced from input :
    
14.      Exchange-correlation= PBE
    
15.                            (   1   4   3   4   0   0   0)
    
16.      Any further DFT definition will be discarded
    
17.      Please, verify this is what you really want
    
18. 
    
19. ……
    
20. 
    
21. 
    
22. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
23.      Error in routine phq_readin (1):
    
24.      no elec. field with metals
    
25. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
26. 
    
27.      stopping ...

复制代码

想在此求助对于含金属元素的有机无机杂化体系，应如何计算其红外光谱？谢谢！

qczgzly

自我解答如下：
首先对于金属元素，epsil字段不能设置为.true.，否则出现如上报错。其次设置epsil=.false.后，用ph.x计算声子未得到dielec相关信息及输出的fildyn文件，原因为声子计算没达到收敛阈值。在pw.x计算中提高收敛精度，例如设置conver_thr=1.0d-10，然后在ph.x计算中设置tr2_ph=1.0d-08，则成功完成计算，得到预期信息。

dr_tao

请问楼主现在还从事IR模拟方面的工作吗？有没有兴趣合作课题？