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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】gaussian中简正坐标的输出
111/2返回列表12下一页
查看: 2007  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tiancj

木虫 (小有名气)

[交流] 【求助】gaussian中简正坐标的输出 已有2人参与

在一篇文章《Vibrational Analysis in Gaussian》中说到“The purpose of this document is to describe how Gaussian calculates the reduced mass, frequencies, force constants, and normal coordinates which are printed out at the end of a frequency calculation.”从这句话看normal coordinates应该在计算频率的时候打印出来了啊，可是不知道哪地方给的是normal coordinates？以H2o为例，计算频率时的部分输出为：
(Enter /home/gau03/g03/l716.exe)
Dipole        =-3.44478646D-17-3.96506834D-17-6.72452773D-01
Polarizability= 4.00613449D-02 3.18281559D-16 5.50808430D+00
                 1.58802971D-16 1.97181875D-16 2.56587227D+00
HyperPolar    =-3.59200065D-17-6.27316867D-18-2.96845704D-16
                -1.10949459D-14 7.62420726D-02-3.00230546D-16
                -7.01237208D+00-1.56024452D-16 2.91713493D-15
                -4.11420102D+00
Full mass-weighted force constant matrix:
Low frequencies ---   -6.0514   -1.0843   -0.0019   -0.0012   -0.0009    4.5788
Low frequencies --- 2170.0372 4140.0019 4391.0737
Diagonal vibrational polarizability:
        0.0000000       0.0418305       0.1108974
Harmonic frequencies (cm**-1), IR intensities (KM/Mole), Raman scattering
activities (A**4/AMU), depolarization ratios for plane and unpolarized
incident light, reduced masses (AMU), force constants (mDyne/A),
and normal coordinates:（这里也显示normal coordinates，但是不知道怎么看，我怎么找不到啊？
                     1                      2                      3
                    A1                     A1                     B2
Frequencies --  2170.0372              4140.0019              4391.0737
Red. masses --     1.0785                 1.0491                 1.0774
Frc consts  --     2.9923                10.5941                12.2392
IR Inten    --     7.2379                44.2862                29.9710
Atom AN      X      Y      Z        X      Y      Z        X      Y      Z
   1   1     0.00   0.45   0.54     0.00  -0.57   0.42     0.00   0.54  -0.45
   2   1     0.00  -0.45   0.54     0.00   0.57   0.42     0.00   0.54   0.45
   3   8     0.00   0.00  -0.07     0.00   0.00  -0.05     0.00  -0.07   0.00
谢谢指点。
回复此楼

» 收录本帖的淘帖专辑推荐

量化 Gaussian

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-08-01 15:53:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
★ ★ ★
tiancj(金币+5): 2010-08-01 16:21:00
zhou2009(金币+3): 2010-09-23 16:43:18
Atom AN      X      Y      Z        X      Y      Z        X      Y      Z
   1   1     0.00   0.45   0.54     0.00  -0.57   0.42     0.00   0.54  -0.45
   2   1     0.00  -0.45   0.54     0.00   0.57   0.42     0.00   0.54   0.45
   3   8     0.00   0.00  -0.07     0.00   0.00  -0.05     0.00  -0.07   0.00
这个就是。例如2号原子的X,Y,Z坐标构成第1种简正坐标的系数为0.00 -0.45 0.54。
你把这里每个振动模式的3*3矩阵排列成9*1就容易看出简正坐标与原先3N坐标的关系了。如第一个简正坐标：
0.00
0.45
0.54
0.00
-0.45
0.54
0.00
0.00
-0.07

[ Last edited by sobereva on 2010-8-1 at 16:18 ]
赞一下(1人) 回复此楼
2楼2010-08-01 16:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiancj

木虫 (小有名气)
引用回帖:
Originally posted by sobereva at 2010-08-01 16:04:47:
Atom AN      X      Y      Z        X      Y      Z        X      Y      Z
   1   1     0.00   0.45   0.54     0.00  -0.57   0.42     0.00   0.54  -0.45
   2   1     0.00  -0.45   0.54     0.00   ...

