24小时热门版块排行榜

返回列表

当前主题已经存档。

sxllgaoli898

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1172
帖子: 102
在线: 1.7小时
虫号: 514052
注册: 2008-02-28

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

计算溶剂效应时，我的输入文件为
%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 ]

回复此楼

» 猜你喜欢

固态锂金属电池迁移数计算求助已经有1人回复
荷兰奈梅亨大学赵文博与Galimberti课题组招收理论计算CSC博士 2000 欧元/月+房补已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有220人回复
26储能博士申请自荐已经有22人回复
全奖博士英国利物浦大学 × 台湾清华大学联合培养已经有0人回复
寻博导已经有11人回复
26申博求博导-第一性原理/计算化学-可提前进组已经有0人回复
上海交通大学--宁波东方理工大学电池方向博士招生已经有10人回复
基金正文30页指的是报告正文还是整个申请书已经有5人回复
海南大学全国重点实验室杨金霖老师招收储能电池方向博士生（2026年3月报名，9月入学）已经有1人回复
海南师范大学2026年博士研究生招收 (在职想提升学历人员可报考) 申请考核制+普通招考已经有0人回复

1楼 2008-07-22 17:00:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jghe

金虫 (著名写手)

应助: 10 (幼儿园)
贵宾: 0.225
金币: 112.2
散金: 5
红花: 3
帖子: 1481
在线: 123小时
虫号: 322338
注册: 2007-03-11
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 量子化学

★ ★ ★
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）忽略不收敛，继续往下做。

赞一下(10人)

回复此楼

2楼2008-07-22 21:35:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sxllgaoli898

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1172
帖子: 102
在线: 1.7小时
虫号: 514052
注册: 2008-02-28

我试过你说的方法了可是还是不行，对了，在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

应助: 0 (幼儿园)
金币: 1490
散金: 18
红花: 3
帖子: 1194
在线: 291.4小时
虫号: 312628
注册: 2006-12-31
专业: 催化化学

★ ★
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

赞一下(10人)

回复此楼

乐观、执着、勤奋