24小时热门版块排行榜     先丰纳米 买纳米材料就点这里

查看: 650  |  回复: 13

zhangzihang

金虫 (初入文坛)

[求助] 关于Gussian扫描Scan问题的求助已有1人参与

大家好,我在练习一个算例,扫描C2H6的势能面,考察旋转角的影响,输入如下
%chk=F:\Guass calculation\Rotation barreier\try\C2H6-scan.chk
# scan b3lyp/6-31g nosymm

Scan of C2H6

0,1
C
C  1  r1
H  1  r2  2  a1
H  1  r2  2  a1  3  120.0
H  1  r2  2  a1  3  240.0
H  2  r2  1  a1  4  d1
H  2  r2  1  a1  3  d1
H  2  r2  1  a1  5  d1

r1 1.5481
r2 1.0956
a1 111.7037
d1 0.0 10 12.0

但是总是提示下面的错误,我少了什么东西呢?求高手指点我应该怎么改才能算出来,谢谢

关于Gussian扫描Scan问题的求助
QQ图片20150405142318.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

小范范1989

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3, 感谢应助 2015-04-05 20:46:55
zhangzihang: 金币+5, ★★★★★最佳答案, 谢谢 2015-04-06 22:14:24
1:首先,你要确定好,你想转动哪个角度?测量Hf随哪个角度的变化而变化的势能面。
2:解决的办法如下
  first:通过高斯view找到你要转动的哪个角度,看看是哪几个原子之间的二面角:比如说为第77号21号1号11号的二面角。-161.69311为初始的二面角。
second:那么表投就这样写:#p opt=modredundant B3LYP/6-31g*  
                在输入文件的末尾加上: D 77 21 1 11 -161.69311 S 9 20.0       意思是: S为scan扫描,9 为步数。20为步长。也就是说,从-161.69311开始,每隔20度,做一次优化,一共做了9次

3:例子:
%chk=test.chk
%mem=2000MB
%nprocs=2
#p opt=modredundant B3LYP/6-31g*  

test

0 1
N              
C                  1    1.42182339
C                  2    1.41187693    1  108.55046268
C                  2    1.39732515    1  131.13987479    3 -177.31353020    0
C                  3    1.39505464    2  120.61294909    1 -179.02190081    0
C                  4    1.39291131    2  118.11991600    1  177.93272519    0
H                  4    1.07960679    2  121.12765296    1   -1.49950297    0
C                  5    1.40063869    3  120.09363495    2    0.72771633    0
H                  5    1.08632591    3  119.86964238    2 -177.99781868    0
H                  6    1.08620200    4  118.36116394    2 -178.30693988    0
C                  1    1.42316809    2  107.95367792    4 -177.66632680    0
C                 11    1.39773374    1  130.75238984    2 -177.47912075    0
C                 11    1.41348402    1  108.39791969    2   -0.19859891    0
C                 12    1.39550606   11  117.90611901    1  177.85827315    0
H                 12    1.07945426   11  121.03899554    1   -2.24183509    0
C                 13    1.39734584   11  120.34077346    1 -178.71538008    0
C                 16    1.39171268   13  119.05405080   11    0.49609107    0
H                 14    1.08665134   12  118.71639877   11  179.89565304    0
H                 16    1.08667774   13  120.25920879   11 -179.95716843    0
H                 17    1.08635764   16  119.96701143   13 -179.99164538    0
C                  1    1.38859868    2  125.95852057    4    2.58039812    0
C                 21    2.26466521    1  149.95522655    2 -161.70024141    0
C                 21    2.26474372    1  149.87452673    2   18.33190511    0
C                  8    1.48517534    5  121.10278398    3 -179.81249205    0
C                 24    1.40065804    8  120.85847741    5  141.90910954    0
C                 24    1.41480731    8  120.53264236    5  -37.75523727    0
C                 25    1.39715394   24  120.16878268    8  179.96330807    0
H                 25    1.08649711   24  119.81771191    8   -1.94007934    0
C                 26    1.39075858   24  122.43301328    8  179.92197397    0
H                 26    1.08612543   24  118.84637431    8   -1.61353609    0
C                 29    1.39715482   26  117.95828052   24    0.37693719    0
H                 29    1.08496473   26  120.66445278   24 -179.98606977    0
N                 31    1.39852652   29  129.81601119   26 -178.84668073    0
C                 33    1.40043451   31  108.29932110   29  178.04834231    0
C                 34    1.41837810   33  108.97988675   31    0.06005001    0
C                 34    1.39753910   33  129.41959510   31  178.30250415    0
C                 35    1.39966721   34  119.52663070   33  179.40243178    0
C                 36    1.39330735   34  117.71693790   33 -179.00524758    0
H                 36    1.08495400   34  121.27112064   33    0.68667418    0
C                 37    1.39184035   35  119.16190218   34   -0.31417692    0
H                 37    1.08693310   35  120.33926028   34  179.92154911    0
H                 38    1.08669867   36  119.07709229   34 -179.89813070    0
H                 40    1.08627549   37  119.84319096   35  179.89500738    0
C                 33    1.42121087   31  125.82842355   29   -1.40761080    0
C                 44    1.40076718   33  120.11550497   31  123.98160334    0
C                 44    1.40069693   33  120.10836674   31  -55.98314800    0
C                 45    1.39461717   44  119.92600860   33  179.56884420    0
H                 45    1.08552833   44  119.41365541   33   -1.08442949    0
C                 46    1.39473885   44  119.92014470   33  179.47035114    0
H                 46    1.08555961   44  119.41460171   33   -1.24143809    0
C                 49    1.39622804   46  120.35207466   44    0.96911948    0
H                 47    1.08657783   45  119.50691503   44 -179.90195456    0
H                 49    1.08656821   46  119.50280175   44 -179.88839640    0
H                 51    1.08631380   49  120.16536404   46  179.51704600    0
C                 22    1.48262325   21  149.97021600    1    1.13350816    0
C                 55    1.40365757   22  120.23797060   21 -172.11281891    0
C                 55    1.40404417   22  120.50921697   21    8.02845219    0
C                 56    1.39264233   55  120.29346431   22  179.69062639    0
H                 56    1.08404915   55  118.85278810   22    0.02044172    0
C                 57    1.39267236   55  120.27967616   22 -179.95674400    0
H                 57    1.08429200   55  119.07035143   22    0.86101545    0
C                 58    1.39706375   56  120.15478512   55    0.36554585    0
H                 58    1.08661813   56  119.76156474   55 -179.64086211    0
H                 60    1.08659700   57  119.75001694   55 -179.68066487    0
H                 62    1.08685556   58  120.07824147   56  179.85561201    0
C                 23    1.48260564   21  149.98066130    1    1.08627660    0
C                 66    1.40401950   23  120.50953694   21    7.79234781    0
C                 66    1.40366355   23  120.23619696   21 -172.37015748    0
C                 67    1.39268231   66  120.28010860   23 -179.98970523    0
H                 67    1.08427907   66  119.07167858   23    0.83730467    0
C                 68    1.39264933   66  120.29049493   23  179.74289390    0
H                 68    1.08405171   66  118.85235902   23    0.07710336    0
C                 71    1.39704887   68  120.15691393   66    0.33777514    0
H                 69    1.08658976   67  119.75229600   66 -179.68664096    0
H                 71    1.08661910   68  119.76094706   66 -179.65183209    0
H                 73    1.08684794   71  120.07648837   68  179.87316490    0
N                 21    1.33969060    1  117.52172574    2   18.01645212    0
N                 21    1.33960525    1  117.59853917    2 -161.95558389    0
N                 22    1.34186134   21   92.14009672    1 -179.96112094    0

