应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 分子结构优化
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 2310  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

rock0927

铁虫 (正式写手)

[求助] 分子结构优化

将两个反应物画在gaussian view里，然后进行结构优化可以吗？我计算的时候总是出现错误，求指导。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-06-12 14:48:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

rock0927(金币+2): 谢谢 2011-06-12 19:10:46
将两个反应物画在gaussian view里，然后进行结构优化是可以的。它是反应的初始起点。

计算的时候总是出现错误，因没有贴出相关文件，可能性太多，很难说是什么原因。

首先两个反应物位置放得不对。
其次两个反应物放在一起，它们有各种作用，如果仅是范德华力，要考虑包含电子相关的计算方法，如G09中的b97d之类。
赞一下 回复此楼
2楼2011-06-12 15:49:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rock0927

铁虫 (正式写手)
引用回帖:
Originally posted by zhou2009 at 2011-06-12 15:49:32:
将两个反应物画在gaussian view里，然后进行结构优化是可以的。它是反应的初始起点。

计算的时候总是出现错误，因没有贴出相关文件，可能性太多，很难说是什么原因。

首先两个反应物位置放得不对。
其次两 ...

非常感谢你的回答，我的反应是丙烷碳正离子与丙烯的反应，碳正离子攻击双键碳原子，要算反应的过渡态，现在是想将反应物进行结构优化，最后出现的错误是 link died, 没有优化成功。
赞一下 回复此楼
3楼2011-06-12 19:10:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rock0927

铁虫 (正式写手)
高人再给我指点一下啊
赞一下 回复此楼
4楼2011-06-13 08:58:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)

独行者

【答案】应助回帖

rock0927(金币+1): 谢谢 2011-06-13 11:10:43
构型貌似没搭错，把输出文件最后的具体错误信息贴出来呀。link died太抽象了。
赞一下 回复此楼
轮回的过往，注定的未来；独孤的行者，孤独的路途……
5楼2011-06-13 10:24:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice625

铜虫 (正式写手)

【答案】应助回帖

不能放在一起吧，反应物应该单独优化，把优化好的两个结构放在一起做的应该是过渡态优化

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
没有一颗心会因为追求自己的梦想而受伤！
6楼2011-06-13 10:40:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rock0927

铁虫 (正式写手)
引用回帖:
Originally posted by 独孤狼 at 2011-06-13 10:24:22:
构型貌似没搭错，把输出文件最后的具体错误信息贴出来呀。link died太抽象了。

我输入的文件是
%chk=E:\gaussian\calculation\c3c3\1.chk
%mem=6MW

#T RHF/6-31G(d) Opt Test

reactant

0 2
C                  0.06735864   -0.76970904   -0.62111435
H                 -0.11569932   -0.62705807   -1.66564318
H                 -0.20598915    0.11615528   -0.08686216
H                  1.10616345   -0.97372046   -0.46566688
C                 -0.77086082   -1.95637869   -0.11042537
H                 -0.71284328   -2.24678131    0.91777708
C                 -1.58172344   -2.63410329   -0.95881104
H                 -2.07919676   -3.41798387   -0.42687915
H                 -2.30845274   -1.96437347   -1.36896530
H                  0.40907944   -3.59327904   -2.37308325
C                 -0.03634820   -3.09853305   -1.53539559
H                 -0.98055444   -2.68684555   -1.82502982
H                  0.60814395   -2.31241579   -1.20140527
C                 -0.24606686   -4.11041246   -0.39359082
H                  0.69813904   -4.52210088   -0.10395680
C                 -1.17365424   -5.24183320   -0.87428693
H                 -1.31781169   -5.94590659   -0.08157054
H                 -2.11852062   -4.83037902   -1.16209415
H                 -0.72912263   -5.73533065   -1.71318616

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 7 2.0
6
7 8 1.0 9 1.0
8
9
10 11 1.0
11 12 1.0 13 1.0 14 1.0
12
13
14 15 1.0 16 1.0
15
16 17 1.0 18 1.0 19 1.0
17
18
19

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 7 2.0
6
7 8 1.0 9 1.0
8
9
10 11 1.0
11 12 1.0 13 1.0 14 1.0
12
13
14 15 1.0 16 1.0
15
16 17 1.0 18 1.0 19 1.0
17
18
19

最后出现的结果是：
Rotational constants (GHZ):      8.2517924      2.2549388      2.1280543
   116 basis functions,   220 primitive gaussians,   116 cartesian basis functions
    25 alpha electrons       24 beta electrons
       nuclear repulsion energy       266.0872330984 Hartrees.
NAtoms=   19 NActive=   19 NUniq=   19 SFac= 1.00D+00 NAtFMM=   60 Big=F
RHF requested, but multiplicity is not 1.
Error termination via Lnk1e in d:\gaussian\l401.exe at Mon Jun 13 11:08:26 2011.
Job cpu time:  0 days  0 hours  0 minutes  1.0 seconds.
File lengths (MBytes):  RWF=     12 Int=      0 D2E=      0 Chk=      4 Scr=      1
赞一下 回复此楼
7楼2011-06-13 11:12:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowerge

禁虫 (小有名气)
★ ★
zhou2009(金币+2): 2011-06-13 12:04:08
rock0927(金币+1): 谢谢 2011-06-13 12:29:19
本帖内容被屏蔽
8楼2011-06-13 11:17:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rock0927

铁虫 (正式写手)
引用回帖:
Originally posted by flowerge at 2011-06-13 11:17:16:
你电荷是0，多重度是2，应该用UHF，不能用RHF

但是你说的碳正离子，应该是电荷1，多重度1吧？

电荷是1，多重度也是1的话，是用RHF吗
赞一下 回复此楼
9楼2011-06-13 12:30:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rock0927 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320分人工智能调剂 +8 振—TZ 2026-04-03 8/400 2026-04-05 22:33 by 范式思维
[考研] 材料专硕322分 +10 哈哈哈吼吼吼哈 2026-04-04 10/500 2026-04-05 21:22 by 学员8dgXkO
[考研] 348求调剂 +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by 猪会飞
[考研] 346分的生物与医药08600求调剂 +5 常雨阳上岸 2026-04-05 6/300 2026-04-05 13:42 by imissbao
[考研] 0703化学 +12 妮妮ninicgb 2026-04-04 13/650 2026-04-05 10:46 by 啊俊！
[考研] 材料工程085601数二英一335求调剂 +6 双马尾痞老板2 2026-03-31 6/300 2026-04-04 22:29 by hemengdong
[考研] 333求调剂 +12 wfh030413@ 2026-04-03 13/650 2026-04-04 21:02 by jj987
[考研] 一志愿安徽某211 0703化学总分339求调剂 +6 晚风不晚 2026-04-04 6/300 2026-04-04 20:11 by dongzh2009
[考研] 306求调剂 +3 hyb上名工 2026-04-02 3/150 2026-04-04 18:12 by 热情沙漠
[考研] 学硕机械工程303求调剂 +6 无名所以叫吴明 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[考研] 366求调剂 +7 sbdnd 2026-04-03 7/350 2026-04-03 12:40 by cymywx
[考研] 285求调剂 +6 FZAC123 2026-03-30 6/300 2026-04-03 12:22 by xingguangj
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 一志愿深大085601材料工程专业（专硕）300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 求调剂22408 288分 +5 new382 2026-04-02 5/250 2026-04-03 09:13 by 醉在风里
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 08工科275求调剂，可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 11:35 by wangjy2002
信息提示
关闭
请填处理意见
关闭