24小时热门版块排行榜    

查看: 422  |  回复: 9
当前主题已经存档。

railic

新虫 (初入文坛)

[交流] 【求助】为什么冻结了的坐标做优化时还是会变化呢....

小弟在做优化时,将一个二面角坐标冻结,但是结果的优化结果还是和冻结的值不是太一样,不知有什么错误?
in文件片段:
# RHF/6-31G(d) Opt=ModRedundant

test

0        1       
C1
O2                   1    B1
C3                   2    B2       1    A1
O4                   2    B3       1    A2       3    D1       0
C5                   4    B4       2    A3       1    D2       0
O6                   4    B5       2    A4       1    D3       0
C7                   6    B6       4    A5       2    D4       0
C8                   1    B7       7    A6       6    D5       0
C9                   1    B8       8    A7       2    D6       0
C10                  8    B9       1    A8       9    D7       0
C11                  9    B10      1    A9       8    D8       0
......

1 8 9 10         F
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

railic

新虫 (初入文坛)

顶一顶,实在不明白为什么....求助~
例如上面的D7设定的是80度,最后结果里D7已经变成70点零几度了。。。
2楼2009-04-05 00:11:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niliu

铁杆木虫 (著名写手)

★ ★
yuhuobuku(金币+2):欢迎参加讨论
你直接在结构描述部分用80.0替换D7,后面不用设。
3楼2009-04-05 17:17:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

railic

新虫 (初入文坛)

这样不相当于没有冻结吗?...不过我还是试了下,确实相当于没有冻结啊....
是不是我误会你的意思了?
4楼2009-04-06 12:07:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★
yuhuobuku(金币+2):欢迎参加讨论
加上geom=modrundant
5楼2009-04-06 17:58:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

railic

新虫 (初入文坛)

指的是opt=modredundant吗?因为geom那个会报错的~
opt=modredundant选项我加了啊.....
6楼2009-04-06 22:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

★ ★
yuhuobuku(金币+2):欢迎参加讨论
我记得3楼的方法可以的啊,我原来还用过。
莫非gauss版本更新后此法不通了?
7楼2009-04-06 23:48:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

想起来了,你应该把最后一列的0都去掉
8楼2009-04-06 23:49:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

railic

新虫 (初入文坛)

在z-marix里面把那列0去掉了,列表里直接就是那个值,还是不成啊好像....

# Opt RHF/6-31G(d)

testcase for TMS single point opt

0        1       
C1
O2                   1    B1
C3                   2    B2       1    A1
O4                   2    B3       1    A2       3    D1      
C5                   4    B4       2    A3       1    D2      
O6                   4    B5       2    A4       1    D3      
C7                   6    B6       4    A5       2    D4      
C8                   1    B7       7    A6       6    D5      
C9                   1    B8       8    A7       2    D6      
C10                  8    B9       1    A8       9    20.      
C11                  9    B10      1    A9       8    D8      
C12                  9    B11      1    A10      8    D9      
C13                  10   B12      8    A11      1    D10      
......
      Variables:
......
D1                  177.04064                  
D2                  117.51173                  
D3                    8.84396                  
D4                 -119.62742                  
D5                   19.77555                  
D6                  -25.8056                                    
D8                   34.55914                  
D9                 -149.06025                  
D10                -150.75611                  
......

我觉得还是要冻结的吧,这样算出来的和没冻结一样...
9楼2009-04-07 23:52:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

railic

新虫 (初入文坛)

sigh...还是不成,为什么呢???
10楼2009-04-08 11:30:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 railic 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见