24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 964  |  回复: 4
当前主题已经存档。

kohn

木虫 (正式写手)

[交流] [求助]gaussian优化结构的时候怎么固定原子阿?

使用gaussian优化结构的时候想固定住几个原子,看了Help,好像就只能用ModRedundant命令.但这个命令的语法实在是没有搞明白.
哪位能否给一个例子?最好是一个gif文件。

谢谢了!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyy291

铁杆木虫 (小有名气)

★ ★ ★
kohn(金币+2,VIP+0):谢谢!
yyx19840628(金币+1,VIP+0):谢谢共享
关键词里用Opt=ModRedundant,分子说明部分之后空行后对冻结做说明,如用内坐标冻结实例:
# HF/6-31G* Opt=ModRedundant Test
Partial optimization
1 1
C
H 1 R1
H 1 R1 2 A1
O 1 R2 2 A2 3 120.0
H 4 R3 3 A3 2 180.0
A1=120.0
...
R3=1.1

4 5             1.3 F
5 4 3 2             F


这个输入固定O-H键和到最后一个氢原子的二面角,用这种方法定义的任何值都不能和前面Z矩阵列出的值相同(这里就是O-H键长);通过调整结构来保证这个限制,限制的值可选,如该例中第二个修改的内坐标默认为Z矩阵的值(180.0)。





或者分子说明部分中可选的第二个位置,设为-1(默认为0),那么几何优化时该原子就被冻结了。输入部分之后不用再做冻结说明了。

C     -1     0.0     0.0    0.0
H      0     0.0     0.0    0.9
...


上例中C被冻结
如果无法忘记,就不要忘记好了,真正的忘记是不需要努力的。
2楼2008-11-27 09:18:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kohn

木虫 (正式写手)

非常感谢yyy291 虫虫!

我的email是kohnzhang@yahoo.com.cn,希望能进一步向你请教,能留个email么?
3楼2008-11-27 10:43:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyy291

铁杆木虫 (小有名气)

★ ★ ★
kohn(金币+3,VIP+0):只有你一个人回复,金币全部给你吧:)能否给我发两个个完整的gif文件。谢谢!
深入我就不太清楚了,呵呵,也是一边学一边尝试着的
如果无法忘记,就不要忘记好了,真正的忘记是不需要努力的。
4楼2008-11-27 14:21:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyy291

铁杆木虫 (小有名气)

★ ★
erylingjet(金币+2,VIP+0):感谢交流
%mem=800mb
%chk=trans-cage.chk
#t opt am1 scf=maxcycle=999

Title

