| 查看: 1860 | 回复: 2 | ||||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||||
我的小名木虫 (小有名气)
|
[求助]
振动模折合质量的计算
|
|||
|
两个质点的折合质量:u=m1×m2/(m1+m2) 在计算分子的振动频率的时候,每一个振动模中都有好几个原子在运动。g09给出的reduced masses是如何计算的? 在我计算的分子中,有一个C=C振动模,g09给的折合质量是6.05,看着就是这么回事儿。 还有一个C=O振动模,g09给的折合质量是10.77。自己按公式算一下6.86 ,问题就来了。 我想这一个C=O振动模,还有其它原子也在运动,因此折合质量不是这么简单计算的。 求大牛指教。 |
» 猜你喜欢
书籍求助:汽车市场营销理论与实务(电子版)——章小平
已经有0人回复
七嗪类物质合成求助
已经有0人回复
物理化学论文润色/翻译怎么收费?
已经有237人回复
东方理-中科大联合博士生招聘
已经有0人回复
钴酸锂半电池小倍率容量上不去
已经有1人回复
吉林大学材料物理本科生求问调剂信息
已经有23人回复
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯!
已经有1人回复
光电子赋能·半导体创芯!
已经有1人回复
光电子赋能·半导体创芯!
已经有0人回复
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
分子动力学模拟合作需求
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐,对您同样有帮助:
某一电子态上各振动能级的能量如何计算?
已经有4人回复
振动能级之间的能量差怎么计算?
已经有7人回复
振动分辨的电子光谱的计算
已经有30人回复
吸收和发射光谱计算后振子强度大于2,这合理吗?
已经有3人回复
求教:声固耦合中薄膜的振动问题
已经有20人回复
Dmol3计算周期性体系振动频率时k点的设置问题
已经有3人回复
如何用拉格朗日方程计算考虑弹簧质量的振动
已经有8人回复
如何用gaussian软件计算分子的振动光谱?
已经有8人回复
如何用VASP计算H原子的振动频率
已经有16人回复
GenAdia:计算键长/键角/二面角的力常数和振动频率的程序
已经有36人回复
【交流】振动模式与红外光谱、拉曼光谱的原理与计算问题的讨论
已经有14人回复
vasp计算振动强度
已经有15人回复
【求助】castep计算声子谱,为什么不能进行振动分析
已经有9人回复
【求助】振动熵的第一原理计算
已经有3人回复
【求助】一维谐振子的振动频率
已经有7人回复
【求助】计算振动模式
已经有14人回复
【求助】ph.x计算振动频率的问题
已经有15人回复
【求助】用MS中的castep计算声子振动模频率,具体怎么算
已经有10人回复

我的小名
木虫 (小有名气)
- 应助: 16 (小学生)
- 金币: 1691.5
- 散金: 941
- 红花: 5
- 帖子: 261
- 在线: 195.3小时
- 虫号: 1744705
- 注册: 2012-04-08
- 性别: GG
- 专业: 理论和计算化学

3楼2014-05-22 09:30:45
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-

专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 12839.4
- 散金: 27307
- 红花: 300
- 帖子: 4667
- 在线: 925.7小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★ ★
gmy1990: 金币+3 2014-05-22 04:42:22
gmy1990: 金币+3 2014-05-22 04:42:22
|
多原子分子的频率计算可以归结为一个数学问题:对角化矩阵方程F * R = M * R * E,其中F是力常数矩阵,M是原子质量矩阵,特征矢量R和特征值是待求的。用Matlab语言描述,就是[R,E]=eig(F,M)。 特征矢量R满足这样的关系:R^T * M * R = I,因此在不做质量加权的情况下,每个特征值Ei对应的特征矢量Ri并不是归一化的,也就无法得出正确的力常数。因此在实际应用中,往往对每个特征矢量Ri重新归一化:Li = Ri/sqrt(Ri^T * Ri)。重新归一化后的Li就是简正振动模式,而1/(Ri^T * Ri)就是该模式的约化质量mi。因此Li = sqrt(mi) * Ri。重新归一化以后,就可以计算力常数:Ki = Li^T * F * Li,或者Ki = Ei * mi |

2楼2014-05-22 02:22:40












回复此楼