| 查看: 2467 | 回复: 2 | ||
shasha87铁虫 (小有名气)
|
[求助]
高斯优化中怎样固定内坐标
|
| 各位大神好，我有一个问题想请教一下。我在高斯03中优化Hg(CH3)Cl分子，想把其内坐标固定，但是Hg-X不固定（英文是：Internal coordinates were held fixed, except Hg-X)。输入文件要怎么写，需要哪些关键字。谢谢。帮帮忙吧大家！ |
回复此楼» 猜你喜欢
散金求好运
已经有315人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有253人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
祈福青基
已经有34人回复
-
牛来！米来！面来！
已经有8人回复
-
2027申博咨询
已经有11人回复
-
寻求CrC2Al的cif文件，感谢
已经有0人回复
-
ORCA 纯新手向入门指南0.7版
已经有0人回复
-
哈尔滨工业大学（深圳）绿色合成化学团队招收27年春季或秋季入学的博士生
已经有5人回复
-
吉林大学化学学院理论化学研究所刘靖尧教授课题组招收2027级申请考核制博士研究生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯优化查看结构，发现H-O键消失，但是H原子还在？
已经有31人回复
- 如何让高斯在优化结构时保持其中某个单键不变
已经有7人回复
- 高斯优化过渡态，有两个虚频，一个想要，想把另一个消掉-----请高手指点
已经有10人回复
- 在进行高斯计算时如何消去虚频（转载）
已经有18人回复
- 求助如何固定原子坐标进行优化，快搞疯了
已经有16人回复
- 高斯考虑溶剂效应后优化出错
已经有8人回复
- 如何 用高斯 固定部分原子 进行优化
已经有4人回复
- 高斯优化求助
已经有7人回复
- 高斯对多分子体系的构型优化
已经有4人回复
- 高斯扫势能面的问题
已经有5人回复
- TiO2（101）晶面的优化
已经有13人回复
- 【求助】Ms如何导入通过高斯优化好的结构
已经有9人回复
- 【讨论】这简直是爆炸反应啊（TNT算毛毛啊），高斯怎能这样优化，有图有真相
已经有8人回复
- 【求助】键长键角等固定后的优化..悬赏30...【已经解决～感谢～】
已经有22人回复
- 【求助】咨询：高斯优化中局域能量最小值和全局最小值的问题
已经有6人回复
- 【求助】高斯中坐标固定和冻结核是一回事吗？
已经有5人回复
- 【求助】高斯优化好后如何看标记碳原子到环糊精XY平面的距离
已经有3人回复
- 【求助】高斯优化过渡态出错
已经有14人回复
- 【讨论】话题：Ir(111)配合物用高斯优化时的自旋多重度
已经有3人回复
zjx187
金虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1489.9
- 帖子: 176
- 在线: 255.1小时
- 虫号: 1876020
- 注册: 2012-07-02
- 性别: MM
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gkf高: 金币+2 2013-09-09 23:03:41
shasha87: 金币+5, ★有帮助 2013-09-10 23:16:06
感谢参与，应助指数 +1
gkf高: 金币+2 2013-09-09 23:03:41
shasha87: 金币+5, ★有帮助 2013-09-10 23:16:06
|
如果你用内坐标表示来表示原子，是不需要加关键词的，比如说下面的CH4分子的内坐标形式 0 1 C H 1 1.07000000 H 1 B2 2 A1 H 1 B3 3 A2 2 D1 0 H 1 B4 3 A3 2 D2 0 B2 1.07000000 B3 1.07000000 B4 1.07000000 A1 109.47120255 A2 109.47125080 A3 109.47121829 D1 -119.99998525 D2 120.00000060 你把B1的数值放到上面就是固定了1号碳与2号H院子之间的距离，一次类推，就可以固定原子之间的键角以及二面角，希望对你有帮助。 |
2楼2013-09-09 20:21:43
shasha87
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 3664.8
- 散金: 5
- 红花: 1
- 帖子: 275
- 在线: 70小时
- 虫号: 2124965
- 注册: 2012-11-13
- 性别: MM
- 专业: 同位素地球化学
|
您好，我这样做过，但是发现固定不了，最后还是变了。我把我的输入文件发给你看一下：谢谢你！ %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 |
3楼2013-09-10 09:17:52
