计算化学公社

标题: 如何根据每个原子的质量加权笛卡尔速度得到整个分子的平动速度 [打印本页]

作者
Author: diaolanxinyu 时间: 2015-3-31 09:17
标题: 如何根据每个原子的质量加权笛卡尔速度得到整个分子的平动速度
在用Gaussian的BOMD模拟Cl 碰撞甲烷的过程，初始时Cl的动能是10ev（7.4km/s），甲烷是静止的，产物还是Cl和甲烷，下面给出的是110fs时的质量加权笛卡尔速度，
我的问题是：如何根据每个原子的质量加权笛卡尔速度得到整个分子的平动速度呢,（我主要是想通过整个分子的平动速度得到整个分子的平动能）？谢谢大家！

Summary information for step 552
Time (fs) 110.000000
EKin =    0.3644458; EPot = -500.0075671; ETot = -499.6431213 A.U.
Angular momentum (instantaneous)
JX =  0.1293873955D-14  JY = -0.3735567519D-02  JZ =  0.1081672049D-01
Jtot =  0.1144359676D-01 H-BAR;  J (Quantum Number) =  0.1309387619D-03
Total energy  -4.996431D+02  Delta-E 1.175109D-07 A.U.
Total angular momentum 1.144360D-02  Delta-A 3.870168D-15 h-bar
Cartesian coordinates: (bohr)
I= 1 X=  -1.576954159062D+01 Y=  -7.195060983615D-04 Z= 2.598979628732D-04
I= 2 X=  -1.369849043678D+01 Y=  -5.807053272838D-03 Z= 5.473875778448D-03
I= 3 X=  -1.635599005944D+01 Y=  -9.937013931357D-01 Z=  -1.723499295775D+00
I= 4 X=  -1.635269189896D+01 Y= 1.989550051811D+00 Z=  -2.363673793412D-03
I= 5 X=  -1.636471889861D+01 Y=  -9.928058858894D-01 Z= 1.721645971130D+00
I= 6 X=  -8.112000781387D+00 Y= 3.265758917910D-04 Z=  -1.254112635563D-04
MW cartesian velocity: (sqrt(amu)*bohr/sec)
C 1 X=  -6.924931935129D+14 Y= 1.536419348832D+10 Z=  -3.585213293840D+10
H 2 X=  -1.968869985520D+14 Y=  -2.987345107320D+11 Z= 2.627037425623D+11
H 3 X=  -1.531223649273D+14 Y=  -5.383258550676D+12 Z=  -9.317799028956D+12
H 4 X=  -1.530372935514D+14 Y= 1.069418920586D+13 Z=  -2.177834484579D+10
H 5 X=  -1.533672060781D+14 Y=  -5.446398917046D+12 Z= 9.407631472662D+12
Cl 6 X=  -3.071614488483D+14 Y= 6.471270915764D+10 Z=  -3.514937941433D+10

作者
Author: sobereva 时间: 2015-4-2 10:09
根据每个原子的质量，把MW cartesian velocity除以sqrt(原子质量)就可以得到每个原子的速度。然后用mv^2/2你可以得到每个原子在每个方向的动能。原子在每个方向的动能加和就是分子的动能（加和时注意运动方向相反的会抵消），然后再除以(m/2)再开根号就是分子每个方向的速度。

由于最近时间紧没时间仔细琢磨，我想应该是这样。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)