| 查看: 1499 | 回复: 15 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
gongxd325至尊木虫 (文坛精英)
|
[交流]
【讨论】基组在关键词行输入和用GEN输入的差别
|
||
|
最近用B3LYP/GEN计算含CNOAg的体系，发现个问题： （1）B3LYP/GEN C N O 0 CEP-121G **** Ag 0 CEP-121G **** 67 basis functions, 110 primitive gaussians, 70 cartesian basis functions Convergence failure -- run terminated. （2）B3LYP/CEP-121G 能顺利结束 70 basis functions, 110 primitive gaussians, 70 cartesian basis functions 按道理两种输入方法都是对所有元素采用了CEP-121G基组，但结果居然大不一样： 一个能正常完成另一个不能，而且基函数数目也不等； 两种输入方法究竟有何差别？ 用gen做计算的不少，究竟结果的可靠度如何？ gen常出现不收敛的问题（Convergence failure或Inaccurate quadrature in CalDSu），如何解决？ [ Last edited by gongxd325 on 2009-10-9 at 10:28 ] |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有58人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
gongxd325
至尊木虫 (文坛精英)
- 应助: 177 (高中生)
- 金币: 37616.6
- 散金: 13187
- 红花: 109
- 帖子: 34003
- 在线: 1844小时
- 虫号: 233962
- 注册: 2006-03-30
- 性别: MM
- 专业: 理论和计算化学
|
这是输入 % chk=1 # B3LYP/GENECP OPT GFinput optcyc=200 Test1 0 1 Ag 0.000000 0.973585 0.000000 O -0.035321 -1.043566 0.000000 N 0.003644 -2.322076 0.000000 C 0.042843 -3.525907 0.000000 C O N 0 CEP-121G **** Ag 0 CEP-121G **** Ag 0 CEP-121G --Link1-- % chk=2 # B3LYP/CEP-121G OPT GFinput optcyc=200 Test2 0 1 Ag 0.000000 0.973585 0.000000 O -0.035321 -1.043566 0.000000 N 0.003644 -2.322076 0.000000 C 0.042843 -3.525907 0.000000 测试了5D，对结果没有影响。 用LanL2DZ，两种输入方式的计算结果完全一样，能量都是-313.697668au。 但CEP-121G的结果就不一样，basis functions分别是67和70，能量是-281.819101和-177.65283au。 似乎是基组的问题？ |
13楼2009-10-11 17:20:48
abbott
金虫 (著名写手)
不要用QQ问我东西
- 应助: 16 (小学生)
- 金币: 1071.6
- 散金: 2787
- 红花: 10
- 帖子: 1015
- 在线: 105.2小时
- 虫号: 452267
- 注册: 2007-11-05
- 性别: GG
- 专业: 计算机硬件技术
2楼2009-10-09 11:42:36
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 553.8
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
3楼2009-10-09 12:53:18
hairan
木虫 (著名写手)
瞌睡木虫
- QC强帖: 3
- 应助: 61 (初中生)
- 贵宾: 1.06
- 金币: 4984.3
- 散金: 8072
- 红花: 44
- 帖子: 2900
- 在线: 733.1小时
- 虫号: 432069
- 注册: 2007-08-12
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 10-9 21:19
gongxd325(金币+2,VIP+0):Thanks 10-10 08:33
lei0736(金币+4,VIP+0):谢谢 10-9 21:19
gongxd325(金币+2,VIP+0):Thanks 10-10 08:33
|
第一种方式，你使用了赝势基组却没有读入赝势，当然不能正常结束 估计优化过程中肯定还出现围绕Ag的键长变得异常的短 正确的使用方式是 B3LYP/genecp 或者 B3LYP/gen pseudo=read C N O 0 CEP-121G **** Ag 0 CEP-121G **** ！此处有一个空行 Ag 0 CEP-121G 第二种方式与上面的写这么多是等效的，所以如果体系内所有的原子都用一种基组，而且是Gaussian自带的基组，那么就不必使用gen相关的关键词 [ Last edited by hairan on 2009-10-9 at 20:56 ] |
5楼2009-10-09 20:51:29