版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

shasha87

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3664.8
散金: 5
紅花: 1
帖子: 275
在線: 70小時(shí)
蟲號: 2124965
注冊: 2012-11-13
性別: MM
專業(yè): 同位素地球化學(xué)

[求助] 高斯中固定部分原子的內(nèi)坐標(biāo)進(jìn)行優(yōu)化

求大神幫幫忙！很是痛苦，我搞了一晚上，依然是錯的。
我想優(yōu)化Hg（CH3)Cl，固定住一部分內(nèi)坐標(biāo)，結(jié)果死掉了。我的輸入文件是這樣的：
%chk=F:\Calculation\1.chk
# RHF/Gen Pseudo=Read Opt=ModRedundant Test

rhf/cc-pvdz-pp Opt of Hg(CH3)Cl

0 1
Cl
Hg 1 R1
C  2 R2    1 180.00
H  3 1.0899 2 109.5901 1 109.5901
H  3 1.0899 2 109.5901 4 109.3521
H  3 1.0899 2 109.5901 5 109.3521

R1=2.3417
R2=2.1040

3 4             F
3 5             F
3 6             F
3 2 1          F
4 3 2          F
5 3 2          F
6 3 2          F
4 3 2 1          F
5 3 2 4          F
6 3 2 5          F

C H Cl 0
cc-pvdz
****
Hg 0
S 7 1.00
   41.4804000             0.0100720
   22.5328000          -0.1147810
   14.0790000             0.3564080
   5.7247400          -0.8484430
   1.4794800             0.9405960
   0.6882660             0.4353290
   0.1641810             0.0131780
S 7 1.00
   41.4804000             0.0020320
   22.5328000          -0.0323140
   14.0790000             0.1153500
   5.7247400          -0.3177630
   1.4794800             0.5135230
   0.6882660             0.2500500
   0.1641810          -0.6488410
S 7 1.00
   41.4804000             0.0218280
   22.5328000          -0.1555640
   14.0790000             0.3860250
   5.7247400          -0.8961040
   1.4794800             2.3321950
   0.6882660          -1.5976600
   0.1641810          -1.0536460
S 1 1.00
   0.0572230             1.0000000
P 6 1.00
   10.5805000             0.1286920
   6.8053900          -0.3608420
   1.7901600             0.5641170
   0.8622610             0.5002740
   0.3752770             0.1137850
   0.1284210             0.0029390
P 6 1.00
   10.5805000          -0.0344430
   6.8053900             0.1009290
   1.7901600          -0.1924390
   0.8622610          -0.1960690
   0.3752770             0.1042040
   0.1284210             0.5873870
P 6 1.00
   10.5805000             0.0714680
   6.8053900          -0.2108590
   1.7901600             0.4452040
   0.8622610             0.3916080
   0.3752770          -0.6670830
   0.1284210          -0.6195740
P 1 1.00
   0.0434470             1.0000000
D 5 1.00
   11.5841000             0.0167290
   7.2487900          -0.0695580
   1.9366200             0.3034370
   0.8987480             0.4548780
   0.3869310             0.3368080
D 5 1.00
   11.5841000          -0.0249840
   7.2487900             0.1041030
   1.9366200          -0.6169860
   0.8987480          -0.3538260
   0.3869310             0.6872800
D 1 1.00
   0.1504660             1.0000000
F 1 1.00
   0.9974000             1.0000000
****

Hg 0
HG-ECP    5    60
h-ul potential
  1
2    1.0000000             0.0000000
s-ul potential
  2
2    12.4130710          275.7747970
2    6.8979130          49.2678980
p-ul potential
  4
2    11.3103200          80.5069840
2    10.2107730          161.0348240
2    5.9398040             9.0834160
2    5.0197550          18.3677730
d-ul potential
  4
2    8.4078950          51.1372560
2    8.2140860          76.7074590
2    4.0126120             6.5618210
2    3.7953980             9.8180700
f-ul potential
  2
2    3.2731060             9.4290010
2    3.2083210          12.4948560
g-ul potential
  2
2    4.4852960          -6.3384140
2    4.5132000          -8.0998630

