| 查看: 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 ] |
回复此楼» 猜你喜欢
化工安全有关的复习提纲
已经有0人回复
-
国家级人才课题组招收2026年入学博士
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有230人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
散金
已经有80人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
已经有32人回复
-
香港中文大学（深圳）2026年诚聘海内外优青
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
gongxd325
至尊木虫 (文坛精英)
- 应助: 177 (高中生)
- 金币: 37616.6
- 散金: 13187
- 红花: 109
- 帖子: 34003
- 在线: 1844小时
- 虫号: 233962
- 注册: 2006-03-30
- 性别: MM
- 专业: 理论和计算化学
★ ★ ★
yjcmwgk(金币+3,VIP+0):问题解决 10-12 13:48
yjcmwgk(金币+3,VIP+0):问题解决 10-12 13:48
|
现在搞清楚了，如下输入，两种方式才等价： # B3LYP/GENECP OPT 6D 10F 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 **** C O N 0 CEP-121G Ag 0 CEP-121G --Link1-- # B3LYP/CEP-121G OPT 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 用GENECP时CON也要读ECP，还要加关键词6D和10F，但LanL2DZ基组不需要这样。 |
16楼2009-10-12 09:43:38
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