| 查看: 2613 | 回复: 1 | ||
shasha87铁虫 (小有名气)
|
[求助]
高斯中固定部分原子的内坐标进行优化
|
|
求大神帮帮忙！很是痛苦，我搞了一晚上，依然是错的。 我想优化Hg（CH3)Cl，固定住一部分内坐标，结果死掉了。我的输入文件是这样的： %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文件里的是本该固定的坐标却变化了。输入文件是这样的：%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 |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
汕头大学化学招调剂，相关内容如下。
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有202人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
高斯优化查看结构，发现H-O键消失，但是H原子还在？
已经有31人回复- 卟啉（100—150个原子）高斯计算的计算机配置问题
已经有14人回复
- 高斯09寻找过渡态，急需帮助，十分感谢
已经有6人回复
- 高斯中双原子分子计算 如何调设置使离解能正常?
已经有4人回复
- gauss输入文件中那个原子X是什么意思
已经有3人回复
- 从高斯激发态优化输出文件中提取原子坐标
已经有9人回复
- 高斯中两个原子的位置重合
已经有6人回复
- 如何 用高斯 固定部分原子 进行优化
已经有4人回复
- 高斯优化求助
已经有7人回复
- 高斯对多分子体系的构型优化
已经有4人回复
- 高斯view中我确定不了原子的位子,想通过两条键长的交点确定,怎么弄
已经有5人回复
- 高斯view中调节键长时怎么使四五个原子组成一个团（group）,和另一个原子之间进行调节
已经有4人回复
- 【求助】Ms如何导入通过高斯优化好的结构
已经有9人回复
- 【讨论】这简直是爆炸反应啊（TNT算毛毛啊），高斯怎能这样优化，有图有真相
已经有8人回复
- 【求助】无法用高斯软件算出含吡啶、呋喃这类分子的能级轨道
已经有10人回复
- 【求助】高斯中坐标固定和冻结核是一回事吗？
已经有5人回复
- 【求助】高斯优化好后如何看标记碳原子到环糊精XY平面的距离
已经有3人回复
- 【求助】请问高斯中怎样得到分子轨道中各原子轨道的组成
已经有4人回复
- 【求助】高斯计算中如何冻结一部分原子
已经有8人回复
29050801
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1065.5
- 散金: 6
- 红花: 1
- 帖子: 163
- 在线: 78.4小时
- 虫号: 1765116
- 注册: 2012-04-19
- 性别: GG
- 专业: 生物化学
2楼2013-11-16 09:12:28