錯誤顯示： Symmetry A" KE= 6.458561633319D+01
2 Symmetry operations used in ECPInt.
ECPInt:  NShTT=    630 NPrTT= 8835 LenC2=    625 LenP2D= 6879.
LDataN:  DoStor=F MaxTD1= 9 Len=  602
LDataN:  DoStor=T MaxTD1= 9 Len=  602
***** Axes restored to original set *****
-------------------------------------------------------------------
Center    Atomic                Forces (Hartrees/Bohr)
Number    Number             X             Y             Z
-------------------------------------------------------------------
1       17          0.000030377 -0.000086492 -0.000104841
2       80       -0.000142390 0.000405421 0.000535743
3       6       -0.010894671 0.031019845 -0.003355223
4       1       -0.003600434 -0.024419441 0.005719256
5       1       -0.003472078 0.009885872 -0.008514192
6       1          0.018079197 -0.016805204 0.005719256
-------------------------------------------------------------------
Cartesian Forces:  Max    0.031019845 RMS    0.011923168

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Internal  Forces:  Max          nan RMS          nan
Search for a local minimum.
Step number 1 out of a maximum of  25
All quantities printed in internal units (Hartrees-Bohrs-Radians)
LAPACK or BLAS routine DLASCL set Info=       -4.
Bad arguments to LAPack or BLAS routine.
Error termination via Lnk1e in C:\G03W\l103.exe at Mon Sep 09 23:23:25 2013.
Job cpu time:  0 days  0 hours  0 minutes 31.0 seconds.
File lengths (MBytes):  RWF=    16 Int=    0 D2E=    0 Chk=    10 Scr=    1
請大神幫幫忙吧！灰常感謝！
我用另外一種方法固定，但是out文件里的是本該固定的坐標(biāo)卻變化了。輸入文件是這樣的：%chk=F:\Calculation\1.chk
# RHF/Gen Pseudo=Read Opt

rhf/cc-pvdz-pp Opt of Hg(CH3)Cl

0 1
Cl
Hg 1 R1
C  2 R2    1 180.00
H  3 1.0899 2 109.5901 1 109.5901
H  3 1.0899 2 109.5901 4 109.3521
H  3 1.0899 2 109.5901 5 109.3521

R1=2.3417
R2=2.1040

C H Cl 0
cc-pvdz
****
Hg 0
S 7 1.00
   41.4804000             0.0100720
   22.5328000          -0.1147810
   14.0790000             0.3564080
   5.7247400          -0.8484430
   1.4794800             0.9405960
   0.6882660             0.4353290
   0.1641810             0.0131780
S 7 1.00
   41.4804000             0.0020320
   22.5328000          -0.0323140
   14.0790000             0.1153500
   5.7247400          -0.3177630
   1.4794800             0.5135230
   0.6882660             0.2500500
   0.1641810          -0.6488410
S 7 1.00
   41.4804000             0.0218280
   22.5328000          -0.1555640
   14.0790000             0.3860250
   5.7247400          -0.8961040
   1.4794800             2.3321950
   0.6882660          -1.5976600
   0.1641810          -1.0536460
S 1 1.00
   0.0572230             1.0000000
P 6 1.00
   10.5805000             0.1286920
   6.8053900          -0.3608420
   1.7901600             0.5641170
   0.8622610             0.5002740
   0.3752770             0.1137850
   0.1284210             0.0029390
P 6 1.00
   10.5805000          -0.0344430
   6.8053900             0.1009290
   1.7901600          -0.1924390
   0.8622610          -0.1960690
   0.3752770             0.1042040
   0.1284210             0.5873870
P 6 1.00
   10.5805000             0.0714680
   6.8053900          -0.2108590
   1.7901600             0.4452040
   0.8622610             0.3916080
   0.3752770          -0.6670830
   0.1284210          -0.6195740
P 1 1.00
   0.0434470             1.0000000
D 5 1.00
   11.5841000             0.0167290
   7.2487900          -0.0695580
   1.9366200             0.3034370
   0.8987480             0.4548780
   0.3869310             0.3368080
D 5 1.00
   11.5841000          -0.0249840
   7.2487900             0.1041030
   1.9366200          -0.6169860
   0.8987480          -0.3538260
   0.3869310             0.6872800
D 1 1.00
   0.1504660             1.0000000
F 1 1.00
   0.9974000             1.0000000
****

Hg 0
HG-ECP    5    60
h-ul potential
  1
2    1.0000000             0.0000000
s-ul potential
  2
2    12.4130710          275.7747970
2    6.8979130          49.2678980
p-ul potential
  4
2    11.3103200          80.5069840
2    10.2107730          161.0348240
2    5.9398040             9.0834160
2    5.0197550          18.3677730
d-ul potential
  4
2    8.4078950          51.1372560
2    8.2140860          76.7074590
2    4.0126120             6.5618210
2    3.7953980             9.8180700
f-ul potential
  2
2    3.2731060             9.4290010
2    3.2083210          12.4948560
g-ul potential
  2
