近红外计算求助，针对3-quanta的应该怎么做？

haos423

想请教一下各位老师，纯有机分子红外/振动分析的计算，针对3-quanta的应该怎么做？Gaussian和ORCA好像都只能做到基频和2-quanta的非谐振计算；

对3-quanta的计算，免费的软件/方法我只找到nwchem的VSCF，参考软件的手册(https://nwchemgit.github.io/VSCF.html)，以水分子为例，附件中的输入文件能计算3倍泛频，但似乎缺少和频的部分

我的输入文件编写得正确吗？或者有其他方法能完整计算3-quanta IR谱（包括所有overtone和combination bands）？谢谢！

(以下是输出文件末尾部分的节选)

  RESULTS OF VIBRATIONAL SCF CALCULATION: Frequencies, cm-1
  Mode      Harmonic     Diagonal      VSCF       PT2-VSCF
   1       11813.18     12371.95     10966.29     10873.87
   2       11489.95     10997.81     10899.68     10668.23
   3        4952.31      4867.02      4706.89      4613.12
......Finished vibrational SCF......


  IR INTENSITIES CALCULATED USING DIPOLE MOMENT:
  mode   frequency, cm-1  intensity, km/mol
   1       10873.87          1.52
   2       10668.23          0.01
   3        4613.12          0.01


  (Note: SCF dipole moments used to calculate intensities)
......Finished calculating IR intensities......

bxyao

我在之前也算过这个，但是忘记具体的关键词了，服务器宕机一次，把我的数据搞丢了，隐约记得是翻看高斯手册，然后这个Spectro=MaxQuanta=quanta: Compute transition integrals to states with up to the specified quanta. The default is 2.，freq还要加个readxxx（具体忘记了）好像，才能计算额外的非谐振部分，你可以试试

haos423

bxyao 发表于 2026-1-30 09:51
我在之前也算过这个，但是忘记具体的关键词了，服务器宕机一次，把我的数据搞丢了，隐约记得是翻看高斯手册 ...

非常感谢老师的指点！readanharm加上您说的Spectro=MaxQuanta=3，确实这样就可以得到了