应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Multiwfn » 转动配分函数的惯性矩
【调剂】北京石油化工学院2024年16个专业接受调剂
51/1返回列表
查看: 1389  |  回复: 4
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lingao

金虫 (小有名气)

[求助] 转动配分函数的惯性矩

采用Multiwfn软件提取高斯输出文件中转动配分函数的惯性矩阵以计算配分函数转动所占的比例值发现在out文件中如下
CH3OCl
Principal axes and moments of inertia in atomic units:
                           1         2         3
     Eigenvalues --    44.62145 294.08298 326.93014
           X            0.99948  -0.03232   0.00000
           Y            0.03232   0.99948   0.00001
           Z            0.00000  -0.00001   1.00000
而在Multiwfn中提取出的信息为
************ Moments of inertia tensor (amu*Angstrom^2) ************
             1             2             3
     1  0.125817D+02 -0.231982D+01  0.384147D-03
     2 -0.231982D+01  0.827644D+02  0.102662D-03
     3  0.384147D-03  0.102662D-03  0.920486D+02
The moments of inertia relative to X,Y,Z axes (amu*Angstrom^2):
  0.12581740E+02  0.82764360E+02  0.92048584E+02

请问这两者之间是如何换算的？在计算中可否直接使用
       X            0.99948  -0.03232   0.00000
           Y            0.03232   0.99948   0.00001
           Z            0.00000  -0.00001   1.00000
来计算det（I）
回复此楼

» 本帖已获得的红花（最新10朵）

Zwaard

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
耐得了寂寞，拥得住繁华。。。
1楼2013-07-05 18:51:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
gkf高: 金币+3 2013-07-05 21:06:10
Multiwfn的结果不是从高斯的输出文件中提取的，而是根据分子坐标和元素由Multiwfn自己计算的。
你所列的两套数据间无法对比。应当用Multiwfn输出的************ Principal axes ************后面的数据来与之对比。Principal moments of inertia (amu*Angstrom^2):下面的三个值与高斯输出的Eigenvalues后面的三个值相对应，转换关系是除以0.5292^2
赞一下 回复此楼
2楼2013-07-05 20:23:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingao

金虫 (小有名气)
引用回帖:
2楼: Originally posted by sobereva at 2013-07-05 20:23:03
Multiwfn的结果不是从高斯的输出文件中提取的，而是根据分子坐标和元素由Multiwfn自己计算的。
你所列的两套数据间无法对比。应当用Multiwfn输出的************ Principal axes ************后面的数据来与之对比。 ...

感谢您的回答，
请问可否详细说明在multiwfn计算出的这些数值所代表含义的区别 ************ Moments of inertia tensor (amu*Angstrom^2) ************
             1             2             3
     1  0.125817D+02 -0.231982D+01  0.384147D-03
     2 -0.231982D+01  0.827644D+02  0.102662D-03
     3  0.384147D-03  0.102662D-03  0.920486D+02
The moments of inertia relative to X,Y,Z axes (amu*Angstrom^2):
  0.12581740E+02  0.82764360E+02  0.92048584E+02
Rotational constant relative to X,Y,Z axes (GHz):
     40.16765550      6.10623950      5.49035071
************ Principal axes ************
             1             2             3
     1  0.999455D+00 -0.330000D-01  0.454441D-05
     2  0.330000D-01  0.999455D+00  0.992219D-05
     3 -0.486936D-05 -0.976682D-05  0.100000D+01
Principal moments of inertia (amu*Angstrom^2):
  0.12505144E+02  0.82840956E+02  0.92048584E+02
Rotational constant relative to principal axes (GHz):
     40.41368817      6.10059359      5.49035071

我完全迷糊了，在计算det（I）-determinant of moment of inertia 时选取哪组数据最佳？  非常感谢。。
赞一下 回复此楼
耐得了寂寞，拥得住繁华。。。
3楼2013-07-06 10:17:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gkf高: 金币+3 2013-07-06 22:49:18
lingao: 金币+20, ★★★★★最佳答案, 非常感谢您耐心的回复，讲解真是太详细了，非常感谢！！ 2013-07-08 08:25:59
我不太清楚你算det(I)的目的，如果有公式涉及到这个量，你再多看看相关的书看看有没有内容具体说明此问题。
Moments of inertia tensor是相对于X,Y,Z笛卡尔坐标轴来计算的转动惯量张量。The moments of inertia relative to X,Y,Z axes显然就是这个矩阵的三个对角元，即相对于X,Y,Z轴的转动惯量。
Principal axes（主轴）是转动惯量张量的本征向量，这个矩阵中每一列对应一个本征向量。Principal moments of inertia也就是转动惯量张量的本征值，物理意义是相对于那三个主轴的转动惯量。

