24小时热门版块排行榜    

北京石油化工学院2025年硕士研究生招生接受调剂公告
查看: 1365  |  回复: 4

oyljw

至尊木虫 (著名写手)

格物致知Professional

[求助] 计算晶体怎样减少键长键角的相对误差?怎样确定中心原子周围键的键组成?

Sample TextSample Text如题: 计算晶体选择啥样的计算模型能尽量减少键长键角的相对误差?怎样确定中心原子周围键的组成?


以前自己做了个晶体,结构一般,想结合量化计算发个很小的杂志。原来用半经验软件MOPAC用PM3计算http://www.muchong.com/bbs/viewthread.php?tid=2117576

投稿快一年了,现在编辑来信说:
The application of the PM3 method to a single molecule of medium size is not acceptable.
I suggest that authors perform ab initio calculations and NBO analysis, in order to provide bond orders, net charges, and compositions of the bonds around Zn. The calculated bond lengths and bond angles should be included in Table 1. The relative errors for bond lengths and bond angles should be calculated. The low values for the errors (< 1%) should support the applicability of the computational model.

我改用高斯从头算,现在的问题是:
1.如何在高斯的log文件中provide compositions of the bonds around Zn.需要找哪个关键词。
2.计算键长键角相对误差,我用b3lyp+自定义机组(输入文件在下面注释中)计算得出结果后,应该对中心原子周围的主要的键长键角进行计算相对误差,以键长相对误差为例:键长相对误差=(晶体结构键长数据-计算的键长数据)的绝对值/晶体结构键长数据  X100%  这样做对吗?请指正。
另外,编辑说The low values for the errors (< 1%) should support the applicability of the computational model.低的误差(< 1%) 表明计算模型适用。我用b3lyp+自定义机组得到的结果相对误差都大于1%(最大的19.51%,最小的1.15%)是不是说明我的方法基组不合适?如果不合适应该换用什么基组?甚至用什么软件更好?





注释:我的输入文件为
%chk=aa.chk
%mem=2500MB
%nprocshared=8
#p opt freq b3lyp/genecp pop=NBOread

CSD ENTRY aa

0 1
Zn
N                  1    2.34247177
N                  1    2.26583009    2   74.67133567
N                  1    2.26470418    3  108.45387328    2  -84.27206505    0
N                  1    2.33272502    4   74.75545221    2  160.99037416    0
N                  3    2.17823069    1  148.25169943    4  -86.55567440    0
H                  6    0.86092684    3  160.41542137    1   -1.01784456    0
N                  4    2.20504647    1  147.50844368    3  -96.89107326    0
H                  8    0.86012208    4  160.80157179    1    8.22173323    0
N                  2    2.20403925    1  144.49312508    4  118.34174000    0
H                 10    0.85948066    2  159.86688321    1   -4.67754538    0
N                  5    2.19585177    1  145.31437552    4   -2.31310391    0
H                 12    0.85905646    5  160.71415979    1    7.78773274    0
O                  1    2.26530550    4  165.05672296    8   73.65400787    0
O                  1    2.32188738    4   87.36376751    8   89.42200324    0
O                 14    2.21802187    1  127.16598495    4   -3.87840996    0
O                 15    2.20426155    1  131.98817650    4   60.68101982    0
C                  8    1.36535783    4   72.37628372    1 -170.55864504    0
H                 18    0.92910710    8  126.28865769    4 -179.25965044    0
C                  4    1.31792299    1  112.56362863   14   80.11781918    0
C                  2    1.37645487    1  143.65164353    4  -72.08145626    0
H                 21    0.92798761    2  125.36748534    1    6.25531045    0
C                 17    2.37046936   15   62.48969697    1   82.92040062    0
C                  5    1.32101855    1  110.41221636    4    1.02756083    0
C                  6    1.37106163    3   73.33731042    1  179.37514784    0
H                 25    0.93157233    6  127.24743374    3  179.81070515    0
C                 10    1.37825578    2   73.88690557    1  173.66505552    0
H                 27    0.92993817   10  127.37067947    2  179.92442106    0
C                 23    1.39411978   17   92.06842273   15 -168.68933121    0
H                 29    0.92997258   23  119.51905308   17    1.41090098    0
C                 29    1.37268604   23  121.03945760   17 -178.64582025    0
H                 31    0.92937237   29  120.44902562   23 -177.30992652    0
C                 25    1.35085826    6  105.44858219    3   -0.11328181    0
H                 33    0.92967682   25  124.98905021    6 -179.74821115    0
C                 17    1.23806098   15   29.12935628    1   84.24421924    0
C                 31    1.39838837   29  119.24370484   23    2.64235067    0
H                 36    0.93014515   31  120.28024937   29 -179.27820822    0
C                 23    1.39360863   17  147.87132833   15   21.01926127    0
C                  2    1.31646990    1  110.27011083    4  113.83472528    0
C                 16    1.22801832   14   28.56196567    1   79.99893105    0
C                  3    1.31382114    1  113.69408285    4  -86.82019555    0
C                  5    1.37841939    1  143.49945056    4 -172.33062941    0
H                 42    0.93078945    5  125.55868867    1   -6.13936993    0
C                 36    1.36194457   31  119.67088988   29    0.77534788    0
H                 44    0.93133560   36  118.94760131   31  177.72557813    0
C                 42    1.34835195    5  108.87826344    1  173.68511212    0
H                 46    0.92933202   42  126.34591429    5 -179.63740854    0
C                 18    1.34070764    8  107.49813937    4    0.72070981    0
H                 48    0.92983708   18  125.26220807    8  178.75375593    0

