24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 670  |  回复: 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的回帖

ccc13


zzgyb(金币+1,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
我找过渡态计算过程中,有L502出错的,手工调过渡态构型效果比较好
8楼2008-11-11 15:14:20
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 土木304求调剂 +4 兔突突突, 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突,
[考研] 环境科学与工程334分求调剂 +7 王一一依依 2026-03-30 9/450 2026-04-02 21:15 by 1104338198
[论文投稿] chinese chemical letters英文版投稿求助 120+4 Yishengeryi 2026-03-30 6/300 2026-04-02 17:19 by Yishengeryi
[考研] 348求调剂 +11 zzzzyk123 2026-04-01 11/550 2026-04-02 16:52 by Wang200018
[考博] 材料工程专业硕士申博 +3 麟正宇 2026-03-30 3/150 2026-04-02 15:04 by greychen00
[考研] 085900土木水利336分求调剂 +4 Zhangjiangj 2026-03-31 6/300 2026-04-02 11:40 by 1753564080
[考研] 301求调剂 +13 A_JiXing 2026-04-01 13/650 2026-04-02 09:01 by sanrepian
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 材料调剂 +14 一样YWY 2026-04-01 14/700 2026-04-01 21:07 by lijunpoly
[考研] 284求调剂 +12 小熊~~ 2026-03-31 12/600 2026-04-01 20:23 by 花??
[考研] 283求调剂 +9 A child 2026-03-28 9/450 2026-04-01 14:20 by Jaylen.
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[考研] 311(085601)求调剂 +12 liziyeyeye 2026-03-28 13/650 2026-04-01 00:34 by fmesaito
[考研] 289求调剂 +7 BrightLL 2026-03-29 7/350 2026-03-31 22:05 by 544594351
[考研] 考研调剂 +9 小蜡新笔 2026-03-29 10/500 2026-03-31 19:52 by Dyhoer
[考研] 276求调剂 +3 赵久华 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 一志愿双一流机械285分求调剂 +4 幸运的三木 2026-03-29 5/250 2026-03-29 14:49 by Miko19
信息提示
请填处理意见