24小时热门版块排行榜    

查看: 826  |  回复: 7
当前主题已经存档。

sxllgaoli898

金虫 (小有名气)

[交流] 【求助】计算溶剂效应时出现的问题

计算溶剂效应时,我的输入文件为
%mem=1000mb
%rwf=loc1,1900mb,loc2,1900mb,loc3,-1
#p b3lyp/6-31+g(d,p) SCF=Tight SCRF=(Dipole) maxdisk=170gb

NH3

0 1
N
H                  1              B1
H                  1              B2    2              A1
H                  1              B3    3              A2    2              D1

   B1             1.00000000
   B2             1.00000000
   B3             1.00000000
   A1           109.47120255
   A2           109.47125080
   D1          -119.99998525


为什么老出现Error termination via Lnk1e in d:\G03-C01\l502.exe ?

[ Last edited by zzgyb on 2008-7-23 at 09:03 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)

★ ★ ★
zzgyb(金币+3,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
自洽场不收敛,可以找找看解决自洽场不收敛的方法,很多的,转一些,如下:

对于SCF不收敛,通常有以下的解决方法:
1. 使用小基组,或低级算法计算,得到scf收敛的波函数,用guess=read读初始波函数。
2. 使用scf=qc,这个计算会慢,而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛,会提示L508错误。
3. 改变键长,一般是缩小一点,有时会有用。
4. 计算相同体系的其他电子态,比如相应的阴离子、阳离子体系或单重态体系,得到的收敛波函数作为初始猜测进行计算。

1,改变初猜,使用GUESS=CORE,HUCKEL,MIX关键词;
2,对于开壳层体系,先对于其相应的闭壳层离子计算到收敛,然后使用GUESS=READ关键词从chk文件中读取分子轨道初猜;
3,先用小基组得到分子轨道,然后同样使用GUESS=READ关键词从chk文件中读取;
4,使用能级提升法(Vshift),将HOMO和LUMO的能量间隔加大;
5,将SCF=SLEAZY收敛限放松,或者用SCFCON关键词改变收敛限;
6,对于使用了弥散函数的情形,需要提高积分精度,加关键词SCF=NoVarAcc;
7,略微改变分子几何结构,得到收敛的波函数后,再使用原来的几何构型读分子轨道计算;
8,考虑使用不同级别的计算方法和基组(小一些的);
9,使用强制收敛方法SCF=QC(经常用到)或者SCF=DM。
10,用iop(5/13=1)忽略不收敛,继续往下做。
2楼2008-07-22 21:35:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxllgaoli898

金虫 (小有名气)

我试过你说的方法了可是还是不行,对了,在Error termination via Lnk1e in d:\G03-C01\l502.exe 之前有这么一句,Epsi is zero in input, but no built-in Epsi for this solvent,什么意思啊 ?
3楼2008-07-24 08:34:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
引用回帖:
Originally posted by sxllgaoli898 at 2008-7-24 08:34:
我试过你说的方法了可是还是不行,对了,在Error termination via Lnk1e in d:\G03-C01\l502.exe 之前有这么一句,Epsi is zero in input, but no built-in Epsi for this solvent,什么意思啊 ?

在输入文件中是要给出溶剂的参数(溶剂分子半径和介电常数)的,你没有给出参数。举例介电常数为80.0,半径为1.6。
建议你用SCRF=(dipole, dielectric=80.0 ,a0=16.0 )
或者:采用你自己的输入方法SCRF=dipole,但是在分子结构信息后面,给出溶剂的参数:1.6 80.0
乐观、执着、勤奋
4楼2008-07-25 16:12:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxllgaoli898

金虫 (小有名气)

溶剂是水也需要输入吗?水的半径是多少?
5楼2008-07-28 16:08:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun


zzgyb(金币+1,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
引用回帖:
Originally posted by sxllgaoli898 at 2008-7-28 16:08:
溶剂是水也需要输入吗?水的半径是多少?

gaussian默认的溶剂是水,水溶剂是不需要输入的吧。
乐观、执着、勤奋
6楼2008-07-28 20:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingkai-xingkai

木虫 (小有名气)


zzgyb(金币+1,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
像NH3这样小的分子,用SCRF=IPCM或SCIPCM就好了,也不用计算分子半径
7楼2008-08-03 08:13:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ccc13


zzgyb(金币+1,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
我找过渡态计算过程中,有L502出错的,手工调过渡态构型效果比较好
8楼2008-11-11 15:14:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sxllgaoli898 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:13 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:01 by Aj1rhIDL5ixY
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:00 by Aj1rhIDL5ixY
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 12:35 by 40ms4Wlo4umh
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 11:13 by 23jxep3nCNZb
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-14 10:59 by 23jxep3nCNZb
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-14 10:56 by 23jxep3nCNZb
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-14 10:44 by 23jxep3nCNZb
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-14 10:32 by 23jxep3nCNZb
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
[有机交流] 科研求助 10+4 我吃一只鱼 2026-09-07 4/200 2026-09-10 08:32 by 0940628
信息提示
请填处理意见