C H N O 0
6-31g(d)
****
Zn  0
Lanl2dz
****

Zn  0
$nbo bndidx $end


对我上面的两个问题 我该怎样处理?麻烦高手出招!!谢谢!!
回复此楼
独立之精神自由之思想
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

涂易

新虫 (小有名气)

【答案】应助回帖

oyljw(金币+4): 谢谢提醒 2011-07-20 09:57:11
你的分子中  有很多共轭体系  最好加上弥散函数和极化函数,你用posmol打开输出文件 可以知道哪些原子间形成键 及键矩等信息 不过键已经形式化了  不能区分单双键
首孝悌,次谨信,泛爱众,而亲人,有余力,则学文
2楼2011-07-20 09:53:48
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

oyljw

至尊木虫 (著名写手)

格物致知Professional

引用回帖:
Originally posted by 涂易 at 2011-07-20 09:53:48:
你的分子中  有很多共轭体系  最好加上弥散函数和极化函数,你用posmol打开输出文件 可以知道哪些原子间形成键 及键矩等信息 不过键已经形式化了  不能区分单双键

请问 极化是加d呢 还是d,p呢
弥散是加+,还是++?
C H N O 0
6-31g(d)
****
Zn  0
Lanl2dz
****

Zn  0


改成
C H N O 0
6-311++g(d,p)
****
Zn  0
Lanl2dz
****

Zn  0

可以吗?还是对氢有特殊的要求?我根基不好请帮助指正 谢谢
独立之精神自由之思想
3楼2011-07-20 10:45:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

涂易

新虫 (小有名气)

【答案】应助回帖

★ ★
oyljw(金币+2): 2011-07-20 14:03:41
gmy1990(金币+2): 谢谢你对论坛的贡献,奖励一下! 2011-07-20 21:55:48
对于有孤电子体系和共轭体系的  分子在选择基组上最好加上弥散和极化函数如:6-31G+(d),是在6-31G(d)的基础上对重原子添加弥散函数;6-31G++(d)是在6-31G+(d)的基础上对H添加弥散函数;6-31G(d,p)的意思是对6-31G的基础上对重原子加一个d函数,对H添加一个p函数,就是分别添加了极化函数;除了H元素外,其他都简称重元素
首孝悌,次谨信,泛爱众,而亲人,有余力,则学文
4楼2011-07-20 11:12:50
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

oyljw

至尊木虫 (著名写手)

格物致知Professional