D 77 21 1 11 -161.69311 S 9 20.0
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
2楼2015-04-05 15:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

柒月小鱼

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 小范范1989 at 2015-04-05 15:17:27
1:首先,你要确定好,你想转动哪个角度?测量Hf随哪个角度的变化而变化的势能面。
2:解决的办法如下
  first:通过高斯view找到你要转动的哪个角度,看看是哪几个原子之间的二面角:比如说为第77号21号1号11号的 ...

您好 我想问一下这个D是什么意思 可以不写吗
5楼2017-08-17 21:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zhangzihang

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by 小范范1989 at 2015-04-05 15:17:27
1:首先,你要确定好,你想转动哪个角度?测量Hf随哪个角度的变化而变化的势能面。
2:解决的办法如下
  first:通过高斯view找到你要转动的哪个角度,看看是哪几个原子之间的二面角:比如说为第77号21号1号11号的 ...

谢谢啦
3楼2015-05-18 10:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZHOUQIAOQIAO

金虫 (初入文坛)

超级玛丽

引用回帖:
2楼: Originally posted by 小范范1989 at 2015-04-05 15:17:27
1:首先,你要确定好,你想转动哪个角度?测量Hf随哪个角度的变化而变化的势能面。
2:解决的办法如下
  first:通过高斯view找到你要转动的哪个角度,看看是哪几个原子之间的二面角:比如说为第77号21号1号11号的 ...

大侠,我能问问 您 双氢键的势能面扫描吗 ???
加油!
4楼2015-11-08 10:32:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

引用回帖:
5楼: Originally posted by 柒月小鱼 at 2017-08-17 21:17:27
您好 我想问一下这个D是什么意思 可以不写吗...

D是二面角 ,后面的四个数字是这个二面角的定义。
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
6楼2017-08-18 09:07:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

柒月小鱼

木虫 (小有名气)

引用回帖:
6楼: Originally posted by 小范范1989 at 2017-08-18 09:07:59
D是二面角 ,后面的四个数字是这个二面角的定义。...

谢谢老师回复 今天我在做的时候 没有用edit  redundant coordinate 载入 然后发现没有写D 直接写这四个原子序号 也成功提交了任务 这样可以么 我等会再计算一个简单分子 对比一下

