| 查看: 923 | 回复: 3 | |||
[交流]
【求助】能量优化、频率分析中的原子限定问题
|
 我想对上述Claisen重排反应过程用高斯进行计算模拟。 首先是对反应物A、产物B进行能量优化及频率分析 # opt freq b3lyp/6-31g(d) geom=connectivity test 然后对过渡态进行能量优化及频率分析。 刚刚接触高斯，由于分子A,B都比较大，我想对图中红色标记刚性的原子进行固定，以节省计算时间。 小弟有几个问题想请较一下各位老手： 1、 这样的原子固定是否可行，会不会对结果有大的影响？ 2、如何在gjf文件中如何对这些原子的进行固定？ 3、该反应经过6元环过渡态TS，如何在GaussView中画过渡态结构，怎么画快要断裂的键，及快要形成的键？ 谢谢！ 以下是A,B两个分子的gjf文件。 分子A： %chk=D:\A.chk %nprocshared=1 # opt freq b3lyp/6-31g(d) geom=connectivity test A 0 1 C 3.58908800 -16.05509900 0.00008000 C 3.54701600 -17.38645700 0.00007000 C 4.79870200 -18.15398600 0.00004300 C 5.96693200 -17.51469000 0.00002800 C 6.00051200 -16.08917400 0.00003800 C 4.88109800 -15.37948600 0.00006400 O 7.13111900 -18.21046200 0.00000200 C 8.39010900 -17.61212100 -0.00003300 C 8.48420000 -16.27227600 -0.00002600 C 7.26949600 -15.43947000 0.00001800 C 9.67596500 -18.46575000 -0.00009200 C 10.87199600 -17.81381400 -0.00014400 C 10.89684400 -16.35481300 -0.00012600 C 9.78873100 -15.63101200 -0.00006900 O 12.20556800 -18.41206600 -0.00022200 F 12.08074400 -15.72779200 -0.00016900 F 9.89177700 -14.29462500 -0.00005800 N 7.33720400 -14.16538600 0.00002600 C 4.97010800 -13.93642800 0.00007400 C 6.17300500 -13.36602300 0.00005600 O 9.64764900 -19.92180500 -0.00009800 C 8.38756900 -20.60039400 -0.00003100 C 12.36806400 -19.83661100 -0.00026700 C 13.81953300 -20.19664200 -0.00036800 C 14.18268900 -21.47623900 -0.00050900 H 2.69148500 -15.49089100 0.00009900 H 2.61367200 -17.88972100 0.00008000 H 4.77354600 -19.21412000 0.00003500 H 4.09139200 -13.34315600 0.00009500 H 6.25713300 -12.30912900 0.00006400 H 8.55856600 -21.68606900 -0.00004300 H 7.81849000 -20.31987800 -0.89861400 H 7.81859000 -20.31988300 0.89861800 H 11.89251000 -20.28139700 0.88556600 H 11.89239800 -20.28135900 -0.88605900 H 14.55673800 -19.43366700 -0.00033400 H 13.44226400 -22.23641000 -0.00054800 H 15.21128300 -21.73628000 -0.00058500 1 2 2.0 6 1.0 26 1.0 2 3 1.0 27 1.0 3 4 2.0 28 1.0 4 5 1.0 7 1.0 5 6 2.0 10 1.0 6 19 1.0 7 8 1.0 8 9 2.0 11 1.0 9 10 1.0 14 1.0 10 18 2.0 11 12 2.0 21 1.0 12 13 1.0 15 1.0 13 14 2.0 16 1.0 14 17 1.0 15 23 1.0 16 17 18 20 1.0 19 20 2.0 29 1.0 20 30 1.0 21 22 1.0 22 31 1.0 32 1.0 33 1.0 23 24 1.0 34 1.0 35 1.0 24 25 2.0 36 1.0 25 37 1.0 38 1.0 26 27 28 29 30 31 32 33 34 35 36 37 38 分子B： %chk=D:\B.chk %nprocshared=1 # opt freq b3lyp/6-31g(d) geom=connectivity test B 0 1 C 34.43129200 -16.17151100 0.41808200 C 34.48192600 -17.50180900 0.37144200 C 35.77667300 -18.17617700 0.21409700 C 36.89139000 -17.45426400 0.11827600 C 36.82613000 -16.02743000 0.16846000 C 35.66578500 -15.40287000 0.31143500 O 38.09310400 -18.06245500 -0.02757500 C 39.28845100 -17.36463200 -0.14875200 C 39.30597400 -16.01887800 -0.09956200 C 38.04118500 -15.28208900 0.06255600 C 40.57316800 -18.18671900 -0.31326000 C 41.81468200 -17.32125700 -0.36197900 C 41.73141200 -15.88023700 -0.31395000 C 40.56581200 -15.26180600 -0.19532400 O 42.92521100 -17.83956500 -0.52011400 O 40.57055500 -18.99454100 -1.56361200 C 40.64817800 -19.17417500 0.90360900 F 42.85444900 -15.15466600 -0.39493700 F 40.55506700 -13.92226700 -0.15595800 N 38.01633600 -14.00735300 0.10726800 C 35.65199100 -13.95766300 0.35937800 C 36.80447800 -13.29879900 0.25919200 C 39.46159400 -19.91799900 -1.79061300 C 41.85218000 -20.07713300 0.84949600 C 42.07742600 -20.93500300 1.83971000 H 33.50149300 -15.67485500 0.53163300 H 33.59073600 -18.07134300 0.44778800 H 35.82475200 -19.23469600 0.17644600 H 40.67856000 -18.60833100 1.84455500 H 39.75367100 -19.80916000 0.94413100 H 34.73825300 -13.43204600 0.47319700 H 36.81226900 -12.23910100 0.29291100 H 39.61371600 -20.42581500 -2.75270400 H 38.51449300 -19.36513200 -1.82565300 H 39.41589600 -20.67228300 -0.99481000 H 42.51487600 -20.03609300 0.02345600 H 41.41052700 -20.97906900 2.66337100 H 42.92267500 -21.57515100 1.80522100 1 2 2.0 6 1.0 26 1.0 2 3 1.0 27 1.0 3 4 2.0 28 1.0 4 5 1.0 7 1.0 5 6 2.0 10 1.0 6 21 1.0 7 8 1.0 8 9 2.0 11 1.0 9 10 1.0 14 1.0 10 20 2.0 11 12 1.0 16 1.0 17 1.0 12 13 1.0 15 2.0 13 14 2.0 18 1.0 14 19 1.0 15 16 23 1.0 17 24 1.0 29 1.0 30 1.0 18 19 20 22 1.0 21 22 2.0 31 1.0 22 32 1.0 23 33 1.0 34 1.0 35 1.0 24 25 2.0 36 1.0 25 37 1.0 38 1.0 26 27 28 29 30 31 32 33 34 35 36 37 38 [ Last edited by chemozhang on 2011-4-9 at 13:20 ] |
