24小时热门版块排行榜    

查看: 406  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

railic

新虫 (初入文坛)

[交流] 【求助】为什么在优化计算中冻结二面角没有作用?

不好意思,已经问过一次了...但是似乎没有结论....
厚着脸再问一次-_-这次我说清楚些
这几天快被这个东东整的不行了,见下:
# rhf/6-31G(d) Opt=Modredundant geom=connectivity

cisTMS opt

0 1
C
O                  1              B1
C                  2              B2    1              A1
O                  2              B3    1              A2    3              D1
C                  4              B4    2              A3    1              D2
O                  4              B5    2              A4    1              D3
C                  6              B6    4              A5    2              D4
C                  1              B7    7              A6    6              D5
C                  1              B8    8              A7    2              D6
C                  8              B9    1              A8    9              D7
C                  9             B10    1              A9    8              D8
...
   D6           -25.70980734
   D7             20.
   D8            34.30284260
...
1 8 9 10 F
...
结果输出最优文件里面:
D1    D(9,1,8,10)            19.6418         -DE/DX =   -0.005
若是D7变化到大一点的值,优化结果能和原先设定差10度以上,为什么啊...我不是冻结了吗?急求帮助!
谢啦谢啦~
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

railic

新虫 (初入文坛)

o....hoho,了解咯~谢谢各位啦~
5楼2009-04-14 00:30:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

quantumor

金虫 (著名写手)

快乐兔子

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 4-20 17:29
你的输入文件里对D1的定义是D(1,8,9,8),但冻结坐标部分对该角的定义却是D(1,8,9,10),定义的对象不同啊。
建议你将坐标里的D7定义为“1 8 9 10” 试试。
愿好运与快乐伴随你!
2楼2009-04-13 07:26:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

railic

新虫 (初入文坛)

没有没有,d7那一行对应的碳原子编号就是10啊~所以d7定义的就是1,8,9,10四个点定义的二面角啊~
3楼2009-04-13 10:44:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

半支烟5591

银虫 (正式写手)

专家



xuefei06(金币+1,VIP+0):感谢你的参与! 4-13 23:43
Opt=Modredundant 本来就是一边优化一边扫描,费那事干吗,推荐阿A RIGID  SCAN !
4楼2009-04-13 13:56:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见