24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tiancj

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1200
散金: 78
红花: 1
帖子: 256
在线: 148.5小时
虫号: 541434
注册: 2008-04-07
专业: 原子和分子物理

[交流] 【求助】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
谢谢指点。

回复此楼

1楼 2010-08-01 15:53:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tiancj

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1200
散金: 78
红花: 1
帖子: 256
在线: 148.5小时
虫号: 541434
注册: 2008-04-07
专业: 原子和分子物理

引用回帖:

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的回帖

查看全部 11 个回答

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫

QC强帖: 17
应助: 226 (大学生)
贵宾: 5.4
金币: 2180.5
红花: 217
帖子: 1241
在线: 366.8小时
虫号: 281271
注册: 2006-09-23
专业: 宇宙学

★ ★ ★
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

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 1200
散金: 78
红花: 1
帖子: 256
在线: 148.5小时
虫号: 541434
注册: 2008-04-07
专业: 原子和分子物理

引用回帖:

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的回帖