24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2036  |  回复: 5
本帖产生 1 个 QC强帖 ,点击这里进行查看

coolrainbow

木虫 (著名写手)

未来国家冻凉

[交流] 【coolrainbow个人文集】小分子实验结构与理论计算的对应关系 已有5人参与

小分子的几何结构,一般可以用光谱方法,比如转动光谱,很精确的测定,比如氮气分子 N2, 我们可以得到极精确的值:1.100105Å。

但是事情没有那么简单。大家都知道分子是在不停的振动的,即使0K也有所谓零点能ZPE存在。因此,所谓 “氮气键长多少” 就成了一个不明确(ill-defined)的问题,在振动时候,N2的两个氮原子的距离可以在(比如) 1.10009~1.0011 Å之间振动,究竟哪个算键长?就不好说了。这就引出了所谓 r0 和 re 结构的问题。

以双原子转动光谱的测定为例:

             若把分子是为刚性转子,则转动光谱项 F(J)=BeJ(J+1) (J为角动量量子数),

Be与键长有关(有确定的公式,这里写作Be=C/r^2),因此只要测出Be,就能够得到键长 re。

              对于大J,需加“离心校正”,则转动光谱项 F(J)=BeJ(J+1)-De(J(J+1))^2,

其中De=4Be^3/w,w为分子的振动频率。只要我们的光谱范围能包含足够多的量子数,就能得到这些值。

但是分子在不停地振动,在简谐近似下,振动可以用 v = 0,1,2,3, ...作为好量子数,在不同的振动下,分子的转动也不同,于是F(J)中的Be与振动转态 v 有关,所以:

             考虑振动情况下,转动光谱项 Fv(J)=BvJ(J+1)-Dv(J(J+1))^2,  Bv = Be-alpha*(v+1/2)

好吧,问题来了。最“真实”的情形,应该是我们说的最后一个公式 Fv ,根据Boltzman分布,v = 0是最多的,于是,光谱中强度最大的就是v=0的项:

             F0(J)=B0J(J+1)-D0(J(J+1))^2,  B0 = Be-alpha*1/2

也就是说,光谱中实际测出来的常数,如果忽略掉 v>0 项,都是B0而不是Be值,如果直接令B0=C/r^2,则算出来的 r 显然是某种对应于零振动态的平均值,这个结构我们称为r0 结构。 如果我们的仪器足够好,可以测出 v=1态的转动光谱,就可以用两个数据解两个未知数,同时求得Be和alpha,再用Be=C/r^2求出 r,这个键长 r 就是对应于势能曲线中的最小值的位置,也就是“平衡结构”,我们称为re 结构。

         那么,r0 和 re 那个是键长?我们看 N2:

                             r0                      re
(14N)2            1.1001             1.0976
(15N)2            1.1000             1.0976
14N+15N       1.1000                 N/A

从 re 的定义来看,re 与 核质量无关,而由于 r0 涉及到振动,因此与核质量有关,也就是说,r0 是同位素相关的。从上面明显的看出,同样分子,r0与同位素有关,re无关。同时,对同一分子,r0和re可以相差仅0.0025Å,可能不明显。不过,个别微波转动光谱的实验中,r0 和 re 差别可达0.01Å。可见,似乎re当键长更合理些。

更加复杂情况是,在大的小分子(如4个原子)结构分析时,由于实验值不足以决定 r0 结构的全部的变量,有时候则需要做些“constraint”。例如,对于分子ZnCH3,用转动光谱计算其结构时,需要估计C-H键长并将这个值固定(因为光谱参数不够。可笑的是,这个固定的值是用量化算出来的),然后再用光谱公式算出Zn-C键长,称之为“实验值“。(见J. Chem. Phys. 1993, 99, 9376–9388)这个"实验值"的可靠性就值得怀疑了。
        在量子化学计算中,基于Born-Oppenheimer近似,我们可以知道,量子化学计算中优化的结构就是 re 结构,因此,我们量化计算的结果应该和实验的 re 结构比较。不幸的是,从技术的角度看,对于大多数分子,测定 re 结构是非常困难的,事实上,很多利用转动光谱测得的小分子的结构都是 r0 结构。因此,量化计算预测的结构和 r0 结构是必然有微小差别的。

       当然,这并不是给量化人员提供了预测结构与实验不符时的借口。对于大多数情况,r0 结构和re 结构的差别不会到0.005Å以上。当你的结构与实验值(比如 r0 结构)相差达0.01Å数量级时,最好还是先考虑下你的方法模型问题,而不是责怪实验人员没有给你 re 结构~~~~

[ Last edited by 375642546 on 2011-4-19 at 15:34 ]
回复此楼
技术博客:http://hi.baidu.com/coolrainbow/blog
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qlihytc

新虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
老板经常告诉我们:要知道化学精度!键长多少范围内,你看有些做量化的表格的小数点范围就能看出他懂不懂计算!
2楼2011-03-21 12:25:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

x7511413

金虫 (正式写手)

科研小菜


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by qlihytc at 2011-03-21 12:25:32:
老板经常告诉我们:要知道化学精度!键长多少范围内,你看有些做量化的表格的小数点范围就能看出他懂不懂计算!

请问一般是多少位呢 谢谢啦
哦也
3楼2011-03-26 13:00:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖交流
(1)键长化学精度多少?
(2)键能(能量)化学精度多少? 0.01meV/atom?
4楼2011-03-26 15:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiechong

金虫 (著名写手)

好贴,要顶一下
5楼2011-03-27 09:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yasistar

禁虫 (正式写手)

本帖内容被屏蔽

6楼2011-04-01 14:59:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 coolrainbow 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 材料工程302分求调剂 +9 zyx上岸! 2026-04-04 9/450 2026-04-05 22:08 by 醉翁wl
[考研] 一志愿同济大学323分(080500)求调剂 +8 yikeniu 2026-04-01 8/400 2026-04-05 18:15 by cql1109
[考研] 080200学硕,机械工程专业277分,求带走! +7 瓶子PZ 2026-03-31 7/350 2026-04-05 17:49 by liucky
[考研] 288求调剂 +7 没有答案_ 2026-04-05 7/350 2026-04-05 17:01 by yulian1987
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 男生,一志愿沪9生物学071000,初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 340求调剂 +4 jhx777 2026-03-29 4/200 2026-04-04 20:08 by 无际的草原
[考研] 一志愿沪9,求生物学调剂,326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[论文投稿] 求文献 5+3 ys879651$ 2026-04-02 3/150 2026-04-04 17:22 by bobvan
[考研] 363求调剂 +7 zh096 2026-04-04 7/350 2026-04-04 17:13 by dongzh2009
[考研] 考研调剂 +4 zybz冲冲冲 2026-04-03 6/300 2026-04-04 13:08 by zybz冲冲冲
[考研] 336求调剂 +8 kiyy 2026-04-01 8/400 2026-04-03 19:41 by lijunpoly
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 一志愿厦门大学材料工程专硕354找调剂!!! +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-02 3/150 2026-04-02 15:06 by cal0306
[考研] 322求调剂 +5 熹僖XX 2026-03-31 6/300 2026-04-02 10:08 by 求调剂zz
[考研] 一志愿北京科技,085601总分305求调剂 +9 半生瓜! 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
信息提示
请填处理意见