| 查看: 2173 | 回复: 7 | |||||||
| 本帖产生 1 个 QC强帖 ，点击这里进行查看 | |||||||
beefly专家顾问 (职业作家)
地沟油冶炼专家
|
[交流]
【beefly个人文集】用Gaussian计算超重元素 已有7人参与
|
||||||
|
用Gaussian计算超重元素，可能涉及到以下几个问题 1）如何指定没有正式命名的元素 以第八周期120号元素Ubn（与镭同一主族）为例。Gaussian不支持用系统命名法命名的元素符号。遇到这种情况，可以输入原子的核电荷数。 2）初始猜测 Gaussian大多数的初始猜测方法都会失效。可以用guess=core，这种方法用单电子哈密顿作初始猜测，不依赖于元素，但是对有些体系可能不好。 3）相对论效应 对于超重元素必须考虑相对论效应。一类方法是相对论ECP，另一类方法是全电子相对论。在Gaussian提供的全电子标量相对论方法中，常用的有二阶Douglas-Kroll-Hess（Int=DKH2）和相对论排除小分量（Int=RESC），以及Gaussian 09加入的四阶Douglas-Kroll-Hess（Int=DKHSO）。参见使用手册的Integral关键词。与相对论计算密切相关的还有原子核的电荷模型。Gaussian对所有的全电子相对论方法默认用Gaussian型电荷分布，但是有时为了和文献作比较，或者基组是用点核模型收缩的，就要用到点核模型。开关是IOP(3/93=1 or 2)。但是，DKHSO只能用Gaussian型电荷分布。 4）基组 对于103号之前的元素，Gaussian提供了非收缩的UGBS基组。103号以后的元素，需要从文献查找基组，或者根据一些公式自己造基组。 如果对未正式命名的元素使用自定义基组，那么在基组输入部分的第一个位置无法用元素符号。可以用原子在坐标中的序号代替，不能用核电荷数。 使用自定义ECP的方法类似。 这些基组有时无法通过Gaussian的积分精度测试，需要加上选项int=noxctest避免检查。 5）限制： 只能用于HF，Post-HF。不能用于半经验方法，因为没有参数化。第八周期以上元素不能做DFT，因为格点生成程序不支持主量子数N>7。需要修改源代码。 例子：120号元素的类氢离子（去掉119个电子）。因为只有1s轨道占据，使用s型基函数就够了。这里用了32个s型Gaussian原函数，指数的计算公式为 a_i = exp(-3.84 + (i-1) * 0.72), i=1,2,...,32 参见文献J. Chem. Phys., 117, 9215, 2002. 标量相对论方法采用DKH2。为了与上面文献的结果进行比较，使用点核模型。输入文件如下： # hf/gen guess=core int(dkh2,noxctest) iop(3/93=1) [120] ^ +119 119 2 120 1 0 s 1 1.0 0.02149360134509 1.0 s 1 1.0 0.04415716841969 1.0 s 1 1.0 0.09071795328941 1.0 s 1 1.0 0.18637397603941 1.0 ... s 1 1.0 25140735.0760292 1.0 s 1 1.0 51649961.0801941 1.0 s 1 1.0 106111395.371615 1.0 **** 这个计算得到的总能量是-9569.1016 a.u.，与上面文献报道的DKH2能量一致。 [ Last edited by beefly on 2011-1-9 at 00:59 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 我的收藏 | 基组的介绍总结 |
» 猜你喜欢
Renlab@monash全奖博士机会，最后三天！
已经有0人回复
-
今年面上没有消息，想来是没了，散个金攒攒人品，为明年祈福
已经有698人回复
-
物理化学论文润色/翻译怎么收费?
已经有83人回复
-
合肥工业大学854仪器技术综合专业课考研
已经有0人回复
-
27年电化学博士毕业，想找海外QS前50的高校读博后求推荐
已经有3人回复
-
27年电池/材料类博士申请
已经有2人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
2027年博士申请或转博科研助理 求推荐
已经有1人回复
-
氯化钠与有机磺酸分离
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 怎么用gaussian计算计算双原子分子势能曲线？
已经有15人回复
- 用gaussian计算丙烷的标准生成热，结果与实验值相差很大，盼高手指导。
已经有7人回复
- 那位大侠指导下如何利用Gaussian计算过渡态
已经有7人回复
- 【求助】用gaussian09进行优化频率的计算出来L103错误，什么原因啊？？
已经有3人回复
- 【求助】用gaussian view打开优化好的结构后,将断裂键连上,是否影响后面计算
已经有9人回复
- 【求助】用Gaussian计算分子的熵时，不同分子的熵的计算基准是一样的吗？谢谢
已经有11人回复
- 【讨论】用gaussian进行SCS-mp2计算
已经有12人回复
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
2楼2011-01-08 13:55:44
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
3楼2011-01-08 14:12:46
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13319.5
- 散金: 27307
- 红花: 301
- 帖子: 4683
- 在线: 928.6小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
4楼2011-01-08 14:15:54
green1986
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 6458.7
- 散金: 13
- 帖子: 254
- 在线: 215.7小时
- 虫号: 1153216
- 注册: 2010-11-22
- 性别: MM
- 专业: 材料物理化学
5楼2011-01-08 20:39:34
qzhaosdu
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1401.6
- 散金: 2052
- 红花: 7
- 帖子: 1969
- 在线: 263.5小时
- 虫号: 493351
- 注册: 2008-01-10
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
6楼2011-01-08 23:23:18
九仙居士
捐助贵宾 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 8414.2
- 散金: 2940
- 红花: 11
- 帖子: 999
- 在线: 472.8小时
- 虫号: 454954
- 注册: 2007-11-10
- 性别: GG
- 专业: 聚合物共混与复合材料
7楼2011-01-09 08:56:20
sculhf
荣誉版主 (知名作家)
断翼天使
-
专家经验: +1 - 应助: 243 (大学生)
- 贵宾: 4.018
- 金币: 43355.9
- 散金: 1634
- 红花: 64
- 沙发: 16
- 帖子: 7688
- 在线: 2008.4小时
- 虫号: 497659
- 注册: 2008-02-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 基金申请
8楼2011-01-09 23:04:25
