计算化学公社

 找回密码 Forget password
 注册 Register

求助,还是非简谐频率计算,尝试多种方法都不行,是结构的问题吗

查看数: 913 | 评论数: 9 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-5-19 21:23

正文摘要:

本帖最后由 xuejunfei1987 于 2023-5-21 16:14 编辑 最近计算一种SH和NH3的复合物,老师要非简谐频率的结果,简谐频率是正常的,但是非简谐计算结果一直算不对,一开始scf不收敛,加了calcall之后正常结束,但是 ...

回复 Reply

xuejunfei1987 发表于 Post on 2023-5-24 09:49:20
wzkchem5 发表于 2023-5-21 17:17
由输出文件显然可以看出,问题在于你调用的是闭壳层的MP2。应该调用开壳层MP2,调用方法见说明书

老师, molpro我不熟悉,看了手册也没找到调用开壳层mp2的方法,您能给我说一下吗  
wzkchem5 发表于 Post on 2023-5-21 17:17:42
xuejunfei1987 发表于 2023-5-21 09:14
好的 我已经上传了   NH3-HS.out

由输出文件显然可以看出,问题在于你调用的是闭壳层的MP2。应该调用开壳层MP2,调用方法见说明书
xuejunfei1987 发表于 Post on 2023-5-21 16:14:36
wzkchem5 发表于 2023-5-21 05:39
上传完整输出文件

好的 我已经上传了   NH3-HS.out
wzkchem5 发表于 Post on 2023-5-21 05:39:01
xuejunfei1987 发表于 2023-5-20 14:50
我提交了一个下面的计算,一直说
NUMBER OF OCCUPIED ORBITALS NOT CONSISTENT WITH NUMBER OF ELECTRO ...

上传完整输出文件
xuejunfei1987 发表于 Post on 2023-5-20 21:50:51
wzkchem5 发表于 2023-5-20 15:47
http://sobereva.com/203

我提交了一个下面的计算,一直说
NUMBER OF OCCUPIED ORBITALS NOT CONSISTENT WITH NUMBER OF ELECTRONS FOR CCSD,不知道改哪里,能请教一下吗

**,HS
memory,4000,m
geometry
  N        -2.126910000000      0.000149000000     -0.000161000000
  H        -2.502044000000     -0.681266000000     -0.647022000000
  H        -2.502915000000      0.901189000000     -0.265086000000
  H        -2.501600000000     -0.221098000000      0.913177000000
  S         1.397022000000      0.000007000000     -0.000016000000
  H         0.042575000000      0.000009000000      0.000316000000
end

basis=vtz
logfile,scratch

{uhf,
occ,14
wf,27,1,1}

uhf
mp2
optg
{frequencies,symm=auto
print,low=50}

label1
{uhf
start,atden}
{mp2
cphf,1}

surf,start1D=label1,dipole=1,sym=auto
vscf,dipole=1,combi=1
vci,dipole=1,version=3,combi=1
xuejunfei1987 发表于 Post on 2023-5-20 18:19:55
wzkchem5 发表于 2023-5-20 15:47
http://sobereva.com/203

谢谢,能不能提供一个molpro提交vscf计算的模版啊  
wzkchem5 发表于 Post on 2023-5-20 15:47:13
xuejunfei1987 发表于 2023-5-20 06:22
谢谢,您说的这些方法我不会用啊,您能不能帮我提交个计算呀

http://sobereva.com/203
xuejunfei1987 发表于 Post on 2023-5-20 13:22:13
wzkchem5 发表于 2023-5-19 23:32
这种非谐性效应很大的体系,高斯的基于VPT2的方法可能不适用,可以考虑用VSCF、VCI,乃至路径积分分子动力 ...

谢谢,您说的这些方法我不会用啊,您能不能帮我提交个计算呀
wzkchem5 发表于 Post on 2023-5-19 23:32:14
这种非谐性效应很大的体系,高斯的基于VPT2的方法可能不适用,可以考虑用VSCF、VCI,乃至路径积分分子动力学的方法计算

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 06:37 , Processed in 0.173021 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list