xzhfood

應(yīng)助: 5 (幼兒園)
貴賓: 22.859
金幣: 712.2
散金: 31
紅花: 23
帖子: 1820
在線: 837小時
蟲號: 6596
注冊: 2003-04-11
專業(yè): 食品科學(xué)基礎(chǔ)

[交流] 計算勢能面3D圖的討論

看到一篇論文，覺得確定化合物絕對構(gòu)型思路不錯，就想計算一下文章中的勢能面圖，化合物結(jié)構(gòu)與勢能面圖如圖所示。
文章見附件，這篇文章曾經(jīng)是論壇中推薦的一篇論文，影響因子5以上：

http://www.gaoyang168.com/bbs/viewthread.php?tid=2037219&fpage=1&target=blank

勢能面計算采用AM1水平
The potential energy surface of compound 1 was scanned at the AM1 level by rotating about the C-8-C-7 and C-7-C-1 bonds, resulting in eight minimum energy points (Fig. 3).

我按照這個思路，算法是：

%mem=1500MB
%nprocshared=2
%chk=C:\Users\xzhfood\Desktop\SHINeng\compound.Conf.M0001-scan.chk
# opt=modredundant am1 geom=connectivity

坐標設(shè)置按照 C-8-C-7 和 C-7-C-1及二者的夾角
剛開始設(shè)置：
B 8 7 S 100 0.05000 （鍵長，C-8-C-7）
B 7 5 S 100 0.05000  （鍵長， C-7-C-1）
A 8 7 5 S 100 0.05000  （鍵角， C-7-C-1）
計算一開始出錯：
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GEDIIS/GDIIS optimizer.
NTrRot= -1 NTRed= 235 NAtoms= 41 NSkip= 118 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in C:\G09W\l103.exe at Sat Aug 02 17:21:31 2014.
Job cpu time:  0 days  0 hours 13 minutes 55.0 seconds.
File lengths (MBytes):  RWF= 113 Int=    0 D2E=    0 Chk= 105 Scr=    1

變換坐標后：

B 8 7 S 40 1.000 （鍵長，C-8-C-7）
B 7 5 S 40 1.000  （鍵長， C-7-C-1）
A 8 7 5 S 40 1.000  （鍵角， C-7-C-1）

高斯計算2小時后，又出現(xiàn)相同的錯誤

最后再變換坐標：
B 8 7 S 5 1.0000
B 7 5 S 5 1.0000
A 8 7 5 S 5 1.0000

計算幾分種后，又出錯。

現(xiàn)在不知道如何設(shè)置了，請各位高手指點一下，多謝。

完整的輸入文件：

%mem=1500MB
%nprocshared=2
%chk=C:\Users\xzhfood\Desktop\SHINeng\compound.Conf.M0001-scan.chk
# opt=modredundant am1 geom=connectivity

M0001