这里是简正坐标吗？这不是振动模式1下振动时原子2在x，y，z方向的变化量吗？感觉这是笛卡尔坐标啊?!谢谢你的帮助!
赞一下 回复此楼
3楼2010-08-01 16:26:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
★
tiancj(金币+5): 2010-08-01 23:21:52
heyo_123(金币+1):鼓励交流！！！ 2010-08-02 08:26:03
简正振动坐标和原始3N个笛卡尔坐标是一个线性变换的关系，前者可写为后者的线性组合。每个振动模式对应一简正坐标，由于简正坐标在每个原子笛卡尔坐标上都有分量，所以振动自然会体现在原子笛卡尔坐标的变化上。
赞一下(1人) 回复此楼
4楼2010-08-01 17:01:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiancj

木虫 (小有名气)
引用回帖:
Originally posted by sobereva at 2010-08-01 17:01:37:
简正振动坐标和原始3N个笛卡尔坐标是一个线性变换的关系，前者可写为后者的线性组合。每个振动模式对应一简正坐标，由于简正坐标在每个原子笛卡尔坐标上都有分量，所以振动自然会体现在原子笛卡尔坐标的变化上。

谢谢！那这里的0.00 -0.45 0.54等指的是前者（简正坐标）还是后者（笛卡尔坐标）啊？输出结果里能不能前者后者都看到呢？或者是看到具体怎么变换的吗？
赞一下 回复此楼
5楼2010-08-01 23:34:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
heyo_123(金币+1):鼓励交流！！！ 2010-08-02 08:26:11
引用回帖:
Originally posted by tiancj at 2010-08-01 23:34:49:



谢谢！那这里的0.00 -0.45 0.54等指的是前者（简正坐标）还是后者（笛卡尔坐标）啊？输出结果里能不能前者后者都看到呢？或者是看到具体怎么变换的吗？

是简正坐标在笛卡尔坐标上的分量
你最好先了解了解基、坐标变换的概念，线数书上一般就有相关知识。
赞一下(2人) 回复此楼
6楼2010-08-01 23:41:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiancj

木虫 (小有名气)
引用回帖:
Originally posted by sobereva at 2010-08-01 23:41:03:


是简正坐标在笛卡尔坐标上的分量
你最好先了解了解基、坐标变换的概念，线数书上一般就有相关知识。

谢谢！你说的是不是这个意思？Q=∑ Lq  Q简正坐标，L变换矩阵，q笛卡尔质量加权坐标，这样的话知道了分量q，还需要知道L才能得到简正坐标对吗？如果是这样L怎么得到啊？
赞一下 回复此楼
7楼2010-08-02 09:59:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
heyo_123(金币+1):鼓励交流 2010-08-02 20:29:27
若不考虑质权、笛卡尔-内坐标转换这些细节问题，L就是Hessian矩阵的本征向量矩阵，也就是Vibrational Analysis in Gaussian的公式7。

关于变换关系的细节问题在Vibrational Analysis in Gaussian里面都详细说明。
赞一下(2人) 回复此楼
8楼2010-08-02 18:24:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiancj

木虫 (小有名气)
引用回帖:
Originally posted by sobereva at 2010-08-02 18:24:03:
若不考虑质权、笛卡尔-内坐标转换这些细节问题，L就是Hessian矩阵的本征向量矩阵，也就是Vibrational Analysis in Gaussian的公式7。

关于变换关系的细节问题在Vibrational Analysis in Gaussian里面都详细说明。

谢谢！
回复此楼
9楼2010-08-07 00:15:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiancj

木虫 (小有名气)
引用回帖:
Originally posted by sobereva at 2010-08-01 16:04:47:
Atom AN      X      Y      Z        X      Y      Z        X      Y      Z
   1   1     0.00   0.45   0.54     0.00  -0.57   0.42     0.00   0.54  -0.45
   2   1     0.00  -0.45   0.54     0.00   ...

麻烦再确认一下这里的数字0.00 0.45等是指的Q=∑ Lq 中的q吧？谢谢！
赞一下 回复此楼
10楼2010-08-17 20:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tiancj 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +11 医学老男孩 2026-08-13 29/1450 2026-08-14 00:04 by 大快人心啦
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
信息提示
关闭
请填处理意见
关闭