非常感谢
独立之精神自由之思想
5楼2011-08-12 20:39:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 oyljw 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 211理学化学0703调剂 +7 陈橙ii 2025-03-30 7/350 2025-04-05 01:22 by 海誓山盟
[考研] 269求调剂 +5 若春白雪 2025-04-04 5/250 2025-04-04 23:51 by 褚迎松小小
[考研] 化学0703调剂自荐 一志愿东北大学312 +16 爱科研的小王同� 2025-04-03 21/1050 2025-04-04 23:12 by 河大薇薇
[考研] 085602化学工程,273分 +18 青阳21st 2025-04-01 19/950 2025-04-04 19:01 by DNDYR
[考研] 求调剂!一志愿南京理工大学材料工程,总分344,有1篇SCI,1个专利 +8 ishimou 2025-04-04 8/400 2025-04-04 18:49 by DNDYR
[文学芳草园] 阿里的牦牛肉(下) +3 南泽洲Siren 2025-03-29 6/300 2025-04-04 17:14 by 南泽洲Siren
[考研] 278求调剂 +3 芸生wfy 2025-04-04 4/200 2025-04-04 09:24 by 川哥v5
[考研] 324材料化工求调剂 5+3 HOPE!! ! 2025-04-03 14/700 2025-04-04 07:34 by liangzixue
[考研] 求调剂 +5 芜湖呢 2025-04-01 5/250 2025-04-03 19:10 by 诗与自由
[考研] 278求调剂 +9 202103847 2025-04-02 9/450 2025-04-03 14:57 by fmesaito
[考研] 283分机械专硕不挑专业 +7 汪某. 2025-03-31 7/350 2025-04-03 11:02 by 郭艺艺加油
[考研] 材料调剂 +11 @Wg1210618 2025-04-01 12/600 2025-04-03 10:07 by 青春木花
[考研] 考研材料求调剂 +5 wwwj12300 2025-04-02 6/300 2025-04-02 21:38 by nanaliuyun
[考研] 中科大有机化学,312分寻求调剂 +5 欢欢想变强 2025-04-01 7/350 2025-04-01 22:29 by 长江白鲟
[考研] 0856材料专硕求调剂331分一志愿西安交通大学 +3 巧克力酱? 2025-04-01 3/150 2025-04-01 18:32 by ximiga713
[考研] 调剂 +4 慕贤笙1 2025-03-31 5/250 2025-04-01 11:28 by 慕贤笙1
[考研] 求化工调剂 +10 清河清 2025-03-30 10/500 2025-04-01 10:12 by houyaoxu
[考研] 调剂 +9 1141257783 2025-03-29 10/500 2025-04-01 10:08 by Cloris-y
[考研] 313求调剂 +6 Pj. 2025-03-31 6/300 2025-03-31 18:49 by DNDYR
[考研] 303求调剂 +8 TuTuXii 2025-03-30 9/450 2025-03-31 18:46 by 发疯的石头
信息提示
请填处理意见
澶嶈瘯璋冨墏澶嶈瘯璋冨墏
澶嶈瘯璋冨墏浜岀淮鐮�鑾峰彇澶嶈瘯 & 璋冨墏鎸囧崡锛�寰俊鎵竴鎵�
鍟嗗姟鍚堜綔鍟嗗姟鍚堜綔
鍟嗗姟鍚堜綔浜岀淮鐮�鍟嗗姟鍚堜綔鍜ㄨ寰俊鎵竴鎵�(澶囨敞锛氬晢鍔″悎浣�)
瀛︽湳骞茶揣瀛︽湳骞茶揣
瀛︽湳骞茶揣浜岀淮鐮�寰俊鎵竴鎵�浜嗚В鏈€鏂板鏈共璐�
璁烘枃杈呭璁烘枃杈呭
璁烘枃杈呭浜岀淮鐮�寰俊鎵竴鎵�娣诲姞鑰佸笀锛屽厤璐瑰挩璇�
甯堣祫鎷涘嫙甯堣祫鎷涘嫙