| 2 | 1/1 | 返回列表 |
| 查看: 2660 | 回復(fù): 1 | ||
shasha87鐵蟲(chóng) (小有名氣)
|
[求助]
高斯中固定部分原子的內(nèi)坐標(biāo)進(jìn)行優(yōu)化
|
|
求大神幫幫忙!很是痛苦,我搞了一晚上,依然是錯(cuò)的。 我想優(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 錯(cuò)誤顯示: 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 請(qǐng)大神幫幫忙吧!灰常感謝! 我用另外一種方法固定,但是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 |
新蟲(chóng) (小有名氣)

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +5 | 深V宿舍吧 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過(guò)六級(jí)。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +4 | 篤行致遠(yuǎn). 2026-03-28 | 4/200 |
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè),初試351分,復(fù)試被刷求調(diào)劑! +3 | 真名有冰 2026-03-29 | 4/200 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研]
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|