0 1
C                1.30809652 -0.10155412 -1.59795358
C                0.46488955 -1.13359132 -2.30423701
C                -0.75973305 -0.32554395 -2.71648227
C                -0.91792950 0.68815622 -1.57608994
C                0.52618255 1.15640147 -1.23510403
O                2.49900702 -0.27607303 -1.36180499
C                0.65557829 1.59857018 0.25963276
C                0.25984831 0.52915068 1.27788930
C                -0.97208647 0.62837022 1.94435317
C                -1.37163744 -0.34422680 2.86566457
C                -0.54803968 -1.43928669 3.14705413
C                0.68466790 -1.52869796 2.49774451
C                1.09134670 -0.55951202 1.57877042
O                1.49824040 -2.59130734 2.77689589
C                -0.97544770 -2.48046881 4.13935119
C                2.02227520 2.20032811 0.63376373
C                0.99503296 2.27370851 -2.18955858
O                -1.74244568 1.77439734 -1.97287445
C                -1.99007969 -1.20109198 -2.91628176
H                0.99118228 -1.53839175 -3.17213966
H                0.22155241 -1.93157414 -1.59600296
H                -0.56616231 0.20021376 -3.66049352
H                -1.40096498 0.19741334 -0.72608194
H                -0.05469351 2.42943139 0.38891428
H                -1.63779076 1.46739165 1.74899207
H                -2.33237204 -0.23769672 3.36442172
H                2.06274061 -0.65015026 1.10022233
H                2.32275574 -2.49008845 2.27248578
H                -1.98541542 -2.28560594 4.51523521
H                -0.98276834 -3.47087151 3.67271003
H                -0.29475806 -2.48809835 4.99668795
H                2.04931338 2.47194976 1.69587529
H                2.22073248 3.11676922 0.06931922
H                2.84950173 1.50833736 0.45199681
H                0.49341030 3.22203931 -1.96870828
H                2.07630623 2.43580348 -2.11827099
H                0.79540133 2.02694117 -3.23817342
H                -2.63823679 1.42936686 -2.13293197
H                -2.85377490 -0.59966488 -3.21678173
H                -1.80916159 -1.94447000 -3.69970222
H                -2.25456400 -1.73677404 -1.99830705

1 2 1.0 5 1.0 6 2.0
2 3 1.0 20 1.0 21 1.0
3 4 1.0 19 1.0 22 1.0
4 5 1.0 18 1.0 23 1.0
5 7 1.0 17 1.0
6
7 8 1.0 16 1.0 24 1.0
8 9 1.0 13 2.0
9 10 2.0 25 1.0
10 11 1.0 26 1.0
11 12 2.0 15 1.0
12 13 1.0 14 1.0
13 27 1.0
14 28 1.0
15 29 1.0 30 1.0 31 1.0
16 32 1.0 33 1.0 34 1.0
17 35 1.0 36 1.0 37 1.0
18 38 1.0
19 39 1.0 40 1.0 41 1.0
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

B 8 7 S 5 1.0000
B 7 5 S 5 1.0000
A 8 7 5 S 5 1.0000

1.JPG

2.JPG

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : Cycloabiesesquine_A,_a_unique_sesquiterpenoid_from_Abies_delavay.pdf

2014-08-02 20:13:50, 734.86 K

» 猜你喜歡

高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有2人回復(fù)
化學(xué)化工學(xué)院招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）已經(jīng)有1人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費? 已經(jīng)有69人回復(fù)
河海大學(xué) （211.雙一流高校）鈣鈦礦光電實驗室丁勇團隊歡迎有志青年！已經(jīng)有3人回復(fù)
黑龍江省原子能研究院化學(xué)課題組招收2026級化學(xué)專業(yè)或化學(xué)相關(guān)專業(yè)碩士生已經(jīng)有0人回復(fù)
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有0人回復(fù)
溫州大學(xué)新能源材料微觀溯源課題組招收2026級博士研究生（材料、化學(xué)方向）已經(jīng)有6人回復(fù)
材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流已經(jīng)有5人回復(fù)
26申博已經(jīng)有2人回復(fù)
巴黎西岱大學(xué)（Université Paris Cité）Fayna MAMMERI 導(dǎo)師尋找CSC公派博士生已經(jīng)有9人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

這種圖是用什么軟件畫的，最好是容易上手的軟件已經(jīng)有16人回復(fù)
團簇吸附氫分子的吸附能居然是正的？已經(jīng)有36人回復(fù)
計算同一個分子在基態(tài)時不同構(gòu)型的能量已經(jīng)有5人回復(fù)
求助Ce3+5d能級劈裂大小以及5d重心位置計算方法已經(jīng)有5人回復(fù)
【討論】請教這幾張3D效果是用什么軟件做的？已經(jīng)有6人回復(fù)
【求助】怎樣判斷優(yōu)化后的構(gòu)型是穩(wěn)定的、合理的？已經(jīng)有10人回復(fù)
【討論】機械制圖基礎(chǔ)—看看你的立體感—2D畫3D 已經(jīng)有28人回復(fù)