2    4.4852960          -6.3384140
2    4.5132000          -8.0998630

回復(fù)此樓

» 猜你喜歡

求助各位大佬，球球了，在這科薩頭已經(jīng)有0人回復(fù)
求標(biāo)準(zhǔn)卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340 已經(jīng)有2人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有110人回復(fù)
量子點(diǎn)電致發(fā)光器件（QLED）：藍(lán)光量子點(diǎn)（QDs）和納米氧化鋅電子傳輸材料已經(jīng)有0人回復(fù)
需要一個(gè)圓偏振光誘導(dǎo)化學(xué)反應(yīng)，請問使用圓偏振光作為光源，應(yīng)該需要哪些設(shè)備搭建已經(jīng)有0人回復(fù)
半透性容器進(jìn)行高濕或者低濕研究已經(jīng)有0人回復(fù)
求ENDNOTE破解軟件安裝包已經(jīng)有4人回復(fù)
河海大學(xué) （211.雙一流高校）鈣鈦礦光電實(shí)驗(yàn)室丁勇團(tuán)隊(duì)歡迎有志青年！已經(jīng)有3人回復(fù)
黑龍江省原子能研究院化學(xué)課題組招收2026級化學(xué)專業(yè)或化學(xué)相關(guān)專業(yè)碩士生已經(jīng)有0人回復(fù)
26博士申請已經(jīng)有0人回復(fù)
招生材料與化工、環(huán)境科學(xué)調(diào)劑研究生已經(jīng)有4人回復(fù)

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

高斯優(yōu)化查看結(jié)構(gòu)，發(fā)現(xiàn)H-O鍵消失，但是H原子還在？已經(jīng)有31人回復(fù)
卟啉（100—150個(gè)原子）高斯計(jì)算的計(jì)算機(jī)配置問題已經(jīng)有14人回復(fù)
高斯09尋找過渡態(tài)，急需幫助，十分感謝已經(jīng)有6人回復(fù)
高斯中雙原子分子計(jì)算如何調(diào)設(shè)置使離解能正常? 已經(jīng)有4人回復(fù)
gauss輸入文件中那個(gè)原子X是什么意思已經(jīng)有3人回復(fù)
從高斯激發(fā)態(tài)優(yōu)化輸出文件中提取原子坐標(biāo) 已經(jīng)有9人回復(fù)
高斯中兩個(gè)原子的位置重合已經(jīng)有6人回復(fù)
如何用高斯固定部分原子進(jìn)行優(yōu)化已經(jīng)有4人回復(fù)
高斯優(yōu)化求助已經(jīng)有7人回復(fù)
高斯對多分子體系的構(gòu)型優(yōu)化已經(jīng)有4人回復(fù)
高斯view中我確定不了原子的位子,想通過兩條鍵長的交點(diǎn)確定,怎么弄已經(jīng)有5人回復(fù)
高斯view中調(diào)節(jié)鍵長時(shí)怎么使四五個(gè)原子組成一個(gè)團(tuán)（group）,和另一個(gè)原子之間進(jìn)行調(diào)節(jié) 已經(jīng)有4人回復(fù)
【求助】Ms如何導(dǎo)入通過高斯優(yōu)化好的結(jié)構(gòu) 已經(jīng)有9人回復(fù)
【討論】這簡直是爆炸反應(yīng)�。═NT算毛毛�。�，高斯怎能這樣優(yōu)化，有圖有真相已經(jīng)有8人回復(fù)
【求助】無法用高斯軟件算出含吡啶、呋喃這類分子的能級軌道已經(jīng)有10人回復(fù)
【求助】高斯中坐標(biāo)固定和凍結(jié)核是一回事嗎？已經(jīng)有5人回復(fù)
【求助】高斯優(yōu)化好后如何看標(biāo)記碳原子到環(huán)糊精XY平面的距離已經(jīng)有3人回復(fù)
【求助】請問高斯中怎樣得到分子軌道中各原子軌道的組成已經(jīng)有4人回復(fù)
【求助】高斯計(jì)算中如何凍結(jié)一部分原子已經(jīng)有8人回復(fù)

1樓 2013-09-09 23:48:48

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

29050801

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1065.5
散金: 6
紅花: 1
帖子: 163
在線: 78.4小時(shí)
蟲號: 1765116
注冊: 2012-04-19
性別: GG
專業(yè): 生物化學(xué)

你用的windows版本吧？

贊一下

回復(fù)此樓

個(gè)人主頁http://cnrdn.com/pXg5

2樓2013-11-16 09:12:28

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 shasha87 的主題更新