| 查看: 633 | 回复: 6 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 majun04 的 36 个金币 | |||
[交流]
【求助】金币换取有离子参加的化学反应过渡态的计算实例
|
|||
|
用Gaussian03做有正价离子（氢离子或者金属离子）催化的化学反应，在寻找反应路径中的过渡态。 目前困难： 游离的氢离子和乙腈作用，氢离子和氮原子成键过渡态找不到。成键后的稳定态已经找到。但是乙腈和氢离子在同一个文件中优化怎么也不收敛，改变方法依然没有找到。因此用QST2方法找，不成功。自己设计过渡态优化结果不收敛，要么502错误，要么氢优化过程中，和氮原子成键（回到定点，并非过渡态）。 求助内容： 提供金属离子（或者氢离子）与有机物分子在同一个文件中优化的方法。在此状况下离子和分子并没有任何相互作用，仅仅是为了读取能量基点（或者方便用于QST2方法的输入文件）。 提供有氢离子或者金属离子（带正电荷的）参与的催化反应中，金属离子（或者氢离子）与反应物成键（或者断键），过渡态的寻找方法。 最好是提供具体的命令参数 调试很久依然没有成功，请高手帮忙！重金答谢！ |
回复此楼» 猜你喜欢
书籍求助：汽车市场营销理论与实务（电子版）——章小平
已经有0人回复
-
东方理-中科大联合博士生招聘
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有266人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
钴酸锂半电池小倍率容量上不去
已经有1人回复
-
吉林大学材料物理本科生求问调剂信息
已经有23人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
分子动力学模拟合作需求
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2008-11-20 23:15:43
 |
3楼2008-11-21 00:30:52
dolphin618
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1206.3
- 散金: 1756
- 帖子: 624
- 在线: 106小时
- 虫号: 211713
- 注册: 2006-03-08
- 性别: GG
- 专业: 能源化工
4楼2008-11-21 04:41:25
5楼2008-11-21 09:16:03
dolphin618
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1206.3
- 散金: 1756
- 帖子: 624
- 在线: 106小时
- 虫号: 211713
- 注册: 2006-03-08
- 性别: GG
- 专业: 能源化工
6楼2008-11-22 04:37:23
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
majun04(金币+10,VIP+0):谢谢你！
lei0736(金币+3,VIP+0):谢谢
majun04(金币+10,VIP+0):谢谢你！
lei0736(金币+3,VIP+0):谢谢
|
帮你扫了一下势能面，H+和N原子的距离从1.0扫到2.0，能量一直上升，感觉这个反应没有过渡态。 楼上说的很有道理，换成H3O+的话，应该就能找到过渡态了。这样的反应用QST2找过渡态应该不成问题的。 附上扫描势能面的输入文件供你参考： 真空中 #p b3lyp/6-31G(d,p) opt=z-matrix nosymm Title Card Required 1 1 H N 1 B1 H 2 B2 1 A1 C 3 B3 2 A2 1 D1 C 4 B4 3 A3 2 D2 H 5 B5 4 A4 3 D3 H 5 B6 4 A5 3 D4 B1 1.00000001 S 10 0.10 B2 3.20611278 B3 2.14826264 B4 1.54000000 B5 1.07000000 B6 1.07000000 A1 161.66026897 A2 9.66804955 A3 28.00778050 A4 109.47125053 A5 109.47121812 D1 -180.00000000 D2 -180.00000000 D3 119.99998499 D4 -120.00000125 溶剂中： #p b3lyp/6-31G(d,p) opt=z-matrix nosymm scrf(pcm,solvent=CH3CN,read) Title Card Required 1 1 H N 1 B1 H 2 B2 1 A1 C 3 B3 2 A2 1 D1 C 4 B4 3 A3 2 D2 H 5 B5 4 A4 3 D3 H 5 B6 4 A5 3 D4 B1 1.00000001 S 10 0.10 B2 3.20611278 B3 2.14826264 B4 1.54000000 B5 1.07000000 B6 1.07000000 A1 161.66026897 A2 9.66804955 A3 28.00778050 A4 109.47125053 A5 109.47121812 D1 -180.00000000 D2 -180.00000000 D3 119.99998499 D4 -120.00000125 radii=uff |
7楼2008-11-22 12:23:45