1樓 2014-08-02 20:12:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xzhfood

榮譽版主 (著名寫手)

應(yīng)助: 5 (幼兒園)
貴賓: 22.859
金幣: 712.2
散金: 31
紅花: 23
帖子: 1820
在線: 837小時
蟲號: 6596
注冊: 2003-04-11
專業(yè): 食品科學(xué)基礎(chǔ)

3D 圖已經(jīng)算出來了，如何得到等高線圖呢？如何從等高線圖等到穩(wěn)定構(gòu)象呢？

贊一下

回復(fù)此樓

2樓2014-08-20 21:40:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xzhfood

榮譽版主 (著名寫手)

應(yīng)助: 5 (幼兒園)
貴賓: 22.859
金幣: 712.2
散金: 31
紅花: 23
帖子: 1820
在線: 837小時
蟲號: 6596
注冊: 2003-04-11
專業(yè): 食品科學(xué)基礎(chǔ)

3D 圖已經(jīng)算出來了，如何得到等高線圖呢？如何從等高線圖等到穩(wěn)定構(gòu)象呢？

贊一下

回復(fù)此樓

3樓2014-08-20 21:40:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

戴世杰

木蟲 (著名寫手)

應(yīng)助: 103 (高中生)
金幣: 4018.4
散金: 1099
紅花: 18
帖子: 1218
在線: 745.7小時
蟲號: 1275484
注冊: 2011-04-23
專業(yè): 食品加工技術(shù)

★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
gmy1990: 金幣+2 2014-08-22 04:10:49

引用回帖:

3樓: Originally posted by xzhfood at 2014-08-20 21:40:26
3D 圖已經(jīng)算出來了，如何得到等高線圖呢？如何從等高線圖等到穩(wěn)定構(gòu)象呢？

樓主，你用的是笛恰爾坐標，我實在沒看懂分子說明后面的變量B和A的意思。只要把坐標和能量的數(shù)據(jù)導(dǎo)出去，就可以用專門的軟件做登高線圖，比如matlab

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

4樓2014-08-21 22:10:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xzhfood

榮譽版主 (著名寫手)

應(yīng)助: 5 (幼兒園)
貴賓: 22.859
金幣: 712.2
散金: 31
紅花: 23
帖子: 1820
在線: 837小時
蟲號: 6596
注冊: 2003-04-11
專業(yè): 食品科學(xué)基礎(chǔ)

B 8 7 S 5 1.0000
B 7 5 S 5 1.0000
A 8 7 5 S 5 1.0000

B:Bond 表示鍵長，8 7 表示坐標， A ：Angle 表示鍵角；

我后來改成計算二面角，就得到三維圖了，打開log文件可以通過Results-scan看到3D圖，但是results--surfaces/contours顯示為灰色，不知道如何激活。

計算坐標設(shè)置：
D 5 4 10 11 S 9 30.000
D 4 10 11 17 S 9 30.000

贊一下

回復(fù)此樓

5樓2014-08-22 11:40:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 265求調(diào)劑 +10	梁梁校校 2026-03-17	10/500	2026-03-23 21:17 by 一切OK
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 材料292調(diào)劑 +3	橘頌思美人 2026-03-23	3/150	2026-03-23 19:36 by pswait
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5	葵梓衛(wèi)隊 2026-03-18	7/350	2026-03-23 16:26 by lingjue
[考研] 08工學(xué)調(diào)劑 +7	用戶573181 2026-03-20	11/550	2026-03-23 15:47 by 我愛學(xué)習(xí)學(xué)習(xí)使?/a>
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 考研調(diào)劑 +4	來好運來來來 2026-03-21	4/200	2026-03-22 12:15 by 星空星月
[基金申請] 山東省面上項目限額評審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學(xué)員8dgXkO
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學(xué)員8dgXkO
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 070300化學(xué)319求調(diào)劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3	123阿標 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