用气象中的过渡态放在水里，做出的irc就只有三个点，求助是什么原因，需要如何修改？

F范小呆

NH3+H3O+->NH4++H2O  在气象中找到过渡态（虚频是-248），加scrf之后虚频-61，而且IRC只有三个点，改了输入文件调整了stepsize,maxpionts之后irc没有任何变化，想请教这是什么原因，我应该怎么修改输入文件？以下是irc的输入文件，输出文件正常输出了，后附有输出文件和irctu%chk=sNH3+H3O02-NH4+H2O02TS+freq+irc.chk%mem=60mw
%nproc=1
#P MP2/6-31G(d,p) irc(calcall,stepsize=5,maxpoints=200) scrf nosymm

irc

0 2
N1
H2 1 r1
H3 1 r2 2 A1
H4 1 r3 3 A2 2 D1
H5 1 r4 4 A3 3 D2
O6 1 r5 5 A4 4 D3
H7 6 r6 5 A5 4 D4
H8 6 r7 7 A6 5 D5

r1=1.0139103
r2=4.09048058
r3=1.01390058
r4=1.01385923
r5=2.87516579
r6=0.96176238
r7=0.98065858
A1=106.39279735
A2=12.44953486
A3=105.82995262
A4=113.3995576
A5=102.56441348
A6=103.33698573
D1=-89.08607599
D2=153.46477915
D3=124.07104159
D4=62.83692976
D5=16.08818363

谢谢各位前辈百忙之中解决我的问题

Accelerator

频率和IRC的计算水平必须和优化完全相同

F范小呆

Accelerator 发表于 2018-4-16 21:41
频率和IRC的计算水平必须和优化完全相同

谢谢您的回复，从优化到算频率到IRC计算水平都是MP2/6-31G（d,p）

liyuanhe211

F范小呆 发表于 2018-4-16 21:48
谢谢您的回复，从优化到算频率到IRC计算水平都是MP2/6-31G（d,p）

包含用的溶剂模型、用的积分格点等所有影响势能面及势能面上某点导数的因素必须都相同

xingyuan

你好，我能问你一下你的过渡态是怎么算出来的吗？我计算的结果很奇怪，不过虚频也在-460多

zjxitcc

楼主你这直接把气相下的过渡态放到溶剂模型里，不对啊。你是否有看过算freq时四个收敛指标？如果前两个指标不是YES，说明力不为零，你这算的频率就算只有1个虚频，也根本没有意义。