| 查看: 461 | 回复: 2 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】过渡态计算后如何计算activation和barrier energy
|
|||
|
简单的SiH2+H2->SiH4的反应。用QST2算法。 运算完后，如何分析输出文件得到activation和barrier energy？ 部分输出文件如下： Berny optimization. Internal Forces: Max 0.000000190 RMS 0.000000120 Search for a local minimum. Step number 3 out of a maximum of 20 All quantities printed in internal units (Hartrees-Bohrs-Radians) Update second derivatives using D2CorX and points 1 2 3 Trust test= 1.05D+00 RLast= 9.35D-03 DXMaxT set to 3.00D-01 The second derivative matrix: R1 R2 R3 R4 A1 R1 0.18612 R2 0.00359 0.18612 R3 0.00359 0.00359 0.18612 R4 0.00359 0.00359 0.00359 0.18612 A1 0.00000 0.00000 0.00000 0.00000 0.16000 A2 0.00000 0.00000 0.00000 0.00000 0.00000 A3 0.00000 0.00000 0.00000 0.00000 0.00000 A4 0.00000 0.00000 0.00000 0.00000 0.00000 A5 0.00000 0.00000 0.00000 0.00000 0.00000 A6 0.00000 0.00000 0.00000 0.00000 0.00000 A2 A3 A4 A5 A6 A2 0.16000 A3 0.00000 0.16000 A4 0.00000 0.00000 0.16000 A5 0.00000 0.00000 0.00000 0.16000 A6 0.00000 0.00000 0.00000 0.00000 0.16000 Eigenvalues --- 0.16000 0.16000 0.16000 0.16000 0.16000 Eigenvalues --- 0.18253 0.18253 0.18253 0.196871000.00000 RFO step: Lambda= 0.00000000D+00. Quartic linear search produced a step of -0.00020. Iteration 1 RMS(Cart)= 0.00000061 RMS(Int)= 0.00000000 Variable Old X -DE/DX Delta X Delta X Delta X New X (Linear) (Quad) (Total) R1 2.81577 0.00000 0.00000 0.00000 0.00000 2.81578 R2 2.81577 0.00000 0.00000 0.00000 0.00000 2.81578 R3 2.81577 0.00000 0.00000 0.00000 0.00000 2.81578 R4 2.81577 0.00000 0.00000 0.00000 0.00000 2.81578 A1 1.91063 0.00000 0.00000 0.00000 0.00000 1.91063 A2 1.91063 0.00000 0.00000 0.00000 0.00000 1.91063 A3 1.91063 0.00000 0.00000 0.00000 0.00000 1.91063 A4 1.91063 0.00000 0.00000 0.00000 0.00000 1.91063 A5 1.91063 0.00000 0.00000 0.00000 0.00000 1.91063 A6 1.91063 0.00000 0.00000 0.00000 0.00000 1.91063 Item Value Threshold Converged? Maximum Force 0.000000 0.000450 YES RMS Force 0.000000 0.000300 YES Maximum Displacement 0.000001 0.001800 YES RMS Displacement 0.000001 0.001200 YES Predicted change in Energy=-3.652736D-13 Optimization completed. -- Stationary point found. ---------------------------- ! Optimized Parameters ! ! (Angstroms and Degrees) ! -------------------------- -------------------------- ! Name Definition Value Derivative Info. ! -------------------------------------------------------------------------------- ! R1 R(1,2) 1.49 -DE/DX = 0.0 ! ! R2 R(1,3) 1.49 -DE/DX = 0.0 ! ! R3 R(1,4) 1.49 -DE/DX = 0.0 ! ! R4 R(1,5) 1.49 -DE/DX = 0.0 ! ! A1 A(2,1,3) 109.4712 -DE/DX = 0.0 ! ! A2 A(2,1,4) 109.4712 -DE/DX = 0.0 ! ! A3 A(2,1,5) 109.4712 -DE/DX = 0.0 ! ! A4 A(3,1,4) 109.4712 -DE/DX = 0.0 ! ! A5 A(3,1,5) 109.4712 -DE/DX = 0.0 ! ! A6 A(4,1,5) 109.4712 -DE/DX = 0.0 ! -------------------------------------------------------------------------------- GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad Input orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 14 0 4.594300 0.447240 3.685398 2 1 0 5.090964 -0.957593 3.685398 3 1 0 5.090990 1.149648 4.902014 4 1 0 5.090990 1.149648 2.468782 5 1 0 3.104256 0.447259 3.685398 --------------------------------------------------------------------- Distance matrix (angstroms): 1 2 3 4 5 1 Si 0.000000 2 H 1.490044 0.000000 3 H 1.490044 2.433232 0.000000 4 H 1.490044 2.433232 2.433232 0.000000 5 H 1.490044 2.433232 2.433232 2.433232 0.000000 Stoichiometry H4Si Framework group TD[O(Si),4C3(H)] Deg. of freedom 1 Full point group TD Largest Abelian subgroup D2 NOp 4 Largest concise Abelian subgroup D2 NOp 4 Standard orientation: --------------------------------------------------------------------- Center Atomic Atomic Coordinates (Angstroms) Number Number Type X Y Z --------------------------------------------------------------------- 1 14 0 0.000000 0.000000 0.000000 2 1 0 0.860277 0.860277 0.860277 3 1 0 -0.860277 -0.860277 0.860277 4 1 0 -0.860277 0.860277 -0.860277 5 1 0 0.860277 -0.860277 -0.860277 --------------------------------------------------------------------- Rotational constants (GHZ): 84.6964445 84.6964445 84.6964445 ********************************************************************** Population analysis using the SCF density. ********************************************************************** Orbital symmetries: Occupied (A1) (A1) (T2) (T2) (T2) (A1) (T2) (T2) (T2) Virtual (A1) (T2) (T2) (T2) (T2) (T2) (T2) (A1) (A1) (T2) (T2) (T2) (T2) (T2) (T2) (A1) The electronic state is 1-A1. Alpha occ. eigenvalues -- -68.78060 -6.14880 -4.25364 -4.25364 -4.25364 Alpha occ. eigenvalues -- -0.74612 -0.48377 -0.48377 -0.48377 Alpha virt. eigenvalues -- 0.06637 0.08818 0.08818 0.08818 0.18421 Alpha virt. eigenvalues -- 0.18421 0.18421 0.31273 0.46991 0.48975 Alpha virt. eigenvalues -- 0.48975 0.48975 1.21628 1.21628 1.21628 Alpha virt. eigenvalues -- 1.21664 Condensed to atoms (all electrons): 1 2 3 4 5 1 Si 11.804787 0.357348 0.357348 0.357348 0.357348 2 H 0.357348 0.894691 -0.020195 -0.020195 -0.020195 3 H 0.357348 -0.020195 0.894691 -0.020195 -0.020195 4 H 0.357348 -0.020195 -0.020195 0.894691 -0.020195 5 H 0.357348 -0.020195 -0.020195 -0.020195 0.894691 Mulliken atomic charges: 1 1 Si 0.765821 2 H -0.191455 3 H -0.191455 4 H -0.191455 5 H -0.191455 Sum of Mulliken charges= 0.00000 Atomic charges with hydrogens summed into heavy atoms: 1 1 Si 0.000000 2 H 0.000000 3 H 0.000000 4 H 0.000000 5 H 0.000000 Sum of Mulliken charges= 0.00000 Electronic spatial extent (au): Charge= 0.0000 electrons Dipole moment (field-independent basis, Debye): X= 0.0000 Y= 0.0000 Z= 0.0000 Tot= 0.0000 Quadrupole moment (field-independent basis, Debye-Ang): XX= -17.5500 YY= -17.5500 ZZ= -17.5500 XY= 0.0000 XZ= 0.0000 YZ= 0.0000 Traceless Quadrupole moment (field-independent basis, Debye-Ang): XX= 0.0000 YY= 0.0000 ZZ= 0.0000 XY= 0.0000 XZ= 0.0000 YZ= 0.0000 Octapole moment (field-independent basis, Debye-Ang**2): XXX= 0.0000 YYY= 0.0000 ZZZ= 0.0000 XYY= 0.0000 XXY= 0.0000 XXZ= 0.0000 XZZ= 0.0000 YZZ= 0.0000 YYZ= 0.0000 XYZ= -2.1944 Hexadecapole moment (field-independent basis, Debye-Ang**3): XXXX= -45.7252 YYYY= -45.7252 ZZZZ= -45.7252 XXXY= 0.0000 XXXZ= 0.0000 YYYX= 0.0000 YYYZ= 0.0000 ZZZX= 0.0000 ZZZY= 0.0000 XXYY= -16.0642 XXZZ= -16.0642 YYZZ= -16.0642 XXYZ= 0.0000 YYXZ= 0.0000 ZZXY= 0.0000 N-N= 2.119282818723D+01 E-N=-7.362834439632D+02 KE= 2.913222622995D+02 Symmetry A KE= 2.131983161207D+02 Symmetry B1 KE= 2.604131539294D+01 Symmetry B2 KE= 2.604131539294D+01 Symmetry B3 KE= 2.604131539294D+01 |
回复此楼» 猜你喜欢
中教金源的流动池
已经有0人回复
-
RSC ADV状态问题
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有97人回复
-
湖南大学招收锂离子电池、离子液体、高分子与凝胶材料、电催化方向博士生多名
已经有8人回复
-
南京大学蔡亮课题组诚招2026年申请-考核制博士生2-3名（电解水制氢，XAFS谱学等）
已经有32人回复
-
求助各位大佬配位化学学习教材！！！
已经有0人回复
-
大连海事大学国家级人才团队2026年博士研究生招生启事
已经有10人回复
-
内蒙古大学青年长江学者王蕾教授课题组2026年博士招聘，2-3人
已经有11人回复
-
多相催化-2026年CSC博士及洪堡合作博士后招生！
已经有0人回复
-
中国地质大学（武汉）—国家级青年人才杨明教授组-招收博士-新能源材料化学及催化材料
已经有28人回复
-
求光催化产过氧化氢的PPT ！！！毕业答辩需要
已经有1人回复
wuhanhgf2002
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 323.8
- 帖子: 435
- 在线: 43.1小时
- 虫号: 121302
- 注册: 2005-12-03
- 性别: GG
- 专业: 有机合成
3楼2009-06-05 20:00:00
冬天里的骄阳
铁杆木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 7179.7
- 散金: 10
- 帖子: 387
- 在线: 318小时
- 虫号: 439463
- 注册: 2007-10-13
- 性别: GG
- 专业: 理论和计算化学
2楼2009-06-05 10:54:07