回复此楼» 猜你喜欢
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
【补录】天大-利物浦福州联合学院 生物与生物医药科学硕士
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有196人回复
-
2026第二届先进功能材料与器件学术研讨会——探功能材料精髓·研先进器件新章！
已经有4人回复
-
中科院西安光机所——聚焦功能材料创新｜诚邀参与 2026 年度学术研讨会。
已经有1人回复
-
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有1人回复
-
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有0人回复
-
分子动力学 Plumed教程求助
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于冻结原子优化后解冻的问题
已经有4人回复
- 关于TD-DFT优化构型并计算频率的问题
已经有12人回复
- MS 32原子金属体系结构优化问题求助
已经有4人回复
- castep优化结果中各原子能量与最终能量的关系？
已经有6人回复
- 【求助】优化结构前后，怎么看原子移动多少
已经有11人回复
- 【求助】关于dmol中部分原子固定的几何优化问题
已经有10人回复
- 【求助】关于晶胞中取代一个原子的优化问题
已经有14人回复
- 【求助】TD-DFT优化激发态构型后的频率问题
已经有11人回复
- 【求助】请问关于wien2k的原子位置优化问题
已经有6人回复
» 抢金币啦！回帖就可以得到:
-
郑州大学急招1名2026级博士生
+1/474
-
⚽️赛场争冠，科研进阶 | 原位光谱，适配钙钛矿 、OPV等全流程表征
+1/86
-
Cu2O纳米线
+5/80
-
欢迎报考2027年浙大宁理化学工程（085602）专业研究生，一志愿高录取稳上岸！
+1/34
-
东南大学有机多孔功能材料团队（国家杰青团队）招收2027级推夏令营免硕士/直博生
+1/33
-
东南大学有机多孔功能材料团队（国家杰青团队） 2027级博士研究生招生
+1/29
-
湖南师范大学(211)—招收2026年入学的“申请-考核”制博士（有微纳光学经验者优先）
+1/25
-
美国普渡大学(Purdue University)生物化学系招聘博士后
+1/18
-
东北师范大学荒漠与草地生态学方向诚聘博士后青年才俊
+1/13
-
哈尔滨工业大学招收硕士研究生（欢迎环境、市政、生物、化学、农业等专业，长期有效）
+1/10
-
四川大学周加境课题组招聘博士后/博士/研究助理（生物质与藻类资源利用/自组装材料）
+1/9
-
密苏里大学生物材料合成生物学博士后招聘
+1/8
-
哈尔滨工业大学招收博士研究生（欢迎环境、生物、市政、农业、化学等专业）长期有效
+1/8
-
求助最新版ISO 5817
+1/7
-
化学识别核酸的综述投稿
+1/6
-
哈尔滨工业大学（深圳）赵怡潞课题组诚招博士后、2027学年博士生
+1/5
-
上海大学微电子学院杨军教授团队招聘带编专任教师
+1/4
-
中科院深圳先进院成会明院士/郭鑫团队招聘钠离子电池方向博士后
+1/3
-
上海大学微电子学院杨军教授团队招聘带编专任教师
+1/2
-
招科研助理，提前招27年博士生
+1/1
★ ★ ★
chemozhang(金币+2): 多谢！还请较如何画过渡态的结构呢？ 2011-04-09 16:33:58
小红豆(金币+3): 感谢交流 常来 2011-04-11 09:50:58
chemozhang(金币+2): 多谢！还请较如何画过渡态的结构呢？ 2011-04-09 16:33:58
小红豆(金币+3): 感谢交流 常来 2011-04-11 09:50:58
|
我们版的置顶帖儿GaussianFAQ会对你有帮助的。 http://muchong.com/bbs/viewthread.php?tid=1566146 你在这里可以看到（1）怎么固定某部分原子 也可以看到（2）什么是过渡态，什么是中间态 但是我不建议你固定原子，与其那样还不如直接模型化，把分子结构调小呢。 |
2楼2011-04-09 16:14:44
3楼2011-04-09 20:06:37
4楼2013-08-20 11:41:04