【求助】采用溶剂模型是否要使用弥散函数(急)
在考虑有溶剂的情况时,进行计算选择基组的时候是否要采用弥散函数?《Exploring Chemistry with Electronic Structure Methods》书上进行溶剂计算的时候,怎么都带弥散函数的基组?
恳请虫友答疑。
[ Last edited by majun04 on 2008-7-14 at 14:59 ]
返回小木虫查看更多
今日热帖
在考虑有溶剂的情况时,进行计算选择基组的时候是否要采用弥散函数?《Exploring Chemistry with Electronic Structure Methods》书上进行溶剂计算的时候,怎么都带弥散函数的基组?
恳请虫友答疑。
[ Last edited by majun04 on 2008-7-14 at 14:59 ]
返回小木虫查看更多
1、在计算溶剂的情况下,选用基组的时候是不是一定要用弥散函数?用了弥散函数是不是一定要用参数SCF=Tight?
这个问题,你先要搞清楚概念,什么是弥散基组:
弥散基组
弥散基组是对劈裂价键基组的另一种扩大。在高斯函数中,变量α对函数形态有极大的作用,当α的取值很大时,函数图像会向原点附近聚集,而当α取值很小的时候,函数的图像会向着远离原点的方向弥散,这种α很小的高斯函数被称为弥散函数。所谓弥散基组就是在劈裂价键基组的基础上添加了弥散函数的基组,这样的基组可以用于非键相互作用体系的计算。
所以要看你计算的体系,如果存在非键相互作用体系,最好使用加弥散函数,这样计算出来的结果会比没加弥散函数的结果准确一些。
附:劈裂价键基组
根据量子化学理论,基组规模越大,量化计算的精度就越高,当基组规模趋于无限大时,量化计算的结果也就逼近真实值,为了提高量子化学计算精度,需要加大基组的规模,即增加基组中基函数的数量,增大基组规模的一个方法是劈裂原子轨道,即使用多于一个基函数来表示一个原子轨道。
劈裂价键基组就是应用上述方法构造的较大型基组,所谓劈裂价键就是将价层电子的原子轨道用两个或以上基函数来表示。常见的劈裂价键基组有3-21G、4-21G、4-31G、6-31G、6-311G等,在这些表示中前一个数字用来表示构成内层电子原子轨道的高斯型函数数目,“-”以后的数字表示构成价层电子原子轨道的高斯型函数数目。如6-31G所代表的基组,每个内层电子轨道是由6个高斯型函数线性组合而成,每个价层电子轨道则会被劈裂成两个基函数,分别由3个和1个高斯型函数线性组合而成。
劈裂价键基组能够比STO-3G基组更好地描述体系波函数,同时计算量也比最小基组有显著的上升需要根据研究的体系不同而选择相应的基组进行计算。
SCF=Tight 是指示使用完全的收敛条件。一般情况下,SCF 和 DFT计算若使用包含扩散函数的基底函数应该使用关键词SCF= Tight要求使用较严格的 SCF 收敛条件。
2、溶剂的介电常数问题,我发现手册上给出的溶剂介电常数与《兰氏化学手册》查到的溶剂的介电常数有差别,请问相信哪一个?
GS手册里给出的可能是计算值,兰氏给出的都是实验值。如果是作实验比较,建议使用兰氏,如果是计算比较,建议使用GS推荐的。
3、还有就是在计算溶剂的影响问题时候,在一样的溶质条件下,是不是单独的更改介电常数值(前提是这个溶剂在gaussian手册中给出了)就可以了?
是
只要感变介电常数,就代表改变溶剂。
4、gaussian手册中没有提供介电常数得溶剂,如果我想算,它得基本参数应该到哪里去找才可靠些?
这个就不好说了,如果没有,就只能近似使用兰氏,或者是其他的一些文献资料,这些只能自己查文献了。,
非常感谢!
问题3,还想问一下:如果不是gaussian中默认的溶剂,是否能仅仅更改介电常数就可以了?还是要有其他的改动?
我是想和实验相比较,打算采用兰氏的数据。在设置输入文件的时候,是否是仅仅更改介电常数就可以了,命令行统一做如此设置:“# opt rb3lyp/6-31+g(d) scf=tight scrf ” 这样做行吗?
[ Last edited by majun04 on 2008-7-16 at 10:37 ]
GS离可选的是溶剂模型,就是用什么方法来模拟溶剂,不同的数学模型,选用的方法不同,结果有所差别,具体选用什么模型,最好有参考文献做依据;
选溶剂是主要通过介电常数。
但要说明一点:
介电常数只是定义溶剂的众多参数之一,改变介电常数值并不一定就是定义一个新溶剂。
[ Last edited by zzgyb on 2008-7-16 at 12:39 ]
在命令行,我设置的scrf命令,表示在溶剂下的,《gaussian03手册》上不是说在缺省的条件下,scrf命令等价于scrf=iefpcm吗?(我理解的对吗?)
对于溶剂,我的意思是:仅仅是在给出物质结构参数后面加上介电常数,就可以了?还是要在scrf命令后面设置溶剂名称?
例如我算乙腈作为溶剂的:
例子一:
# opt freq=noraman b3lyp/6-31g(d) scrf=iefpcm
结构后面给出乙腈的介电常数
例子二:
# opt freq=noraman b3lyp/6-31g(d) scrf=(iefpcm,solvent=acetonitrile)
(结构后面不再给出乙腈的介电常数)
(两个例子给的分子结构式是一致的)
上面两个例子,对于高斯计算是否是一样任务?
计算出结果发现不一样:
例子一结果:
After PCM corrections, the SCF energy is -269.576720510 a.u.
Dipole Moment 2.9672 Debye
例子二结果:
After PCM corrections, the SCF energy is -269.575977323 a.u.
Dipole Moment 2.9409 Debye
请问为什么会不样?
我需要计算不同溶剂对反应物的影响,在更换溶剂的时候:
对于例子1是否是仅仅更换介电常数值就可以了,要是对于gaussian手册上没有提供介电常数值的溶剂,我从兰氏手册上查到介电常数,输入了(仅仅是更改了介电常数,其他设置不变),gaussian程序能够“认可”,并正确计算吗?
对于例子二,因为gaussview中提供的溶剂较少,我需要的在选项中没有,怎么办?是自己设置参数吗?要设置的参数,除了介电常数还有那些?怎么设置才正确。
课题需要,在自学gaussian,课题组没有可以讨论的人,所以问题较多,也比较初级,恳请答疑,麻烦你了。
[ Last edited by majun04 on 2008-7-16 at 12:41 ]
scrf(Self-Consistent ReactionField)是一个自洽反应场的通用模型方法,这个方法将溶剂描述为连续的,均一的。是目前最常用的,不要把他理解为溶剂反应模型。还有其他的溶剂反应模型,他们描述溶剂的方法是不一样的。
就像我们计算DFT,你选择了DFT并不能计算,还要制定其他的,如果不制定,就按程序默认的基组方法计算。
SCRF也是一样,使用它也要选择具体的模型方法
那就是Onsager等,如果你不选择,系统默认就是Onsager。
从道理上说,你的例1、2,这两个计算结果应该是一样的。
对你你改变介电常数,是否改变溶剂。这样要看你是用的是什么方法,如果有的方法只要求输入介电常数,那么你改变介电常数是可以的。
GS手册没有,你查到也可以使用,GS会认可的,只不过计算比较起来会有误差。
GSV只是一部分,不必全通过他设置,可以自己修改参数设置。不必顾及GSV。
因为别人没有他照样可以设置计算。
谢谢你的帮助。
我是在用gaussian03计算。
刚刚看了一下手册,P226中的gaussian03和gaussian98功能的差别:SCRF和SCRF=PCM现在都默认为SCRF=IEFPCM。刚刚我也算了一下,在命令行用scrf和用scrf=iefpcm的结果是一样的。