| 查看: 2312 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 我的收藏 |
» 猜你喜欢
散金求好运
已经有315人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有91人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
祈福青基
已经有34人回复
-
瑞典-博后-固态电池-粘结剂-锂离子电池方向
已经有53人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
-
牛来！米来！面来！
已经有8人回复
-
广东省面上有什么动态么，啥时候公布呀？
已经有1人回复
-
2027申博咨询
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 求问一般优化小分子结构用什么方法和机组
已经有10人回复
- 【求助】计算荧光光谱的疑惑，求教！
已经有18人回复
2楼2011-03-21 12:25:32
3楼2011-03-26 13:00:35
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
4楼2011-03-26 15:03:20
5楼2011-03-27 09:37:11
|
本帖内容被屏蔽 |
6楼2011-04-01 14:59:33
10