| 查看: 880 | 回复: 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 ] |
回复此楼» 猜你喜欢
求助文献
已经有0人回复
-
锂电池正极材料前驱体法制备正极材料
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有133人回复
-
招计算化学博士，要求本或研的学位计算机相关，英语六级，第一作者SCI或第一专利授权
已经有0人回复
-
深圳理工大学-湖南大学招收电化学方向联培博士生1名
已经有1人回复
-
求助文献
已经有3人回复
-
双一流高效天津工业大学博士第二批招生
已经有11人回复
-
请问一下大家，那里有LaSrCoO3（钙钛矿相，纯相）样品呀，想要购买几公斤
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于冻结原子优化后解冻的问题
已经有4人回复
- 关于TD-DFT优化构型并计算频率的问题
已经有12人回复
- MS 32原子金属体系结构优化问题求助
已经有4人回复
- castep优化结果中各原子能量与最终能量的关系？
已经有6人回复
- 【求助】优化结构前后，怎么看原子移动多少
已经有11人回复
- 【求助】关于dmol中部分原子固定的几何优化问题
已经有10人回复
- 【求助】关于晶胞中取代一个原子的优化问题
已经有14人回复
- 【求助】TD-DFT优化激发态构型后的频率问题
已经有11人回复
- 【求助】请问关于wien2k的原子位置优化问题
已经有6人回复
» 抢金币啦！回帖就可以得到:
-
坐标南京，诚租女友
+5/9321
-
上海嘉定周边俩小时车程
+5/170
-
东南大学废弃碳资源利用团队（杰青/海优团队）27级博士研究生和科研助理招生
+1/127
-
坐标广州，征女友
+2/124
-
2026年江西师范大学流动化学与工程研究院陈芬儿院士团队化学博士招生
+3/68
-
希望论文顺利接受
+1/29
-
智合健物课题组2026年博士生招生(5月6日—2026年5月13日)
+1/22
-
坐标南京
+1/20
-
南京邮电大学李巍教授招收2026博士生（5月12日前有效）
+1/19
-
中国农业大学资源与环境学院 - 招聘杰出人才、青年研究员、优秀人才和博士后
+1/18
-
类器官/器官芯片 有编制研究员
+1/11
-
吉林大学任雷课题组招机器人和触觉传感方向2026年入学博士生
+1/9
-
香港理工大学蒲俊宏课题组诚聘博士生、研究助理及博士后
+1/7
-
重庆医科大学感染病教育部重点实验室唐开福课题组2026年拟招学术型博士研究生一名
+1/5
-
港科大广州 - 量子计算全额奖学金博士招聘（申请制，2026和2027均有名额）
+1/4
-
成都中医药大学药学院 谢恬教授招收2026级药物化学，有机合成方向的博士生
+1/4
-
UIO-66-NH2/RGO如何碳化，求助！！！！！
+1/3
-
美国Vanderbilt University博士后招聘
+1/3
-
双一流高校南林化学工程学院-国家级青年人才团队招2026级博士（最后一波）
+1/1
-
东北林业大学（211，双一流）招收高分子、化学、复合材料、生物质材料等方向博士生
+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