0 1
O   -1   7.033600  -3.946100  -2.047400
O   -1  -0.913200  -9.215600  -2.011800
O   -1  -6.186200  -1.274300  -2.055000   
O   -1   1.759900   3.998700  -2.102500   
O   -1   1.350900  -7.728700  -4.740400   
O   -1  -0.913200  -9.215600   2.011800   
O   -1   4.769000  -5.403900   4.778900   
O   -1   3.261400   1.713800  -4.805600   
O   -1   1.759900   3.998700   2.102500   
O   -1   5.531600  -1.635300   4.753900   
O   -1  -4.686600  -3.560600  -4.758000   
O   -1  -0.499800   2.545300   4.726700   
O   -1   1.350900  -7.728700   4.740400   
O   -1   3.261400   1.713800   4.805600   
O   -1   4.769000  -5.403900  -4.778900   
O   -1   5.531600  -1.635300  -4.753900   
O   -1  -2.412100  -6.909300   4.791900   
O   -1  -3.923900   0.207800  -4.787900   
O   -1  -6.186200  -1.274300   2.055000   
O   -1   7.033600  -3.946100   2.047400   
O   -1  -0.499800   2.545300  -4.726700   
O   -1  -3.923900   0.207800   4.787900   
O   -1  -4.686600  -3.560600   4.758000   
O   -1  -2.412100  -6.909300  -4.791900   
O   -1   5.705000  -6.190800  -2.404700   
O   -1  -3.155800  -7.889100  -2.390600   
O   -1  -4.857900   0.965100  -2.446200   
O   -1   4.002200   2.667700  -2.466600   
O   -1   1.669100  -8.865300  -2.387500   
O   -1   1.669100  -8.865300   2.387500   
O   -1   5.705000  -6.190800   2.404700   
O   -1   1.099400   0.717900  -5.931300   
O   -1  -0.818800   3.678500   2.372300   
O   -1   3.781900  -3.235100   5.899300   
O   -1  -5.833900  -3.859500  -2.407700   
O   -1   1.099400   0.717900   5.931300   
O   -1  -0.250100  -5.915400   5.919100   
O   -1   4.002200   2.667700   2.466600   
O   -1   3.781900  -3.235100  -5.899300   
O   -1   6.681800  -1.366000  -2.436500   
O   -1  -3.155800  -7.889100   2.390600   
O   -1  -2.938400  -1.962000  -5.907000   
O   -1  -5.833900  -3.859500   2.407700   
O   -1   6.681800  -1.366000   2.436500   
O   -1  -0.818800   3.678500  -2.372300   
O   -1  -4.857900   0.965100   2.446200   
O   -1  -2.938400  -1.962000   5.907000   
O   -1  -0.250100  -5.915400  -5.919100   
O   -1   6.426300  -5.461500   0.000000   
O   -1  -2.441300  -8.595400   0.000000   
O   -1  -5.576200   0.270500   0.000000   
O   -1   3.290499   3.406000   0.000000   
O   -1   3.577700  -7.370100  -3.392300   
O   -1   0.732900  -9.238000   0.000000   
O   -1   3.577700  -7.370100   3.392300   
O   -1   3.267700  -0.747100  -5.725700   
O   -1   0.116300   4.048000   0.000000   
O   -1   3.267700  -0.747100   5.725700   
O   -1  -4.334100  -5.776300  -3.389800   
O   -1  -1.458400   0.292600   5.681700   
O   -1   2.314600  -5.403800   5.717300   
O   -1   5.182800   0.584900   3.384100   
O   -1   2.314600  -5.403800  -5.717300   
O   -1   5.182800   0.584900  -3.384100   
O   -1  -4.334100  -5.776300   3.389800   
O   -1  -1.458400   0.292600  -5.681700   
O   -1  -6.217900  -2.903600   0.000000   
O   -1   7.068800  -2.287500   0.000000   
O   -1  -2.733900   2.132800  -3.452700   
O   -1  -2.733900   2.132800   3.452700   
O   -1  -2.417800  -4.449600   5.715300   
O   -1  -2.417800  -4.449600  -5.715300   
Si  -1   6.852200  -5.463700  -1.541500   
Si  -1  -2.433800  -9.031200  -1.516200   
Si  -1  -6.004100   0.250400  -1.571000   
Si  -1   3.281200   3.821100  -1.606400   
Si  -1   2.409300  -8.421500  -3.744600   
Si  -1   0.560500  -9.618300   1.591400   
Si  -1   5.001100  -6.656800   3.794900   
Si  -1   2.707900   0.744100  -5.966100   
Si  -1   0.290900   4.448600   1.497500   
Si  -1   4.431300  -1.763100   5.923000   
Si  -1  -5.383700  -4.610900  -3.756600   
Si  -1  -0.362800   1.452700   5.902300   
Si  -1   1.221600  -6.564700   5.955400   
Si  -1   4.514500   2.008100   3.733600   
Si  -1   3.806600  -4.921900  -5.927400   
Si  -1   6.229800  -0.633700  -3.795500   
Si  -1  -3.659200  -7.147500   3.801800   
Si  -1  -2.964500  -0.353800  -5.952800   
Si  -1  -6.605500  -2.718500   1.558800   
Si  -1   7.458000  -2.450200   1.544300   
Si  -1  -1.564900   3.180800  -3.820100   
Si  -1  -4.145500   1.519700   3.745200   
Si  -1  -3.571800  -3.342900   5.942100   
Si  -1  -1.865700  -6.012100  -5.916900   
Si  -1  -6.004100   0.250400   1.571000   
Si  -1   3.281200   3.821100   1.606400   
Si  -1   6.852200  -5.463700   1.541500   
Si  -1  -2.433800  -9.031200   1.516200   
Si  -1  -1.564900   3.180800   3.820100   
Si  -1   0.290900   4.448600  -1.497500   
Si  -1  -4.145500   1.519700  -3.745200   
Si  -1  -1.865700  -6.012100   5.916900   
Si  -1   0.560500  -9.618300  -1.591400   
Si  -1  -3.571800  -3.342900  -5.942100   
Si  -1   6.229800  -0.633700   3.795500   
Si  -1   1.221600  -6.564700  -5.955400   
Si  -1  -0.362800   1.452700  -5.902300   
Si  -1  -3.659200  -7.147500  -3.801800   
Si  -1  -2.964500  -0.353800   5.952800   
Si  -1  -5.383700  -4.610900   3.756600   
Si  -1   4.514500   2.008100  -3.733600   
Si  -1   3.806600  -4.921900   5.927400   
Si  -1   7.458000  -2.450200  -1.544300   
Si  -1  -6.605500  -2.718500  -1.558800   
Si  -1   2.409300  -8.421500   3.744600   
Si  -1   5.001100  -6.656800  -3.794900   
Si  -1   4.431300  -1.763100  -5.923000   
Si  -1   2.707900   0.744100   5.966100   
H   -1   8.087000  -6.166400  -1.703900   
H   -1  -3.135600 -10.267000  -1.675400   
H   -1  -7.239200   0.950700  -1.741800   
H   -1   3.982500   5.054800  -1.782600   
H   -1   2.969800  -9.574300  -4.378700   
H   -1   0.726900 -11.027700   1.766500   
H   -1   5.841800  -7.615900   4.441400   
H   -1   3.147800   1.228800  -7.237500   
H   -1   0.124400   5.860200   1.654100   
H   -1   5.057900  -1.540600   7.189000   
H   -1  -6.533600  -5.176700  -4.391300   
H   -1  -0.576000   2.088400   7.165300   
H   -1   1.436500  -7.183700   7.226400   
H   -1   5.478500   2.853500   4.367200   
H   -1   4.300300  -5.371200  -7.191900   
H   -1   7.378800  -0.077500  -4.439800   
H   -1  -4.622300  -7.984100   4.447800   
H   -1  -3.460300   0.077400  -7.222900   
H   -1  -8.015700  -2.883900   1.728800   
H   -1   8.868600  -2.282900   1.709800   
H   -1  -2.126300   4.325000  -4.468600   
H   -1  -4.985200   2.488000   4.379500   
H   -1  -4.196600  -3.547600   7.211900   
H   -1  -2.307300  -6.513700  -7.181000   
H   -1  -7.239200   0.950700   1.741800   
H   -1   3.982500   5.054800   1.782600   
H   -1   8.087000  -6.166400   1.703900   
H   -1  -3.135600 -10.267000   1.675400   
H   -1  -2.126300   4.325000   4.468600   
H   -1   0.124400   5.860200  -1.654100   
H   -1  -4.985200   2.488000  -4.379500   
H   -1  -2.307300  -6.513700   7.181000   
H   -1   0.726900 -11.027700  -1.766500   
H   -1  -4.196600  -3.547600  -7.211900   
H   -1   7.378800  -0.077500   4.439800   
H   -1   1.436500  -7.183700  -7.226400   
H   -1  -0.576000   2.088400  -7.165300   
H   -1  -4.622300  -7.984100  -4.447800   
H   -1  -3.460300   0.077400   7.222900   
H   -1  -6.533600  -5.176700   4.391300   
H   -1   5.478500   2.853500  -4.367200   
H   -1   4.300300  -5.371200   7.191900   
H   -1   8.868600  -2.282900  -1.709800   
H   -1  -8.015700  -2.883900  -1.728800   
H   -1   2.969800  -9.574300   4.378700   
H   -1   5.841800  -7.615900  -4.441400   
H   -1   5.057900  -1.540600  -7.189000   
H   -1   3.147800   1.228800   7.237500   
C    0   0.000000   0.000000   0.000000   
H    0   1.083600   0.000000   0.000000   
C    0  -0.728500   1.183200   0.000000   
H    0  -0.208000   2.136600   0.000000   
C    0  -2.129000   1.154400   0.000000   
C    0  -2.802200  -0.068000   0.000000   
C    0  -2.077700  -1.258500   0.000000   
H    0  -2.572200  -2.224600   0.000000   
C    0  -0.677000  -1.231700   0.000000   
N    0  -0.045400  -2.501900   0.000000   
N    0   1.214700  -2.465700   0.000000   
C    0   1.846300  -3.735900   0.000000   
C    0   1.169100  -4.967800   0.000000   
H    0   0.085400  -4.967600   0.000000   
C    0   1.897400  -6.151200   0.000000   
H    0   1.376800  -7.104400   0.000000   
C    0   3.297800  -6.122600   0.000000   
C    0   3.971300  -4.900100   0.000000   
C    0   3.247100  -3.709600   0.000000   
H    0   3.741600  -2.743500   0.000000   
H    0  -3.887700  -0.093500   0.000000   
H    0   5.056700  -4.875000   0.000000   
H    0  -2.690100   2.084200   0.000000   
H    0   3.858900  -7.052500   0.000000   


