应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】溶剂化，输入文件错误在哪里？
111/2返回列表12下一页
查看: 739  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

majun04

木虫 (著名写手)

majun

[交流] 【求助】溶剂化，输入文件错误在哪里？

在做溶剂化影响，溶剂为乙腈，自定义参数，输入文件如下：
－－－－－－－－－
%chk=BPGs_acetonitrile
%mem=400MW
%nproc=2
# opt rb3lyp/6-31+g(d) guess=save scf=tight SCRF=(iefpcm,read)

The opt of PG in acetonitrile

0,1
O,-0.7433319673,-0.0863902621,-1.2378623616
H,0.0461751547,0.0521206253,-1.7464405408
O,1.8853082878,0.2119830342,-0.3836961661
H,2.0865488667,1.120587132,-0.2059516254
C,0.8866721394,-0.2341253658,0.50449513
H,0.8937484307,-1.3146415464,0.4465503013
H,1.1168359262,0.0548795481,1.5270224637
C,-0.4794162669,0.2898370909,0.0857491385
H,-0.4627896387,1.379965466,0.1599450413
C,-1.6007247038,-0.2474780049,0.9569667588
H,-1.6268705924,-1.3312100913,0.913224437
H,-1.474331297,0.0580812904,1.9912890701
H,-2.5543144265,0.1260730782,0.6035629106

Eps=36.64
Rsolv=3.17
Density=0.7875

－－－－－－－－－－－－－－
计算提示错误，输出文件如下：
－－－－－－－－－－－－－－
Error termination via Lnk1e in /home/majun/gaussian03/g03/l301.exe at Thu Jul 17 16:20:00 2008.
－－－－－－－－－－－－－－
查手册是产生基组信息错误，可是实在是看不出哪里有错误啊，请高手帮忙指正！
回复此楼

» 猜你喜欢
乐观、执着、勤奋
1楼 2008-07-17 16:28:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)
用的也不好，呵呵。问一下，guess=save 这个选项是干啥用的，

把它去掉可不可以算？

[ Last edited by wxjbuilder on 2008-7-17 at 16:44 ]
赞一下 回复此楼
悠悠天地，奈何独立苍茫？
2楼2008-07-17 16:43:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun
引用回帖:
问一下，guess=save 这个选项是干啥用的，

把它去掉可不可以算？


这个是保存轨道信息到chk文件中，我去掉试试。

刚刚试过了，还是不行，可能问题不是出在这。还有没有别的情况？

[ Last edited by majun04 on 2008-7-17 at 16:49 ]
赞一下 回复此楼
乐观、执着、勤奋
3楼2008-07-17 16:47:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)
保存轨道到chk文件好像要调用chk文件，具体命令我看一下
赞一下 回复此楼
悠悠天地，奈何独立苍茫？
4楼2008-07-17 16:55:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)
★ ★
suntao1982(金币+2,VIP+0):xiexie
中文手册150页
--Link1--
%Chk=name
# Guess=(Read,Save,Only,NaturalOrbitals) Geom=AllCheck

你可以把你原来算出来的chk文件先做个备份，然后用
%Chk=name
# Guess=(Read,Save,Only,NaturalOrbitals) Geom=AllCheck

标题

0  1

eof
试一下，看可不可以
赞一下(10人) 回复此楼
悠悠天地，奈何独立苍茫？
5楼2008-07-17 16:58:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)
也可能是基组的问题，换个基组再试一下
赞一下 回复此楼
悠悠天地，奈何独立苍茫？
6楼2008-07-17 16:59:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun
你的这个输入文件是不是解决了轨道信息的问题？
我的输入文件的目的是在溶剂中优化结构式，后面的工作是要算频率。目前这个文件就是第一个文件。
现在是在溶剂中，自定义参数情况下，不能优化结构式。
赞一下 回复此楼
乐观、执着、勤奋
7楼2008-07-17 17:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun
引用回帖:
Originally posted by wxjbuilder at 2008-7-17 16:59:
也可能是基组的问题，换个基组再试一下

你看这个，就是直接用gaussian提供的溶剂参数的命令：
# opt rb3lyp/6-31+g(d) guess=save scf=tight SCRF=(iefpcm,solvent=acetonitrile)
这个文件我计算成功了。
但是，后来为了和实验值相比较，自己定义参数计算就没有成功。
这样看的话，是基组问题的可能性比较小吧。
赞一下 回复此楼
乐观、执着、勤奋
8楼2008-07-17 17:06:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)
★
suntao1982(金币+1,VIP+0):辛苦了，欢迎常来！
那就应该是# opt rb3lyp/6-31+g(d) guess=save scf=tight SCRF=(iefpcm,read，solvent=acetonitrile)，后面的不变
你这样再试一下吧，呵呵，我是黔驴技穷了
赞一下(10人) 回复此楼
悠悠天地，奈何独立苍茫？
9楼2008-07-17 17:13:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun
引用回帖:
Originally posted by wxjbuilder at 2008-7-17 17:13:
那就应该是# opt rb3lyp/6-31+g(d) guess=save scf=tight SCRF=(iefpcm,read，solvent=acetonitrile)，后面的不变

尝试了，还是不行，都是一个错误提示。
呵呵，谢谢你的帮助啊，祝你课题顺利。我再想想办法。
赞一下 回复此楼
乐观、执着、勤奋
10楼2008-07-17 17:17:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 majun04 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