24小时热门版块排行榜    

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

xingkai-xingkai

木虫 (小有名气)


zzgyb(金币+1,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
像NH3这样小的分子,用SCRF=IPCM或SCIPCM就好了,也不用计算分子半径
7楼2008-08-03 08:13:29
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 277求调剂 +3 12A3 2026-04-02 3/150 2026-04-02 21:30 by swb0218
[考研] 310求调剂 +17 争取九点睡 2026-03-30 17/850 2026-04-02 16:40 by guanxin1001
[考博] 材料工程专业硕士申博 +3 麟正宇 2026-03-30 3/150 2026-04-02 15:04 by greychen00
[考研] 0805求调剂 +8 是水分 2026-03-31 8/400 2026-04-02 10:46 by guanxin1001
[考研] 304求调剂 +12 素年祭语 2026-03-31 15/750 2026-04-01 22:41 by peike
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 省双一流重点一本大学招收调剂 +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[考研] 086000生物与医药 初试274求调剂 +5 小叮当来了 2026-03-30 6/300 2026-04-01 14:25 by mxj1203
[考研] 085600,321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[考研] 求调剂,一志愿北林食品与营养095500,301分,已过六级,有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 277跪求调剂 +8 1915668 2026-03-27 13/650 2026-03-31 14:58 by 王亮_大连医科大
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
[考研] 328求调剂 +8 嗯滴的基本都 2026-03-27 8/400 2026-03-30 17:20 by Wang200018
[考研] 复试调剂 +3 raojunqi0129 2026-03-28 3/150 2026-03-28 15:27 by 落睿可思
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
信息提示
请填处理意见