24小时热门版块排行榜    

北京石油化工学院2025年硕士研究生招生接受调剂公告
查看: 1368  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 351求调剂 +5 nbnbnbb 2025-04-07 6/300 2025-04-07 14:48 by lijunpoly
[考研] 292求调剂 +12 书馆明公 2025-04-05 13/650 2025-04-07 14:35 by lijunpoly
[考研] 化学学硕339分求调剂 +13 灯塔的水母 2025-04-07 14/700 2025-04-07 13:36 by DNDYR
[考研] 277求调剂,求老师收留 +3 hxx43 2025-04-06 3/150 2025-04-07 08:56 by 学员3lybWG
[考研] 材料与化工求调剂 +13 cs0106 2025-04-03 15/750 2025-04-07 07:32 by houyaoxu
[考研] 材料与化工284调剂 +7 老田的暖手宝 2025-04-03 7/350 2025-04-07 06:56 by houyaoxu
[考研] 材料化工调剂 +11 α,β,γ 2025-04-04 11/550 2025-04-07 06:52 by houyaoxu
[考研] 294求调剂 +6 微笑男孩~ 2025-04-06 6/300 2025-04-06 22:32 by pphmily
[考研] 机械专硕085501,已过四级,本科双一流求调剂 +5 黑空海 2025-04-06 6/300 2025-04-06 16:12 by zsg5560
[考研] 化学调剂 +4 不知道起什么名� 2025-04-06 4/200 2025-04-06 15:36 by calebby
[考研] 283求调剂 +3 reht 2025-04-04 3/150 2025-04-05 16:53 by yulian1987
[考研] 263求调剂 +10 36木木 2025-04-03 10/500 2025-04-04 21:34 by Ade.
[考研] 一志愿中科院 315 材料与化工,求调剂 +14 假老练1 2025-04-01 16/800 2025-04-04 20:12 by asdhh1991
[考研] 安徽理工大学材料工程调剂总分302 +14 _梦_ 2025-04-01 16/800 2025-04-04 18:15 by 考研人15000
[考研] 求调剂,不挑专业 +3 2246197870 2025-04-03 3/150 2025-04-04 08:07 by Ldh2580
[考研] 调剂!一志愿郑大283材料与化工专硕 +6 不吃香菜菜菜 2025-03-31 9/450 2025-04-03 19:35 by 13137570862
[考研] 山东中医药大学青岛中医药科学院考研调剂 +3 是个甜豆豆呀 2025-04-02 6/300 2025-04-02 23:44 by 是个甜豆豆呀
[考研] 材料调剂,302一志愿湖南大学专硕,强烈意愿读博 +6 高歌奋进 2025-04-02 6/300 2025-04-02 23:44 by lixin_bift
[考研] 材料与化工310求调剂院校 +18 五毛冰淇淋a 2025-04-01 20/1000 2025-04-02 23:19 by 星空星月
[考研] 化工一志愿北京理工大学292分求调剂 +5 化工方向 2025-04-01 5/250 2025-04-01 23:04 by cqjd0301
信息提示
请填处理意见
婢跺秷鐦拫鍐ㄥ婢跺秷鐦拫鍐ㄥ
婢跺秷鐦拫鍐ㄥ娴滃瞼娣惍锟�閼惧嘲褰囨径宥堢槸 & 鐠嬪啫澧忛幐鍥у础閿涳拷瀵邦喕淇婇幍顐$閹碉拷
閸熷棗濮熼崥鍫滅稊閸熷棗濮熼崥鍫滅稊
閸熷棗濮熼崥鍫滅稊娴滃瞼娣惍锟�閸熷棗濮熼崥鍫滅稊閸溿劏顕�瀵邦喕淇婇幍顐$閹碉拷(婢跺洦鏁為敍姘櫌閸斺€虫値娴o拷)
鐎涳附婀抽獮鑼舵彛鐎涳附婀抽獮鑼舵彛
鐎涳附婀抽獮鑼舵彛娴滃瞼娣惍锟�瀵邦喕淇婇幍顐$閹碉拷娴滃棜袙閺堚偓閺傛澘顒熼張顖氬叡鐠愶拷
鐠佺儤鏋冩潏鍛嚤鐠佺儤鏋冩潏鍛嚤
鐠佺儤鏋冩潏鍛嚤娴滃瞼娣惍锟�瀵邦喕淇婇幍顐$閹碉拷濞h濮為懓浣哥瑎閿涘苯鍘ょ拹鐟版尒鐠囷拷
鐢牐绁幏娑樺珯鐢牐绁幏娑樺珯