24小时热门版块排行榜    

查看: 1045  |  回复: 4

laohengzi

木虫 (正式写手)

[求助] 欲进行势能面扫描,但只进行了一次计算,求高手指教!

欲进行二面角势能面扫描,输入文件如下
%chk=D:\anion\22\fsi-scan.chk
# b3lyp/6-31g(d) Scan

LiFSI

-1 1
N
S                  1              B1
S                  1              B2    2              A1
O                  2              B3    1              A2    3              D1
O                  3              B4    1              A3    2              D2
O                  2              B5    1              A4    3              D3
O                  3              B6    1              A5    2              D4
F                  3              B7    1              A6    2              D5
F                  2              B8    1              A7    3              D6

   B1             1.61028751
   B2             1.60924050
   B3             1.45823976
   B4             1.45955908
   B5             1.45925626
   B6             1.45839243
   B7             1.66604238
   B8             1.66603985
   A1           123.57284359
   A2           107.68033351
   A3           117.89251141
   A4           117.18416602
   A5           108.15454154
   A6           101.54938199
   A7           103.05984767
   D1           179.51051695
   D2           -68.27507510
   D3           -42.67842645
   D4            71.75662252
   D5           180.00000000
   D6            70.00000000

1 2 1.0 3 1.0
2 4 2.0 6 2.0 9 1.0
3 5 2.0 7 2.0 8 1.0
4
5
6
7
8
9


D5 180.0 S 72 5.000000
D 9 2 1 3 F




输出文件如下:
%chk=D:\anion\22\fsi-scan.chk
Default route:  MaxDisk=2000MB
---------------------
# scan b3lyp/6-31g(d)
---------------------
1/38=1,60=1/1,8;
2/17=6,18=5,29=3,40=1/2;
3/5=1,6=6,7=1,11=2,16=1,25=1,30=1,74=-5/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,18=1,28=1/1;
1/60=1/8(1);
99/9=1/99;
2/29=3/2;
3/5=1,6=6,7=1,11=2,16=1,25=1,30=1,74=-5/1,2,3;
4/5=5,16=3/1;
5/5=2,38=5/2;
1/60=1/8(-4);
99/9=1/99;
-----
LiFSI
-----
Symbolic Z-matrix:
Charge = -1 Multiplicity = 1
N
S                    1    B1
S                    1    B2       2    A1
O                    2    B3       1    A2       3    D1       0
O                    3    B4       1    A3       2    D2       0
O                    2    B5       1    A4       3    D3       0
O                    3    B6       1    A5       2    D4       0
F                    3    B7       1    A6       2    D5       0
F                    2    B8       1    A7       3    D6       0
       Variables:
  B1                    1.61029                  
  B2                    1.60924                  
  B3                    1.45824                  
  B4                    1.45956                  
  B5                    1.45926                  
  B6                    1.45839                  
  B7                    1.66604                  
  B8                    1.66604                  
  A1                  123.57284                  
  A2                  107.68033                  
  A3                  117.89251                  
  A4                  117.18417                  
  A5                  108.15454                  
  A6                  101.54938                  
  A7                  103.05985                  
  D1                  179.51052                  
  D2                  -68.27508                  
  D3                  -42.67843                  
  D4                   71.75662                  
  D5                  180.                       
  D6                   70.                       

Scan the potential surface.
Variable   Value     No. Steps Step-Size
-------- ----------- --------- ---------
A total of      1 points will be computed.


为什么没有按照设置的72步 进行呢?请高手指教!

[ Last edited by laohengzi on 2012-4-14 at 20:36 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laohengzi

木虫 (正式写手)

木有高手帮忙
2楼2012-04-15 10:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+3 2012-04-15 13:46:48
你希望D5扫描72步

然后你却在输入文件里把D5规定死了
“D5           180.00000000”

高斯无法理解你的要求
3楼2012-04-15 10:29:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laohengzi

木虫 (正式写手)

引用回帖:
3楼: Originally posted by yjcmwgk at 2012-04-15 10:29:44:
你希望D5扫描72步

然后你却在输入文件里把D5规定死了
“D5           180.00000000”

高斯无法理解你的要求

谢谢您的回复,我在输入文件中规定的D5值时D5的初始值,如果不给出初始值,高斯就无法计算了,那么请问我应该如何修改输入文件?谢谢!
4楼2012-04-15 10:44:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laohengzi

木虫 (正式写手)

问题已解决,谢谢版主
5楼2012-04-17 21:50:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 laohengzi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[硕博家园] 极限拉扯,如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见