| 查看: 619 | 回复: 2 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】请教高斯计算中能量校正的问题
|
|||
|
初学高斯，想进行能量校正，请问我的操作是否正确？ 下面是H2O计算out文件的部分内容： 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: 1 2 3 A1 A1 B2 Frequencies -- 1613.2798 3881.7625 3990.2135 Red. masses -- 1.0826 1.0452 1.0823 Frc consts -- 1.6602 9.2791 10.1530 IR Inten -- 67.8281 12.7620 62.8282 Atom AN X Y Z X Y Z X Y Z 1 8 0.00 0.00 0.07 0.00 0.00 0.05 0.00 0.07 0.00 2 1 0.00 -0.43 -0.56 0.00 0.58 -0.40 0.00 -0.56 0.43 3 1 0.00 0.43 -0.56 0.00 -0.58 -0.40 0.00 -0.56 -0.43 ------------------- - Thermochemistry - ------------------- Temperature 298.150 Kelvin. Pressure 1.00000 Atm. Atom 1 has atomic number 8 and mass 15.99491 Atom 2 has atomic number 1 and mass 1.00783 Atom 3 has atomic number 1 and mass 1.00783 Molecular mass: 18.01056 amu. Principal axes and moments of inertia in atomic units: 1 2 3 EIGENVALUES -- 2.17514 4.15411 6.32926 X 0.00000 0.00000 1.00000 Y 1.00000 0.00000 0.00000 Z 0.00000 1.00000 0.00000 This molecule is an asymmetric top. Rotational symmetry number 2. Rotational temperatures (Kelvin) 39.81982 20.85013 13.68467 Rotational constants (GHZ): 829.71138 434.44674 285.14265 Zero-point vibrational energy 56734.4 (Joules/Mol) 13.55985 (Kcal/Mol) Vibrational temperatures: 2321.15 5584.98 5741.02 (Kelvin) Zero-point correction= 0.021609 (Hartree/Particle) Thermal correction to Energy= 0.024445 Thermal correction to Enthalpy= 0.025389 Thermal correction to Gibbs Free Energy= 0.003979 Sum of electronic and zero-point Energies= -76.413073 Sum of electronic and thermal Energies= -76.410237 Sum of electronic and thermal Enthalpies= -76.409293 Sum of electronic and thermal Free Energies= -76.430703 E (Thermal) CV S KCal/Mol Cal/Mol-Kelvin Cal/Mol-Kelvin Total 15.339 6.012 45.062 Electronic 0.000 0.000 0.000 Translational 0.889 2.981 34.608 Rotational 0.889 2.981 10.446 Vibrational 13.562 0.050 0.007 Q Log10(Q) Ln(Q) Total Bot 0.147901D-01 -1.830028 -4.213794 Total V=0 0.128648D+09 8.109403 18.672590 Vib (Bot) 0.115014D-09 -9.939250 -22.885968 Vib (V=0) 0.100042D+01 0.000181 0.000416 Electronic 0.100000D+01 0.000000 0.000000 Translational 0.300432D+07 6.477746 14.915562 Rotational 0.428032D+02 1.631476 3.756612 ***** Axes restored to original set ***** ------------------------------------------------------------------- Center Atomic Forces (Hartrees/Bohr) Number Number X Y Z ------------------------------------------------------------------- 1 8 -0.000000193 0.000000000 -0.000000136 2 1 0.000000105 0.000000000 0.000000056 3 1 0.000000088 0.000000000 0.000000080 ------------------------------------------------------------------- Cartesian Forces: Max 0.000000193 RMS 0.000000097 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Berny optimization. Internal Forces: Max 0.000000186 RMS 0.000000118 Search for a local minimum. Step number 1 out of a maximum of 2 All quantities printed in internal units (Hartrees-Bohrs-Radians) Second derivative matrix not updated -- analytic derivatives used. The second derivative matrix: R1 R2 A1 R1 0.55672 R2 -0.00588 0.55672 A1 0.02929 0.02929 0.15377 Eigenvalues --- 0.14950 0.55511 0.56260 Angle between quadratic step and forces= 90.00 degrees. Linear search not attempted -- first point. Iteration 1 RMS(Cart)= 0.00000067 RMS(Int)= 0.00000000 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 1.80999 0.00000 0.00000 0.00000 0.00000 1.80999 R2 1.80999 0.00000 0.00000 0.00000 0.00000 1.80999 A1 1.83194 0.00000 0.00000 0.00000 0.00000 1.83193 Item Value Threshold Converged? Maximum Force 0.000000 0.000015 YES RMS Force 0.000000 0.000010 YES Maximum Displacement 0.000001 0.000060 YES RMS Displacement 0.000001 0.000040 YES Predicted change in Energy=-1.293968D-13 Optimization completed. -- Stationary point found. ---------------------------- ! Optimized Parameters ! ! (Angstroms and Degrees) ! -------------------------- -------------------------- ! Name Definition Value Derivative Info. ! -------------------------------------------------------------------------------- ! R1 R(1,2) 0.9578 -DE/DX = 0.0 ! ! R2 R(1,3) 0.9578 -DE/DX = 0.0 ! ! A1 A(2,1,3) 104.9622 -DE/DX = 0.0 ! -------------------------------------------------------------------------------- GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad 1|1|UNPC-UNK|Freq|RmPW1PW91|6-311+G(d,p)|H2O1|PCUSER|14-Sep-2009|0||#N GEOM=ALLCHECK GUESS=READ SCRF=CHECK GENCHK RMPW1PW91/6-311+G(D,P) FRE Q||Title Card Required||0,1|O,-0.0952734031,0.,-0.0673325842|H,-0.0573 541316,0.,0.8897195459|H,0.8195413565,0.,-0.3510588722||Version=x86-Wi n32-G03RevB.03|State=1-A1|HF=-76.4346818|RMSD=1.110e-010|RMSF=9.664e-0 08|Dipole=0.6952636,0.,0.4913637|DipoleDeriv=-0.430751,0.,0.0776326,0. ,-0.7723616,0.,0.0776326,0.,-0.4857332,0.2777636,0.,-0.0837002,0.,0.38 61808,0.,-0.038118,0.,0.1804785,0.1529874,0.,0.0060676,0.,0.3861808,0. ,-0.0395146,0.,0.3052547|Polar=7.0700986,0.,6.5393675,-0.681324,0.,7.5 526376|PG=C02V [C2(O1),SGV(H2)]|NImag=0||0.54534068,0.,-0.00012431,-0. 11480464,0.,0.62664959,-0.04445742,0.,-0.05313401,0.04643966,0.,0.0000 6215,0.,0.,-0.00009869,0.00630978,0.,-0.54153772,0.00406631,0.,0.55722 138,-0.50088326,0.,0.16793865,-0.00198225,0.,-0.01037609,0.50286551,0. ,0.00006215,0.,0.,0.00003654,0.,0.,-0.00009869,0.10849486,0.,-0.085111 87,0.04906770,0.,-0.01568366,-0.15756256,0.,0.10079554||0.00000019,0., 0.00000014,-0.00000011,0.,-0.00000006,-0.00000009,0.,-0.00000008|||@ 进行能量校正我按以下步骤操作，请问对吗？请高手指导！谢谢。 水分子的能量校正步骤： 1.计算结果中，读取能量：HF=-76.4346818 2.读取零点能：Zero-point correction= 0.021609 3.如果不考虑校正因子：校正后的能量E=-76.4346818+0.021609=-76.413073 当考虑ZPE校正因子时：校正后的能量E=-76.4346818+0.021609*0.9580=-76.413980 [ Last edited by fegg7502 on 2009-9-15 at 23:34 ] |
回复此楼» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
物理化学论文润色/翻译怎么收费?
已经有298人回复
-
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
zeoliters
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.902
- 金币: 2088.4
- 散金: 16
- 帖子: 445
- 在线: 60.6小时
- 虫号: 772181
- 注册: 2009-05-16
- 专业: 计算化学，量子化学
2楼2009-09-15 23:28:24
3楼2009-09-16 08:03:02