24小时热门版块排行榜    

查看: 1138  |  回复: 7
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 cysqxr 的 1 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cysqxr

木虫 (小有名气)

[交流] 【求助】有人能解释一下局部优化吗?

局部优化在高斯中有几个固定的关键字,但是有哪位高手给指点一下具体应该怎么办?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abbott

金虫 (著名写手)

不要用QQ问我东西

★ ★
gwdavid(金币+2,VIP+0):谢谢交流!欢迎再次光临模拟版!
为什么进行局部优化:
冻结坐标有时是有用的,举几个例子:
1,如优化溶质分子和nH2O的supermolecule,周围的水分子,可能结构变化不大,这时就可以固定水的键长或键角,减少工作量。
2,如用cluster模拟大块固体,这时也许就要按照晶格参数固定原子键长和键角。
3,又如不固定NH3-H-H2O中N-O键长是得不到质子传递势垒。
4,还有在优化稳定构型或过渡态得时候也有时用到部分优化得到接近得结果,再放开优化,这样逐步得到你想要得结果。
关于部分优化的实例,我目前就是采用下面的方法做的。

opt=z-matrix的部分优化。

%Chk=reac-hf3
#p RHF/STO-3G opt(z-matrix,maxcyc=500)
freq=noraman scf(maxcyc=200)

reactant fopt

1 1
N -1 -1.75389 -0.48105 0.31424
H 0  x2 y2 z2
H 0  x3 y3 z3
H 0  x4 y4 z4
O -1 2.16296 0.18089 -0.15503
H 0  x6 y6 z6
H 0  x7 y7 z7
H 0  x8 y8 z8

x2=-2.2088
y2=0.39705
z2=0.46255
x3=-2.14589
y3=-0.9279
z3=-0.4899
x4=-1.88451
y4=-1.06386
z4=1.11629
x6=1.36256
y6=0.28882
z6=-0.67401
x7=2.32973
y7=0.98241
z7=0.34634
x8=2.90378
y8=0.01009
z8=-0.74115

***********************************************
opt=modredundant的部分优化。
分子坐标部分可以为具体的数值,而不用为参数形式。

%Chk=reac-hf3
#p RHF/STO-3G opt(modredundant,maxcyc=500)
freq=noraman scf(maxcyc=200)

reactant fopt

1 1
N -1 -1.75389 -0.48105 0.31424
H 0  
H 0  
H 0  
O -1 2.16296 0.18089 -0.15503
H 0  
H 0  
H 0  

1 5 F


希望这些信息对你有用。

以上的东西转自厦大论坛

[ Last edited by abbott on 2008-9-29 at 22:08 ]
Chemistry[]==[]Chem[]is[]try!!!
6楼2008-09-29 22:07:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

虚谦

至尊木虫 (正式写手)

★ ★
gwdavid(金币+2,VIP+0):谢谢交流!欢迎再次光临模拟版!
用gen分别优化
或者是有金属的用过popt局部优化
大家都说我是个好人,其实,我也是这么想的~!
2楼2008-09-28 12:08:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

朱三娥

铜虫 (小有名气)


gwdavid(金币+1,VIP+0):谢谢你的回答,欢迎再次光临计算模拟版!
好象是用FREZON固定不要优化的一部分
3楼2008-09-28 16:44:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jysgg1015

铁杆木虫 (著名写手)

★ ★
gwdavid(金币+2,VIP+0):谢谢你的回答,欢迎再次光临计算模拟版!
你可以固定一部分参数,利用opt=modredundant.
男儿当自强!!!
4楼2008-09-29 21:36:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见