24小时热门版块排行榜    

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

majun04

木虫 (著名写手)

majun

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

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

[ Last edited by majun04 on 2008-7-14 at 14:59 ]
回复此楼
乐观、执着、勤奋
已阅   回复此楼   关注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的回帖
查看全部 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的回帖

majun04

木虫 (著名写手)

majun

引用回帖:
Originally posted by zzgyb at 2008-7-16 11:18:
GS离默认的都是水
就是介电常数默认为80

在命令行,我设置的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 ]
乐观、执着、勤奋
5楼2008-07-16 11:32:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工调剂求导师收留!一志愿失利,踏实肯干,有植物提取科研经历 +7 yzyzx 2026-04-09 7/350 2026-04-09 23:21 by may_新宇
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 267求调剂 +5 再忙也要吃饭啊 2026-04-09 5/250 2026-04-09 18:47 by stone_128
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 一志愿211,0703化学305分求调剂 +20 严西西戏 2026-04-06 26/1300 2026-04-08 22:12 by 凯凯要变帅
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 化工学硕 285求调剂 +26 Wisjxn 2026-04-07 26/1300 2026-04-08 14:42 by screening
[考研] 288环境专硕,求调材料方向 +35 lllllos 2026-04-04 39/1950 2026-04-07 23:24 by 一只好果子?
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 材料调剂 +17 小刘同学吖吖 2026-04-06 18/900 2026-04-07 11:41 by 诗与自由
[考研] 334分控制工程求调剂 +4 姜尚真sadasd 2026-04-03 4/200 2026-04-07 09:26 by 蓝云思雨
[考研] 085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干 +3 1021075758 2026-04-06 4/200 2026-04-07 09:23 by 蓝云思雨
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 083200 333求调剂 +3 十二!! 2026-04-04 3/150 2026-04-05 08:28 by barlinike
[考研] 085400电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-03 6/300 2026-04-04 21:50 by hemengdong
信息提示
请填处理意见