| 查看: 1151 | 回复: 6 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】采用溶剂模型是否要使用弥散函数（急）
|
|||
|
在考虑有溶剂的情况时，进行计算选择基组的时候是否要采用弥散函数？《Exploring Chemistry with Electronic Structure Methods》书上进行溶剂计算的时候，怎么都带弥散函数的基组？ 恳请虫友答疑。 [ Last edited by majun04 on 2008-7-14 at 14:59 ] |
回复此楼» 猜你喜欢
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
-
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有17人回复
-
物理化学论文润色/翻译怎么收费?
已经有72人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
070300总分298求调剂
已经有2人回复
-
26博士申请
已经有0人回复
-
一志愿厦门大学化学学硕307求调剂
已经有19人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有2人回复
-
硕士研究生调剂
已经有0人回复
-
投稿期刊求助
已经有1人回复
-
江西师大化学工程学院调剂今晚0点开放！ 缺额100个，速冲！
已经有0人回复
|
在命令行，我设置的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
zzgyb
荣誉版主 (文坛精英)
小木虫警察局局长
- 应助: 0 (幼儿园)
- 贵宾: 22.58
- 金币: 42798.1
- 红花: 19
- 帖子: 10334
- 在线: 199.3小时
- 虫号: 136380
- 注册: 2005-12-16
- 性别: GG
- 专业: 催化化学
- 管辖: 竞技体育
★ ★ ★ ★ ★
suntao1982(金币+5,VIP+0):谢谢交流!
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手册中没有提供介电常数得溶剂，如果我想算，它得基本参数应该到哪里去找才可靠些？ 这个就不好说了，如果没有，就只能近似使用兰氏，或者是其他的一些文献资料，这些只能自己查文献了。 |
2楼2008-07-16 10:09:57
3楼2008-07-16 10:14:50
zzgyb
荣誉版主 (文坛精英)
小木虫警察局局长
- 应助: 0 (幼儿园)
- 贵宾: 22.58
- 金币: 42798.1
- 红花: 19
- 帖子: 10334
- 在线: 199.3小时
- 虫号: 136380
- 注册: 2005-12-16
- 性别: GG
- 专业: 催化化学
- 管辖: 竞技体育
4楼2008-07-16 11:18:06
