24小时热门版块排行榜    

查看: 2234  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gougaozhan

至尊木虫 (职业作家)

busy

[求助] PCM模型求助 已有4人参与

计算用PCM模型,溶剂用乙酸酐,但Gaussian 09默认的溶剂里面没有乙酸酐,因此照说明书里的进行自定义,输入文件如下:
%mem=1500mb                           
%chk=L4-PCM2                           
#p  b3lyp/6-311++g** opt freq 5D SCRF=(SMD,Solvent=Generic,PCM=Read) nosymm

L4-PCM2

0   1
6           0       -4.086788    1.164831   -0.001087
6           0       -4.000176   -0.249678   -0.006797
7           0       -2.831932   -0.889750   -0.007595
6           0       -1.713552   -0.162361   -0.002777
6           0       -1.661624    1.252053    0.003164
6           0       -2.920554    1.892355    0.003773
7           0       -5.125477   -1.024625   -0.011694
7           0       -0.512554   -0.858557   -0.003727
6           0        0.721810   -0.211941    0.000938
6           0        0.759503    1.206613    0.006901
6           0       -0.392222    1.942626    0.008090
6           0       -0.320638   -2.234223   -0.009061
6           0        1.028195   -2.464786   -0.007804
6           0        1.712850   -1.207903   -0.001522
6           0        3.154449   -0.967357    0.001861
8           0        3.631530    0.162924    0.007240
6           0        4.068017   -2.183960   -0.001425
6           0       -0.348792    3.448367    0.014450
1           0       -5.052909    1.655368   -0.000589
1           0       -2.969196    2.974461    0.008126
1           0       -5.026303   -2.025444   -0.015939
1           0       -6.044377   -0.620813   -0.012039
1           0        1.736499    1.669112    0.010425
1           0       -1.166074   -2.898116   -0.013174
1           0        1.484174   -3.442671   -0.011091
1           0        5.104004   -1.848064    0.001718
1           0        3.885948   -2.809168    0.877958
1           0        3.889014   -2.802117   -0.886404
1           0        0.683525    3.799507    0.017543
1           0       -0.847278    3.862946    0.897337
1           0       -0.844473    3.870368   -0.866500

stoichiometry=C4H6O3
solventname=Ethanoyl ethanoate
eps=21.00
epsinf=1.9329

运行后报错如下:QPErr --- A syntax error was detected in the input line.
D,Solvent=Generic,PCM=Read) nosymm
                      '
Last state="SCR2"
TCursr=63935 LCursr=   61
Error termination via Lnk1e in F:\G09W\l1.exe at Sun Apr 27 14:52:47 2014.
Job cpu time:  0 days  0 hours  0 minutes  0.0 seconds.
File lengths (MBytes):  RWF=      1 Int=      0 D2E=      0 Chk=      1 Scr=      1
不知道是什么原因,我实在看不出来,求高手帮帮忙,看看错在哪里?或者告诉我用乙酸酐作溶剂要怎么写?非常感谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

^_^不要生气||要争气,不要看破||要突破,不要嫉妒||要欣赏,不要托延||要积极,不要心动||要行动^_^
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gougaozhan

至尊木虫 (职业作家)

busy

引用回帖:
5楼: Originally posted by 枪下游魂 at 2014-04-28 08:08:54
2楼的方法是可以算的。仔细查看说明书,SMD和Read都是PCM的选项,他们的优先级应该一样,那为何SMD可以写在scrf下面,而read却要再写在pcm下面呢?可能是说明书的一个bug。至于你说的还是一算就错,不知道具体出错 ...

Gaussian官网的例子:
# B3LYP/6-31G(d) 5D SCRF(SMD,Solvent=Generic,Read)

Water, solvation by methanol, re-defined as generic solvent.

0 1
O
H,1,0.94
H,1,0.94,2,104.5

stoichiometry=C1H4O1  Input section for PCM keywords
solventname=methanol
eps=32.63
epsinf=1.758

                      Blank line terminates PCM input
我的就照这个写,不知道错在那了
^_^不要生气||要争气,不要看破||要突破,不要嫉妒||要欣赏,不要托延||要积极,不要心动||要行动^_^
6楼2014-04-28 10:29:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

alexwpch

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
gougaozhan: 金币+5 2014-04-27 21:10:27
官网上现在例子的写法改成SCRF=(SMD,Solvent=Generic,Read),这样就可以了。
貌似是因为说明书是老旧的,G09的语法有所改变。
2楼2014-04-27 21:01:48
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gougaozhan

至尊木虫 (职业作家)

busy

引用回帖:
2楼: Originally posted by alexwpch at 2014-04-27 21:01:48
官网上现在例子的写法改成SCRF=(SMD,Solvent=Generic,Read),这样就可以了。
貌似是因为说明书是老旧的,G09的语法有所改变。

我的及尺这样写的,但是一算就错,不知道什么原因
^_^不要生气||要争气,不要看破||要突破,不要嫉妒||要欣赏,不要托延||要积极,不要心动||要行动^_^
3楼2014-04-27 21:10:19
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
3楼: Originally posted by gougaozhan at 2014-04-27 21:10:19
我的及尺这样写的,但是一算就错,不知道什么原因...

打错字,看不懂。
你的那个写法多了"PCM=",请看仔细了。可能G09的写法改变了,但是说明书还没来及修改,所以你用旧的写法会报错。
4楼2014-04-28 07:47:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见