24小时热门版块排行榜    

查看: 620  |  回复: 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的回帖

sxllgaoli898

金虫 (小有名气)

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

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的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见