高斯输出的那个并不是转动惯量张量，而是笛卡尔坐标下的转动惯量张量的三个本征向量（主轴）和三个本征值（相对于主轴的转动惯量），和Multiwfn输出的************ Principal axes ************
后面的内容是相对应的。
赞一下 回复此楼
4楼2013-07-06 22:44:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Zwaard

铜虫 (初入文坛)
送红花一朵
请问虫友们与各位大神，
1. 楼上Sob提到的X,Y,Z笛卡尔坐标是高斯输出中的哪一套笛卡尔坐标系？ Standard Orientation ?
2. 楼主所提到的
************ Principal axes ************
             1             2             3
     1  0.999455D+00 -0.330000D-01  0.454441D-05
     2  0.330000D-01  0.999455D+00  0.992219D-05
     3 -0.486936D-05 -0.976682D-05  0.100000D+01
是否指的是三个主轴与原有的X,Y,Z笛卡尔坐标的转换矩阵 ？
如果是，为什么这个矩阵总是非常非常的接近单位矩阵 ？
赞一下 回复此楼
5楼2015-07-13 10:16:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lingao 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 厦门大学，课题组新增一个博士名额，085803核能工程 +5 wanghx_xmu 2024-04-17 9/450 2024-04-19 16:58 by 1532341449
[教师之家] 美国教授开除中国研究生：我就不该录取你！ +7 苏东坡二世 2024-04-15 7/350 2024-04-19 16:20 by feng6531
[找工作] 杭州国企和浙江高校如何选择？ +15 restart2024 2024-04-15 21/1050 2024-04-19 16:05 by l419110028
[高分子] 聚酰胺650与环氧树脂e44固化 +3 yindingxin 2024-04-15 3/150 2024-04-19 13:55 by weilingdun
[基金申请] 特别资助审核状态 +7 pantray 2024-04-17 7/350 2024-04-19 13:51 by Breezy恋
[考研] 298求调剂 +8 孙大大@ 2024-04-17 8/400 2024-04-19 10:25 by juncaixing
[基金申请] 下雨了 +13 zju2000 2024-04-16 19/950 2024-04-19 09:24 by duxin_30
[论文投稿] 求助 100+3 guhangyu 2024-04-15 3/150 2024-04-19 08:48 by bobvan
[考研] 267求调剂 +3 工科材料267 2024-04-15 3/150 2024-04-18 19:45 by lature00
[高分子] 聚氨酯丙烯酸酯预聚体合成 +3 15692385990 2024-04-14 6/300 2024-04-18 19:43 by zhang1991
[考研] 275求调剂 +8 调剂0856 2024-04-14 10/500 2024-04-17 23:23 by 华北刘兵
[无机/物化] H3PO2在130度附近分解为磷化氢（剧毒），怎样防护，求指点。 1(金币+5) +3 大力2010 2024-04-13 9/450 2024-04-17 21:55 by 大力2010
[考博] 25年申博求助，射频方向 +4 SherlockAH 2024-04-15 10/500 2024-04-16 18:06 by ljysudaee
[考研] 294求调剂 +3 694062003 2024-04-15 4/200 2024-04-16 15:01 by 邹邹哈哈
[考研] 312求调剂 +3 Lauhalo 2024-04-15 3/150 2024-04-16 10:16 by 19862091
[考研] 食品工程考数学283还能有什么调剂选择 +7 宇宙无敌美少女 2024-04-12 8/400 2024-04-16 10:11 by cyclodextrin726
[考博] 博士申请 +5 啦啦啦3737 2024-04-13 9/450 2024-04-16 07:07 by 1147wzzzz
[论文投稿] with efitor 越久是不是越容易拒稿。我的已经一个多月了 +5 lizhengke06 2024-04-14 5/250 2024-04-15 18:33 by jonewore
[考研] 334求调剂 +4 学药救人 2024-04-14 4/200 2024-04-15 15:05 by hunanzang
[考研] 专硕调剂招生 +3 电致发光 2024-04-15 4/200 2024-04-15 07:34 by ashorewmj
信息提示
关闭
请填处理意见
关闭