24小时热门版块排行榜    

查看: 1065  |  回复: 6
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

majun04

木虫 (著名写手)

majun

[交流] 【求助】采用溶剂模型是否要使用弥散函数(急)

在考虑有溶剂的情况时,进行计算选择基组的时候是否要采用弥散函数?《Exploring Chemistry with Electronic Structure Methods》书上进行溶剂计算的时候,怎么都带弥散函数的基组?
恳请虫友答疑。

[ Last edited by majun04 on 2008-7-14 at 14:59 ]
回复此楼
乐观、执着、勤奋
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

谢谢你的帮助。
我是在用gaussian03计算。
刚刚看了一下手册,P226中的gaussian03和gaussian98功能的差别:SCRF和SCRF=PCM现在都默认为SCRF=IEFPCM。刚刚我也算了一下,在命令行用scrf和用scrf=iefpcm的结果是一样的。
乐观、执着、勤奋
7楼2008-07-16 13:05:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

zzgyb

荣誉版主 (文坛精英)

小木虫警察局局长

优秀版主

★ ★ ★ ★ ★
suntao1982(金币+5,VIP+0):谢谢交流!
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《郭靖,这就是盲目扩招的后果
2楼2008-07-16 10:09:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

非常感谢!

问题3,还想问一下:如果不是gaussian中默认的溶剂,是否能仅仅更改介电常数就可以了?还是要有其他的改动?
我是想和实验相比较,打算采用兰氏的数据。在设置输入文件的时候,是否是仅仅更改介电常数就可以了,命令行统一做如此设置:“# opt rb3lyp/6-31+g(d) scf=tight scrf ” 这样做行吗?

[ Last edited by majun04 on 2008-7-16 at 10:37 ]
乐观、执着、勤奋
3楼2008-07-16 10:14:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzgyb

荣誉版主 (文坛精英)

小木虫警察局局长

优秀版主

★ ★
suntao1982(金币+2,VIP+0):辛苦了,兄弟!
GS离可选的是溶剂模型,就是用什么方法来模拟溶剂,不同的数学模型,选用的方法不同,结果有所差别,具体选用什么模型,最好有参考文献做依据;
选溶剂是主要通过介电常数。
但要说明一点:
介电常数只是定义溶剂的众多参数之一,改变介电常数值并不一定就是定义一个新溶剂。

[ Last edited by zzgyb on 2008-7-16 at 12:39 ]
江南七怪+马玉+洪七公+一灯大师+周伯通教一个弱智儿童=郭靖,天下第一王重阳教七个天资不错的小孩=全真七子,全真七子×3《郭靖,这就是盲目扩招的后果
4楼2008-07-16 11:18:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见