发自小木虫IOS客户端
7楼2017-08-18 12:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

引用回帖:
7楼: Originally posted by 柒月小鱼 at 2017-08-18 12:28:16
谢谢老师回复 今天我在做的时候 没有用edit  redundant coordinate 载入 然后发现没有写D 直接写这四个原子序号 也成功提交了任务 这样可以么 我等会再计算一个简单分子 对比一下
...

我没测试过,你测试一下看看,可能后面高斯更加智能了,应该成为默认了。
等你测试了,期待你的回复。
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
8楼2017-08-18 20:10:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

柒月小鱼

木虫 (小有名气)

引用回帖:
8楼: Originally posted by 小范范1989 at 2017-08-18 20:10:12
我没测试过,你测试一下看看,可能后面高斯更加智能了,应该成为默认了。
等你测试了,期待你的回复。...

老师您好 我确实是想早些回复您的 只是这两天在做键长扫描的时候 总是无法读取输入文件 后来才晓得新的高斯不允许设置初始键长B 21 58 s 50 0.0300  应该是这个模式   然后要回复您的是 确实新的高斯很智能 我没有用edit redundant coordinate 载入也能成功读取 即 21 58 s 50 0.0300 也能成功扫描 再次感谢您的回复
9楼2017-08-21 09:10:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

引用回帖:
9楼: Originally posted by 柒月小鱼 at 2017-08-21 09:10:51
老师您好 我确实是想早些回复您的 只是这两天在做键长扫描的时候 总是无法读取输入文件 后来才晓得新的高斯不允许设置初始键长B 21 58 s 50 0.0300  应该是这个模式   然后要回复您的是 确实新的高斯很智能 我没有 ...

好的,谢谢你的回复。
这个初始数值从高斯09C还是哪个版本开始,就不允许添加 初始数值了。你是高斯09D吧?
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
10楼2017-08-21 20:55:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangzihang 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 现在知道基金结果的同志是通过什么渠道得到消息的啊?求助 +6 hanzybit 2018-08-14 8/400 2018-08-14 20:44 by gaoyt
[公派出国] 郑大药学硕士小愿望 +5 15639724387 2018-08-14 7/350 2018-08-14 20:32 by Edwin123
[论文投稿] 今年中科院分区材料类会单独作为大类吗? +4 maohui1988 2018-08-13 8/400 2018-08-14 08:16 by maohui1988
[基金申请] 青基已挂,刚刚才知道 +44 wilfred123 2018-08-10 50/2500 2018-08-14 07:50 by LENA_310
[基金申请] 不中的还在反复刷代码,中的已经准备好几篇论文可以标注了 +9 枢纽 2018-08-13 10/500 2018-08-14 01:40 by flyresearch
[公派出国] 关于2019csc申请,雅思要求/语言证明 +9 Felix_1014 2018-08-11 13/650 2018-08-13 20:33 by Xu_Zhou
[教师之家] 官科其实不比民科强多少——学术界伪科学盛行 +5 取万军首级 2018-08-12 11/550 2018-08-13 19:19 by 取万军首级
[精细化工] 甘草酸二钾和卡波姆 +4 xuzhengnjut 2018-08-12 8/400 2018-08-13 12:19 by 魔兽面膜
[论文投稿] 求写sci论文的大神带飞 +9 Chenmeng 2018-08-12 11/550 2018-08-13 09:15 by P.F.Yin
[育儿交流] 有儿初长成,让他知道成功唯有坚持! (金币+1) +3 wukongwang 2018-08-11 3/150 2018-08-13 00:02 by ymygdd
[考研] 宿舍 +5 40770 2018-08-08 12/600 2018-08-12 19:50 by 禾子一起
[电化学] XPS +4 张玉燕 2018-08-10 7/350 2018-08-12 15:48 by 东交梦厦
[基金申请] 刚发现原来8月17号是七夕呀! +18 独自的奋斗 2018-08-11 22/1100 2018-08-12 13:42 by 独自的奋斗
[考博] 中科院词汇备考问题 +4 liruogan09 2018-08-09 8/400 2018-08-12 11:58 by 小-精-灵
[文学芳草园] 30而立 +11 yunhai29 2018-08-10 11/550 2018-08-11 12:05 by 悬崖挺拔松
[留学DIY] 申博推荐信的推荐人可冒名顶替吗? +6 lm_1123 2018-08-08 14/700 2018-08-10 21:43 by chengmuran
[硕博家园] 抉择 +11 相信爱情123 2018-08-09 12/600 2018-08-10 20:06 by ylc6351
[微米和纳米] 哪里有做固体核磁的? +4 leiyang0302 2018-08-10 4/200 2018-08-10 19:47 by 多啦A梦2833
[硕博家园] word2007版 +6 chenyanning 2018-08-09 8/400 2018-08-09 17:19 by 柳絮虫
[基金申请] 申请基金,成果用写合作作者的文章吗? +8 125402404 2018-08-08 8/400 2018-08-08 10:42 by chenyc828
信息提示
请填处理意见