G09和G16计算NMR结果不一样，连IOp都不一样

trevormx

想算一个化合物的NICS，同时用G09和G16计算，初始gjf文件一模一样
# nmr=giao scrf=(smd,solvent=dichloromethane) def2tzvp m062x

但是计算的时候G09和G16的计算出来的NMR位移值不太一致，影响到小数点后一位了。查看输出文件，发现两个IOp都不一样。
G09的是
1/38=1/1;
2/12=2,17=6,18=5,40=1/2;
3/5=44,7=101,11=2,16=1,25=1,30=1,70=32201,72=9,74=-55/1,2,8,3;
4//1;
5/5=2,38=5,53=9/2;
8/6=1,10=90,11=11/1;
10/13=100,45=16/2;
6/7=2,8=2,9=2,10=2,28=1/1;
99/9=1/99;


G16的是
1/38=1,172=1/1;
2/12=2,17=6,18=5,40=1/2;
3/5=44,7=101,11=2,25=1,30=1,70=32201,72=9,74=-55/1,2,3,8;
4//1;
5/5=2,38=5,53=9/2;
8/6=1,10=90,11=11/1;
10/13=100,45=16/2;
6/7=2,8=2,9=2,10=2,28=1/1;
99/9=1/99;


不知道是怎么回事，是不是IOp不一样影响了输出结果，各位高手帮忙分析分析啊。

sobereva

补充一句，不仅是int设定不一样，g09和g16在溶剂模型实现上的细节也不同。g16默认条件下的结果原理上比g09更准，非要和g09对比就加上g09default

zjxitcc

trevormx 发表于 2019-9-16 12:47
查了查IOp表，这几个不同貌似不影响结果的，应该是和你说的是精度影响的吧，不过差的值也是挺大的。谢谢 ...

见高斯官网http://gaussian.com/g09def/
如果不想写int(acc2e=10,fine)，直接写G09Defaults更省事。

trevormx

zjxitcc 发表于 2019-9-16 12:42
“初始gjf文件一模一样”这句话就是原因啊。众所周知的区别是，G09默认int(acc2e=10,fine)，而G16默认是int ...

查了查IOp表，这几个不同貌似不影响结果的，应该是和你说的是精度影响的吧，不过差的值也是挺大的。谢谢啦！

zjxitcc

“初始gjf文件一模一样”这句话就是原因啊。众所周知的区别是，G09默认int(acc2e=10,fine)，而G16默认是int(acc2e=12,ultrafine)。acc2e=10或12一般不会影响结果，但是fine与ultrafine差别还是挺大的。