计算化学公社

标题: 求助QE模拟红外光谱出错 [打印本页]
作者
Author: qczgzly    时间: 2023-2-6 07:49
标题: 求助QE模拟红外光谱出错
    我有一个氢氧化铜吸附有机小分子体系，首先用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:

  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

  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

  19. ……


  22. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  23.      Error in routine phq_readin (1):
  24.      no elec. field with metals
  25. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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

作者
Author: qczgzly    时间: 2023-2-13 14:21
自我解答如下：
首先对于金属元素，epsil字段不能设置为.true.，否则出现如上报错。其次设置epsil=.false.后，用ph.x计算声子未得到dielec相关信息及输出的fildyn文件，原因为声子计算没达到收敛阈值。在pw.x计算中提高收敛精度，例如设置conver_thr=1.0d-10，然后在ph.x计算中设置tr2_ph=1.0d-08，则成功完成计算，得到预期信息。
作者
Author: dr_tao    时间: 2023-6-28 09:00
请问楼主现在还从事IR模拟方面的工作吗？有没有兴趣合作课题？



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