24小时热门版块排行榜    

查看: 457  |  回复: 2
当前主题已经存档。

su571

[交流] [求助]scan扫描设置

扫描AlCd的势能曲线,在-1价,多重度为1时,求能量最低时对应的键长
输入
%chk=AlCd-.chk
# bpw91/gen pseudo=read scan test nosymm  scf=qc  iop(1/8=1)
AlCd pes scan
-1 1

Al
Cd 1 r

r  2.0

Al 0
6-311++G(2df)   
****
Cd 0
S 9 1.00
0.29293500E+02 0.66660000E-02
0.16682700E+02 -0.49734000E-01
0.10423500E+02 0.13155900E+00
0.30952600E+01 -0.52438200E+00
0.70797100E+00 0.78637900E+00
0.32256700E+00 0.46395800E+00
0.69085000E-01 0.20486000E-01
0.27415000E-01 -0.33670000E-02
0.29840000E+00 0.00000000E+00
S 9 1.00
0.29293500E+02 -0.13580000E-02
0.16682700E+02 0.12392000E-01
0.10423500E+02 -0.37572000E-01
0.30952600E+01 0.17410900E+00
0.70797100E+00 -0.33850900E+00
0.32256700E+00 -0.33541700E+00
0.69085000E-01 0.65814700E+00
0.27415000E-01 0.54051000E+00
0.29840000E+00 0.00000000E+00
S 9 1.00
0.29293500E+02 0.11242000E-01
0.16682700E+02 -0.57401000E-01
0.10423500E+02 0.11064800E+00
0.30952600E+01 -0.37694300E+00
0.70797100E+00 0.12026360E+01
0.32256700E+00 -0.48758200E+00
0.69085000E-01 -0.18636390E+01
0.27415000E-01 0.17964270E+01
0.29840000E+00 0.00000000E+00
S 1 1.00
0.27415000E-01 0.10000000E+01
S 1 1.00
0.29840000E+00 0.10000000E+01
P 8 1.00
0.90271500E+01 0.20396000E-01
0.40242700E+01 -0.15366200E+00
0.99067400E+00 0.43208200E+00
0.47279700E+00 0.49834800E+00
0.21973400E+00 0.19574200E+00
0.80088000E-01 0.14083000E-01
0.30640000E-01 0.60000000E-04
0.33270000E+00 0.00000000E+00
P 8 1.00
0.90271500E+01 -0.56190000E-02
0.40242700E+01 0.46511000E-01
0.99067400E+00 -0.14952700E+00
0.47279700E+00 -0.22433700E+00
0.21973400E+00 -0.11470000E-01
0.80088000E-01 0.63649700E+00
0.30640000E-01 0.47462500E+00
0.33270000E+00 0.00000000E+00
P 8 1.00
0.90271500E+01 0.11131000E-01
0.40242700E+01 -0.96301000E-01
0.99067400E+00 0.29868500E+00
0.47279700E+00 0.88921900E+00
0.21973400E+00 -0.14854780E+01
0.80088000E-01 -0.11311200E+00
0.30640000E-01 0.89298500E+00
0.33270000E+00 0.00000000E+00
P 1 1.00
0.30640000E-01 0.10000000E+01
P 1 1.00
0.33270000E+00 0.10000000E+01
D 6 1.00
0.43045000E+01 -0.43910000E-02
0.13201900E+01 0.11968200E+00
0.56140500E+00 0.29627400E+00
0.22670200E+00 0.39087300E+00
0.87949000E-01 0.35396100E+00
0.32305000E-01 0.15016200E+00
D 6 1.00
0.43045000E+01 0.44330000E-02
0.13201900E+01 -0.14991000E+00
0.56140500E+00 -0.37329300E+00
0.22670200E+00 -0.31298700E+00
0.87949000E-01 0.42189700E+00
0.32305000E-01 0.59366000E+00
D 1 1.00
0.32305000E-01 0.10000000E+01
D 1 1.00
0.98300000E+00 0.10000000E+01
F 1 1.00
0.84990000E+00 0.10000000E+01
F 1 1.00
0.20280000E+00 0.10000000E+01
****

Cd 0
ECP28MDF,4,28
G-Komponente
1
2 1.000000 0.000000
S-G
2
2 7.858275 135.134974
2 3.382128 15.411632
P-G
4
2 6.849791 29.251437
2 6.710092 58.508363
2 3.042159 3.780243
2 2.937330 7.676547
D-G
4
2 5.416315 11.849911
2 5.333416 17.778103
2 1.976212 2.062383
2 1.961111 3.075654
F-G
2
2 5.028590 -6.928078
2 5.005582 -9.155099

以2.0(初值) 5(步数) 1(步长)具体步长步数怎么设啊
iop(1/8=1)为步长?步数又怎么设?

[ Last edited by su571 on 2008-4-1 at 14:58 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuihang

铁杆木虫 (正式写手)

VASPer


lei0736(金币+1,VIP+0):谢谢
看个大致趋势吧,步长再密也不太可能正好把最低能量键长Scan出来
聪明人知道自己是傻子,傻子以为自己是聪明人
2楼2008-04-01 11:04:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

su571

谢谢,已经可以算下去了,就是想看个大概的趋势,还要详细优化呢
3楼2008-04-01 11:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 su571 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见