24小时热门版块排行榜    

汕头大学海洋科学、生物学、生物与医药等3个专业接受调剂
查看: 1363  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂 +5 Summer0701 2025-04-03 6/300 2025-04-03 14:25 by siwuhen
[考研] 材料与化工(085600)求调剂 +10 想上个好学校啊 2025-04-03 10/500 2025-04-03 14:21 by 木卜心
[考研] 268求调剂 +5 千千运气 2025-03-31 5/250 2025-04-03 13:52 by ZZU2011
[考研] 一志愿东南大学材料学硕,总分342分求调剂 +11 好想读研啊啊 2025-03-31 11/550 2025-04-03 13:24 by 诗与自由
[考研] 304分求调剂,一志愿厦门大学085602化学工程 +13 wang.yao. 2025-04-03 13/650 2025-04-03 13:16 by 诗与自由
[考研] 323农学调剂 +4 00@ 2025-03-29 8/400 2025-04-03 12:51 by cbzhang
[考研] 材料与化工286求调剂 +7 Panuu 2025-04-02 8/400 2025-04-03 10:48 by Panuu
[考研] 0805材料科学与工程328本2求调剂 +6 修。。。。 2025-04-01 6/300 2025-04-03 06:57 by microscopy
[考研] 资源与环境初试274求调剂 +11 程123-- 2025-03-30 16/800 2025-04-02 23:48 by zllcz
[考研] 0703化学求调剂 +8 13349027219 2025-03-29 11/550 2025-04-02 23:29 by siwuhen
[考研] 一志愿西南交大,材料085600代码,总分292求调剂 +12 小白小白… 2025-03-30 12/600 2025-04-02 21:49 by luosha500
[考研] 求调剂 +12 Ther蟹黄堡 2025-03-28 16/800 2025-04-02 13:54 by sqc0802
[考研] 求调剂 +11 速回大呼小叫 2025-04-01 11/550 2025-04-02 10:27 by panzt
[考研] 085600 求调剂 +8 cs0106 2025-04-01 8/400 2025-04-02 09:11 by Ade.
[考研] 085500机械,261分,求调剂 +4 一只小肥羊呀 2025-03-29 8/400 2025-04-01 18:18 by 杨光于青云
[考研] 086001,294求调剂 +8 侯hr 2025-03-29 13/650 2025-04-01 09:43 by zhfms
[考研] 机械专硕297调剂 +4 李lfff 2025-03-30 6/300 2025-03-30 20:12 by 阿健1234
[考研] 294求调剂 +3 侯hr 2025-03-30 3/150 2025-03-30 15:11 by 耀之银魂
[考研] 一志愿上海交通大学 085800能源动力专硕 英一数一 本科北京科技大学 求调剂! +4 统计系陈威宁 2025-03-28 4/200 2025-03-29 14:36 by 褚迎松小小
[考研] 求调剂 +3 往你家跑 2025-03-27 3/150 2025-03-29 07:45 by sy天使2009
信息提示
请填处理意见
复试调剂复试调剂
复试调剂二维码获取复试 & 调剂指南!微信扫一扫
商务合作商务合作
商务合作二维码商务合作咨询微信扫一扫(备注:商务合作)
学术干货学术干货
学术干货二维码微信扫一扫了解最新学术干货
论文辅导论文辅导
论文辅导二维码微信扫一扫添加老师,免费咨询
师资招募师资招募