谢谢给金币了,呵呵,这是我跟别人要的,也还在试着算
如果无法忘记,就不要忘记好了,真正的忘记是不需要努力的。
5楼2008-11-28 08:32:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kohn 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工301分求调剂院校 +18 刺痛jk 2026-04-06 19/950 2026-04-07 00:02 by longlong0822
[考研] 材料调剂 +14 小刘同学吖吖 2026-04-06 15/750 2026-04-06 22:37 by qlm5820
[考研] 华南理工0703化学,总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 347材料专硕求调剂 +7 zj8215216 2026-04-06 7/350 2026-04-06 21:57 by dongzh2009
[考研] 一志愿哈工大,初试329,求环境科学与工程调剂! +6 余未辛 2026-04-06 6/300 2026-04-06 20:51 by lbsjt
[考研] 土木水利专硕276分求调剂 +3 我想上学!!6 2026-04-05 6/300 2026-04-06 18:51 by 我想上学!!6
[考研] 0854求调剂 +7 亨氏番茄沙司 2026-04-06 8/400 2026-04-06 17:27 by 土木硕士招生
[考研] 301求调剂 +7 细胞相关蛋白 2026-04-03 7/350 2026-04-06 11:47 by lijunpoly
[考研] 一志愿211生物学280分 求调剂 +5 李rien 2026-04-05 5/250 2026-04-06 10:30 by zhyzzh
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 10/500 2026-04-04 23:12 by 无际的草原
[考研] 413求调剂 +4 柯某某 2026-03-31 4/200 2026-04-04 22:18 by 学员6BFVa3
[考研] 265求调剂 +17 林深温澜 2026-04-01 20/1000 2026-04-04 01:09 by userper
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 338求调剂,一志愿能源动力,外语是日语203 +5 zzz,,r 2026-04-02 5/250 2026-04-03 09:45 by 蓝云思雨
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 0805求调剂 +8 是水分 2026-03-31 8/400 2026-04-02 10:46 by guanxin1001
[考研] 08生物与医药专硕初试346找调剂 +6 dianeeee 2026-04-01 7/350 2026-04-02 08:23 by guoweigw
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
信息提示
请填处理意见