| 查看: 2614 | 回复: 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 |
回复此楼» 猜你喜欢
北京师范大学胡德平课题组以及袁萌伟课题组考研招生
已经有0人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有290人回复
-
祈福论文小修
已经有21人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有26人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有40人回复
» 本主题相关价值贴推荐，对您同样有帮助:
高斯优化查看结构，发现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